检索结果-内蒙古大学图书馆

A diffuse-interface lattice Boltzmann method for the dendritic growth with thermosolutal convection

学校读者我要写书评

暂无评论

arXiv 2022年

作者： Zhan, Chengjie Chai, Zhenhua Shi, Baochang Jiang, Ping Geng, Shaoning Sun, Dongke School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan430074 China Institute of Interdisciplinary Research for Mathematics and Applied Science Huazhong University of Science and Technology Wuhan430074 China Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan430074 China The State Key Laboratory of Digital Manufacturing Equipment and Technology School of Mechanical Science and Engineering Huazhong University of Science and Technology Wuhan430074 China School of Mechanical Engineering Southeast University Nanjing211189 China

In this work, we proposed a diffuse interface model for the dendritic growth with thermosolutal convection. In this model, the sharp boundary between the fluid and solid dendrite is replaced by a thin but nonzero thickness diffuse interface, which is described by the order parameter governed by the phase-field equation for the dendritic growth. The governing equations for solute and heat transfer are modified such that the previous special treatments for source term can be avoided. To solve the model for the dendritic growth with thermosolutal convection, we also developed a diffuse-interface multi-relaxation-time lattice Boltzmann (LB) method. In this method, the order parameter in the phase-field equation is combined into the force caused by the fluid-solid interaction, and the treatment on the complex fluid-solid interface can be avoided. In addition, four LB models are developed for the phase-field equation, concentration equation, temperature equation and the Navier-Stokes equations in a unified framework. Finally, to test the present diffuse-interface LB method, we performed some simulations of the dendritic growth, and found that the numerical results are in good agreements with some previous works. © 2022, CC BY-NC-ND.

关键词： Navier Stokes equations

A high-order discrete energy decay and maximum-principle preserving scheme for time fractional Allen-Cahn equation

学校读者我要写书评

暂无评论

arXiv 2023年

作者： Zhang, Guoyu Huang, Chengming Alikhanov, Anatoly A. Yin, Baoli School of Mathematical Sciences Inner Mongolia University Hohhot010021 China School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan430074 China Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan430074 China North-Caucasus Center for Mathematical Research North-Caucasus Federal University Stavropol355017 Russia

The shifted fractional trapezoidal rule (SFTR) with a special shift is adopted to construct a finite difference scheme for the time-fractional Allen-Cahn (tFAC) equation. Some essential key properties of the weights of SFTR are explored for the first time. Based on these properties, we rigorously demonstrate the discrete energy decay property and maximum-principle preservation for the scheme. Numerical investigations show that the scheme can resolve the intrinsic initial singularity of such nonlinear fractional equations as tFAC equation on uniform meshes without any correction. Comparison with the classic fractional BDF2 and L2-1σ method further validates the superiority of SFTR in solving the tFAC equation. Experiments concerning both discrete energy decay and discrete maximum-principle also verify the correctness of the theoretical results. © 2023, CC BY-NC-SA.

关键词： Maximum principle

Asymptotic behavior of the principal eigenvalue of a linear second order elliptic operator with small/large diffusion coefficient and its application

学校读者我要写书评

暂无评论

arXiv 2018年

作者： Peng, Rui Zhang, Guanghui Zhou, Maolin School of Mathematics and Statistics Jiangsu Normal University Xuzhou Jiangsu Province221116 China School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan430074 China Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan430074 China School of Science and Technology University of New England ArmidaleNSW2351 Australia

In this article, we are concerned with the following eigenvalue problem of a linear second order elliptic operator: − D∆φ − 2α∇m(x) · ∇φ + V (x)φ = λφ in Ω, complemented by a general boundary condition including Dirichlet boundary condition and Robin boundary condition: ∂n∂φ + β(x)φ = 0 on ∂Ω, where β ∈ C(∂Ω) allows to be positive, sign-changing or negative, and n(x) is the unit exterior normal to ∂Ω at x. The domain Ω ⊂ RN is bounded and smooth, the constants D > 0 and α > 0 are, respectively, the diffusive and advection coefficients, and m ∈ C2(Ω) ¯ , V ∈ C(Ω) ¯ are given functions. We aim to investigate the asymptotic behavior of the principal eigenvalue of the above eigenvalue problem as the diffusive coefficient D → 0 or D → ∞. Our results, together with those of [4, 5, 10] where the Nuemann boundary case (i.e., β = 0 on ∂Ω) and Dirichlet boundary case were studied, reveal the important effect of advection and boundary conditions on the asymptotic behavior of the principal eigenvalue. We also apply our results to a reaction-diffusion-advection equation which is used to describe the evolution of a single species living in a heterogeneous stream environment and show some interesting behaviors of the species persistence and extinction caused by the buffer zone and small/large diffusion rate. Copyright © 2018, The Authors. All rights reserved.

关键词： Boundary conditions

