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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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;L^(2)_(loc)(Ω)).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.展开更多
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.展开更多
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.展开更多
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{...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/n(2 + ∑i^n=1 mi).展开更多
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 mediu...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.展开更多
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.展开更多
基金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.
基金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.
文摘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.
文摘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(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 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;L^(2)_(loc)(Ω)).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.
基金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.
基金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 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/n(2 + ∑i^n=1 mi).
文摘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.
文摘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.