From Hölder Continuous Solutions of 3D Incompressible Navier-Stokes Equations to No-Finite Time Blowup on [ 0,∞ ]

This article gives a general model using specific periodic special functions, that is, degenerate elliptic Weierstrass P functions composed with the LambertW function, whose presence in the governing equations through the forcing terms simplify the periodic Navier Stokes equations (PNS) at the centers of arbitrary r balls of the 3-Torus. The continuity equation is satisfied together with spatially periodic boundary conditions. The yicomponent forcing terms consist of a function F as part of its expression that is arbitrarily small in an r ball where it is associated with a singular forcing expression both for inviscid and viscous cases. As a result, a significant simplification occurs with a v3(vifor all velocity components) only governing PDE resulting. The extension of three restricted subspaces in each of the principal directions in the Cartesian plane is shown as the Cartesian product ℋ=Jx,t×Jy,t×Jz,t. On each of these subspaces vi,i=1,2,3is continuous and there exists a linear independent subspace associated with the argument of the W function. Here the 3-Torus is built up from each compact segment of length 2R on each of the axes on the 3 principal directions x, y, and z. The form of the scaled velocities for non zero scaled δis related to the definition of the W function such that e−W(ξ)=W(ξ)ξwhere ξdepends on t and proportional to δ→0for infinite time t. The ratio Wξis equal to 1, making the limit δ→0finite and well defined. Considering r balls where the function F=(x−ai)2+(y−bi)2+(z−ci)2−ηset equal to −1e+rwhere r>0. is such that the forcing is singular at every distance r of centres of cubes each containing an r-ball. At the centre of the balls, the forcing is infinite. The main idea is that a system of singular initial value problems with infinite forcing is to be solved for where the velocities are shown to be locally Hölder continuous. It is proven that the limit of these singular problems shifts the finite time blowup time ti∗for first and higher derivatives to t=∞thereby indicating that there is no finite time blowup. Results in the literature can provide a systematic approach to study both large space and time behaviour for singular solutions to the Navier Stokes equations. Among the references, it has been shown that mathematical tools can be applied to study the asymptotic properties of solutions.

Two regularity criteria for 3D Navier-Stokes equations in a bounded domain

We prove two new regularity criteria for the 3D incompressible Navier-Stokes equations in a bounded domain. Our results also hold for the 3D Boussinesq system with zero heat conductivity.

THE NAVIER-STOKES EQUATIONS IN STREAM LAYER AND ON STREAM SURFACE AND A DIMENSION SPLIT METHODS

In this paper,we proposal stream surface and stream layer.By using classical tensor calculus,we derive 3-D Navier-Stokes Equations(NSE)in the stream layer under semigeodesic coordinate system,Navier-Stokes equation on the stream surface and 2-D Navier-Stokes equations on a two dimensional manifold. After introducing stream function on the stream surface,a nonlinear initial-boundary value problem satisfies by stream function is obtained,existence and uniqueness of its solution are proven.Based this theory we proposal a new method called"dimension split method"to solve 3D NSE.

New algorithm for solving 3D incompressible viscous equations based on projection method

A new algorithm based on the projection method with the implicit finite difference technique was established to calculate the velocity fields and pressure.The calculation region can be divided into different regions according to Reynolds number.In the far-wall region,the thermal melt flow was calculated as Newtonian flow.In the near-wall region,the thermal melt flow was calculated as non-Newtonian flow.It was proved that the new algorithm based on the projection method with the implicit technique was correct through nonparametric statistics method and experiment.The simulation results show that the new algorithm based on the projection method with the implicit technique calculates more quickly than the solution algorithm-volume of fluid method using the explicit difference method.

An adaptive finite-difference method for seismic traveltime modeling based on 3D eikonal equation

3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.

Numerical simulation of standing wave with 3D predictor-corrector finite difference method for potential flow equations

A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.

An explicit method for numerical simulation of wave equations: 3D wave motion

In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.

Computation of the Eigenvalues of 3<i>D</i>“Charged” Integral Equations

The Rayleigh-Ritz and the inverse iteration methods are used in order to compute the eigenvalues of 3D Fredholm-Stieltjes integral equations, i.e. 3D Fredholm equations with respect to suitable Stieltjes-type measures. Some applications are shown, relevant to the problem of computing the eigenvalues of a body charged by a finite number of masses concentrated on points, curves or surfaces lying in.

THEORETICAL RESULTS ON THE EXISTENCE,REGULARITY AND ASYMPTOTIC STABILITY OF ENHANCED PULLBACK ATTRACTORS:APPLICATIONS TO 3D PRIMITIVE EQUATIONS

