The dynamics of a solid spherical body in an oscillating liquid flow in a vertical axisymmetric channel of variable cross section is experimentally studied.It is shown that the oscillating liquid leads to the generati...The dynamics of a solid spherical body in an oscillating liquid flow in a vertical axisymmetric channel of variable cross section is experimentally studied.It is shown that the oscillating liquid leads to the generation of intense averaged flows in each of the channel segments.The intensity and direction of these flows depend on the dimensionless oscillating frequency.In the region of studied frequencies,the dynamics of the considered body is examined when the primary vortices emerging in the flow occupy the whole region in each segment.For a fixed frequency,an increase in the oscillation amplitude leads to a phase-inclusion holding effect,i.e.,the body occupies a quasi-stationary position in one of the cells of the vertical channel,while oscillating around its average position.It is also shown that the oscillating motion of a liquid column generates an averaged force acting on the body,the magnitude of which depends on the properties of the body and its position in the channel.The quasi-stationary position is determined by the relative density and size of the body,as well as the dimensionless frequency.The behavior of the body as a function of the amplitude and frequency of fluid oscillation and relative size is discussed in detail.Such findings may be used in the future to control the position of a phase inclusion and/or to strengthen mass transfer effects in a channel of variable cross section by means of fluid oscillations.展开更多
This investigation aims to analyze the effects of heat transport characteristics in the unsteady flow of nanofluids over a moving plate caused by a moving slot factor.The BRS variable is utilized for the purpose of an...This investigation aims to analyze the effects of heat transport characteristics in the unsteady flow of nanofluids over a moving plate caused by a moving slot factor.The BRS variable is utilized for the purpose of analyzing these characteristics.The process of mathematical computation involves converting the governing partial differential equations into ordinary differential equations that have suitable similarity components.The Keller-Box technique is employed to solve the ordinary differential equations(ODEs)and derive the corresponding mathematical outcomes.Figures and tables present the relationship between growth characteristics and various parameters such as temperature,velocity,skin friction coefficient,concentration,Sherwood number,and Nusselt number.The results are assessed by comparing them to previous findings.The observation reveals that higher dimensionless reference temperature and variable values of the moving slot parameter have a suppressing effect on the velocity and temperature patterns of nanofluids.Higher values of the dimensionless reference temperature and moving slot parameter lead to enhancements in the Sherwood number,skin friction coefficient,and Nusselt number.The conductivity of the nanofluid is ultimately affected by these enhancements.展开更多
A variable coefficient viscoelastic equation with a time-varying delay in the boundary feedback and acoustic boundary conditions and nonlinear source term is considered.Under suitable assumptions, general decay result...A variable coefficient viscoelastic equation with a time-varying delay in the boundary feedback and acoustic boundary conditions and nonlinear source term is considered.Under suitable assumptions, general decay results of the energy are established via suitable Lyapunov functionals and some properties of the convex functions. Our result is obtained without imposing any restrictive growth assumption on the damping term and the elements of the matrix A and the kernel function g.展开更多
The gas temperature within hypersonic boundary layer flow is so high that the specific heat of gas is no longer a constant but relates to temperature. How variable specific heat influences on boundary layer flow stabi...The gas temperature within hypersonic boundary layer flow is so high that the specific heat of gas is no longer a constant but relates to temperature. How variable specific heat influences on boundary layer flow stability is worth researching. The effect of the variable specific heat on the stability of hypersonic boundary layer flows is studied and compared with the case of constant specific heat based on the linear stability theory. It is found that the variable specific heat indeed has some effects on the neutral curves of both the first-mode and the second-mode waves and on the maximum rate of growth also. Therefore, the relationship between specific heat and temperature should be considered in the study of the stability of the boundary layer.展开更多
A new type of displacement pile, the X-section cast-in-place concrete （XCC） pile, has recently been developed in China. Extensive field tests and laboratory experi- ments are undertaken to evaluate its performance a...A new type of displacement pile, the X-section cast-in-place concrete （XCC） pile, has recently been developed in China. Extensive field tests and laboratory experi- ments are undertaken to evaluate its performance and quantify the non-uniform deforma- tion effect （NUDE） of the X-shaped cross section during installation. This paper develops a simplified theoretical model that attempts to capture the NUDE. Based on the theory of complex variable plane elasticity, closed-form solutions of the stress and displacement for the X-shaped cavity boundary value problem are given. Subsequently, the analytical solution is used to evaluate the NUDE, the concrete filling index （CFI）, and the perimeter reduction coefficient of the XCC pile cross section. The computed results are compared with field test results, showing reasonable agreement. The present simplified theoretical model reveals the deformation mechanism of the X-shaped cavity and facilitates applica- tion of the newly developed XCC pile technique in geotechnical engineering.展开更多
The inherent compliance of continuum robots holds great promise in the fields of soft manipulation and safe human–robot interaction.This compliance reduces the risk of damage to the manipulated object and its surroun...The inherent compliance of continuum robots holds great promise in the fields of soft manipulation and safe human–robot interaction.This compliance reduces the risk of damage to the manipulated object and its surroundings.However,continuum robots possess theoretically infinite degrees of freedom,and this high flexibility usually leads to complex deformations when subjected to external forces and positional constraints.Describing these complex deformations is the main challenge in modeling continuum robots.In this study,we investigated a novel variable curvature modeling method for continuum robots,considering external forces and positional constraints.The robot configuration curve is described using the developed mechanical model,and then the robot is fitted to the curve.A ten-section continuum robot prototype with a length of 1 m was developed in order to validate the model.The feasibility and accuracy of the model were verified by the ability of the robot to reach target points and track complex trajectories with a load.This work was able to serve as a new perspective for the design analysis and motion control of continuum robots.展开更多
The boundary layer flow of viscous incompressible fluid over a stretching cylinder has been considered to study flow field and temperature field. Due to non-linearity, a numerical approach called Keller-box technique ...The boundary layer flow of viscous incompressible fluid over a stretching cylinder has been considered to study flow field and temperature field. Due to non-linearity, a numerical approach called Keller-box technique has been used to compute the values of velocity function f and temperature field at different points of dynamic region. The expressions for skin friction and Nusselt number have also been obtained. The dependence of velocity profile and temperature profile on the dimensionless parameter of practical interest has been analyzed in detail by graphs. The dependence of Skin friction and Nusselt number has been seen through tables.展开更多
A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable ...A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable coe?cients. The equations describe the thermoelastic behaviors of nonhomogeneous anisotropic materials with properties that vary smoothly from point to point in space. No restriction is imposed on the spatial variations of the thermoelastic coe?cients as long as all the requirements of the laws of physics are satis?ed. To check the validity and accuracy of the proposed numerical method, some speci?c test problems with known solutions are solved.展开更多
In the present paper, the boundary layer flow of Walters Liquid B Model over a stretching plate has been considered to solve heat flow problem with variable conductivity. First, using similarity transformation, the ve...In the present paper, the boundary layer flow of Walters Liquid B Model over a stretching plate has been considered to solve heat flow problem with variable conductivity. First, using similarity transformation, the velocity components of velocity have been obtained. Then, the heat flow problem has been considered in two ways: 1) prescribed surface temperature (PST), and 2) prescribed stretching plate heat flux (PHF) in case of variable conductivity. Due to variable conductivity, temperature profile has its two part- one mean tempera-ture and other temperature profile induced due to variable conductivity. The related results have been dis-cussed with the help of graphs.展开更多
When the air temperature reaches 600 K or higher, vibration is excited. The specific heat is not a constant but a function of temperature. Under this condition, the transition position of hypersonic sharp wedge bounda...When the air temperature reaches 600 K or higher, vibration is excited. The specific heat is not a constant but a function of temperature. Under this condition, the transition position of hypersonic sharp wedge boundary layer is predicted by using the improved eN method considering variable specific heat. The transition positions with different Mach numbers of oncoming flow, half wedge angles, and wall conditions are computed condition, the nearer to the Mach number The results show that for the same oncoming flow condition and wall transition positions of hypersonic sharp wedge boundary layer move much leading edge than those of the flat plate. The greater the oncoming flow the closer the transition position to the leading edge.展开更多
We consider a strongly non-linear degenerate parabolic-hyperbolic problem with p(x)-Laplacian diffusion flux function. We propose an entropy formulation and prove the existence of an entropy solution.
In this paper, we consider a numerical approximation for the boundary optimal control problem with the control constraint governed by a heat equation defined in a variable domain. For this variable domain problem, the...In this paper, we consider a numerical approximation for the boundary optimal control problem with the control constraint governed by a heat equation defined in a variable domain. For this variable domain problem, the boundary of the domain is moving and the shape of theboundary is defined by a known time-dependent function. By making use of the Galerkin finite element method, we first project the original optimal control problem into a semi-discrete optimal control problem governed by a system of ordinary differential equations. Then, based on the aforementioned semi-discrete problem, we apply the control parameterization method to obtain an optimal parameter selection problem governed by a lumped parameter system, which can be solved as a nonlinear optimization problem by a Sequential Quadratic Programming （SQP） algorithm. The numerical simulation is given to illustrate the effectiveness of our numerical approximation for the variable domain problem with the finite element method and the control parameterization method.展开更多
By the separation of singularity, a special Fourier series solution of the boundary value problem for plane is obtained, which can satisfy all boundary conditions and converges rapidly. II is proved that the solution ...By the separation of singularity, a special Fourier series solution of the boundary value problem for plane is obtained, which can satisfy all boundary conditions and converges rapidly. II is proved that the solution is equal to the result of separation of variables. As a result, the non-linear characteristic equations resulting from the method of separation of variables are transformed into polynomial equations that can provide a foundation for approximate computation and asymptotic analysis.展开更多
The Complex Variable Boundary Element Method (CVBEM) procedure is extended to modeling applications of the two-dimensional linear diffusion partial differential equation (PDE) on a rectangular domain. The methodology ...The Complex Variable Boundary Element Method (CVBEM) procedure is extended to modeling applications of the two-dimensional linear diffusion partial differential equation (PDE) on a rectangular domain. The methodology in this work is suitable for modeling diffusion problems with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The underpinning of the modeling approach is to decompose the global initial-boundary value problem into a steady-state component and a transient component. The steady-state component is governed by the Laplace PDE and is modeled using the Complex Variable Boundary Element Method. The transient component is governed by the linear diffusion PDE and is modeled by a linear combination of basis functions that are the products of a two-dimensional Fourier sine series and an exponential function. The global approximation function is the sum of the approximate solutions from the two components. The boundary conditions of the steady-state problem are specified to match the boundary conditions from the global problem so that the CVBEM approximation function satisfies the global boundary conditions. Consequently, the boundary conditions of the transient problem are specified to be continuously zero. The initial condition of the transient component is specified as the difference between the initial condition of the global initial-boundary value problem and the CVBEM approximation of the steady-state solution. Therefore, when the approximate solutions from the two components are summed, the resulting global approximation function approximately satisfies the global initial condition. In this work, it will be demonstrated that the coupled global approximation function satisfies the governing diffusion PDE. Lastly, a procedure for developing streamlines at arbitrary model time is discussed.展开更多
In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier ser...In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.展开更多
In this paper, we study higher order elliptic partial differential equations with variable growth, and obtain the existence of solutions in the setting of Wm,p（x） spaces by means of an abstract result for variationa...In this paper, we study higher order elliptic partial differential equations with variable growth, and obtain the existence of solutions in the setting of Wm,p（x） spaces by means of an abstract result for variational inequalities obtained by Gossez and Mustonen. Our result generalizes the corresponding one of Kováik and Rákosník.展开更多
Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with...Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.展开更多
To investigate the interaction between the tropical Pacific and China seas a variable-grid global ocean circulation model with fine grid covering the area from 20°S to 50°N and from 99° to 150°E is...To investigate the interaction between the tropical Pacific and China seas a variable-grid global ocean circulation model with fine grid covering the area from 20°S to 50°N and from 99° to 150°E is developed. Numerical computation of the annually cyclic circulation fields is performed. The results of the annual mean zonal currents and deep to abyssal western boundary currents in the equatorial Pacific Ocean are reported. The North Equatorial Current,the North Equatorial Countercurrent, the South Equatorial Current and the Equatorial Undercurrent are fairly well simulated. The model well reproduces the northward flowing abyssal western boundary current.From the model results a lower deep western boundary current east of the Bismarck-Solomon-New Hebrides Island chain at depths around 2 000 m has been found. The model results also show that the currents in the equatorial Pacific Ocean have multi-layer structures both in zonal currents and western boundary currents, indicating that the global ocean overturning thermohaline circulation appears of multi-layer pattern.展开更多
In this paper,we consider the Cahn-Hilliard-Hele-Shaw(CHHS)system with the dynamic boundary conditions,in which both the bulk and surface energy parts play important roles.The scalar auxiliary variable approach is int...In this paper,we consider the Cahn-Hilliard-Hele-Shaw(CHHS)system with the dynamic boundary conditions,in which both the bulk and surface energy parts play important roles.The scalar auxiliary variable approach is introduced for the physical system;the mass conservation and energy dissipation is proved for the CHHS system.Subsequently,a fully discrete SAV finite element scheme is proposed,with the mass conservation and energy dissipation laws established at a theoretical level.In addition,the convergence analysis and error estimate is provided for the proposed SAV numerical scheme.展开更多
Based on the precise integration method （PIM）, a coupling technique of the high order multiplication perturbation method （HOMPM） and the reduction method is proposed to solve variable coefficient singularly pertur...Based on the precise integration method （PIM）, a coupling technique of the high order multiplication perturbation method （HOMPM） and the reduction method is proposed to solve variable coefficient singularly perturbed two-point boundary value prob lems （TPBVPs） with one boundary layer. First, the inhomogeneous ordinary differential equations （ODEs） are transformed into the homogeneous ODEs by variable coefficient dimensional expansion. Then, the whole interval is divided evenly, and the transfer ma trix in each sub-interval is worked out through the HOMPM. Finally, a group of algebraic equations are given based on the relationship between the neighboring sub-intervals, which are solved by the reduction method. Numerical results show that the present method is highly efficient.展开更多
文摘The dynamics of a solid spherical body in an oscillating liquid flow in a vertical axisymmetric channel of variable cross section is experimentally studied.It is shown that the oscillating liquid leads to the generation of intense averaged flows in each of the channel segments.The intensity and direction of these flows depend on the dimensionless oscillating frequency.In the region of studied frequencies,the dynamics of the considered body is examined when the primary vortices emerging in the flow occupy the whole region in each segment.For a fixed frequency,an increase in the oscillation amplitude leads to a phase-inclusion holding effect,i.e.,the body occupies a quasi-stationary position in one of the cells of the vertical channel,while oscillating around its average position.It is also shown that the oscillating motion of a liquid column generates an averaged force acting on the body,the magnitude of which depends on the properties of the body and its position in the channel.The quasi-stationary position is determined by the relative density and size of the body,as well as the dimensionless frequency.The behavior of the body as a function of the amplitude and frequency of fluid oscillation and relative size is discussed in detail.Such findings may be used in the future to control the position of a phase inclusion and/or to strengthen mass transfer effects in a channel of variable cross section by means of fluid oscillations.
文摘This investigation aims to analyze the effects of heat transport characteristics in the unsteady flow of nanofluids over a moving plate caused by a moving slot factor.The BRS variable is utilized for the purpose of analyzing these characteristics.The process of mathematical computation involves converting the governing partial differential equations into ordinary differential equations that have suitable similarity components.The Keller-Box technique is employed to solve the ordinary differential equations(ODEs)and derive the corresponding mathematical outcomes.Figures and tables present the relationship between growth characteristics and various parameters such as temperature,velocity,skin friction coefficient,concentration,Sherwood number,and Nusselt number.The results are assessed by comparing them to previous findings.The observation reveals that higher dimensionless reference temperature and variable values of the moving slot parameter have a suppressing effect on the velocity and temperature patterns of nanofluids.Higher values of the dimensionless reference temperature and moving slot parameter lead to enhancements in the Sherwood number,skin friction coefficient,and Nusselt number.The conductivity of the nanofluid is ultimately affected by these enhancements.
文摘A variable coefficient viscoelastic equation with a time-varying delay in the boundary feedback and acoustic boundary conditions and nonlinear source term is considered.Under suitable assumptions, general decay results of the energy are established via suitable Lyapunov functionals and some properties of the convex functions. Our result is obtained without imposing any restrictive growth assumption on the damping term and the elements of the matrix A and the kernel function g.
基金Project supported by the National Natural Science Foundation of China (Nos. 10772134 and90716007)
文摘The gas temperature within hypersonic boundary layer flow is so high that the specific heat of gas is no longer a constant but relates to temperature. How variable specific heat influences on boundary layer flow stability is worth researching. The effect of the variable specific heat on the stability of hypersonic boundary layer flows is studied and compared with the case of constant specific heat based on the linear stability theory. It is found that the variable specific heat indeed has some effects on the neutral curves of both the first-mode and the second-mode waves and on the maximum rate of growth also. Therefore, the relationship between specific heat and temperature should be considered in the study of the stability of the boundary layer.
基金supported by the National Natural Science Foundation of China(No.51420105013)the State Key Laboratory for Geomechanics and Deep Underground Engineering,China University of Mining and Technology(No.SKLGDUEK1713)the Fundamental Research Funds for the Central Universities(Nos.106112017CDJXY200003 and 106112017CDJPT200001)
文摘A new type of displacement pile, the X-section cast-in-place concrete （XCC） pile, has recently been developed in China. Extensive field tests and laboratory experi- ments are undertaken to evaluate its performance and quantify the non-uniform deforma- tion effect （NUDE） of the X-shaped cross section during installation. This paper develops a simplified theoretical model that attempts to capture the NUDE. Based on the theory of complex variable plane elasticity, closed-form solutions of the stress and displacement for the X-shaped cavity boundary value problem are given. Subsequently, the analytical solution is used to evaluate the NUDE, the concrete filling index （CFI）, and the perimeter reduction coefficient of the XCC pile cross section. The computed results are compared with field test results, showing reasonable agreement. The present simplified theoretical model reveals the deformation mechanism of the X-shaped cavity and facilitates applica- tion of the newly developed XCC pile technique in geotechnical engineering.
基金Supported by National Natural Science Foundation of China(Grant Nos.51975566,61821005,U1908214)Key Research Program of Frontier Sciences,CAS,China(Grant No.ZDBS-LY-JSC011).
文摘The inherent compliance of continuum robots holds great promise in the fields of soft manipulation and safe human–robot interaction.This compliance reduces the risk of damage to the manipulated object and its surroundings.However,continuum robots possess theoretically infinite degrees of freedom,and this high flexibility usually leads to complex deformations when subjected to external forces and positional constraints.Describing these complex deformations is the main challenge in modeling continuum robots.In this study,we investigated a novel variable curvature modeling method for continuum robots,considering external forces and positional constraints.The robot configuration curve is described using the developed mechanical model,and then the robot is fitted to the curve.A ten-section continuum robot prototype with a length of 1 m was developed in order to validate the model.The feasibility and accuracy of the model were verified by the ability of the robot to reach target points and track complex trajectories with a load.This work was able to serve as a new perspective for the design analysis and motion control of continuum robots.
文摘The boundary layer flow of viscous incompressible fluid over a stretching cylinder has been considered to study flow field and temperature field. Due to non-linearity, a numerical approach called Keller-box technique has been used to compute the values of velocity function f and temperature field at different points of dynamic region. The expressions for skin friction and Nusselt number have also been obtained. The dependence of velocity profile and temperature profile on the dimensionless parameter of practical interest has been analyzed in detail by graphs. The dependence of Skin friction and Nusselt number has been seen through tables.
文摘A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable coe?cients. The equations describe the thermoelastic behaviors of nonhomogeneous anisotropic materials with properties that vary smoothly from point to point in space. No restriction is imposed on the spatial variations of the thermoelastic coe?cients as long as all the requirements of the laws of physics are satis?ed. To check the validity and accuracy of the proposed numerical method, some speci?c test problems with known solutions are solved.
文摘In the present paper, the boundary layer flow of Walters Liquid B Model over a stretching plate has been considered to solve heat flow problem with variable conductivity. First, using similarity transformation, the velocity components of velocity have been obtained. Then, the heat flow problem has been considered in two ways: 1) prescribed surface temperature (PST), and 2) prescribed stretching plate heat flux (PHF) in case of variable conductivity. Due to variable conductivity, temperature profile has its two part- one mean tempera-ture and other temperature profile induced due to variable conductivity. The related results have been dis-cussed with the help of graphs.
基金supported by the National Natural Science Foundation of China(Nos.11172203 and91216111)the National Basic Research Program of China(No.2009CB724103)
文摘When the air temperature reaches 600 K or higher, vibration is excited. The specific heat is not a constant but a function of temperature. Under this condition, the transition position of hypersonic sharp wedge boundary layer is predicted by using the improved eN method considering variable specific heat. The transition positions with different Mach numbers of oncoming flow, half wedge angles, and wall conditions are computed condition, the nearer to the Mach number The results show that for the same oncoming flow condition and wall transition positions of hypersonic sharp wedge boundary layer move much leading edge than those of the flat plate. The greater the oncoming flow the closer the transition position to the leading edge.
文摘We consider a strongly non-linear degenerate parabolic-hyperbolic problem with p(x)-Laplacian diffusion flux function. We propose an entropy formulation and prove the existence of an entropy solution.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61374096 and 61104048)the Natural Science Foundation of Zhejiang Province of China(Grant No.Y6110751)
文摘In this paper, we consider a numerical approximation for the boundary optimal control problem with the control constraint governed by a heat equation defined in a variable domain. For this variable domain problem, the boundary of the domain is moving and the shape of theboundary is defined by a known time-dependent function. By making use of the Galerkin finite element method, we first project the original optimal control problem into a semi-discrete optimal control problem governed by a system of ordinary differential equations. Then, based on the aforementioned semi-discrete problem, we apply the control parameterization method to obtain an optimal parameter selection problem governed by a lumped parameter system, which can be solved as a nonlinear optimization problem by a Sequential Quadratic Programming （SQP） algorithm. The numerical simulation is given to illustrate the effectiveness of our numerical approximation for the variable domain problem with the finite element method and the control parameterization method.
基金Supported by the National Natural Science Foundation of Chinathe Doctoral Training of the State Education Commission of China
文摘By the separation of singularity, a special Fourier series solution of the boundary value problem for plane is obtained, which can satisfy all boundary conditions and converges rapidly. II is proved that the solution is equal to the result of separation of variables. As a result, the non-linear characteristic equations resulting from the method of separation of variables are transformed into polynomial equations that can provide a foundation for approximate computation and asymptotic analysis.
文摘The Complex Variable Boundary Element Method (CVBEM) procedure is extended to modeling applications of the two-dimensional linear diffusion partial differential equation (PDE) on a rectangular domain. The methodology in this work is suitable for modeling diffusion problems with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The underpinning of the modeling approach is to decompose the global initial-boundary value problem into a steady-state component and a transient component. The steady-state component is governed by the Laplace PDE and is modeled using the Complex Variable Boundary Element Method. The transient component is governed by the linear diffusion PDE and is modeled by a linear combination of basis functions that are the products of a two-dimensional Fourier sine series and an exponential function. The global approximation function is the sum of the approximate solutions from the two components. The boundary conditions of the steady-state problem are specified to match the boundary conditions from the global problem so that the CVBEM approximation function satisfies the global boundary conditions. Consequently, the boundary conditions of the transient problem are specified to be continuously zero. The initial condition of the transient component is specified as the difference between the initial condition of the global initial-boundary value problem and the CVBEM approximation of the steady-state solution. Therefore, when the approximate solutions from the two components are summed, the resulting global approximation function approximately satisfies the global initial condition. In this work, it will be demonstrated that the coupled global approximation function satisfies the governing diffusion PDE. Lastly, a procedure for developing streamlines at arbitrary model time is discussed.
文摘In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.
文摘In this paper, we study higher order elliptic partial differential equations with variable growth, and obtain the existence of solutions in the setting of Wm,p（x） spaces by means of an abstract result for variational inequalities obtained by Gossez and Mustonen. Our result generalizes the corresponding one of Kováik and Rákosník.
文摘Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.
基金This study is supported by the National Natural Sci-ence Foundation of China under contract No.40136010the Major State Basic Research Program of China under contract No.G1999043808the Youth Fund of National“863”Project of China under contract No.2002AA639350.
文摘To investigate the interaction between the tropical Pacific and China seas a variable-grid global ocean circulation model with fine grid covering the area from 20°S to 50°N and from 99° to 150°E is developed. Numerical computation of the annually cyclic circulation fields is performed. The results of the annual mean zonal currents and deep to abyssal western boundary currents in the equatorial Pacific Ocean are reported. The North Equatorial Current,the North Equatorial Countercurrent, the South Equatorial Current and the Equatorial Undercurrent are fairly well simulated. The model well reproduces the northward flowing abyssal western boundary current.From the model results a lower deep western boundary current east of the Bismarck-Solomon-New Hebrides Island chain at depths around 2 000 m has been found. The model results also show that the currents in the equatorial Pacific Ocean have multi-layer structures both in zonal currents and western boundary currents, indicating that the global ocean overturning thermohaline circulation appears of multi-layer pattern.
基金supported by NSFC(Grant No.11871441)supported by NSF(Grant No.DMS-2012669).
文摘In this paper,we consider the Cahn-Hilliard-Hele-Shaw(CHHS)system with the dynamic boundary conditions,in which both the bulk and surface energy parts play important roles.The scalar auxiliary variable approach is introduced for the physical system;the mass conservation and energy dissipation is proved for the CHHS system.Subsequently,a fully discrete SAV finite element scheme is proposed,with the mass conservation and energy dissipation laws established at a theoretical level.In addition,the convergence analysis and error estimate is provided for the proposed SAV numerical scheme.
基金Project supported by the National Natural Science Foundation of China(Key Program)(Nos.11132004 and 51078145)
文摘Based on the precise integration method （PIM）, a coupling technique of the high order multiplication perturbation method （HOMPM） and the reduction method is proposed to solve variable coefficient singularly perturbed two-point boundary value prob lems （TPBVPs） with one boundary layer. First, the inhomogeneous ordinary differential equations （ODEs） are transformed into the homogeneous ODEs by variable coefficient dimensional expansion. Then, the whole interval is divided evenly, and the transfer ma trix in each sub-interval is worked out through the HOMPM. Finally, a group of algebraic equations are given based on the relationship between the neighboring sub-intervals, which are solved by the reduction method. Numerical results show that the present method is highly efficient.
