A memory reduction technique is proposed for solving stationary kinetic model equations. As implied by an integral solution of the stationary kinetic equation, a velocity distribution function can be reconstructed fro...A memory reduction technique is proposed for solving stationary kinetic model equations. As implied by an integral solution of the stationary kinetic equation, a velocity distribution function can be reconstructed from given macroscopic variables. Based on this fact, we propose a technique to reconstruct distribution function at discrete level, and employ it to develop an implicit numerical method for kinetic equations. The new implicit method only stores the macroscopic quantities which appear in the collision term, and does not store the distribution functions. As a result, enormous memory requirement for solving kinetic equations is totally relieved. Several boundary conditions, such as, inlet, outlet and isothermal boundaries, are discussed. Some numerical tests demonstrate the validity and efficiency of the technique.The new implicit solver provides nearly identical solution as the explicit kinetic solver, while the memory requirement is on the same order as the Navier–Stokes solver.展开更多
This paper presents the effect of mooring diameters, fairlead slopes and pretensions on the dynamic responses of a truss spar platform in intact and damaged line conditions. The platform is modelled as a rigid body wi...This paper presents the effect of mooring diameters, fairlead slopes and pretensions on the dynamic responses of a truss spar platform in intact and damaged line conditions. The platform is modelled as a rigid body with three degrees-of-freedom and its motions are analysed in time-domain using the implicit Newmark Beta technique. The mooring restoring force-excursion relationship is evaluated using quasi-static approach. MATLAB codes DATSpar and QSAML, are developed to compute the dynamic responses of truss spar platform and to determine the mooring system stiffness. To eliminate the conventional trial and error approach in the mooring system design, a numerical tool is also developed and described in this paper for optimising the mooring configuration. It has a graphical user interface and includes regrouping particle swarm optimisation technique combined with DATSpar and QSAML. A case study of truss spar platform with ten mooring lines is analysed using this numerical tool. The results show that optimum mooring system design benefits the oil and gas industry to economise the project cost in terms of material, weight, structural load onto the platform as well as manpower requirements. This tool is useful especially for the preliminary design of truss spar platforms and its mooring system.展开更多
Corrected explicit-implicit domain decomposition(CEIDD) algorithms are studied for parallel approximation of semilinear parabolic problems on distributed memory processors. It is natural to divide the spatial domain i...Corrected explicit-implicit domain decomposition(CEIDD) algorithms are studied for parallel approximation of semilinear parabolic problems on distributed memory processors. It is natural to divide the spatial domain into some smaller parallel strips and cells using the simplest straightline interface(SI) . By using the Leray-Schauder fixed-point theorem and the discrete energy method,it is shown that the resulting CEIDD-SI algorithm is uniquely solvable,unconditionally stable and convergent. The CEIDD-SI method always suffers from the globalization of data communication when interior boundaries cross into each other inside the domain. To overcome this disadvantage,a composite interface(CI) that consists of straight segments and zigzag fractions is suggested. The corresponding CEIDD-CI algorithm is proven to be solvable,stable and convergent. Numerical experiments are presented to support the theoretical results.展开更多
A group of asymmetric difference schemes to approach the Korteweg-de Vries （KdV） equation is given here. According to such schemes, the full explicit difference scheme and the full implicit one, an alternating segme...A group of asymmetric difference schemes to approach the Korteweg-de Vries （KdV） equation is given here. According to such schemes, the full explicit difference scheme and the full implicit one, an alternating segment explicit-implicit difference scheme for solving the KdV equation is constructed. The scheme is linear unconditionally stable by the analysis of linearization procedure, and is used directly on the parallel computer. The numerical experiments show that the method has high accuracy.展开更多
In this work, the effects of boundaries conditions and truncation errors in the distribution of minority carriers in the semiconductor are studied. It is a one-dimensional digital study of a polycrystalline silicon so...In this work, the effects of boundaries conditions and truncation errors in the distribution of minority carriers in the semiconductor are studied. It is a one-dimensional digital study of a polycrystalline silicon solar cell under polychromatic illumination in a dynamic state. Starting from the Boltzmann equation of semiconductors, the author establishes the general equation of particle transport. Assumptions made on the latter allow it to give the equation of distribution of minority carriers in a general way in its case to be studied. This dimensioned distribution equation reveals the parameters of influences on the distribution of carriers. It obtains a partial derivative equation for the carrier distribution function. The boundary conditions are then discretized to order one and then to order two. By considering boundary conditions and the nature of the carriers, the author numerically resolves the discretized general equation by assessing the influence of the nature of the boundary conditions and truncation errors and the influence of the discretization step on the density of the charge carriers by setting certain parameters and varying others. The work ends with a conclusion and logical follow-up to this work.展开更多
In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two ...In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two nonlinear schemes and two linearized schemes are presented. All resulting schemes will be analyzed for accuracy and stability. The exact solution and the conserved quantities are used to highlight the efficiency and the robustness of the proposed schemes. Interaction of two and three solitons will be also conducted. The numerical results show that the interaction behavior is elastic and the conserved quantities are conserved exactly, and this is a good indication of the reliability of the schemes which we derived. A comparison with some existing is presented as well.展开更多
A parabolic equation （PE） based method for analyzing composite scattering under an electromagnetic wave incidence at low grazing angle, which composes of three-dimensional （3-D） electrically large targets and roug...A parabolic equation （PE） based method for analyzing composite scattering under an electromagnetic wave incidence at low grazing angle, which composes of three-dimensional （3-D） electrically large targets and rough surface, is presented and discussed. A superior high-order PE version is used to improve the accuracy at wider paraxial angles, and along with the alternating direction implicit （ADI） differential technique, the computational efficiency is further improved. The formula of bistatic normalized radar cross section is derived by definition and near-far field transformation. Numerical examples are given to show the validity and accuracy of the proposed approach, in which the results are compared with those of Kirchhoff approximation （KA） and moment of method （MoM）. Furthermore, the bistatic scattering properties of composite model in which the 3-D PEC targets on or above the two-dimensional Gaussian rough surfaces under the tapered wave incidence are analyzed.展开更多
This paper introduces a three dielectric layer hybrid solvation model for treating electrostatic interactions of biomolecules in solvents using the PoissonBoltzmann equation.In this model,an interior spherical cavity ...This paper introduces a three dielectric layer hybrid solvation model for treating electrostatic interactions of biomolecules in solvents using the PoissonBoltzmann equation.In this model,an interior spherical cavity will contain the solute and some explicit solvent molecules,and an intermediate buffer layer and an exterior layer contain the bulk solvent.A special dielectric permittivity profile is used to achieve a continuous dielectric transition from the interior cavity to the exterior layer.The selection of this special profile using a harmonic interpolation allows an analytical solution of the model by generalizing the classical Kirkwood series expansion.Discrete image charges are used to speed up calculations for the electrostatic potential within the interior and buffer layer regions.Semi-analytical and least squares methods are used to construct an accurate discrete image approximation for the reaction field due to solvent with or without salt effects.In particular,the image charges obtained by the least squares method provide accurate approximations to the reaction field independent of the ionic concentration of the solvent.Numerical results are presented to validate the accuracy and effectiveness of the image charge methods.展开更多
It is a very common practice to use semi-implicit schemes in various computations,which treat selected linear terms implicitly and the nonlinear terms explicitly.For phase-field equations,the principal elliptic operat...It is a very common practice to use semi-implicit schemes in various computations,which treat selected linear terms implicitly and the nonlinear terms explicitly.For phase-field equations,the principal elliptic operator is treated implicitly to reduce the associated stability constraints while the nonlinear terms are still treated explicitly to avoid the expensive process of solving nonlinear equations at each time step.However,very few recent numerical analysis is relevant to semi-implicit schemes,while”stabilized”schemes have become very popular.In this work,we will consider semiimplicit schemes for the Allen-Cahn equation with general potential function.It will be demonstrated that the maximum principle is valid and the energy stability also holds for the numerical solutions.This paper extends the result of Tang&Yang(J.Comput.Math.,34(5)(2016),pp.471–481),which studies the semi-implicit scheme for the Allen-Cahn equation with polynomial potentials.展开更多
In this paper, an implicit symmetry constraint is calculated and its associated binary nonlinearization of the Lax pairs and the adjoint Lax pairs is carried out for the modified Korteweg-de Vries (mKdV) equation. Aft...In this paper, an implicit symmetry constraint is calculated and its associated binary nonlinearization of the Lax pairs and the adjoint Lax pairs is carried out for the modified Korteweg-de Vries (mKdV) equation. After introducing two new inde-pendent variables, we find that under the implicit symmetry constraint, the spatial part and the temporal part of the mKdV equation are decomposed into two finite-dimensional systems. Furthermore we prove that the obtained finite-dimensional systems are Hamiltonian systems and completely integrable in the Liouville sense.展开更多
In this paper we have made a numerical study on the control of vortex shedding and drag reduction of a cylinder by attaching thin splitter plates. The wake structure of the cylinder of square cross-section with attach...In this paper we have made a numerical study on the control of vortex shedding and drag reduction of a cylinder by attaching thin splitter plates. The wake structure of the cylinder of square cross-section with attached splitter plates is analyzed for a range of Reynolds number, based on the incident stream and height of the cylinder, in the laminar range. The Navier-Stokes equations governing the flow are solved by the control volume method over a staggered grid arrangement. We have used the semi-implicit method for pressure-linked equation （SIMPLE） algorithm for computation. Our results show that the presence of a splitter plate upstream of the cylinder reduces the drag, but it has a small impact on the vortex shedding frequency when the plate length is beyond 1.5 time the height of the cylinder. The presence of a downstream splitter plate dampens the vortex shedding frequency. The entrainment of fluid into the inner side of the separated shear layers is obstructed by the downstream splitter plate. Our results suggest that by attaching in-line splitter plates both upstream and downstream of the cylinder, the vortex shedding can be suppressed, as well as a reduction in drag be obtained. We made a parametric study to determine the optimal length of these splitter plates so as to achieve low drag and low vortex shedding frequency.展开更多
基金supported by the National Natural Science Foundation of China(11602091 and 91530319)the National Key Research and Development Plan(2016YFB0600805)
文摘A memory reduction technique is proposed for solving stationary kinetic model equations. As implied by an integral solution of the stationary kinetic equation, a velocity distribution function can be reconstructed from given macroscopic variables. Based on this fact, we propose a technique to reconstruct distribution function at discrete level, and employ it to develop an implicit numerical method for kinetic equations. The new implicit method only stores the macroscopic quantities which appear in the collision term, and does not store the distribution functions. As a result, enormous memory requirement for solving kinetic equations is totally relieved. Several boundary conditions, such as, inlet, outlet and isothermal boundaries, are discussed. Some numerical tests demonstrate the validity and efficiency of the technique.The new implicit solver provides nearly identical solution as the explicit kinetic solver, while the memory requirement is on the same order as the Navier–Stokes solver.
基金partially supported by YUTP-FRG funded by PETRONAS
文摘This paper presents the effect of mooring diameters, fairlead slopes and pretensions on the dynamic responses of a truss spar platform in intact and damaged line conditions. The platform is modelled as a rigid body with three degrees-of-freedom and its motions are analysed in time-domain using the implicit Newmark Beta technique. The mooring restoring force-excursion relationship is evaluated using quasi-static approach. MATLAB codes DATSpar and QSAML, are developed to compute the dynamic responses of truss spar platform and to determine the mooring system stiffness. To eliminate the conventional trial and error approach in the mooring system design, a numerical tool is also developed and described in this paper for optimising the mooring configuration. It has a graphical user interface and includes regrouping particle swarm optimisation technique combined with DATSpar and QSAML. A case study of truss spar platform with ten mooring lines is analysed using this numerical tool. The results show that optimum mooring system design benefits the oil and gas industry to economise the project cost in terms of material, weight, structural load onto the platform as well as manpower requirements. This tool is useful especially for the preliminary design of truss spar platforms and its mooring system.
基金supported by National Natural Science Foundation of China (Grant No. 10871044)
文摘Corrected explicit-implicit domain decomposition(CEIDD) algorithms are studied for parallel approximation of semilinear parabolic problems on distributed memory processors. It is natural to divide the spatial domain into some smaller parallel strips and cells using the simplest straightline interface(SI) . By using the Leray-Schauder fixed-point theorem and the discrete energy method,it is shown that the resulting CEIDD-SI algorithm is uniquely solvable,unconditionally stable and convergent. The CEIDD-SI method always suffers from the globalization of data communication when interior boundaries cross into each other inside the domain. To overcome this disadvantage,a composite interface(CI) that consists of straight segments and zigzag fractions is suggested. The corresponding CEIDD-CI algorithm is proven to be solvable,stable and convergent. Numerical experiments are presented to support the theoretical results.
基金Project supported by the National Natural Science Foundation of China(No.10671113)the Natural Science Foundation of Shandong Province of China(No.Y2003A04)
文摘A group of asymmetric difference schemes to approach the Korteweg-de Vries （KdV） equation is given here. According to such schemes, the full explicit difference scheme and the full implicit one, an alternating segment explicit-implicit difference scheme for solving the KdV equation is constructed. The scheme is linear unconditionally stable by the analysis of linearization procedure, and is used directly on the parallel computer. The numerical experiments show that the method has high accuracy.
文摘In this work, the effects of boundaries conditions and truncation errors in the distribution of minority carriers in the semiconductor are studied. It is a one-dimensional digital study of a polycrystalline silicon solar cell under polychromatic illumination in a dynamic state. Starting from the Boltzmann equation of semiconductors, the author establishes the general equation of particle transport. Assumptions made on the latter allow it to give the equation of distribution of minority carriers in a general way in its case to be studied. This dimensioned distribution equation reveals the parameters of influences on the distribution of carriers. It obtains a partial derivative equation for the carrier distribution function. The boundary conditions are then discretized to order one and then to order two. By considering boundary conditions and the nature of the carriers, the author numerically resolves the discretized general equation by assessing the influence of the nature of the boundary conditions and truncation errors and the influence of the discretization step on the density of the charge carriers by setting certain parameters and varying others. The work ends with a conclusion and logical follow-up to this work.
文摘In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two nonlinear schemes and two linearized schemes are presented. All resulting schemes will be analyzed for accuracy and stability. The exact solution and the conserved quantities are used to highlight the efficiency and the robustness of the proposed schemes. Interaction of two and three solitons will be also conducted. The numerical results show that the interaction behavior is elastic and the conserved quantities are conserved exactly, and this is a good indication of the reliability of the schemes which we derived. A comparison with some existing is presented as well.
基金Project supported by the National Natural Science Foundation of China(Grant No.61771407)
文摘A parabolic equation （PE） based method for analyzing composite scattering under an electromagnetic wave incidence at low grazing angle, which composes of three-dimensional （3-D） electrically large targets and rough surface, is presented and discussed. A superior high-order PE version is used to improve the accuracy at wider paraxial angles, and along with the alternating direction implicit （ADI） differential technique, the computational efficiency is further improved. The formula of bistatic normalized radar cross section is derived by definition and near-far field transformation. Numerical examples are given to show the validity and accuracy of the proposed approach, in which the results are compared with those of Kirchhoff approximation （KA） and moment of method （MoM）. Furthermore, the bistatic scattering properties of composite model in which the 3-D PEC targets on or above the two-dimensional Gaussian rough surfaces under the tapered wave incidence are analyzed.
基金The authors would like to thank the financial support provided by the National Institutes of Health(Grant No.1R01GM083600-02)Z.Xu is also partially supported by the Charlotte Research Institute through a Duke Postdoctoral FellowshipW.Cai and Z.Xu are also partially supported by the Department of Energy(Grant No.DEFG0205ER25678).
文摘This paper introduces a three dielectric layer hybrid solvation model for treating electrostatic interactions of biomolecules in solvents using the PoissonBoltzmann equation.In this model,an interior spherical cavity will contain the solute and some explicit solvent molecules,and an intermediate buffer layer and an exterior layer contain the bulk solvent.A special dielectric permittivity profile is used to achieve a continuous dielectric transition from the interior cavity to the exterior layer.The selection of this special profile using a harmonic interpolation allows an analytical solution of the model by generalizing the classical Kirkwood series expansion.Discrete image charges are used to speed up calculations for the electrostatic potential within the interior and buffer layer regions.Semi-analytical and least squares methods are used to construct an accurate discrete image approximation for the reaction field due to solvent with or without salt effects.In particular,the image charges obtained by the least squares method provide accurate approximations to the reaction field independent of the ionic concentration of the solvent.Numerical results are presented to validate the accuracy and effectiveness of the image charge methods.
文摘It is a very common practice to use semi-implicit schemes in various computations,which treat selected linear terms implicitly and the nonlinear terms explicitly.For phase-field equations,the principal elliptic operator is treated implicitly to reduce the associated stability constraints while the nonlinear terms are still treated explicitly to avoid the expensive process of solving nonlinear equations at each time step.However,very few recent numerical analysis is relevant to semi-implicit schemes,while”stabilized”schemes have become very popular.In this work,we will consider semiimplicit schemes for the Allen-Cahn equation with general potential function.It will be demonstrated that the maximum principle is valid and the energy stability also holds for the numerical solutions.This paper extends the result of Tang&Yang(J.Comput.Math.,34(5)(2016),pp.471–481),which studies the semi-implicit scheme for the Allen-Cahn equation with polynomial potentials.
文摘In this paper, an implicit symmetry constraint is calculated and its associated binary nonlinearization of the Lax pairs and the adjoint Lax pairs is carried out for the modified Korteweg-de Vries (mKdV) equation. After introducing two new inde-pendent variables, we find that under the implicit symmetry constraint, the spatial part and the temporal part of the mKdV equation are decomposed into two finite-dimensional systems. Furthermore we prove that the obtained finite-dimensional systems are Hamiltonian systems and completely integrable in the Liouville sense.
基金国家重点基础研究发展规划(973)(the National Grand Fundamental Research 973 Program of China under Grant No.2006CB504801)国家自然科学基金(the National Natural Science Foundation of China under Grant No.60521002)
文摘In this paper we have made a numerical study on the control of vortex shedding and drag reduction of a cylinder by attaching thin splitter plates. The wake structure of the cylinder of square cross-section with attached splitter plates is analyzed for a range of Reynolds number, based on the incident stream and height of the cylinder, in the laminar range. The Navier-Stokes equations governing the flow are solved by the control volume method over a staggered grid arrangement. We have used the semi-implicit method for pressure-linked equation （SIMPLE） algorithm for computation. Our results show that the presence of a splitter plate upstream of the cylinder reduces the drag, but it has a small impact on the vortex shedding frequency when the plate length is beyond 1.5 time the height of the cylinder. The presence of a downstream splitter plate dampens the vortex shedding frequency. The entrainment of fluid into the inner side of the separated shear layers is obstructed by the downstream splitter plate. Our results suggest that by attaching in-line splitter plates both upstream and downstream of the cylinder, the vortex shedding can be suppressed, as well as a reduction in drag be obtained. We made a parametric study to determine the optimal length of these splitter plates so as to achieve low drag and low vortex shedding frequency.
