Mathematical modelling has been extensively used to measure intervention strategies for the control of contagious conditions.Alignment between different models is pivotal for furnishing strong substantiation for polic...Mathematical modelling has been extensively used to measure intervention strategies for the control of contagious conditions.Alignment between different models is pivotal for furnishing strong substantiation for policymakers because the differences in model features can impact their prognostications.Mathematical modelling has been widely used in order to better understand the transmission,treatment,and prevention of infectious diseases.Herein,we study the dynamics of a human immunodeficiency virus(HIV)infection model with four variables:S(t),I(t),C(t),and A(t)the susceptible individuals;HIV infected individuals(with no clinical symptoms of AIDS);HIV infected individuals(under ART with a viral load remaining low),and HIV infected individuals with two different incidence functions(bilinear and saturated incidence functions).A novel numerical scheme called the continuous Galerkin-Petrov method is implemented for the solution of themodel.The influence of different clinical parameters on the dynamical behavior of S(t),I(t),C(t)and A(t)is described and analyzed.All the results are depicted graphically.On the other hand,we explore the time-dependent movement of nanofluid in porous media on an extending sheet under the influence of thermal radiation,heat flux,hall impact,variable heat source,and nanomaterial.The flow is considered to be 2D,boundary layer,viscous,incompressible,laminar,and unsteady.Sufficient transformations turn governing connected PDEs intoODEs,which are solved using the proposed scheme.To justify the envisaged problem,a comparison of the current work with previous literature is presented.展开更多
In this paper,we apply local discontinuous Galerkin(LDG)methods for pattern formation dynamical model in polymerizing actin focks.There are two main dificulties in designing effective numerical solvers.First of all,th...In this paper,we apply local discontinuous Galerkin(LDG)methods for pattern formation dynamical model in polymerizing actin focks.There are two main dificulties in designing effective numerical solvers.First of all,the density function is non-negative,and zero is an unstable equilibrium solution.Therefore,negative density values may yield blow-up solutions.To obtain positive numerical approximations,we apply the positivitypreserving(PP)techniques.Secondly,the model may contain stif source.The most commonly used time integration for the PP technique is the strong-stability-preserving Runge-Kutta method.However,for problems with stiff source,such time discretizations may require strictly limited time step sizes,leading to large computational cost.Moreover,the stiff source any trigger spurious filament polarization,leading to wrong numerical approximations on coarse meshes.In this paper,we combine the PP LDG methods with the semi-implicit Runge-Kutta methods.Numerical experiments demonstrate that the proposed method can yield accurate numerical approximations with relatively large time steps.展开更多
Anti-plane electroelastic problems are studied by the Trefftz boundary element method （BEM） in this paper. The Trefftz BEM is based on a weighted residual formulation and indirect boundary approach. In particular th...Anti-plane electroelastic problems are studied by the Trefftz boundary element method （BEM） in this paper. The Trefftz BEM is based on a weighted residual formulation and indirect boundary approach. In particular the point-collocation and Galerkin techniques, in which the basic unknowns are the retained expansion coefficients in the system of complete equations, are considered. Furthermore, special Trefftz functions and auxiliary functions which satisfy exactly the specified boundary conditions along the slit boundaries are also used to derive a special purpose element with local defects. The path-independent integral is evaluated at the tip of a crack to determine the energy release rate for a mode Ⅲ fracture problem. In final, the accuracy and efficiency of the Trefftz boundary element method are illustrated by an example and the comparison is made with other methods.展开更多
This paper investigated the buoyancy and surface tension-driven ferro-thermal-convection (FTC) in a ferrofluid (FF) layer due to influence of general boundary conditions. The lower surface is rigid with insulating to ...This paper investigated the buoyancy and surface tension-driven ferro-thermal-convection (FTC) in a ferrofluid (FF) layer due to influence of general boundary conditions. The lower surface is rigid with insulating to temperature perturbations, while the upper surface is stress-free and subjected to general thermal boundary condition. The numerically Galerkin technique (GT) and analytically regular perturbation technique (RPT) are applied for solving the problem of eigenvalue. It is analyzed that increasing Biot number, decreases the magnetic and Marangoni number is to postponement the onset. Additionally, magnetization nonlinearity parameter has no effect on FTC in the non-existence of Biot number. The results under the limiting cases are found to be in good agreement with those available in the literature.展开更多
In order to obtain an expected numerical solution, a fully discrete discontinuous Galerkin method is applied to a kind of reactive transport problems in two dimension. That is to say, the space variable is discretized...In order to obtain an expected numerical solution, a fully discrete discontinuous Galerkin method is applied to a kind of reactive transport problems in two dimension. That is to say, the space variable is discretized with the symmetric interior penalty Calerkin method （SIPG）, and the time variable is done with the backward Euler method. Making use of the duality technique, hp approximation properties and the interpolation theory, a residual-type a posteriori error estimation is achieved, which can be used for adaptivity. Compared with the analyses of semi-discretization, the current presentation is more challenging and more significant.展开更多
A discontinuous Galerkin (DG) finite element method is presented to solve the thermoelastic coupling problems caused by temperature and pressure dependent thermal contact resistance (TCR).The whole analysis is made up...A discontinuous Galerkin (DG) finite element method is presented to solve the thermoelastic coupling problems caused by temperature and pressure dependent thermal contact resistance (TCR).The whole analysis is made up of two parts,thermal and mechanical analysis.In thermal analysis,the DG method is employed to simulate the temperature jump phenomenon,which satisfies the imperfect thermal contact condition in a straightforward manner.In mechanical analysis,the impenetrability condition is fulfilled through a DG approach with penalty functions.The Picard iteration procedure with a relaxation technique is also adopted to accelerate the rate of convergence and avoid numerical instability.Numerical examples show that the present method is an attractive approach for solving thermoelastic coupling problems caused by TCR.The methodology can also be expanded to solve problems with friction finite deformation contact,node-to-segment contact and node-to-surface contact,etc.in a straightforward manner.展开更多
A finite element solution for the Navier-Stokes equations for steady flow under the porosity effects through a double branched two-dimensional section of a three-dimensional model of a canine aorta was obtained. The n...A finite element solution for the Navier-Stokes equations for steady flow under the porosity effects through a double branched two-dimensional section of a three-dimensional model of a canine aorta was obtained. The numerical solution involves transforming a physical coordinates to a curvilinear boundary fitted coordinate system.The steady flow,branch flow and shear stress under the porous effects were discussed in detail. The shear stress at the wall was calculated for Reynolds number of 1 000 with branch to main aortic flow rate ratio as a parameter. The results are compared with earlier works involving experimental data and it has been observed that our results are very close to the exact solutions.This work is in fact an improvement of the work of Sharma et al. (2001) in the sense that computational technique is economic and (Reynolds) number is large.展开更多
In this paper, a brand-new wavelet-homotopy Galerkin technique is developed to solve nonlinear ordinary or partial differential equations. Before this investigation,few studies have been done for handling nonlinear pr...In this paper, a brand-new wavelet-homotopy Galerkin technique is developed to solve nonlinear ordinary or partial differential equations. Before this investigation,few studies have been done for handling nonlinear problems with non-uniform boundary conditions by means of the wavelet Galerkin technique, especially in the field of fluid mechanics and heat transfer. The lid-driven cavity flow and heat transfer are illustrated as a typical example to verify the validity and correctness of this proposed technique. The cavity is subject to the upper and lower walls’ motions in the same or opposite directions.The inclined angle of the square cavity is from 0 to π/2. Four different modes including uniform, linear, exponential, and sinusoidal heating are considered on the top and bottom walls, respectively, while the left and right walls are thermally isolated and stationary.A parametric analysis of heating distribution between upper and lower walls including the amplitude ratio from 0 to 1 and the phase deviation from 0 to 2π is conducted. The governing equations are non-dimensionalized in terms of the stream function-vorticity formulation and the temperature distribution function and then solved analytically subject to various boundary conditions. Comparisons with previous publications are given,showing high efficiency and great feasibility of the proposed technique.展开更多
The effects of electric field and non-uniform basic temperature gradient on the onset of Rayleigh-Bénard convection in a micropolar fluid are studied using the Galerkin technique. The eigenvalues are obtained for...The effects of electric field and non-uniform basic temperature gradient on the onset of Rayleigh-Bénard convection in a micropolar fluid are studied using the Galerkin technique. The eigenvalues are obtained for free-free, rigid-free and rigid-rigid velocity boundary combinations and for isothermal and/or adiabatic temperature boundaries. The microrotation is assumed to vanish at the boundaries. A linear stability analysis is performed. The influence of various micropolar fluid parameters and electric Rayleigh number on the onset of convection has been analyzed. One linear and five non-uniform temperature profiles are considered and their comparative influence on onset is discussed.展开更多
The number and distribution of the singular points of streamlines in the cross-section of steady flow through a curved tube ate discussed by using the method of topological structure analysis. And a theoretical criter...The number and distribution of the singular points of streamlines in the cross-section of steady flow through a curved tube ate discussed by using the method of topological structure analysis. And a theoretical criterion is obtained for the bifurcation of flow vortexes for the secondary flow turning from two-vortex structure into four-vortex structure. Furthermore, the critical Dean number for bifurcation and the semi-analytical expressions of stream function and axial velocity are given by using Galerkin technique. The result of calculation is consistent with the theoretical criterion.展开更多
The effects of electric field and non-uniform basic temperature gradient on the onset of Rayleigh-Bénard-Marangoni convection in a micropolar fluid are studied using the Galerkin technique. The eigenvalues are ob...The effects of electric field and non-uniform basic temperature gradient on the onset of Rayleigh-Bénard-Marangoni convection in a micropolar fluid are studied using the Galerkin technique. The eigenvalues are obtained for an upper free/adiabatic and lower rigid/isothermal boundaries. The microrotation is assumed to vanish at the boundaries. A linear stability analysis is performed. The influence of various micropolar fluid parameters and electric Rayleigh number on the onset of convection has been analysed. Six different non-uniform temperature profiles are considered and their comparative influence on onset is discussed.展开更多
The effect of viscosity depending exponentially on temperature on the onset of penetrative ferro-thermal-convection (FTC) in a saturated horizontal porous layer in the presence of vertical magnetic field is investigat...The effect of viscosity depending exponentially on temperature on the onset of penetrative ferro-thermal-convection (FTC) in a saturated horizontal porous layer in the presence of vertical magnetic field is investigated. The bounding surface of the ferrofluid layer is considered to be rigid-rigid and insulated to temperature perturbations. The resulting eigenvalue problem is solved numerically using the Galerkin technique and also analytically by a regular perturbation technique with wave number as a perturbation parameter. The analytical and numerical results are found to be concurrence. The characteristics of stability of the system are strongly dependent on the viscosity parameter B. The effect of B on the onset of ferroconvection in a porous layer is dual in nature depending on the choices of physical parameters and a sublayer starts to form at higher values of B. Whereas, increase in magnetic number M1 and the Darcy number Da is to advance the onset of ferroconvection in a porous layer. The nonlinearity of fluid magnetization M3?is found to have no influence on the onset of ferroconvection.展开更多
Composite cylindrical shells,as key components,are widely employed in large rotating machines.However,due to the frequency bifurcations and dense frequency spectra caused by rotation,the nonlinear vibration usually ha...Composite cylindrical shells,as key components,are widely employed in large rotating machines.However,due to the frequency bifurcations and dense frequency spectra caused by rotation,the nonlinear vibration usually has the behavior of complex multiple internal resonances.In addition,the varying temperature fields make the responses of the system further difficult to obtain.Therefore,the multiple internal resonances of composite cylindrical shells with porosities induced by rotation with varying temperature fields are studied in this paper.Three different types of the temperature fields,the Coriolis forces,and the centrifugal force are considered here.The Hamilton principle and the modified Donnell nonlinear shell theory are used to obtain the equilibrium equations of the system,which are transformed into the ordinary differential equations(ODEs)by the multi-mode Galerkin technique.Thereafter,the pseudo-arclength continuation method,which can identify the regions of instability,is introduced to obtain the numerical results.The detailed parametric analysis of the rotating composite shells is performed.Multiple internal resonances caused by the interaction between backward and forward wave modes and the energy transfer phenomenon are detected.Besides,the nonlinear amplitude-frequency response curves are different under different temperature fields.展开更多
文摘Mathematical modelling has been extensively used to measure intervention strategies for the control of contagious conditions.Alignment between different models is pivotal for furnishing strong substantiation for policymakers because the differences in model features can impact their prognostications.Mathematical modelling has been widely used in order to better understand the transmission,treatment,and prevention of infectious diseases.Herein,we study the dynamics of a human immunodeficiency virus(HIV)infection model with four variables:S(t),I(t),C(t),and A(t)the susceptible individuals;HIV infected individuals(with no clinical symptoms of AIDS);HIV infected individuals(under ART with a viral load remaining low),and HIV infected individuals with two different incidence functions(bilinear and saturated incidence functions).A novel numerical scheme called the continuous Galerkin-Petrov method is implemented for the solution of themodel.The influence of different clinical parameters on the dynamical behavior of S(t),I(t),C(t)and A(t)is described and analyzed.All the results are depicted graphically.On the other hand,we explore the time-dependent movement of nanofluid in porous media on an extending sheet under the influence of thermal radiation,heat flux,hall impact,variable heat source,and nanomaterial.The flow is considered to be 2D,boundary layer,viscous,incompressible,laminar,and unsteady.Sufficient transformations turn governing connected PDEs intoODEs,which are solved using the proposed scheme.To justify the envisaged problem,a comparison of the current work with previous literature is presented.
基金supported by the Natural Science Foundation of Shandong Province(ZR2021MA001)the Fundamental Research Funds for the Central Universities(20CX05011A)+1 种基金supported by National Natural Science Foundation of China Grant 11801569supported by NSF grant DMS-1818467 and Simons Foundation 961585.
文摘In this paper,we apply local discontinuous Galerkin(LDG)methods for pattern formation dynamical model in polymerizing actin focks.There are two main dificulties in designing effective numerical solvers.First of all,the density function is non-negative,and zero is an unstable equilibrium solution.Therefore,negative density values may yield blow-up solutions.To obtain positive numerical approximations,we apply the positivitypreserving(PP)techniques.Secondly,the model may contain stif source.The most commonly used time integration for the PP technique is the strong-stability-preserving Runge-Kutta method.However,for problems with stiff source,such time discretizations may require strictly limited time step sizes,leading to large computational cost.Moreover,the stiff source any trigger spurious filament polarization,leading to wrong numerical approximations on coarse meshes.In this paper,we combine the PP LDG methods with the semi-implicit Runge-Kutta methods.Numerical experiments demonstrate that the proposed method can yield accurate numerical approximations with relatively large time steps.
基金Project supported by the National Natural Science Foundation of China (No. 10472086).
文摘Anti-plane electroelastic problems are studied by the Trefftz boundary element method （BEM） in this paper. The Trefftz BEM is based on a weighted residual formulation and indirect boundary approach. In particular the point-collocation and Galerkin techniques, in which the basic unknowns are the retained expansion coefficients in the system of complete equations, are considered. Furthermore, special Trefftz functions and auxiliary functions which satisfy exactly the specified boundary conditions along the slit boundaries are also used to derive a special purpose element with local defects. The path-independent integral is evaluated at the tip of a crack to determine the energy release rate for a mode Ⅲ fracture problem. In final, the accuracy and efficiency of the Trefftz boundary element method are illustrated by an example and the comparison is made with other methods.
文摘This paper investigated the buoyancy and surface tension-driven ferro-thermal-convection (FTC) in a ferrofluid (FF) layer due to influence of general boundary conditions. The lower surface is rigid with insulating to temperature perturbations, while the upper surface is stress-free and subjected to general thermal boundary condition. The numerically Galerkin technique (GT) and analytically regular perturbation technique (RPT) are applied for solving the problem of eigenvalue. It is analyzed that increasing Biot number, decreases the magnetic and Marangoni number is to postponement the onset. Additionally, magnetization nonlinearity parameter has no effect on FTC in the non-existence of Biot number. The results under the limiting cases are found to be in good agreement with those available in the literature.
基金supported by Hunan Provincial Natural Science Foundation of China under Grant No. 10JJ3021Scientific Research Fund of Hunan Provincial Education Department under Grant No.11B032the Planned Science and Technology Project of Hunan Province and Aid program for Science and Technology Innovative Research Team in Higher Educational Institutions of Hunan Province
文摘In order to obtain an expected numerical solution, a fully discrete discontinuous Galerkin method is applied to a kind of reactive transport problems in two dimension. That is to say, the space variable is discretized with the symmetric interior penalty Calerkin method （SIPG）, and the time variable is done with the backward Euler method. Making use of the duality technique, hp approximation properties and the interpolation theory, a residual-type a posteriori error estimation is achieved, which can be used for adaptivity. Compared with the analyses of semi-discretization, the current presentation is more challenging and more significant.
基金supported by the National Natural Science Foundation of China(Grant No. 10872104)the Fundamental Research Funds for the Central Universities(Grant No. FRF-BR-10.007A)
文摘A discontinuous Galerkin (DG) finite element method is presented to solve the thermoelastic coupling problems caused by temperature and pressure dependent thermal contact resistance (TCR).The whole analysis is made up of two parts,thermal and mechanical analysis.In thermal analysis,the DG method is employed to simulate the temperature jump phenomenon,which satisfies the imperfect thermal contact condition in a straightforward manner.In mechanical analysis,the impenetrability condition is fulfilled through a DG approach with penalty functions.The Picard iteration procedure with a relaxation technique is also adopted to accelerate the rate of convergence and avoid numerical instability.Numerical examples show that the present method is an attractive approach for solving thermoelastic coupling problems caused by TCR.The methodology can also be expanded to solve problems with friction finite deformation contact,node-to-segment contact and node-to-surface contact,etc.in a straightforward manner.
文摘A finite element solution for the Navier-Stokes equations for steady flow under the porosity effects through a double branched two-dimensional section of a three-dimensional model of a canine aorta was obtained. The numerical solution involves transforming a physical coordinates to a curvilinear boundary fitted coordinate system.The steady flow,branch flow and shear stress under the porous effects were discussed in detail. The shear stress at the wall was calculated for Reynolds number of 1 000 with branch to main aortic flow rate ratio as a parameter. The results are compared with earlier works involving experimental data and it has been observed that our results are very close to the exact solutions.This work is in fact an improvement of the work of Sharma et al. (2001) in the sense that computational technique is economic and (Reynolds) number is large.
基金Project supported by the National Natural Science Foundation of China(Nos.11272209,11432009,and 11872241)
文摘In this paper, a brand-new wavelet-homotopy Galerkin technique is developed to solve nonlinear ordinary or partial differential equations. Before this investigation,few studies have been done for handling nonlinear problems with non-uniform boundary conditions by means of the wavelet Galerkin technique, especially in the field of fluid mechanics and heat transfer. The lid-driven cavity flow and heat transfer are illustrated as a typical example to verify the validity and correctness of this proposed technique. The cavity is subject to the upper and lower walls’ motions in the same or opposite directions.The inclined angle of the square cavity is from 0 to π/2. Four different modes including uniform, linear, exponential, and sinusoidal heating are considered on the top and bottom walls, respectively, while the left and right walls are thermally isolated and stationary.A parametric analysis of heating distribution between upper and lower walls including the amplitude ratio from 0 to 1 and the phase deviation from 0 to 2π is conducted. The governing equations are non-dimensionalized in terms of the stream function-vorticity formulation and the temperature distribution function and then solved analytically subject to various boundary conditions. Comparisons with previous publications are given,showing high efficiency and great feasibility of the proposed technique.
文摘The effects of electric field and non-uniform basic temperature gradient on the onset of Rayleigh-Bénard convection in a micropolar fluid are studied using the Galerkin technique. The eigenvalues are obtained for free-free, rigid-free and rigid-rigid velocity boundary combinations and for isothermal and/or adiabatic temperature boundaries. The microrotation is assumed to vanish at the boundaries. A linear stability analysis is performed. The influence of various micropolar fluid parameters and electric Rayleigh number on the onset of convection has been analyzed. One linear and five non-uniform temperature profiles are considered and their comparative influence on onset is discussed.
文摘The number and distribution of the singular points of streamlines in the cross-section of steady flow through a curved tube ate discussed by using the method of topological structure analysis. And a theoretical criterion is obtained for the bifurcation of flow vortexes for the secondary flow turning from two-vortex structure into four-vortex structure. Furthermore, the critical Dean number for bifurcation and the semi-analytical expressions of stream function and axial velocity are given by using Galerkin technique. The result of calculation is consistent with the theoretical criterion.
文摘The effects of electric field and non-uniform basic temperature gradient on the onset of Rayleigh-Bénard-Marangoni convection in a micropolar fluid are studied using the Galerkin technique. The eigenvalues are obtained for an upper free/adiabatic and lower rigid/isothermal boundaries. The microrotation is assumed to vanish at the boundaries. A linear stability analysis is performed. The influence of various micropolar fluid parameters and electric Rayleigh number on the onset of convection has been analysed. Six different non-uniform temperature profiles are considered and their comparative influence on onset is discussed.
文摘The effect of viscosity depending exponentially on temperature on the onset of penetrative ferro-thermal-convection (FTC) in a saturated horizontal porous layer in the presence of vertical magnetic field is investigated. The bounding surface of the ferrofluid layer is considered to be rigid-rigid and insulated to temperature perturbations. The resulting eigenvalue problem is solved numerically using the Galerkin technique and also analytically by a regular perturbation technique with wave number as a perturbation parameter. The analytical and numerical results are found to be concurrence. The characteristics of stability of the system are strongly dependent on the viscosity parameter B. The effect of B on the onset of ferroconvection in a porous layer is dual in nature depending on the choices of physical parameters and a sublayer starts to form at higher values of B. Whereas, increase in magnetic number M1 and the Darcy number Da is to advance the onset of ferroconvection in a porous layer. The nonlinearity of fluid magnetization M3?is found to have no influence on the onset of ferroconvection.
基金supported by the National Natural Science Foundation of China(No.11972204)。
文摘Composite cylindrical shells,as key components,are widely employed in large rotating machines.However,due to the frequency bifurcations and dense frequency spectra caused by rotation,the nonlinear vibration usually has the behavior of complex multiple internal resonances.In addition,the varying temperature fields make the responses of the system further difficult to obtain.Therefore,the multiple internal resonances of composite cylindrical shells with porosities induced by rotation with varying temperature fields are studied in this paper.Three different types of the temperature fields,the Coriolis forces,and the centrifugal force are considered here.The Hamilton principle and the modified Donnell nonlinear shell theory are used to obtain the equilibrium equations of the system,which are transformed into the ordinary differential equations(ODEs)by the multi-mode Galerkin technique.Thereafter,the pseudo-arclength continuation method,which can identify the regions of instability,is introduced to obtain the numerical results.The detailed parametric analysis of the rotating composite shells is performed.Multiple internal resonances caused by the interaction between backward and forward wave modes and the energy transfer phenomenon are detected.Besides,the nonlinear amplitude-frequency response curves are different under different temperature fields.
