A supersonic turbulent flowfield involving the pseudo-shock waves in an isolator of a supersonic combustion ramjet is computed using two different LES codes which are a high-order upwind finite volume scheme, and a si...A supersonic turbulent flowfield involving the pseudo-shock waves in an isolator of a supersonic combustion ramjet is computed using two different LES codes which are a high-order upwind finite volume scheme, and a sixth order compact differencing scheme utilizing the localized artificial diffusivity method for stabilizing shock waves and employing a wall model to enable the use of coarse mesh. In the validation study where a supersonic turbulent boundary layer flow over a flat plate is examined, both LES codes are well validated using velocity profile in the boundary layer given by the hot-wire anemometry and normal stress given by the laser Doppler anemometry. In particular, the sixth order compact differencing scheme gives closer agreements with these experimental data. Then, the validated LES codes are applied to solve the Mach number 2.5 supersonic turbulent flowfield involving the pseudo-shock waves. It is shown that typical features of unsteady flowfield of the pseudo-shock waves are well obtained by both schemes. Again, it is indicated that the sixth order compact differencing scheme gives closer agreements with the existing velocity data obtained by particle image velocimetry and pressure fluctuation data on the wall surface. Besides, the computational cost of the compact differencing scheme is found to be 1/7 of that for the upwind finite volume scheme, even though a wall model is solved at each grid point on the wall surface. Therefore, the obtained results in the present study allow recommending the sixth order compact differencing scheme with a wall model for simulating supersonic turbulent flowfield in an isolator involving the pseudo-shock waves.展开更多
This paper presents fabrication and testing of a multiple-evaporator and multiple-condenser loop heat pipe (MLHP) with polytetrafluoroethylene (PTFE) porous media as wicks. The MLHP has two evaporators and two condens...This paper presents fabrication and testing of a multiple-evaporator and multiple-condenser loop heat pipe (MLHP) with polytetrafluoroethylene (PTFE) porous media as wicks. The MLHP has two evaporators and two condensers in a loop heat pipe in order to adapt to various changes of thermal condition in spacecraft. The PTFE porous media was used as the primary wicks to reduce heat leak from evaporators to compensation chambers. The tests were conducted under an atmospheric condition. In the tests that heat loads are applied to both evaporators, the MLHP was stably operated as with a LHP with a single evaporator and a single condenser. The relation between the sink temperature and the thermal resistance was experimentally evaluated. In the test with the heat load to one evaporator, the heat transfer from the heated evaporator to the unheated evaporator was confirmed. In the heat load switching test, in which the heat load is switched from one evaporator to another evaporator repeatedly, the MLHP could be stably operated. The loop operation with the large temperature difference between the heat sinks was also tested. From this result, the stable operation of the MLHP in the various conditions was demonstrated. It was also found that a flow regulator which prevents the uncondensed vapor from the condensers is required at the inlet of the common liquid line when one condenser has higher temperature and cannot condense the vapor in it.展开更多
A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two metho...A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two methods, i.e., the finite element method (FEM) and the boundary element method (BEM). The original problem was restored with continuity and equilibrium conditions being satisfied on the interface of the two sub-regions using an iterative algorithm. To speed up the convergence rate of the iterative algorithm, a dynamically changing relaxation parameter during iteration was introduced. An advantage of the proposed algorithm is that the locations of the nodes on the interface of the two sub-domains can be inconsistent. The validity of the algorithm is demonstrated by the consistence of the results of a numerical example obtained by the proposed method and those by the FEM, the BEM and a present finite element-boundary element (FE-BE) coupling method.展开更多
Based on analysis of deformation in an infinite isotropic elastic matrix containing an embedded elliptic crack, subject to far field triaxial compressive stress, the energy release rate and a mixed fracture criterion ...Based on analysis of deformation in an infinite isotropic elastic matrix containing an embedded elliptic crack, subject to far field triaxial compressive stress, the energy release rate and a mixed fracture criterion are obtained by using an energy balance approach. The additional compliance tensor induced by a single closed elliptic microcrack in a representative volume element and its in-plane growth is derived. The additional compliance tensor induced by the kinked growth of the elliptic microcrack is also obtained. The effect of the microcracks, randomly distributed both in geometric characteristics and orientations, is analyzed with the Taylor's scheme by introducing an appropriate probability density function. A micromechanical damage model for rocks and concretes under triaxial compression is obtained and experimentally verified.展开更多
Correction to:Oki,Y.,Yoshikawa,K.,Takeuchi,H.et al.Orbit insertion strategy of Hayabusa2's rover with large release uncertainty around the asteroid Ryugu.Astrodynamics 2020,4(4):309-329 https://doi.org/10.1007/s42...Correction to:Oki,Y.,Yoshikawa,K.,Takeuchi,H.et al.Orbit insertion strategy of Hayabusa2's rover with large release uncertainty around the asteroid Ryugu.Astrodynamics 2020,4(4):309-329 https://doi.org/10.1007/s42064-020-0080-y The article“Orbit insertion strategy of Hayabusa2’s rover with large release uncertainty around the asteroid Ryugu”written by Yusuke Oki,Kent Yoshikawa,Hiroshi Takeuchi et al.,was originally published electronically on the publisher’s internet portal(currently SpringerLink)on 05 November 2020 without open access.After publication in Volume 4,Issue 4,page 309–329,the author(s)decided to opt for Open Choice and to make the article an open access publication.展开更多
An efficient implicit lower-upper symmetric Gauss-Seidel(LU-SGS)solution approach has been applied to a high order spectral volume(SV)method for unstructured tetrahedral grids.The LU-SGS solver is preconditioned by th...An efficient implicit lower-upper symmetric Gauss-Seidel(LU-SGS)solution approach has been applied to a high order spectral volume(SV)method for unstructured tetrahedral grids.The LU-SGS solver is preconditioned by the block element matrix,and the system of equations is then solved with a LU decomposition.The compact feature of SV reconstruction facilitates the efficient solution algorithm even for high order discretizations.The developed implicit solver has shown more than an order of magnitude of speed-up relative to the Runge-Kutta explicit scheme for typical inviscid and viscous problems.A convergence to a high order solution for high Reynolds number transonic flow over a 3D wing with a one equation turbulence model is also indicated.展开更多
文摘A supersonic turbulent flowfield involving the pseudo-shock waves in an isolator of a supersonic combustion ramjet is computed using two different LES codes which are a high-order upwind finite volume scheme, and a sixth order compact differencing scheme utilizing the localized artificial diffusivity method for stabilizing shock waves and employing a wall model to enable the use of coarse mesh. In the validation study where a supersonic turbulent boundary layer flow over a flat plate is examined, both LES codes are well validated using velocity profile in the boundary layer given by the hot-wire anemometry and normal stress given by the laser Doppler anemometry. In particular, the sixth order compact differencing scheme gives closer agreements with these experimental data. Then, the validated LES codes are applied to solve the Mach number 2.5 supersonic turbulent flowfield involving the pseudo-shock waves. It is shown that typical features of unsteady flowfield of the pseudo-shock waves are well obtained by both schemes. Again, it is indicated that the sixth order compact differencing scheme gives closer agreements with the existing velocity data obtained by particle image velocimetry and pressure fluctuation data on the wall surface. Besides, the computational cost of the compact differencing scheme is found to be 1/7 of that for the upwind finite volume scheme, even though a wall model is solved at each grid point on the wall surface. Therefore, the obtained results in the present study allow recommending the sixth order compact differencing scheme with a wall model for simulating supersonic turbulent flowfield in an isolator involving the pseudo-shock waves.
文摘This paper presents fabrication and testing of a multiple-evaporator and multiple-condenser loop heat pipe (MLHP) with polytetrafluoroethylene (PTFE) porous media as wicks. The MLHP has two evaporators and two condensers in a loop heat pipe in order to adapt to various changes of thermal condition in spacecraft. The PTFE porous media was used as the primary wicks to reduce heat leak from evaporators to compensation chambers. The tests were conducted under an atmospheric condition. In the tests that heat loads are applied to both evaporators, the MLHP was stably operated as with a LHP with a single evaporator and a single condenser. The relation between the sink temperature and the thermal resistance was experimentally evaluated. In the test with the heat load to one evaporator, the heat transfer from the heated evaporator to the unheated evaporator was confirmed. In the heat load switching test, in which the heat load is switched from one evaporator to another evaporator repeatedly, the MLHP could be stably operated. The loop operation with the large temperature difference between the heat sinks was also tested. From this result, the stable operation of the MLHP in the various conditions was demonstrated. It was also found that a flow regulator which prevents the uncondensed vapor from the condensers is required at the inlet of the common liquid line when one condenser has higher temperature and cannot condense the vapor in it.
基金Project supported by China Postdoctoral Science Foundation (No.2004036145)
文摘A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two methods, i.e., the finite element method (FEM) and the boundary element method (BEM). The original problem was restored with continuity and equilibrium conditions being satisfied on the interface of the two sub-regions using an iterative algorithm. To speed up the convergence rate of the iterative algorithm, a dynamically changing relaxation parameter during iteration was introduced. An advantage of the proposed algorithm is that the locations of the nodes on the interface of the two sub-domains can be inconsistent. The validity of the algorithm is demonstrated by the consistence of the results of a numerical example obtained by the proposed method and those by the FEM, the BEM and a present finite element-boundary element (FE-BE) coupling method.
基金supported by the National Natural Science Foundation of China (Nos. 10872220 and 50725414)Japan Society for the Promotion of Science JSPS (No. L08538)
文摘Based on analysis of deformation in an infinite isotropic elastic matrix containing an embedded elliptic crack, subject to far field triaxial compressive stress, the energy release rate and a mixed fracture criterion are obtained by using an energy balance approach. The additional compliance tensor induced by a single closed elliptic microcrack in a representative volume element and its in-plane growth is derived. The additional compliance tensor induced by the kinked growth of the elliptic microcrack is also obtained. The effect of the microcracks, randomly distributed both in geometric characteristics and orientations, is analyzed with the Taylor's scheme by introducing an appropriate probability density function. A micromechanical damage model for rocks and concretes under triaxial compression is obtained and experimentally verified.
文摘Correction to:Oki,Y.,Yoshikawa,K.,Takeuchi,H.et al.Orbit insertion strategy of Hayabusa2's rover with large release uncertainty around the asteroid Ryugu.Astrodynamics 2020,4(4):309-329 https://doi.org/10.1007/s42064-020-0080-y The article“Orbit insertion strategy of Hayabusa2’s rover with large release uncertainty around the asteroid Ryugu”written by Yusuke Oki,Kent Yoshikawa,Hiroshi Takeuchi et al.,was originally published electronically on the publisher’s internet portal(currently SpringerLink)on 05 November 2020 without open access.After publication in Volume 4,Issue 4,page 309–329,the author(s)decided to opt for Open Choice and to make the article an open access publication.
文摘An efficient implicit lower-upper symmetric Gauss-Seidel(LU-SGS)solution approach has been applied to a high order spectral volume(SV)method for unstructured tetrahedral grids.The LU-SGS solver is preconditioned by the block element matrix,and the system of equations is then solved with a LU decomposition.The compact feature of SV reconstruction facilitates the efficient solution algorithm even for high order discretizations.The developed implicit solver has shown more than an order of magnitude of speed-up relative to the Runge-Kutta explicit scheme for typical inviscid and viscous problems.A convergence to a high order solution for high Reynolds number transonic flow over a 3D wing with a one equation turbulence model is also indicated.