Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ (Biot-Squirt) model where the two mechanisms are treated simultaneously, the elastic wa...Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ (Biot-Squirt) model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components (2D2C) or two-dimensions and three-components (2D3C). There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.展开更多
A 2.5D finite-difference(FD)algorithm for the modeling of the electromagnetic(EM)logging-whiledrilling(LWD)tool in anisotropic media is presented.The FD algorithm is based on the Lebedev grid,which allows for the disc...A 2.5D finite-difference(FD)algorithm for the modeling of the electromagnetic(EM)logging-whiledrilling(LWD)tool in anisotropic media is presented.The FD algorithm is based on the Lebedev grid,which allows for the discretization of the frequency-domain Maxwell's equations in the anisotropic media in 2.5D scenarios without interpolation.This leads to a system of linear equations that is solved using the multifrontal direct solver which enables the simulation of multi-sources at nearly the cost of simulating a single source for each frequency.In addition,near-optimal quadrature derived from an optimized integration path in the complex plane is employed to implement the fast inverse Fourier Transform(IFT).The algorithm is then validated by both analytic and 3D solutions.Numerical results show that two Lebedev subgrid sets are sufficient for TI medium,which is common in geosteering environments.The number of quadrature points is greatly reduced by using the near-optimal quadrature method.展开更多
We consider the interior transmission eigenvalue problem corresponding to the scattering for an anisotropic medium of the scalar Helmholtz equation in the case where the boundary?Ωis split into two disjoint parts and...We consider the interior transmission eigenvalue problem corresponding to the scattering for an anisotropic medium of the scalar Helmholtz equation in the case where the boundary?Ωis split into two disjoint parts and possesses different transmission conditions.Using the variational method,we obtain the well posedness of the interior transmission problem,which plays an important role in the proof of the discreteness of eigenvalues.Then we achieve the existence of an infinite discrete set of transmission eigenvalues provided that n≡1,where a fourth order differential operator is applied.In the case of n■1,we show the discreteness of the transmission eigenvalues under restrictive assumptions by the analytic Fredholm theory and the T-coercive method.展开更多
The thermoelastic wave propagation in a tetragonal syngony anisotropic medium of classes 4, 4/m having heterogeneity along z axis has been investigated by employing matrizant method. This medium has an axis of second-...The thermoelastic wave propagation in a tetragonal syngony anisotropic medium of classes 4, 4/m having heterogeneity along z axis has been investigated by employing matrizant method. This medium has an axis of second-order symmetry parallel to z axis. In the case of the fourth-order matrix coefficients, the problems of wave refraction and reflection on the interface of homogeneous anisotropic thermoelastic mediums are solved analytically.展开更多
In the present paper,from the second order partial differential equations for solving the magnetotelluric(MT) fields of general anisotropic medium,we first obtained the second order partial differential equations for ...In the present paper,from the second order partial differential equations for solving the magnetotelluric(MT) fields of general anisotropic medium,we first obtained the second order partial differential equations for some anisotropic media with special conductivity(e.g.diagonal anisotropy,transverse anisotropy,azimuthal anisotropy,etc.) by simplifying the electrical conductivity tensor of anisotropic medium.And then we obtained the analytic solutions to MT fields for the case of transverse and azimuthal anisotropy through converting the conductivity parameter based on that of diagonal anisotropy.We further discussed the influence of the selection of integral limit and step length on precision in solving the analytic solutions for MT fields of isotropic medium.Finally,we presented the MT responses of two transverse and azimuthal anisotropic media as well as some applications of the analytic solutions to MT fields of anisotropic medium.展开更多
Analytical expressions of electric fields inside and outside an anisotropic dielectric sphere are presented by transforming an anisotropic medium into an isotropic one based on the multi-scale transformation of electr...Analytical expressions of electric fields inside and outside an anisotropic dielectric sphere are presented by transforming an anisotropic medium into an isotropic one based on the multi-scale transformation of electromagnetic theory. The theoretical expressions are consistent with those in the literature. The inside electric field, the outside electric field and the angle between their directions are derived in detail. Numerical simulations show that the direction of the outside field influences the magnitude of the inside field, while the dielectric constant tensor greatly affects its direction.展开更多
Non-dimensionalized equations and boundary conditions are presented for the torsion problem of an anisotropic body. The error of the fundamental solution cited in some boundary element books Is pointed out after an ex...Non-dimensionalized equations and boundary conditions are presented for the torsion problem of an anisotropic body. The error of the fundamental solution cited in some boundary element books Is pointed out after an examination of the fundamental solution. Furthermore,a necessary and sufficient boundary integral equation is given for the problem and compared with the conventional boundary integral equation. Numerical results show that great errors of the boundary shear stresses obtained by the conventional boundary integral equation appear with a small error of torsion stiffness. Meanwhile,the necessary and sufficient;boundary integral equation always gives accurate results.展开更多
Through active manipulation of wavelengths,a structure exposed to a water-wave field can achieve a target hydrodynamic performance.Based on the form invariance of the governing equation for shallow water waves,wavelen...Through active manipulation of wavelengths,a structure exposed to a water-wave field can achieve a target hydrodynamic performance.Based on the form invariance of the governing equation for shallow water waves,wavelength modulators have been proposed using the space transformation method,which enables wavelength manipulation by distributing an anisotropic medium that incorporates water depth and gravitational acceleration within the modulation space.First,annular wavelength modulators were designed using the space transformation method to reduce or amplify the wavelength of shallow water waves.The control method of wavelength scaling ratios was investigated.In addition to plane waves,the wavelength modulator was applied to manipulate the wavelength of cylindrical waves.Furthermore,the interactions between a vertical cylinder and modulated water waves were studied.Results indicate that the wavelength can be arbitrarily reduced or amplified by adjusting the dimensional parameters of the modulator.Additionally,the modulator is effective for plane waves and cylindrical waves.This wavelength modulator can enable the structure to achieve the desired scattering characteristics at the target wavelength.展开更多
Fractures significantly impact oil and gas production by altering the mechanical properties of subsurface rocks.This study investigates the forward and inverse problems associated with monoclinic anisotropic media,whi...Fractures significantly impact oil and gas production by altering the mechanical properties of subsurface rocks.This study investigates the forward and inverse problems associated with monoclinic anisotropic media,which arise from the presence of two sets of inclined fractures within an isotropic matrix.We first derive the stiffness matrix for this monoclinic medium by applying Schoenberg's linear slip theory.Subsequently,we employ the point spread function and the steady phase method to derive and simplify the reflection coefficient equation for seismic waves interacting with this complex medium.Finally,we verify the rationality of the method by using synthetic seismic records and actual seismic records.The proposed equations and inversion techniques offer a more accurate framework for characterizing monoclinic anisotropic media affected by inclined fractures,thereby enhancing the interpretation of subsurface structures in hydrocarbon exploration.展开更多
Theoretical model and solutions on power line harmonic radiation (PLHR) propagating in the ground, air, and anisotropic homogeneous ionosphere are presented, The theoretical model is verified by the PLHR events obse...Theoretical model and solutions on power line harmonic radiation (PLHR) propagating in the ground, air, and anisotropic homogeneous ionosphere are presented, The theoretical model is verified by the PLHR events observed by the DEMETER satellite. Some propagation characteristics of PLHR based on the model are obtained. This paper is bene- ficial to quantitatively interpret the formation mechanism of PLHR phenomenon.展开更多
We present a novel efficient approach in calculating induced transmembrane voltage(ITV) on cells based on transformation optics. As cell membrane is much thinner than the dimension of a typical cell, discretizing th...We present a novel efficient approach in calculating induced transmembrane voltage(ITV) on cells based on transformation optics. As cell membrane is much thinner than the dimension of a typical cell, discretizing the membrane needs numerous meshes. Using an anisotropic medium based on transformation optics, the thickness of the membrane can be exaggerated by at least one order, which eliminates rigorous mesh refinement and reduces unknowns greatly. The accuracy and efficiency of the proposed method are verified by a cylindrical cell model. Moreover, the influence on ITV with bound water(BW) layers is also studied. The results show that when cells are exposed to nanosecond electric field, BW layers should be rigorously considered in calculating ITV.展开更多
We establish the existence of fundamental solutions for the anisotropic porous medium equation, ut = ∑n i=1(u^mi)xixi in R^n × (O,∞), where m1,m2,..., and mn, are positive constants satisfying min1≤i≤n{mi...We establish the existence of fundamental solutions for the anisotropic porous medium equation, ut = ∑n i=1(u^mi)xixi in R^n × (O,∞), where m1,m2,..., and mn, are positive constants satisfying min1≤i≤n{mi}≤ 1, ∑i^n=1 mi 〉 n - 2, and max1≤i≤n{mi} ≤1(2 + ∑i^n=1 mi).展开更多
The initial-boundary value problem of an anisotropic porous medium equation■is studied.Compared with the usual porous medium equation,there are two different characteristics in this equation.One lies in its anisotrop...The initial-boundary value problem of an anisotropic porous medium equation■is studied.Compared with the usual porous medium equation,there are two different characteristics in this equation.One lies in its anisotropic property,another one is that there is a nonnegative variable diffusion coefficient a(x,t)additionally.Since a(x,t)may be degenerate on the parabolic boundary∂Ω×(0,T),instead of the boundedness of the gradient|∇u|for the usual porous medium,we can only show that∇u∈L∞(0,T;L2loc(Ω)).Based on this property,the partial boundary value conditions matching up with the anisotropic porous medium equation are discovered and two stability theorems of weak solutions can be proved naturally.展开更多
To investigate the natural convective process in a hydrodynamically and thermally anisotropic porous medium at the representative elementary volume(REV)scale,the present work presented a multiplerelaxation-time lattic...To investigate the natural convective process in a hydrodynamically and thermally anisotropic porous medium at the representative elementary volume(REV)scale,the present work presented a multiplerelaxation-time lattice Boltzmann method(MRT-LBM)based on the assumption of local thermal non-equilibrium conditions(LTNE).Three sets of distribution function were used to solve the coupled momentum and heat transfer equations.One set was used to compute the flow field based on the generalized non-Darcy model;the other two sets were used to solve the temperature fields of fluid and solid under the LTNE.To describe the anisotropy of flow field of the porous media,a permeability tensor and a Forchheimer coefficient tensor were introduced into the model.Additionally,a heat conductivity tensor and a special relaxation matrix with some off-diagonal elements were selected for the thermal anisotropy.Furthermore,by selecting an appropriate equilibrium moments and discrete source terms accounting for the local thermal non-equilibrium effect,as well as choosing an off-diagonal relaxation matrix with some specific elements,the presented model can recover the exact governing equations for natural convection under LTNE with anisotropic permeability and thermal conductivity with no deviation terms through the Chapman-Enskog procedure.Finally,the proposed model was adopted to simulate several benchmark problems.Good agreements with results in the available literatures can be achieved,which indicate the wide practicability and the good accuracy of the present model.展开更多
Self-similar steady natural convection thermal boundary layer flow from a rotating vertical cone to anisotropic Darcian porous medium is investigated theoretically and numerically. The transformed non-dimensional two-...Self-similar steady natural convection thermal boundary layer flow from a rotating vertical cone to anisotropic Darcian porous medium is investigated theoretically and numerically. The transformed non-dimensional two-point boundary value problem is reduced to a system of coupled, highly nonlinear ordinary differential equations, which are solved subject to robust surface and free stream boundary conditions with the MAPLE 17 numerical quadrature software. Validation with earlier non-rotating studies is included, and also further verification of rotating solutions is achieved with a variational finite element method (FEM). The rotational (spin) parameter emerges as an inverse function of the Grashof number. The influence of this parameter, primary Darey number, secondary Darcy number and Prandtl number on tangential velocity and swirl velocity, temperature and heat transfer rate are studied in detail. It is found that the dimensionless tangential velocity increases whilst the dimensionless swirl velocity and temperature decrease with the swirl Darcy number, tangential Darcy number and the rotational parameters. The model finds applications in chemical engineering filtration processing, liquid coating and spinning cone distillation columns.展开更多
The effects of hydrodynamic anisotropy on the mixed-convection in a vertical porous channel heated on its plates with a thermal radiation are investigated analytically for fully developed flow regime.The porous medium...The effects of hydrodynamic anisotropy on the mixed-convection in a vertical porous channel heated on its plates with a thermal radiation are investigated analytically for fully developed flow regime.The porous medium is anisotropic in permeability whose principal axes are oriented in a direction that is oblique to the gravity.The generalized Brinkman-extended Darcy model which allows the no-slip boundary-condition on solid wall is used in the formulation of the problem.The flow reversal,the thermal radiation influence for natural,and forced convection are considered in the limiting cases for low and high porosity media.It was found that the anisotropic permeability ratio,the orientation angle of the principal axes of permeability and the radiation parameter affected significantly the flow regime and the heat transfer.展开更多
Pure qP-wave equations in transversely isotropic(TI)media involve complex spatial-wavenumber coupling terms,but solving them by adopting low-rank approximation and pseudo-spectral methods comes at a high numerical sim...Pure qP-wave equations in transversely isotropic(TI)media involve complex spatial-wavenumber coupling terms,but solving them by adopting low-rank approximation and pseudo-spectral methods comes at a high numerical simulation cost.Scalar operators are introduced to normalize the wavenumber vector into a unit vector aligned with the wavefield propagation direction,and the asymptotic approximation finite-difference method is combined for solution,which can transform the spatial-wavenumber coupling terms into the spatial domain for processing,thereby improving numerical simulation efficiency.However,serving as a rough estimate of the Poynting vector,asymptotic approximation yields direction vectors with alternating orientations and singularities.Additionally,the Poynting vector has obvious null points,which will affect the stability of wavefield extrapolation processes.Especially for some media with partial strong anisotropic characteristics,the results of asymptotic approximation wavefield simulation will produce numerical noise and amplitude instability.To this end,this paper adopts the stabilized Poynting vector approximation and optical flow vector approximation to replace the traditional asymptotic approximation for obtaining more stable and accurate wavefield directions.Meanwhile,scalar operators are combined to solve the pure qP-wave equation in TI media.Numerical experiments demonstrate that in some media with strong anisotropic characteristics,the method based on stable direction vectors can better suppress numerical noise in the wavefield,improve simulation accuracy,and yield higher-resolution reverse time migration(RTM)images.This study expands the application scope of this efficient finite-difference wavefield simulation method in strongly anisotropic media.展开更多
In this paper we extend the standard Ultra Weak Variational Formulation (UWVF) of Maxwell's equations in an isotropic medium to the case of an anisotropic medium. We verify that the underlying theoretical framework...In this paper we extend the standard Ultra Weak Variational Formulation (UWVF) of Maxwell's equations in an isotropic medium to the case of an anisotropic medium. We verify that the underlying theoretical framework carries over to anisotropic media (however error estimates are not yet available) and completely describe the new scheme. We then consider TM mode scattering, show how this results in a Helmholtz equation in two dimensions with an anisotropic coefficient and demonstrate how to formulate the UWVF for it. In one special case, convergence can be proved. We then show some numerical results that suggest that the UWVF can successfully simulate wave propagation in anisotropic media.展开更多
基金National Natural Science Foundation (Project number 40604013).
摘要Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ (Biot-Squirt) model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components (2D2C) or two-dimensions and three-components (2D3C). There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.
摘要A 2.5D finite-difference(FD)algorithm for the modeling of the electromagnetic(EM)logging-whiledrilling(LWD)tool in anisotropic media is presented.The FD algorithm is based on the Lebedev grid,which allows for the discretization of the frequency-domain Maxwell's equations in the anisotropic media in 2.5D scenarios without interpolation.This leads to a system of linear equations that is solved using the multifrontal direct solver which enables the simulation of multi-sources at nearly the cost of simulating a single source for each frequency.In addition,near-optimal quadrature derived from an optimized integration path in the complex plane is employed to implement the fast inverse Fourier Transform(IFT).The algorithm is then validated by both analytic and 3D solutions.Numerical results show that two Lebedev subgrid sets are sufficient for TI medium,which is common in geosteering environments.The number of quadrature points is greatly reduced by using the near-optimal quadrature method.
基金supported by the National Natural Science Foundation of China(11571132,12301542)the Natural Science Foundation of Hubei(2022CFB725)the Natural Science Foundation of Yichang(A23-2-027)。
摘要We consider the interior transmission eigenvalue problem corresponding to the scattering for an anisotropic medium of the scalar Helmholtz equation in the case where the boundary?Ωis split into two disjoint parts and possesses different transmission conditions.Using the variational method,we obtain the well posedness of the interior transmission problem,which plays an important role in the proof of the discreteness of eigenvalues.Then we achieve the existence of an infinite discrete set of transmission eigenvalues provided that n≡1,where a fourth order differential operator is applied.In the case of n■1,we show the discreteness of the transmission eigenvalues under restrictive assumptions by the analytic Fredholm theory and the T-coercive method.
摘要The thermoelastic wave propagation in a tetragonal syngony anisotropic medium of classes 4, 4/m having heterogeneity along z axis has been investigated by employing matrizant method. This medium has an axis of second-order symmetry parallel to z axis. In the case of the fourth-order matrix coefficients, the problems of wave refraction and reflection on the interface of homogeneous anisotropic thermoelastic mediums are solved analytically.
基金supported by the National Natural Science Foundation of China(Grant No.40774035)
摘要In the present paper,from the second order partial differential equations for solving the magnetotelluric(MT) fields of general anisotropic medium,we first obtained the second order partial differential equations for some anisotropic media with special conductivity(e.g.diagonal anisotropy,transverse anisotropy,azimuthal anisotropy,etc.) by simplifying the electrical conductivity tensor of anisotropic medium.And then we obtained the analytic solutions to MT fields for the case of transverse and azimuthal anisotropy through converting the conductivity parameter based on that of diagonal anisotropy.We further discussed the influence of the selection of integral limit and step length on precision in solving the analytic solutions for MT fields of isotropic medium.Finally,we presented the MT responses of two transverse and azimuthal anisotropic media as well as some applications of the analytic solutions to MT fields of anisotropic medium.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 60741003 and 60871047)
摘要Analytical expressions of electric fields inside and outside an anisotropic dielectric sphere are presented by transforming an anisotropic medium into an isotropic one based on the multi-scale transformation of electromagnetic theory. The theoretical expressions are consistent with those in the literature. The inside electric field, the outside electric field and the angle between their directions are derived in detail. Numerical simulations show that the direction of the outside field influences the magnitude of the inside field, while the dielectric constant tensor greatly affects its direction.
摘要Non-dimensionalized equations and boundary conditions are presented for the torsion problem of an anisotropic body. The error of the fundamental solution cited in some boundary element books Is pointed out after an examination of the fundamental solution. Furthermore,a necessary and sufficient boundary integral equation is given for the problem and compared with the conventional boundary integral equation. Numerical results show that great errors of the boundary shear stresses obtained by the conventional boundary integral equation appear with a small error of torsion stiffness. Meanwhile,the necessary and sufficient;boundary integral equation always gives accurate results.
基金supported by National Key R&D Program of China(2022YFC3104200)the China Postdoctoral Science Foundation(2023M742157)+3 种基金the Shandong Province Taishan Scholars Project(tsqn201909172)the Shandong Provincial Natural Science Foundation(ZR2023QA097)the Open Research Fund Program of State Key Laboratory of Coastal and Offshore Engineering,Dalian University of Technology(LP2309)the Innovation Project of Qingdao Post-Doctoral(QDBSH20230202121).
摘要Through active manipulation of wavelengths,a structure exposed to a water-wave field can achieve a target hydrodynamic performance.Based on the form invariance of the governing equation for shallow water waves,wavelength modulators have been proposed using the space transformation method,which enables wavelength manipulation by distributing an anisotropic medium that incorporates water depth and gravitational acceleration within the modulation space.First,annular wavelength modulators were designed using the space transformation method to reduce or amplify the wavelength of shallow water waves.The control method of wavelength scaling ratios was investigated.In addition to plane waves,the wavelength modulator was applied to manipulate the wavelength of cylindrical waves.Furthermore,the interactions between a vertical cylinder and modulated water waves were studied.Results indicate that the wavelength can be arbitrarily reduced or amplified by adjusting the dimensional parameters of the modulator.Additionally,the modulator is effective for plane waves and cylindrical waves.This wavelength modulator can enable the structure to achieve the desired scattering characteristics at the target wavelength.
基金National Natural Science Foundation of China under the Special Fund for Basic Research Operations of the Central Universities(U23B6010)the CNOOC“14th Five Year Plan”Major Science and Technology Project(KJGG2022-0304).
摘要Fractures significantly impact oil and gas production by altering the mechanical properties of subsurface rocks.This study investigates the forward and inverse problems associated with monoclinic anisotropic media,which arise from the presence of two sets of inclined fractures within an isotropic matrix.We first derive the stiffness matrix for this monoclinic medium by applying Schoenberg's linear slip theory.Subsequently,we employ the point spread function and the steady phase method to derive and simplify the reflection coefficient equation for seismic waves interacting with this complex medium.Finally,we verify the rationality of the method by using synthetic seismic records and actual seismic records.The proposed equations and inversion techniques offer a more accurate framework for characterizing monoclinic anisotropic media affected by inclined fractures,thereby enhancing the interpretation of subsurface structures in hydrocarbon exploration.
基金Project supported by the National Natural Science Foundation of China(Grant No.51207006)the Natural Science Foundation of Beijing,China(Grant No.3123038)
摘要Theoretical model and solutions on power line harmonic radiation (PLHR) propagating in the ground, air, and anisotropic homogeneous ionosphere are presented, The theoretical model is verified by the PLHR events observed by the DEMETER satellite. Some propagation characteristics of PLHR based on the model are obtained. This paper is bene- ficial to quantitatively interpret the formation mechanism of PLHR phenomenon.
基金Project supported by the National Key Basic Research Program of China(Grant Nos.2013CB328900 and 2013CB328905)
摘要We present a novel efficient approach in calculating induced transmembrane voltage(ITV) on cells based on transformation optics. As cell membrane is much thinner than the dimension of a typical cell, discretizing the membrane needs numerous meshes. Using an anisotropic medium based on transformation optics, the thickness of the membrane can be exaggerated by at least one order, which eliminates rigorous mesh refinement and reduces unknowns greatly. The accuracy and efficiency of the proposed method are verified by a cylindrical cell model. Moreover, the influence on ITV with bound water(BW) layers is also studied. The results show that when cells are exposed to nanosecond electric field, BW layers should be rigorously considered in calculating ITV.
基金The research is supported by National 973-Project the Trans-Century Training Programme Foundation for the Talents by the Ministry of Education
摘要We establish the existence of fundamental solutions for the anisotropic porous medium equation, ut = ∑n i=1(u^mi)xixi in R^n × (O,∞), where m1,m2,..., and mn, are positive constants satisfying min1≤i≤n{mi}≤ 1, ∑i^n=1 mi 〉 n - 2, and max1≤i≤n{mi} ≤1(2 + ∑i^n=1 mi).
基金supported by Natural Science Foundation of Fujian Province(No.2022J011242),China。
摘要The initial-boundary value problem of an anisotropic porous medium equation■is studied.Compared with the usual porous medium equation,there are two different characteristics in this equation.One lies in its anisotropic property,another one is that there is a nonnegative variable diffusion coefficient a(x,t)additionally.Since a(x,t)may be degenerate on the parabolic boundary∂Ω×(0,T),instead of the boundedness of the gradient|∇u|for the usual porous medium,we can only show that∇u∈L∞(0,T;L2loc(Ω)).Based on this property,the partial boundary value conditions matching up with the anisotropic porous medium equation are discovered and two stability theorems of weak solutions can be proved naturally.
基金supported by the National Natural Science Foundation of China(Grant No.51806067)China Postdoctoral Science Foundation(Granted No.2015M572310)+2 种基金Fundamental Research Funds for the Central Universities(Granted No.2017MS018)Guangdong Province Science and Technology projects(Grante 2017A040402005)Guangdong Bureau of Quality and Technical Supervision Science and Technology projects(Granted No.2016CT23)。
摘要To investigate the natural convective process in a hydrodynamically and thermally anisotropic porous medium at the representative elementary volume(REV)scale,the present work presented a multiplerelaxation-time lattice Boltzmann method(MRT-LBM)based on the assumption of local thermal non-equilibrium conditions(LTNE).Three sets of distribution function were used to solve the coupled momentum and heat transfer equations.One set was used to compute the flow field based on the generalized non-Darcy model;the other two sets were used to solve the temperature fields of fluid and solid under the LTNE.To describe the anisotropy of flow field of the porous media,a permeability tensor and a Forchheimer coefficient tensor were introduced into the model.Additionally,a heat conductivity tensor and a special relaxation matrix with some off-diagonal elements were selected for the thermal anisotropy.Furthermore,by selecting an appropriate equilibrium moments and discrete source terms accounting for the local thermal non-equilibrium effect,as well as choosing an off-diagonal relaxation matrix with some specific elements,the presented model can recover the exact governing equations for natural convection under LTNE with anisotropic permeability and thermal conductivity with no deviation terms through the Chapman-Enskog procedure.Finally,the proposed model was adopted to simulate several benchmark problems.Good agreements with results in the available literatures can be achieved,which indicate the wide practicability and the good accuracy of the present model.
摘要Self-similar steady natural convection thermal boundary layer flow from a rotating vertical cone to anisotropic Darcian porous medium is investigated theoretically and numerically. The transformed non-dimensional two-point boundary value problem is reduced to a system of coupled, highly nonlinear ordinary differential equations, which are solved subject to robust surface and free stream boundary conditions with the MAPLE 17 numerical quadrature software. Validation with earlier non-rotating studies is included, and also further verification of rotating solutions is achieved with a variational finite element method (FEM). The rotational (spin) parameter emerges as an inverse function of the Grashof number. The influence of this parameter, primary Darey number, secondary Darcy number and Prandtl number on tangential velocity and swirl velocity, temperature and heat transfer rate are studied in detail. It is found that the dimensionless tangential velocity increases whilst the dimensionless swirl velocity and temperature decrease with the swirl Darcy number, tangential Darcy number and the rotational parameters. The model finds applications in chemical engineering filtration processing, liquid coating and spinning cone distillation columns.
摘要The effects of hydrodynamic anisotropy on the mixed-convection in a vertical porous channel heated on its plates with a thermal radiation are investigated analytically for fully developed flow regime.The porous medium is anisotropic in permeability whose principal axes are oriented in a direction that is oblique to the gravity.The generalized Brinkman-extended Darcy model which allows the no-slip boundary-condition on solid wall is used in the formulation of the problem.The flow reversal,the thermal radiation influence for natural,and forced convection are considered in the limiting cases for low and high porosity media.It was found that the anisotropic permeability ratio,the orientation angle of the principal axes of permeability and the radiation parameter affected significantly the flow regime and the heat transfer.
基金funded by the Taishan Scholars Research Program of Shandong Provincethe Chinese Academy of Sciences Start-up Funds+4 种基金the National Natural Science Foundation of China(42374164)the Changqing Oilfield Projects“High Precision Imaging Study of Small-scale and High-angle Structures in the Western Ordos Basin”(2024D2ZZ01)“the Imaging Study of Q Least Squares Migration”(671024115010)“the Research on Image Deconvolution Technology Based on Regularization”(202418018212)“the Research on the masked autoencoder method based on vit Network”(30200020-24-ZC0613-0044)
摘要Pure qP-wave equations in transversely isotropic(TI)media involve complex spatial-wavenumber coupling terms,but solving them by adopting low-rank approximation and pseudo-spectral methods comes at a high numerical simulation cost.Scalar operators are introduced to normalize the wavenumber vector into a unit vector aligned with the wavefield propagation direction,and the asymptotic approximation finite-difference method is combined for solution,which can transform the spatial-wavenumber coupling terms into the spatial domain for processing,thereby improving numerical simulation efficiency.However,serving as a rough estimate of the Poynting vector,asymptotic approximation yields direction vectors with alternating orientations and singularities.Additionally,the Poynting vector has obvious null points,which will affect the stability of wavefield extrapolation processes.Especially for some media with partial strong anisotropic characteristics,the results of asymptotic approximation wavefield simulation will produce numerical noise and amplitude instability.To this end,this paper adopts the stabilized Poynting vector approximation and optical flow vector approximation to replace the traditional asymptotic approximation for obtaining more stable and accurate wavefield directions.Meanwhile,scalar operators are combined to solve the pure qP-wave equation in TI media.Numerical experiments demonstrate that in some media with strong anisotropic characteristics,the method based on stable direction vectors can better suppress numerical noise in the wavefield,improve simulation accuracy,and yield higher-resolution reverse time migration(RTM)images.This study expands the application scope of this efficient finite-difference wavefield simulation method in strongly anisotropic media.
摘要In this paper we extend the standard Ultra Weak Variational Formulation (UWVF) of Maxwell's equations in an isotropic medium to the case of an anisotropic medium. We verify that the underlying theoretical framework carries over to anisotropic media (however error estimates are not yet available) and completely describe the new scheme. We then consider TM mode scattering, show how this results in a Helmholtz equation in two dimensions with an anisotropic coefficient and demonstrate how to formulate the UWVF for it. In one special case, convergence can be proved. We then show some numerical results that suggest that the UWVF can successfully simulate wave propagation in anisotropic media.