Convergence analysis of one-point large deviations rate functions of numerical discretizations for stochastic wave equations with small noise

学校读者我要写书评

暂无评论

arXiv 2022年

作者： Jin, Diancong Hong, Jialin Sheng, Derui School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan430074 China Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan430074 China Academy of Mathematics and Systems Science Chinese Academy of Sciences Beijing100190 China School of Mathematical Sciences University of Chinese Academy of Sciences Beijing100049 China

In this work, we present the convergence analysis of one-point large deviations rate functions (LDRFs) of the spatial finite difference method (FDM) for stochastic wave equations with small noise, which is essentially about the asymptotical limit of minimization problems and not a trivial task for the nonlinear cases. In order to overcome the difficulty that objective functions for the original equation and the spatial FDM have different effective domains, we propose a new technical route for analyzing the pointwise convergence of the one-point LDRFs of the spatial FDM, based on the Γ -convergence of objective functions. Based on the new technical route, the intractable convergence analysis of one-point LDRFs boils down to the qualitative analysis of skeleton equations of the original equation and its numerical discretizations. Copyright © 2022, The Authors. All rights reserved.

关键词： Finite difference method

Phase-field-based lattice Boltzmann method for the transport of insoluble surfactant in two-phase flows

学校读者我要写书评

暂无评论

arXiv 2024年

作者： Zhan, Chengjie Liang, Hong Chai, Zhenhua Shi, Baochang School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan430074 China Institute of Interdisciplinary Research for Mathematics and Applied Science Huazhong University of Science and Technology Wuhan430074 China Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan430074 China Department of Physics Hangzhou Dianzi University Hangzhou310018 China

In this work, we present a general second-order phase-field model for the transport of insoluble surfactant in incompressible two-phase flows. In this model, the second-order local Allen-Cahn equation is applied for interface capturing, a general form of the simple scalar transport equation [S. S. Jain, J. Comput. Phys. 515, 113277 (2024)] is adopted for interface-confined surfactant, and the consistent and conservative Navier-Stokes equations with the Marangoni force is used for fluid flows. To solve this model, we further developed a mesoscopic lattice Boltzmann (LB) method, in which the LB model for surfactant transport equation is proposed under the general LB framework for the convection-diffusion type equation, and it can correctly recover the governing equation for surfactant transport. The accuracy of the present LB method is tested by several benchmark problems, and the numerical results show it has a good performance for the transport of the insoluble surfactant in two-phase flows. © 2024, CC BY-NC-ND.

关键词： Navier Stokes equations

Deepstsnet: Reconstructing the Quantum State-Resolved Thermochemical Nonequilibrium Flowfield Using Deep Neural Operator Learning With Scarce Data

学校读者我要写书评

暂无评论

SSRN

SSRN 2023年

作者： Lv, Jiaqi Hong, Qizhen Wang, Xiaoyong Mao, Zhiping Sun, Quanhua State Key Laboratory of High Temperature Gas Dynamics Institute of Mechanics Chinese Academy of Sciences Beijing100190 China School of Engineering Science University of Chinese Academy of Sciences Beijing100049 China School of Mathematical Sciences Fujian Provincial Key Laboratory of Mathematical Modeling and High-Performance Scientific Computing Xiamen University Xiamen361005 China

The hypersonic flow is in a thermochemical nonequilibrium state due to the high-temperature caused by the strong shock compression. In a thermochemical nonequilibrium flow, the distribution of molecular internal energy levels strongly deviates from the equilibrium distribution (i.e., the Boltzmann distribution). It is intractable to directly obtain the microscopic nonequilibrium distribution from existed experimental measurements usually described by macroscopic field variables such as temperature or velocity. Motivated by the idea of deep multi-scale multi-physics neural network (DeepMMNet) proposed in [Mao et al, J. Comput. Phys., 2021], we develop in this paper a data assimilation framework called DeepStSNet to accurately reconstruct the quantum state-resolved thermochemical nonequilibrium flowfield by using sparse experimental measurements of vibrational temperature and pre-trained deep neural operator networks (DeepONets). In particular, we first construct several DeepONets to express the coupled dynamics between field variables in the thermochemical nonequilibrium flow and to approximate the state-to-state (StS) approach, which traces the variation of each vibrational level of molecule accurately. These proposed DeepONets are then trained by using the numerical simulation data, and would later be served as building blocks for the DeepStSNet. We demonstrate the effectiveness and accuracy of DeepONets with different test cases showing that the density and energy of vibrational groups as well as the temperature and velocity fields are predicted with high accuracy. We then extend the architectures of DeepMMNet by considering a simplified thermochemical nonequilibrium model, i.e., the 2T model, showing that the entire thermochemical nonequilibrium flowfield is well predicted by using scattered measurements of full or even partial field variables. We next consider a more accurate and complex thermochemical nonequilibrium model, i.e., the StS-CGM model, and develop a