Several new concepts of enhanced pullback attractors for nonautonomous dynamical systems are introduced here by uniformly enhancing the compactness and attraction of the usual pullback attractors over an infinite forward time-interval under strong and weak topologies.Then we provide some theoretical results for the existence,regularity and asymptotic stability of these enhanced pullback attractors under general theoretical frameworks which can be applied to a large class of PDEs.The existence of these enhanced attractors is harder to obtain than the backward case[33],since it is difficult to uniformly control the long-time pullback behavior of the systems over the forward time-interval.As applications of our theoretical results,we consider the famous 3D primitive equations modelling the large-scale ocean and atmosphere dynamics,and prove the existence,regularity and asymptotic stability of the enhanced pullback attractors in V×V and H^(2)×H^(2) for the time-dependent forces which satisfy some weak conditions.

Modeling Microbial Decomposition in Real 3D Soil Structures Using Partial Differential Equations

Partial Differential Equations (PDEs) have been already widely used to simulate various complex phenomena in porous media. This paper is one of the first attempts to apply PDEs for simulating in real 3D structures. We apply this scheme to the specific case study of the microbial decomposition of organic matter in soil pore space. We got a 3D geometrical representation of the pore space relating to a network of volume primitives. A mesh of the pore space is then created by using the network. PDEs system is solved by free finite elements solver Freefem3d in the particular mesh. We validate our PDEs model to experimental data with 3D Computed Tomography (CT) images of soil samples. Regarding the current state of art on soil organic matter decay models, our approach allows taking into account precise 3D spatialization of the decomposition process by a pore space geometry description.

关于3D不可压Navier-Stokes方程H1正则性的注记

这篇论文主要研究了3D不可压Navier-Stokes方程解的H1正则性. 首先, 本论文给出井详细证明 了3D不可压Navier-Stokes方程解的局部适定性引理. 其次, 应用上述解的局部适定性引理, 第一, 可以严格证明解在小初始数据情形时的全局正则性. 第二, 对于最大可能的所有U0和最大可能的所 有F , 证明出3D不可压Navier-Stokes方程解的H1 正则性. 本论文强调不仅对最大可能U0和某一 固定F 这一情形, 解具有H1正则性, 而且对最大可能U0和最大可能的F 这一情形, 解同样具有H1正则性。

THREE HIGH-ORDER SPLITTING SCHEMES FOR 3D TRANSPORT EQUATION

Two high-order splitting schemes based on the idea of the operators splitting method are given. The three-dimensional advection-diffusion equation was split into several one-dimensional equations that were solved by these two schemes, only three computational grid points were needed in each direction but the accuracy reaches the spatial fourth-order. The third scheme proposed is based on the classical ADI scheme and the accuracy of the advection term of it can reach the spatial fourth-order. Finally, two typical numerical experiments show that the solutions of these three schemes compare well with that given by the analytical solution when the Peclet number is not bigger than 5.

Coupled thermo-hydro-mechanical simulation of CO2 enhanced gas recovery with an extended equation of state module for TOUGH2MP-FLAC3D

As one of the most important ways to reduce the greenhouse gas emission,carbon dioxide（CO2）enhanced gas recovery（CO2-EGR） is attractive since the gas recovery can be enhanced simultaneously with CO2sequestration.Based on the existing equation of state（EOS） module of TOUGH2 MP,extEOS7C is developed to calculate the phase partition of H2O-CO2-CH4-NaCl mixtures accurately with consideration of dissolved NaCI and brine properties at high pressure and temperature conditions.Verifications show that it can be applied up to the pressure of 100 MPa and temperature of 150℃.The module was implemented in the linked simulator TOUGH2MP-FLAC3 D for the coupled hydro-mechanical simulations.A simplified three-dimensional（3D）1/4 model（2.2 km×1 km×1 km） which consists of the whole reservoir,caprock and baserock was generated based on the geological conditions of a gas field in the North German Basin.The simulation results show that,under an injection rate of 200,000 t/yr and production rate of 200,000 sm3/d,CO2breakthrough occurred in the case with the initial reservoir pressure of 5 MPa but did not occur in the case of 42 MPa.Under low pressure conditions,the pressure driven horizontal transport is the dominant process;while under high pressure conditions,the density driven vertical flow is dominant.Under the considered conditions,the CO2-EGR caused only small pressure changes.The largest pore pressure increase（2 MPa） and uplift（7 mm） occurred at the caprock bottom induced by only CO2injection.The caprock had still the primary stress state and its integrity was not affected.The formation water salinity and temperature variations of ±20℃ had small influences on the CO2-EGR process.In order to slow down the breakthrough,it is suggested that CO2-EGR should be carried out before the reservoir pressure drops below the critical pressure of CO2.

Solving Cauchy Issues of Highly Nonlinear Elliptic Equations Using a Meshless Method

