To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direc...To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direct current method, we propose a new mesh refinement and recoarsement method for a two-dimensional point source. We introduce the mesh refinement and mesh recoarsement into the traditional structured mesh subdivision. By refining the horizontal grids, the singularity owing to the point source is minimized and the topography is simulated. By recoarsening the horizontal grids, the number of grid cells is reduced significantly and computational efficiency is improved. Model tests show that the proposed method solves the singularity problem and reduces the number of grid cells by 80% compared to the uniform grid refinement.展开更多
In this paper we develop several new refinement relations of Z for multiple viewpoints oriented requirements method (MVORM). The original motivation is that we found the standard Z refinement relation is not adequate ...In this paper we develop several new refinement relations of Z for multiple viewpoints oriented requirements method (MVORM). The original motivation is that we found the standard Z refinement relation is not adequate or correct when considering specifications that have temporal relationships of operations. The concept of temporal state variables is introduced into Z. Then new implementation relations are defined and new refinement relations are deduced, mainly for temporal state variables to process temporal relationships of operations. We use state transition systems to abstract the temporal state transitions. A simple example is used to show the procedures of MVORM. Finally some directions of further work are forwarded.展开更多
V-Ti-Fe master alloys were prepared by metaUothermic reduction method, and the influences of the mass ratio of V205 to TiO2, Al and Al-Mg alloy addition amounts on the metal recovery rates and alloy compositions were ...V-Ti-Fe master alloys were prepared by metaUothermic reduction method, and the influences of the mass ratio of V205 to TiO2, Al and Al-Mg alloy addition amounts on the metal recovery rates and alloy compositions were investigated. The results show that appropriate technological parameters are： the mass ratio of V205 to TiO2 is 0.5：1, Al addition represents 95% of the theoretical value, and the Al-Mg alloy addition amount is one third that of the Al addition. The results from energy spectrum analysis show that V and Fe distribute uniformly in the prepared alloy, while the segregation for Ti, i.e. Ti-rich phase is detected. A spray refming process was carried out to reduce the impurity contents of Al and O in the prepared alloys. The Al content drops from 4.27% to 1.86%, and the O content drops from 2.10% to 0.91% after the refining process.展开更多
A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the ...A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the Reynolds number based on the hull length is 1.0x 105, An immersed boundary method based on the moving-least-squares reconstruction is used to handle the complex geometric boundaries. The adaptive mesh refinement is utilized to resolve the flows near the hull, The parallel scalabilities of the flow solver are tested on meshes with the number of cells varying from 50 million to 3.2 billion, The parallel solver reaches nearly linear scalability for the flows around the underwater vehicle model, The present simulation captures the essential features of the vortex structures near the hull and in the wake, Both of the time-averaged pressure coefficients and srreamwise velocity profiles obtained from the LES are consistent with the characteristics of the flows pass an appended axisymmetric body. The code efficiency and its correct predictions on flow features allow us to perform the full-scale simulations on tens of thousands of cores with billions of grid points for higher-Reynolds-number flows around the underwater vehicles.展开更多
China experienced a decline of water use intensity in the 11th Five Year Plan,but the water use intensity actually increased in 2009.To the best of our knowledge,the index decomposition analysis method was rarely used...China experienced a decline of water use intensity in the 11th Five Year Plan,but the water use intensity actually increased in 2009.To the best of our knowledge,the index decomposition analysis method was rarely used to analyze changes in water use,and no decomposition analysis has investigated the role of regional economy in the decline of water use intensity.In this paper,we use logarithmic mean Divisia index(LMDI)techniques to decompose the change of water use intensity in the period 2006-2010.We find that the change of industrial water use intensity is confirmed as the dominant contributor to the decline in the overall water use intensity;the regional structure effect and the industrial structure effect is positive to the decline of overall water use intensity;the decline of China's water use intensity is mainly attributed to the effect of developed eastern provinces;meanwhile,the effect of central and undeveloped western is also positive to the decline of overall water use intensity;at least one out of three effects is positive to the decline of water use intensity in the different provinces;the intensity effect is positive and the industrial structure effect is positive to the declines of China's water use intensity based on chaining approach except the period 2008-2009,individually;and the deviation of regional structure effect and industrial structure effect between with regional economy and without regional economy in LMDI is 0.9 and2.3 m^3/10~4 RMB,respectively.展开更多
We present a high performance modularly-built open-source software-OpenIFEM.OpenIFEM is a C++implementation of the modified immersed finite element method(mIFEM)to solve fluid-structure interaction(FSI)problems.This s...We present a high performance modularly-built open-source software-OpenIFEM.OpenIFEM is a C++implementation of the modified immersed finite element method(mIFEM)to solve fluid-structure interaction(FSI)problems.This software is modularly built to perform multiple tasks including fluid dynamics(incompressible and slightly compressible fluid models),linear and nonlinear solid mechanics,and fully coupled fluid-structure interactions.Most of open-source software packages are restricted to certain discretization methods;some are under-tested,under-documented,and lack modularity as well as extensibility.OpenIFEM is designed and built to include a set of generic classes for users to adapt so that any fluid and solid solvers can be coupled through the FSI algorithm.In addition,the package utilizes well-developed and tested libraries.It also comes with standard test cases that serve as software and algorithm validation.The software can be built on cross-platform,i.e.,Linux,Windows,and Mac OS,using CMake.Efficient parallelization is also implemented for high-performance computing for large-sized problems.OpenIFEM is documented using Doxygen and publicly available to download on GitHub.It is expected to benefit the future development of FSI algorithms and be applied to a variety of FSI applications.展开更多
This paper presents a refined method for estimating the annual extreme wave heights at a coastal or offshore project site on the basis of the data acquired at some nearby routine hydrographic stations. This method is ...This paper presents a refined method for estimating the annual extreme wave heights at a coastal or offshore project site on the basis of the data acquired at some nearby routine hydrographic stations. This method is based on the orthogonality principle in linear mean square estimation of stochastic processes. The error of the method is analyzed and compared with that of the conventional method. It is found that the method is able to effectively reduce the error so long as some feasible measures are adopted. A simulated test of the method has been conducted in a large scale wind wave flume. The test results are in good agreement with those given by theoretical error analysis. A scheme to implement the method is proposed on the basis of error analysis. The scheme is so designed as to reduce the estimation error as far as possible. This method is also suitable to utilizing satellite wave data for the estimation.展开更多
The adaptive mesh refinement （AMR） method is applied in the 2-D Euler multi-component elasticplastic hydrodynamics code （MEPH2Y）. It is applied on detonation. Firstly, the AMR method is described, including a cons...The adaptive mesh refinement （AMR） method is applied in the 2-D Euler multi-component elasticplastic hydrodynamics code （MEPH2Y）. It is applied on detonation. Firstly, the AMR method is described, including a conservative spatial interpolation, the time integration methodology with the adapitve time increment and an adaptive computational region method. The advantage of AMR technique is exhibited by numerical examples, including the 1-D C-J detonation and the 2-D implosion ignited from a single point. Results show that AMR can promote the computational efficiency, keeping the accuracy in interesting regions.展开更多
Britholites with the general formula Ca10-xLnx（PO4）6-x（SiO4）xF2, （0≤x≤6） are considered to be promising matrices for the confinement of the by-products in the nuclear industry. A thermodynamic study showed tha...Britholites with the general formula Ca10-xLnx（PO4）6-x（SiO4）xF2, （0≤x≤6） are considered to be promising matrices for the confinement of the by-products in the nuclear industry. A thermodynamic study showed that the stability of these compounds de-creased as the substitution rate increased. The present work was an attempt to gain more information about the structural changes in-duced by the substitution, in order to understand the observed stability decrease. The samples were successfully synthesized as a sin-gle-phase apatite by a solid-state reaction between 1200 and 1400 °C. The structural refinement indicated that the La3＋ions preferen-tially occupied the 6（h） sites. A progressive shift of F-along the c-axis outside its ideal position occurred as a result of the substitution increase. This might be the cause of the observed stability decrease, especially as the energies of the La-O, La-F and Si-O bonds are higher than those of Ca-O, Ca-F and P-O. The distribution of La3＋ between the two non-equivalent sites was confirmed by the charge distribution method.展开更多
The Carboniferous volcanic reservoir in the Chepaizi Uplift became an exploration hot target in recent years for its substantial amount of oils discovered. However, most of the Carboniferous heavy oils were biodegrade...The Carboniferous volcanic reservoir in the Chepaizi Uplift became an exploration hot target in recent years for its substantial amount of oils discovered. However, most of the Carboniferous heavy oils were biodegraded to PM7 or higher with orders of magnitude variation in oil viscosities. Two oil groups (I and II) exactly corresponding to the western and eastern Chepaizi Uplift were distinguished according to their source diagnose. Furthermore, three oil families (II1, II2 and II3), with the biodegradation level of PM7, PM8–8+, PM9+, respectively, were classified based on molecular compositions and parameter-stripping method of strongly bioresistant parameters. Allowing for this extremely high biodegradation case, more biodegradation refractory compound class were added to establish a refined Manco scale to quantitatively evaluate the biodegradation extent. Refined Manco number (RMN2) positively correlated with the oil density, NSO contents, and absolute concentrations of diasteranes and gammacerane, negatively correlated with the absolute concentrations of diahopane, summed tricyclic terpanes and pentacyclic terpanes. This refined scale showed higher resolution than the PM one to differentiate the biodegradation extent of Carboniferous heavy oils from the Chepaizi Uplift, especially those with same PM values but different oil viscosities.展开更多
Hybrid molecule/cluster statistical thermodynamics （HMCST） method is an efficient tool to simulate nano-scale systems under quasi-static loading at finite temperature. In this paper, a self-adaptive algorithm is dev...Hybrid molecule/cluster statistical thermodynamics （HMCST） method is an efficient tool to simulate nano-scale systems under quasi-static loading at finite temperature. In this paper, a self-adaptive algorithm is developed for this method. Explicit refinement criterion based on the gradient of slip shear deformation and a switching criterion based on generalized Einstein approximation is proposed respectively. Results show that this self-adaptive method can accurately find clusters to be refined or transferred to molecules, and efficiently refine or trans- fer the clusters. Furthermore, compared with fully atomistic simulation, the high computational efficiency of the self-adaptive method appears very attractive.展开更多
An essentially conservative adaptive space time conservation element and solution element （CE/SE） method is pro- posed for the effective simulation of shock-induced instability with low computational cost. Its imple...An essentially conservative adaptive space time conservation element and solution element （CE/SE） method is pro- posed for the effective simulation of shock-induced instability with low computational cost. Its implementation is based on redefined conservation elements （CEs） and solution elements （SEs）, optimized interpolations and a Courant number insensitive CE/SE scheme. This approach is used in two applications, the Woodward double Mach reflection and a two- component Richtmyer-Meshkov instability experiment. This scheme reveals the essential features of the investigated cases, captures small unstable structures, and yields a solution that is consistent with the results from experiments or other high order methods.展开更多
The present study provides a three-dimensional volume-of-fluid method based on the adaptive mesh refinement technique.The projection method on the adaptive mesh is introduced for solving the incompressible Navier-Stok...The present study provides a three-dimensional volume-of-fluid method based on the adaptive mesh refinement technique.The projection method on the adaptive mesh is introduced for solving the incompressible Navier-Stokes equations.The octree structure mesh is employed to solve the flow velocities and the pressure.The developed solver is applied to simulate the deformation of the cubic droplet driven by the surface tension without the effect of the gravity.The numerical results well predict the shape evolution of the droplet.展开更多
We propose a high-order enriched partition of unity finite element method for linear and nonlinear time-dependent diffusion problems.The solution of this class of problems often exhibits non-smooth features such as st...We propose a high-order enriched partition of unity finite element method for linear and nonlinear time-dependent diffusion problems.The solution of this class of problems often exhibits non-smooth features such as steep gradients and boundary layers which can be very challenging to recover using the conventional low-order finite element methods.A class of steady-state exponential functions has been widely used for enrichment and its performance to numerically solve these challenges has been demonstrated.However,these enrichment functions have been used only in context of the standard h-version refinement or the so-called q-version refinement.In this paper we demonstrate that the p-version refinement can also be a very attractive option in terms of the efficiency and the accuracy in the enriched partition of unity finite element method.First,the transient diffusion problem is integrated in time using a semi-implicit scheme and the semi-discrete problem is then integrated in space using the p-version enriched finite elements.Numerical results are presented for three test examples of timedependent diffusion problems in both homogeneous and heterogeneous media.The computed results show the significant improvement when using the p-version refined enriched approximations in the finite element analysis.In addition,these results support our expectations for a robust and high-order accurate enriched partition of unity finite element method.展开更多
The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh disto...The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh distortions. In a previous work, we constructed an 8-node quadrilateral spline element L8 using the triangular area coordinates and the B- net method, which can be insensitive to mesh distortions and possess the second order completeness in the Cartesian co- ordinates. In this paper, a thin plate spline element is devel- oped based on the spline element L8 and the refined tech- nique. Numerical examples show that the present element indeed possesses higher accuracy than the DKQ element for distorted meshes.展开更多
Straight-and curved-bar refining plates are two important types of plates commonly used in disc refiners in the papermaking industry.Theoretically,the curved-bar refining plate has a relatively uniform bar interaction...Straight-and curved-bar refining plates are two important types of plates commonly used in disc refiners in the papermaking industry.Theoretically,the curved-bar refining plate has a relatively uniform bar interaction angle,which indicates uniform refining effects.The bar angle of the curved bar was proposed and two typical curved-bar plates,the three-stage radial curved-bar plate and isometric curved-bar plate,were designed in this paper.The arc equations of the curved-bar center line and curved-bar edges were established and finally,the specific edge load(SEL)of the curved-bar plate was derived.The determination of bar parameters was discussed,which provides a theoretical basis for the design of curved-bar plates.展开更多
基金financially supported by the National Natural Science Foundation of China(No.41574127 and 41174104)the National Key Technology R&D Program for the 13th five-year plan(No.2016ZX05018006-006)
文摘To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direct current method, we propose a new mesh refinement and recoarsement method for a two-dimensional point source. We introduce the mesh refinement and mesh recoarsement into the traditional structured mesh subdivision. By refining the horizontal grids, the singularity owing to the point source is minimized and the topography is simulated. By recoarsening the horizontal grids, the number of grid cells is reduced significantly and computational efficiency is improved. Model tests show that the proposed method solves the singularity problem and reduces the number of grid cells by 80% compared to the uniform grid refinement.
基金Supported by Natural Science Foundation of Hubei Province (98J0 75 ) Ziqiang Technical Innovation Foundation ofWuhan Universi
文摘In this paper we develop several new refinement relations of Z for multiple viewpoints oriented requirements method (MVORM). The original motivation is that we found the standard Z refinement relation is not adequate or correct when considering specifications that have temporal relationships of operations. The concept of temporal state variables is introduced into Z. Then new implementation relations are defined and new refinement relations are deduced, mainly for temporal state variables to process temporal relationships of operations. We use state transition systems to abstract the temporal state transitions. A simple example is used to show the procedures of MVORM. Finally some directions of further work are forwarded.
基金Project (2006AA068128) supported by the High-tech Research and Development Program of China
文摘V-Ti-Fe master alloys were prepared by metaUothermic reduction method, and the influences of the mass ratio of V205 to TiO2, Al and Al-Mg alloy addition amounts on the metal recovery rates and alloy compositions were investigated. The results show that appropriate technological parameters are： the mass ratio of V205 to TiO2 is 0.5：1, Al addition represents 95% of the theoretical value, and the Al-Mg alloy addition amount is one third that of the Al addition. The results from energy spectrum analysis show that V and Fe distribute uniformly in the prepared alloy, while the segregation for Ti, i.e. Ti-rich phase is detected. A spray refming process was carried out to reduce the impurity contents of Al and O in the prepared alloys. The Al content drops from 4.27% to 1.86%, and the O content drops from 2.10% to 0.91% after the refining process.
基金supported by the National Natural Science Foundation of China (11302238, 11232011. and 11572331)support from the Strategic Priority Research Program (XDB22040104)+1 种基金the Key Research Program of Frontier Sciences of the Chinese Academy of Sciences (QYZDJ-SSW-SYS002)the National Basic Research Program of China (973 Program 2013CB834100: Nonlinear science)
文摘A large eddy simulation （LES） of the flows around an underwater vehicle model at intermediate Reynolds numbers is performed. The underwater vehicle model is taken as the DARPA SUBOFF with full appendages, where the Reynolds number based on the hull length is 1.0x 105, An immersed boundary method based on the moving-least-squares reconstruction is used to handle the complex geometric boundaries. The adaptive mesh refinement is utilized to resolve the flows near the hull, The parallel scalabilities of the flow solver are tested on meshes with the number of cells varying from 50 million to 3.2 billion, The parallel solver reaches nearly linear scalability for the flows around the underwater vehicle model, The present simulation captures the essential features of the vortex structures near the hull and in the wake, Both of the time-averaged pressure coefficients and srreamwise velocity profiles obtained from the LES are consistent with the characteristics of the flows pass an appended axisymmetric body. The code efficiency and its correct predictions on flow features allow us to perform the full-scale simulations on tens of thousands of cores with billions of grid points for higher-Reynolds-number flows around the underwater vehicles.
基金subsidized by the Central Project of Water Resource Fees[grant number 1261320212020]
文摘China experienced a decline of water use intensity in the 11th Five Year Plan,but the water use intensity actually increased in 2009.To the best of our knowledge,the index decomposition analysis method was rarely used to analyze changes in water use,and no decomposition analysis has investigated the role of regional economy in the decline of water use intensity.In this paper,we use logarithmic mean Divisia index(LMDI)techniques to decompose the change of water use intensity in the period 2006-2010.We find that the change of industrial water use intensity is confirmed as the dominant contributor to the decline in the overall water use intensity;the regional structure effect and the industrial structure effect is positive to the decline of overall water use intensity;the decline of China's water use intensity is mainly attributed to the effect of developed eastern provinces;meanwhile,the effect of central and undeveloped western is also positive to the decline of overall water use intensity;at least one out of three effects is positive to the decline of water use intensity in the different provinces;the intensity effect is positive and the industrial structure effect is positive to the declines of China's water use intensity based on chaining approach except the period 2008-2009,individually;and the deviation of regional structure effect and industrial structure effect between with regional economy and without regional economy in LMDI is 0.9 and2.3 m^3/10~4 RMB,respectively.
文摘We present a high performance modularly-built open-source software-OpenIFEM.OpenIFEM is a C++implementation of the modified immersed finite element method(mIFEM)to solve fluid-structure interaction(FSI)problems.This software is modularly built to perform multiple tasks including fluid dynamics(incompressible and slightly compressible fluid models),linear and nonlinear solid mechanics,and fully coupled fluid-structure interactions.Most of open-source software packages are restricted to certain discretization methods;some are under-tested,under-documented,and lack modularity as well as extensibility.OpenIFEM is designed and built to include a set of generic classes for users to adapt so that any fluid and solid solvers can be coupled through the FSI algorithm.In addition,the package utilizes well-developed and tested libraries.It also comes with standard test cases that serve as software and algorithm validation.The software can be built on cross-platform,i.e.,Linux,Windows,and Mac OS,using CMake.Efficient parallelization is also implemented for high-performance computing for large-sized problems.OpenIFEM is documented using Doxygen and publicly available to download on GitHub.It is expected to benefit the future development of FSI algorithms and be applied to a variety of FSI applications.
文摘This paper presents a refined method for estimating the annual extreme wave heights at a coastal or offshore project site on the basis of the data acquired at some nearby routine hydrographic stations. This method is based on the orthogonality principle in linear mean square estimation of stochastic processes. The error of the method is analyzed and compared with that of the conventional method. It is found that the method is able to effectively reduce the error so long as some feasible measures are adopted. A simulated test of the method has been conducted in a large scale wind wave flume. The test results are in good agreement with those given by theoretical error analysis. A scheme to implement the method is proposed on the basis of error analysis. The scheme is so designed as to reduce the estimation error as far as possible. This method is also suitable to utilizing satellite wave data for the estimation.
基金Sponsored by the National Natural Science Foundation of China(10676120)Laboratory of Computational Physics Foundation(9140C690101070C69)
文摘The adaptive mesh refinement （AMR） method is applied in the 2-D Euler multi-component elasticplastic hydrodynamics code （MEPH2Y）. It is applied on detonation. Firstly, the AMR method is described, including a conservative spatial interpolation, the time integration methodology with the adapitve time increment and an adaptive computational region method. The advantage of AMR technique is exhibited by numerical examples, including the 1-D C-J detonation and the 2-D implosion ignited from a single point. Results show that AMR can promote the computational efficiency, keeping the accuracy in interesting regions.
文摘Britholites with the general formula Ca10-xLnx（PO4）6-x（SiO4）xF2, （0≤x≤6） are considered to be promising matrices for the confinement of the by-products in the nuclear industry. A thermodynamic study showed that the stability of these compounds de-creased as the substitution rate increased. The present work was an attempt to gain more information about the structural changes in-duced by the substitution, in order to understand the observed stability decrease. The samples were successfully synthesized as a sin-gle-phase apatite by a solid-state reaction between 1200 and 1400 °C. The structural refinement indicated that the La3＋ions preferen-tially occupied the 6（h） sites. A progressive shift of F-along the c-axis outside its ideal position occurred as a result of the substitution increase. This might be the cause of the observed stability decrease, especially as the energies of the La-O, La-F and Si-O bonds are higher than those of Ca-O, Ca-F and P-O. The distribution of La3＋ between the two non-equivalent sites was confirmed by the charge distribution method.
基金This work was funded by the National Natural Science Foundation of China(Grant No.42072172,41772120)the Shandong Province Natural Science Fund for Distinguished Young Scholars(Grant No.JQ201311)+1 种基金the SDUST Research Fund(Grant No.2015TDJH101)Associated editor Jie Hao and seven anonymous reviewers were deeply acknowledged for their critical comments and helpful suggestions,which greatly improved the early version of this manuscript.
文摘The Carboniferous volcanic reservoir in the Chepaizi Uplift became an exploration hot target in recent years for its substantial amount of oils discovered. However, most of the Carboniferous heavy oils were biodegraded to PM7 or higher with orders of magnitude variation in oil viscosities. Two oil groups (I and II) exactly corresponding to the western and eastern Chepaizi Uplift were distinguished according to their source diagnose. Furthermore, three oil families (II1, II2 and II3), with the biodegradation level of PM7, PM8–8+, PM9+, respectively, were classified based on molecular compositions and parameter-stripping method of strongly bioresistant parameters. Allowing for this extremely high biodegradation case, more biodegradation refractory compound class were added to establish a refined Manco scale to quantitatively evaluate the biodegradation extent. Refined Manco number (RMN2) positively correlated with the oil density, NSO contents, and absolute concentrations of diasteranes and gammacerane, negatively correlated with the absolute concentrations of diahopane, summed tricyclic terpanes and pentacyclic terpanes. This refined scale showed higher resolution than the PM one to differentiate the biodegradation extent of Carboniferous heavy oils from the Chepaizi Uplift, especially those with same PM values but different oil viscosities.
基金supported by the National Natural Science Foundation of China (Nos.10932011,10772181,10732090,10772012 and 11021262)the National Basic Research Program of China (No.2007CB814803)
文摘Hybrid molecule/cluster statistical thermodynamics （HMCST） method is an efficient tool to simulate nano-scale systems under quasi-static loading at finite temperature. In this paper, a self-adaptive algorithm is developed for this method. Explicit refinement criterion based on the gradient of slip shear deformation and a switching criterion based on generalized Einstein approximation is proposed respectively. Results show that this self-adaptive method can accurately find clusters to be refined or transferred to molecules, and efficiently refine or trans- fer the clusters. Furthermore, compared with fully atomistic simulation, the high computational efficiency of the self-adaptive method appears very attractive.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.10732010,10972010,and 11028206)the Opening Project of State Key Laboratory of Explosion Science and Technology,China (Grant No.KFJJ13-5M)
文摘An essentially conservative adaptive space time conservation element and solution element （CE/SE） method is pro- posed for the effective simulation of shock-induced instability with low computational cost. Its implementation is based on redefined conservation elements （CEs） and solution elements （SEs）, optimized interpolations and a Courant number insensitive CE/SE scheme. This approach is used in two applications, the Woodward double Mach reflection and a two- component Richtmyer-Meshkov instability experiment. This scheme reveals the essential features of the investigated cases, captures small unstable structures, and yields a solution that is consistent with the results from experiments or other high order methods.
基金This work was supported by the National Natural Science Foun-dation of China(No.41776194).
文摘The present study provides a three-dimensional volume-of-fluid method based on the adaptive mesh refinement technique.The projection method on the adaptive mesh is introduced for solving the incompressible Navier-Stokes equations.The octree structure mesh is employed to solve the flow velocities and the pressure.The developed solver is applied to simulate the deformation of the cubic droplet driven by the surface tension without the effect of the gravity.The numerical results well predict the shape evolution of the droplet.
文摘We propose a high-order enriched partition of unity finite element method for linear and nonlinear time-dependent diffusion problems.The solution of this class of problems often exhibits non-smooth features such as steep gradients and boundary layers which can be very challenging to recover using the conventional low-order finite element methods.A class of steady-state exponential functions has been widely used for enrichment and its performance to numerically solve these challenges has been demonstrated.However,these enrichment functions have been used only in context of the standard h-version refinement or the so-called q-version refinement.In this paper we demonstrate that the p-version refinement can also be a very attractive option in terms of the efficiency and the accuracy in the enriched partition of unity finite element method.First,the transient diffusion problem is integrated in time using a semi-implicit scheme and the semi-discrete problem is then integrated in space using the p-version enriched finite elements.Numerical results are presented for three test examples of timedependent diffusion problems in both homogeneous and heterogeneous media.The computed results show the significant improvement when using the p-version refined enriched approximations in the finite element analysis.In addition,these results support our expectations for a robust and high-order accurate enriched partition of unity finite element method.
基金supported by the National Natural Science Foundation of China(11001037,11102037,11290143)the Fundamental Research Funds for the Central Universities(DUT13LK07)
文摘The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh distortions. In a previous work, we constructed an 8-node quadrilateral spline element L8 using the triangular area coordinates and the B- net method, which can be insensitive to mesh distortions and possess the second order completeness in the Cartesian co- ordinates. In this paper, a thin plate spline element is devel- oped based on the spline element L8 and the refined tech- nique. Numerical examples show that the present element indeed possesses higher accuracy than the DKQ element for distorted meshes.
基金funding by the National Natural Science Foundation (Grant No. 50745048)
文摘Straight-and curved-bar refining plates are two important types of plates commonly used in disc refiners in the papermaking industry.Theoretically,the curved-bar refining plate has a relatively uniform bar interaction angle,which indicates uniform refining effects.The bar angle of the curved bar was proposed and two typical curved-bar plates,the three-stage radial curved-bar plate and isometric curved-bar plate,were designed in this paper.The arc equations of the curved-bar center line and curved-bar edges were established and finally,the specific edge load(SEL)of the curved-bar plate was derived.The determination of bar parameters was discussed,which provides a theoretical basis for the design of curved-bar plates.