关键词： Multiphysics

Phase retrieval for sub-Gaussian measurements

学校读者我要写书评

暂无评论

arXiv 2019年

作者： Gao, Bing Liu, Haixia Wang, Yang School of Mathematical Sciences Nankai University Tianjin300071 China School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan430074 China Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan430074 Department of Mathematics Hong Kong University of Science and Technology Clear Water Bay Kowloon Hong Kong

Generally, phase retrieval problem can be viewed as the reconstruction of a function/signal from only the magnitude of the linear measurements. These measurements can be, for example, the Fourier transform of the density function. Computationally the phase retrieval problem is very challenging. Many algorithms for phase retrieval are based on i.i.d. Gaussian random measurements. However, Gaussian random measurements remain one of the very few classes of measurements. In this paper, we develop an efficient phase retrieval algorithm for sub-gaussian random frames. We provide a general condition for measurements and develop a modified spectral initialization. In the algorithm, we first obtain a good approximation of the solution through the initialization, and from there we use Wirtinger Flow to solve for the solution. We prove that the algorithm converges to the global minimizer linearly. Copyright © 2019, The Authors. All rights reserved.

关键词： Gaussian distribution

An adjoint-free algorithm for conditional nonlinear optimal perturbations (CNOPs) via sampling

学校读者我要写书评

暂无评论

arXiv 2022年

作者： Shi, Bin Sun, Guodong State Key Laboratory of Scientific and Engineering Computing Academy of Mathematics and Systems Science Chinese Academy of Sciences Beijing100190 China State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics Institute of Atmospheric Physics Chinese Academy of Sciences Beijing100029 China University of Chinese Academy of Sciences Beijing100049 China

In this paper, we propose a sampling algorithm based on state-of-the-art statistical machine learning techniques to obtain conditional nonlinear optimal perturbations (CNOPs), which is different from traditional (deterministic) optimization methods.1 Specifically, the traditional approach is unavailable in practice, which requires numerically computing the gradient (first-order information) such that the computation cost is expensive since it needs a large number of times to run numerical models. However, the sampling approach directly reduces the gradient to the objective function value (zeroth-order information), which also avoids using the adjoint technique that is unusable for many atmosphere and ocean models and requires large amounts of storage. We show an intuitive analysis for the sampling algorithm from the law of large numbers and further present a Chernoff-type concentration inequality to rigorously characterize the degree to which the sample average probabilistically approximates the exact gradient. The experiments are implemented to obtain the CNOPs for two numerical models, the Burgers equation with small viscosity and the Lorenz-96 model. We demonstrate the CNOPs obtained with their spatial patterns, objective values, computation times and nonlinear error growth. Compared with the performance of the three approaches, all the characters for quantifying the CNOPs are nearly consistent, while the computation time using the sampling approach with fewer samples is extremely shorter. In other words, the new sampling algorithm shortens the computation time to the utmost at the cost of losing little accuracy. © 2022, CC BY-NC-SA.

关键词： Numerical models

A Consistent and Conservative Axisymmetric Lattice Boltzmann Method for Multiphase Electro-Hydrodynamic Flows

学校读者我要写书评

暂无评论

SSRN

SSRN 2023年

作者： Liu, Xi Chai, Zhenhua Shi, Baochang Yuan, Xiaolei School of Mathematics and Statistics Huazhong University of Science and Technology Wuhan430074 China Institute of Interdisciplinary Research for Mathematics and Applied Science Huazhong University of Science and Technology Wuhan430074 China Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan430074 China College of Mathematics and Information Science Hebei University Baoding071002 China

In this paper, a consistent and conservative mathematical model is first developed for multiphase electro-hydrodynamic (EHD) flows, which has some distinct features in the volume conservation, the consistency of reduction, and the consistency of mass-momentum transport. Then a two-relaxation-time axisymmetric lattice Boltzmann (LB) method is proposed for multiphase EHD flows, and the present LB method consists of three LB models, which are used to solve the axisymmetric reduction-consistent and conservative Allen-Cahn equation for the multiphase field, the axisymmetric Poisson equation for the electric potential, and the axisymmetric consistent and conservative Navier-Stokes equations for the flow field. Owing to the reduction-consistent property, the present LB method can be used to handle a set of M (1≤M≤N) EHD fluids in N-phase system. To test the capacity of the LB method, the deformations of a single two-phase droplet and a ternary-phase compound droplet under a uniform electric field are considered. It is found that different from a single droplet with two deformation modes, the compound droplet has four deformation modes, and additionally, the effects of the electric strength, the conductivity ratio and the permittivity ratio on the compound droplet are also investigated. Finally, the compound droplet in the quaternary-phase system is also explored, and there are eight deformation modes that have not been reported in the available literature. © 2023, The Authors. All rights reserved.

关键词： Navier Stokes equations