In this paper,we address 3D inverse Cauchy issues of highly nonlinear elliptic equations in large cuboids by utilizing the new 3D homogenization functions of different orders to adapt all the specified boundary data.We also add the average classification as an approximate solution to the nonlinear operator part,without requiring to cope with nonlinear equations to resolve the weighting coefficients because these constructions are owned many conditions about the true solution.The unknown boundary conditions and the result can be retrieved straightway by coping with a small-scale linear system when the outcome is described by a new 3D homogenization function,which is right to find the numerical solutions with the errors smaller than the level of noise being put on the over-specified Neumann conditions on the bottom of the cuboid.Besides,note that the new homogenization functions method(HFM)does not require dealing with the regularization and highly nonlinear equations.The robustness and accuracy of the HFM are verified by comparing the recovered results of several numerical experiments to the exact solutions in the entire region,even though a very large level of noise 50%is imposed on the over specified Neumann conditions.The numerical errors of our scheme are in the order of O(10^(−1))-O(10^(−4)).

Mathematical Wave Functions and 3D Finite Element Modelling of the Electron and Positron

The wave/particle duality of particles in Physics is well known. Particles have properties that uniquely characterize them from one another, such as mass, charge and spin. Charged particles have associated Electric and Magnetic fields. Also, every moving particle has a De Broglie wavelength determined by its mass and velocity. This paper shows that all of these properties of a particle can be derived from a single wave function equation for that particle. Wave functions for the Electron and the Positron are presented and principles are provided that can be used to calculate the wave functions of all the fundamental particles in Physics. Fundamental particles such as electrons and positrons are considered to be point particles in the Standard Model of Physics and are not considered to have a structure. This paper demonstrates that they do indeed have structure and that this structure extends into the space around the particle’s center (in fact, they have infinite extent), but with rapidly diminishing energy density with the distance from that center. The particles are formed from Electromagnetic standing waves, which are stable solutions to the Schrödinger and Classical wave equations. This stable structure therefore accounts for both the wave and particle nature of these particles. In fact, all of their properties such as mass, spin and electric charge, can be accounted for from this structure. These particle properties appear to originate from a single point at the center of the wave function structure, in the same sort of way that the Shell theorem of gravity causes the gravity of a body to appear to all originate from a central point. This paper represents the first two fully characterized fundamental particles, with a complete description of their structure and properties, built up from the underlying Electromagnetic waves that comprise these and all fundamental particles.

Saint-Venant and Navier-Stokes Models for Tsunami Simulations

Into the frame of the French TANDEM project (Tsunamis in the Atlantic and the English ChaNnel: Definition of the Effects through numerical Modelling) Principia has been working on the development and qualification of two in-house CFD software: the 2D EOLE-SV (Saint-Venant) model for simulation of large scale tsunami propagation from the source up to coastal scale and the 3D EOLE-NS (Navier-Stokes) model dedicated to tsunami coastal impact modelling. This paper presents a large range of test cases carried out into the frame of the project and dedicated to the validation of numerical codes in various tsunami wave conditions. The main aspects of phenomena such as wave generation, propagation and coastal impact are investigated on academic situations. A real case simulation is concerned as well, the devastating 2011 Tohoku event which is compared with in-situ data.

Three-dimensional acoustic wave equation modeling based on the optimal finite-difference scheme

Generally, FD coefficients can be obtained by using Taylor series expansion （TE） or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares （LS） can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional （2D） forward modeling to three-dimensional （3D） and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.

3-D整体变分视频去隔行方法

提出了一种基于3D整体变分模型的视频去隔行方法。首先根据整体变分修复模型,建立3D灰度空间的光滑性能量函数,形成对隔行区域插值的约束;然后利用变分法将最小化能量函数的问题,转化为对Euler-La-grange方程的求解;最后根据隔行视频的特性,给出了一种非迭代的离散实现方法。实验结果表明,和常见的几种去隔行方法相比,本文方法能够较好地抑制边沿锯齿和运动虚像。

自选酿酒酵母(KDLYS9-3)产β-D-葡萄糖苷酶动力学研究

自选高产β-D-葡萄糖苷酶的酿酒酵母(KDLYS9-3)在葡萄酒酿造过程中具有增强香气的效果。依据菌体生长和产酶试验,利用Logistic方程、Dose Resp方程和Nelder方程建立了菌体生长和产酶,以及菌体生长速率与酶生成速率之间关系的动力学模型,通过Origin8.0软件进行非线性拟合,并利用Lineweaver-Burk法作图测定了该菌所产β-D-葡萄糖苷酶的动力学参数K_m值和V_(max)值。结果显示:自选酿酒酵母KDLYS9-3的菌体生长与产酶的相关性为部分偶联型,动力学模型与试验值吻合度好,方程能够反映菌体生长与产酶的变化规律。菌体生长8 h后开始产酶,菌体进入对数生长期时酶大量生成,到39 h菌体生长进入稳定期,随着菌体生长进入衰亡期后酶也随之停止产生;利用Lineweaver-Burk法求得酶动力学参数K_m=8.492 579 mmol/L,V_(max)=1.030 715(μmol/L)/min。研究结果为该菌株的理论研究和实际应用奠定了基础。