This paper focuses on the research of a semi-submersible platform equipped with a DP-assisted mooring system. Based on the working principles of the DP-assisted mooring system and the model of the platform motion, a t...This paper focuses on the research of a semi-submersible platform equipped with a DP-assisted mooring system. Based on the working principles of the DP-assisted mooring system and the model of the platform motion, a time domain simulation program is applied to analyze the impact, in the case of one line failure, on the platform motion, power consumption of the thrusters and the tension of the mooring lines. The results show that, under the 10-year wind dominant, a one line failure will have little impact on the tension of the mooring lines. When the failure line is windward, the power consumption will increase greatly with a weakened position of accuracy. However when the failure line is leeward, the power consumption will be reduced with a partly strengthened oosition of accuracy.展开更多
Direct time-domain simulation of floating structures has advantages:it can calculate wave pressure fields and forces directly; and it is useful for coupled analysis of floating structures with a mooring system. A time...Direct time-domain simulation of floating structures has advantages:it can calculate wave pressure fields and forces directly; and it is useful for coupled analysis of floating structures with a mooring system. A time-domain boundary integral equation method is presented to simulate three-dimensional water wave radiation problems. A stable form of the integration free-surface boundary condition (IFBC) is used to update velocity potentials on the free surface. A multi-transmitting formula (MTF) method with an artificial speed is introduced to the artificial radiation boundary (ARB). The method was applied to simulate a semi-spherical liquefied natural gas (LNG) carrier and a semi-submersible undergoing specified harmonic motion. Numerical parameters such as the form of the ARB, and the time and space discretization related to this method are discussed. It was found that a good agreement can be obtained when artificial speed is between 0.6 and 1.6 times the phase velocity of water waves in the MTF method. A simulation can be done for a long period of time by this method without problems of instability, and the method is also accurate and computationally efficient.展开更多
The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems wit...The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems within a time domain framework,the free water surface needs to simultaneously satisfy both the kinematic and dynamic boundary conditions of the free water surface.This provides conditions for adding artificial damping layers.Using the Runge−Kutta method to solve equations related to time.An upwind differential scheme is used in the present method to deal with the convection terms on the free surface to prevent waves upstream.Through the comparison with the available experimental data and other numerical methods,the present method is proved to have good mesh convergence,and satisfactory results can be obtained.The constant panel method is applied to calculate the hydrodynamic interaction responses of two parallel ships advancing in head waves.Numerical simulations are conducted on the effects of forward speed,different longitudinal and lateral distances on the motion response of two modified Wigley ships in head waves.Then further investigations are conducted on the effects of different ship types on the motion response.展开更多
Based on a simplified domain rotation model, the rotation path of internal domains and cor- responding magnetostriction of <112>, <111> oriented single Tedenol-D crystals under com- pressive prestresses ha...Based on a simplified domain rotation model, the rotation path of internal domains and cor- responding magnetostriction of <112>, <111> oriented single Tedenol-D crystals under com- pressive prestresses have been simulated. Comparisons with results of experiment and other calculation have been made. Results of simulation showed that the <111> oriented single crystal has better low-field magnetostriction properties than the <112> oriented one. Under a compressive prestress of 10 MPa, up to 2300×10-6 saturation magnetostriction of <111> oriented crystal has been obtained at 800 Oe, while for <112> oriented. 1600×10-6 maximum magnetostriction has been reached at 1000 Oe.展开更多
The three-bucket jacket foundation is a new type of foundation for offshore wind turbine that has the advantages of fast construction speed and suitability for deep water. The study of the hoisting and launching proce...The three-bucket jacket foundation is a new type of foundation for offshore wind turbine that has the advantages of fast construction speed and suitability for deep water. The study of the hoisting and launching process is of great significance to ensure construction safety in actual projects. In this paper, a new launching technology is proposed that is based on the foundation of the three-bucket jacket for offshore wind turbine. A complete time domain simulation of the launching process of three-bucket jacket foundation is carried out by a theoretical analysis combined with hydrodynamic software Moses. At the same time, the effects of different initial air storage and sea conditions on the motion response of the structure and the hoisting cable tension are studied. The results show that the motion response of the structure is the highest when it is lowered to 1.5 times the bucket height. The natural period of each degree of freedom of the structure increases with the increase of the lowering depth. The structural motion response and the hoisting cable tension vary greatly in the early phases of Stages Ⅰ and Ⅲ, smaller in Stage Ⅱ, and gradually stabilize in the middle and late phases of Stage Ⅲ.展开更多
Tank sloshing in ship cargo is excited by ship motions, which induces impact load on tank wall and then affects the ship motion. Wave forces acting on ship hull and the retardation function are solved by using three-d...Tank sloshing in ship cargo is excited by ship motions, which induces impact load on tank wall and then affects the ship motion. Wave forces acting on ship hull and the retardation function are solved by using three-dimensional frequency domain theory and an impulse response function method based on the potential flow theory, and global ship motion is examined coupling with nonlinear tank sloshing which is simulated by viscous flow theory. Based on the open source Computational Fluid Dynamics (CFD) development platform Open Field Operation and Manipulation (OpenFOAM), numerical calculation of ship motion coupled with tank sloshing is achieved and the corresponding numerical simulation and validation are carried out. With this method, the interactions of wave, ship body and tank sloshing are completely taken into consideration. This method has quite high efficiency for it takes advantage of potential flow theory for outer flow field and viscous flow theory for inside tank sloshing respectively. The numerical and experimental results of the ship motion agree well with each other.展开更多
The real-time computer-controlled actuators are used to connect the truncated parts of moorings and risers in the active hybrid model testing system. This must be able to work in model-scale real time, based on feedba...The real-time computer-controlled actuators are used to connect the truncated parts of moorings and risers in the active hybrid model testing system. This must be able to work in model-scale real time, based on feedback input from the floater motions. Thus, mooring line dynamics and damping effects are artificially simulated in real time, based on a computer-based model of the problem. In consideration of the nonlinear characteristics of the sea platform catenary mooring line, the equations of the mooring line motion are formulated by using the lumped-mass method and the dynamic response of several points on the mooring line is investigated by the time and frequency domain analysis method. The dynamic response of the representative point on the mooring line is analyzed under the condition of two different corresponding upper endpoint movements namely sine wave excitation and random wave excitation. The corresponding laws of the dynamic response between the equivalent water depth truncated points at different locations and the upper endpoint are obtained, which can provide technical support for further study of the active hybrid model test.展开更多
The cross-flow(CF)vortex-induced vibration(VIV)of a deepwater steep wave riser(SWR)subjected to uniform or shear flow loads is investigated numerically.The model is based on a three-dimensional(3D)nonlinear elastic ro...The cross-flow(CF)vortex-induced vibration(VIV)of a deepwater steep wave riser(SWR)subjected to uniform or shear flow loads is investigated numerically.The model is based on a three-dimensional(3D)nonlinear elastic rod theory coupled with a wake oscillator model.In this numerical simulation,the nonlinear motion equations of the riser with large deformation features are established in a global coordinate system to avoid the transformation between global and local coordinate systems,and are discretized with the time-domain finite element method(FEM).A wakeoscillator model is employed to study the vortex shedding,and the lift force generated by the wake flow is described in a van der Pol equation.A Newmark-βiterative scheme is used to solve their coupling equation for the VIV response of the SWR.The developed model is validated against the existing experimental results for the VIV response of the top-tension riser(TTR).Then,the numerical simulations are executed to determine VIV characteristics of the SWR.The effects of both flow velocity and the spanwise length of the flow field on the drag coefficient in the inline(IL)direction and the lift coefficient in the CF direction are investigated systematically.The results indicate that compared with TTR,the low frequency and multi-modal vibration are the main components of the SWR due to the large deformation and flexible characteristics.For shear flow,the multi-frequency resonance dominates the VIV response of the SWR,especially at the hang-off segment.展开更多
By utilizing a Fabry–Perot (FP) nanocavity adjacent to T-shaped gap waveguide ports, spectrally selective filtering is realized. When the wavelength of incident light corresponds to the resonance wavelength of the ...By utilizing a Fabry–Perot (FP) nanocavity adjacent to T-shaped gap waveguide ports, spectrally selective filtering is realized. When the wavelength of incident light corresponds to the resonance wavelength of the FP nanocavity, the surface plasmons are captured inside the nanocavity, and light is highly reflected from this port. The resonance wavelength is determined by using Fabry–Perot resonance condition for the nanocavity. For any desired filtering frequency the dimension of the nanocavity can be tailored. The numerical results are based on the two-dimensional finite difference time domain simulation under a perfectly matched layer absorbing boundary condition. The analytical and simulation results indicate that the proposed structure can be utilized for filtering and splitting applications.展开更多
Aiming at the dynamic response of reticulated shell structures under wind load,systematic parameter analyses on wind-induced responses of Kiewitt6-6 type single-layer spherical reticulated shell structures and three-w...Aiming at the dynamic response of reticulated shell structures under wind load,systematic parameter analyses on wind-induced responses of Kiewitt6-6 type single-layer spherical reticulated shell structures and three-way grid single-layer cylindrical reticulated shell structures were performed with the random simulation method in time domain,including geometric parameters,structural parameters and aerodynamic parameters.Moreover,a wind-induced vibration coefficient was obtained,which can be a reference to the wind-resistance design of reticulated shell structures.The results indicate that the geometric parameters are the most important factor influencing wind-induced responses of the reticulated shell structures;the wind-induced vibration coeffi-cient is 3.0-3.2 for the spherical reticulated shell structures and that is 2.8-3.0 for the cylindrical reticula-ted shell structures,which shows that the wind-induced vibration coefficients of these two kinds of space frames are well-proportioned.展开更多
Choanoid fluidized bed bioreactors (CFBBs) are newly developed core devices used in bioartificial liver- support systems to detoxify blood plasma of patients with microencapsulated liver cells. Direct numerical simu...Choanoid fluidized bed bioreactors (CFBBs) are newly developed core devices used in bioartificial liver- support systems to detoxify blood plasma of patients with microencapsulated liver cells. Direct numerical simulations (DNS) with a direct-forcing/fictitious domain (DF/FD) method were conducted to study the hydrodynamic performance of a CFBB. The effects of particle-fluid density ratio, particle number, and fil- ter screens preventing particles flowing out of the reactor were investigated. Depending on density ratio, two flow patterns are evident: the circulation mode in which the suspension rises along one sidewall and descends along the other sidewall, and the non-circulation mode in which the whole suspension roughly flows upward. The circulation mode takes place under non-neutral-buoyancy where the particle sedimentation dominates, whereas the non-circulation mode occurs under pure or near-neutral buoy- ancy with particle-fluid density ratios of unity or near unity. With particle-fluid density ratio of 1.01, the bioartificial liver reactor performs optimally as the significant particle accumulation existing in the non-circulation mode and the large shear forces on particles in the circulation mode are avoided. At higher particle volume fractions, more particles accumulate at the filter screens and a secondary counter circulation to the primary flow is observed at the top of the bed. Modelled as porous media, the filter screens play a negative role on particle fluidization velocities; without screens, particles are fluidized faster because of the higher fluid velocities in the jet center region. This work extends the DF/FD-based DNS to a fluidized bed and accounts for effects from inclined side walls and porous media, providing some hydrodynamics insight that is important for CFBB design and operation optimization.展开更多
The Taylor model arithmetic is introduced to deal with uncertainty.The uncertainty of model parameters is described by Taylor models and each variable in functions is replaced with the Taylor model(TM).Thus,time domai...The Taylor model arithmetic is introduced to deal with uncertainty.The uncertainty of model parameters is described by Taylor models and each variable in functions is replaced with the Taylor model(TM).Thus,time domain simulation under uncertainty is transformed to the integration of TM-based differential equations.In this paper,the Taylor series method is employed to compute differential equations;moreover,power system time domain simulation under uncertainty based on Taylor model method is presented.This method allows a rigorous estimation of the influence of either form of uncertainty and only needs one simulation.It is computationally fast compared with the Monte Carlo method,which is another technique for uncertainty analysis.The proposed method has been tested on the 39-bus New England system.The test results illustrate the effectiveness and practical value of the approach by comparing with the results of Monte Carlo simulation and traditional time domain simulation.展开更多
Based on characteristic functions of variants, we developed an unconventional phase field modeling for investigating domains formation and evolution in tetragonal ferroelectrics. In order to develop this computational...Based on characteristic functions of variants, we developed an unconventional phase field modeling for investigating domains formation and evolution in tetragonal ferroelectrics. In order to develop this computational approach, we constructed the anisotropy energy of tetragonal variants, which is used instead of Landau-Devonshire potential in the conventional phase field method, resulting in that much fewer parameters are needed for simulations. This approach is advantageous in simulations of emerging ferroelectric materials. We employ it to study the formation and evolution of domains in tetragonal barium titanate single crystal, as well as the nonlinear behaviors under cyclical stress and electric field loading. A multi-rank laminated ferroelectric domain pattern, 90° domain switching accompanied by polarization rotation, and 180° domain switching accompanied by move of domain wall are predicted. It is found that the speed of 90° domain switching is slower than that of 180° domain switching, due to both polarization and transformation strain changed in 90° domain switching. It also suggests that large strain actuation can be generated in single crystal ferroelectrics via combined electromechanical loading inducing 90° domain switching. The good agreement between simulation results and experimental measurements is observed.展开更多
The unique structural features represented by micro-nanoneedle tip structure reflect wonderful physical and chemical properties.The tip effect includes the concentration of energy such as electrons,photons and magneti...The unique structural features represented by micro-nanoneedle tip structure reflect wonderful physical and chemical properties.The tip effect includes the concentration of energy such as electrons,photons and magnetism in the tip region,which has promising applications in the fields of energy conversion,water capture,environmental restoration and so on.In this review,a comprehensive and systematic summary of the latest advances in the application of the tip effect in different fields is provided.Utilizing advanced Finite Difference Time Domain simulation,we further propose our understanding of the fundamental mechanism of the tip effect induced by micro-nanostructure.However,we need to forge the present study to further reveal the essential law of the tip effect from the perspective of theoretical calculations.This review would provide a solid foundation for further development and application of the tip effect.展开更多
In this paper, we combine the direct-forcing fictitious domain (DF/FD) method and the sharp interface method to resolve the problem of particle dielectrophoresis in two dimensions. The flow field and the motion of p...In this paper, we combine the direct-forcing fictitious domain (DF/FD) method and the sharp interface method to resolve the problem of particle dielectrophoresis in two dimensions. The flow field and the motion of particles are solved with the DF/FD method, the electric field is solved with the sharp inter- face method, and the electrostatic force on the particles is computed using the Maxwell stress tensor method. The proposed method is validated via three problems: effective conductivity of particle compos- ite between two planar plates, cell trapping in a channel, and motion of particles due to both conventional and traveling wave dielectrophoretic forces.展开更多
基金Suppirted by the Programme of Introducing Talents of Discipline to Universities(B07019)
文摘This paper focuses on the research of a semi-submersible platform equipped with a DP-assisted mooring system. Based on the working principles of the DP-assisted mooring system and the model of the platform motion, a time domain simulation program is applied to analyze the impact, in the case of one line failure, on the platform motion, power consumption of the thrusters and the tension of the mooring lines. The results show that, under the 10-year wind dominant, a one line failure will have little impact on the tension of the mooring lines. When the failure line is windward, the power consumption will increase greatly with a weakened position of accuracy. However when the failure line is leeward, the power consumption will be reduced with a partly strengthened oosition of accuracy.
基金Supported by the National Natural Science Foundation of China under Grant No.10572041,50779008 and the 111 Project
文摘Direct time-domain simulation of floating structures has advantages:it can calculate wave pressure fields and forces directly; and it is useful for coupled analysis of floating structures with a mooring system. A time-domain boundary integral equation method is presented to simulate three-dimensional water wave radiation problems. A stable form of the integration free-surface boundary condition (IFBC) is used to update velocity potentials on the free surface. A multi-transmitting formula (MTF) method with an artificial speed is introduced to the artificial radiation boundary (ARB). The method was applied to simulate a semi-spherical liquefied natural gas (LNG) carrier and a semi-submersible undergoing specified harmonic motion. Numerical parameters such as the form of the ARB, and the time and space discretization related to this method are discussed. It was found that a good agreement can be obtained when artificial speed is between 0.6 and 1.6 times the phase velocity of water waves in the MTF method. A simulation can be done for a long period of time by this method without problems of instability, and the method is also accurate and computationally efficient.
基金supported by the National Natural Science Foundation of China(Grant Nos.52271278 and 52111530137)the Natural Science Found of Jiangsu Province(Grant No.BK20221389)the Newton Advanced Fellowships(Grant No.NAF\R1\180304)by the Royal Society.
文摘The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems within a time domain framework,the free water surface needs to simultaneously satisfy both the kinematic and dynamic boundary conditions of the free water surface.This provides conditions for adding artificial damping layers.Using the Runge−Kutta method to solve equations related to time.An upwind differential scheme is used in the present method to deal with the convection terms on the free surface to prevent waves upstream.Through the comparison with the available experimental data and other numerical methods,the present method is proved to have good mesh convergence,and satisfactory results can be obtained.The constant panel method is applied to calculate the hydrodynamic interaction responses of two parallel ships advancing in head waves.Numerical simulations are conducted on the effects of forward speed,different longitudinal and lateral distances on the motion response of two modified Wigley ships in head waves.Then further investigations are conducted on the effects of different ship types on the motion response.
文摘Based on a simplified domain rotation model, the rotation path of internal domains and cor- responding magnetostriction of <112>, <111> oriented single Tedenol-D crystals under com- pressive prestresses have been simulated. Comparisons with results of experiment and other calculation have been made. Results of simulation showed that the <111> oriented single crystal has better low-field magnetostriction properties than the <112> oriented one. Under a compressive prestress of 10 MPa, up to 2300×10-6 saturation magnetostriction of <111> oriented crystal has been obtained at 800 Oe, while for <112> oriented. 1600×10-6 maximum magnetostriction has been reached at 1000 Oe.
基金financially supported by the National Natural Science Foundation of China (Grant No.52171274)。
文摘The three-bucket jacket foundation is a new type of foundation for offshore wind turbine that has the advantages of fast construction speed and suitability for deep water. The study of the hoisting and launching process is of great significance to ensure construction safety in actual projects. In this paper, a new launching technology is proposed that is based on the foundation of the three-bucket jacket for offshore wind turbine. A complete time domain simulation of the launching process of three-bucket jacket foundation is carried out by a theoretical analysis combined with hydrodynamic software Moses. At the same time, the effects of different initial air storage and sea conditions on the motion response of the structure and the hoisting cable tension are studied. The results show that the motion response of the structure is the highest when it is lowered to 1.5 times the bucket height. The natural period of each degree of freedom of the structure increases with the increase of the lowering depth. The structural motion response and the hoisting cable tension vary greatly in the early phases of Stages Ⅰ and Ⅲ, smaller in Stage Ⅱ, and gradually stabilize in the middle and late phases of Stage Ⅲ.
文摘Tank sloshing in ship cargo is excited by ship motions, which induces impact load on tank wall and then affects the ship motion. Wave forces acting on ship hull and the retardation function are solved by using three-dimensional frequency domain theory and an impulse response function method based on the potential flow theory, and global ship motion is examined coupling with nonlinear tank sloshing which is simulated by viscous flow theory. Based on the open source Computational Fluid Dynamics (CFD) development platform Open Field Operation and Manipulation (OpenFOAM), numerical calculation of ship motion coupled with tank sloshing is achieved and the corresponding numerical simulation and validation are carried out. With this method, the interactions of wave, ship body and tank sloshing are completely taken into consideration. This method has quite high efficiency for it takes advantage of potential flow theory for outer flow field and viscous flow theory for inside tank sloshing respectively. The numerical and experimental results of the ship motion agree well with each other.
基金financially supported by the Natural Science Foundation of Zhejiang Province(Grant Nos.Y14E090034 and Y13F020140)the Young Scientist Training Program in Zhejiang Province(Grant No.2013R60G7160040)+1 种基金the State Key Laboratory of Ocean Engineering of Shanghai Jiao Tong University for the Open Fund Project(Grant No.1516)the Open Fund Project of Second Institute of Oceanography(Grant No.SOED1706)
文摘The real-time computer-controlled actuators are used to connect the truncated parts of moorings and risers in the active hybrid model testing system. This must be able to work in model-scale real time, based on feedback input from the floater motions. Thus, mooring line dynamics and damping effects are artificially simulated in real time, based on a computer-based model of the problem. In consideration of the nonlinear characteristics of the sea platform catenary mooring line, the equations of the mooring line motion are formulated by using the lumped-mass method and the dynamic response of several points on the mooring line is investigated by the time and frequency domain analysis method. The dynamic response of the representative point on the mooring line is analyzed under the condition of two different corresponding upper endpoint movements namely sine wave excitation and random wave excitation. The corresponding laws of the dynamic response between the equivalent water depth truncated points at different locations and the upper endpoint are obtained, which can provide technical support for further study of the active hybrid model test.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52111530137 and 52025112)the Natural Science Found of Jiangsu Province(Grant No.BK20160556)the Jiangsu Provincial Higher Education Natural Science Research Major Project(Grant No.18KJA580003)。
文摘The cross-flow(CF)vortex-induced vibration(VIV)of a deepwater steep wave riser(SWR)subjected to uniform or shear flow loads is investigated numerically.The model is based on a three-dimensional(3D)nonlinear elastic rod theory coupled with a wake oscillator model.In this numerical simulation,the nonlinear motion equations of the riser with large deformation features are established in a global coordinate system to avoid the transformation between global and local coordinate systems,and are discretized with the time-domain finite element method(FEM).A wakeoscillator model is employed to study the vortex shedding,and the lift force generated by the wake flow is described in a van der Pol equation.A Newmark-βiterative scheme is used to solve their coupling equation for the VIV response of the SWR.The developed model is validated against the existing experimental results for the VIV response of the top-tension riser(TTR).Then,the numerical simulations are executed to determine VIV characteristics of the SWR.The effects of both flow velocity and the spanwise length of the flow field on the drag coefficient in the inline(IL)direction and the lift coefficient in the CF direction are investigated systematically.The results indicate that compared with TTR,the low frequency and multi-modal vibration are the main components of the SWR due to the large deformation and flexible characteristics.For shear flow,the multi-frequency resonance dominates the VIV response of the SWR,especially at the hang-off segment.
基金Project supported by the National Key Basic Research Program of China (Grant No. 2013CB328702)
文摘By utilizing a Fabry–Perot (FP) nanocavity adjacent to T-shaped gap waveguide ports, spectrally selective filtering is realized. When the wavelength of incident light corresponds to the resonance wavelength of the FP nanocavity, the surface plasmons are captured inside the nanocavity, and light is highly reflected from this port. The resonance wavelength is determined by using Fabry–Perot resonance condition for the nanocavity. For any desired filtering frequency the dimension of the nanocavity can be tailored. The numerical results are based on the two-dimensional finite difference time domain simulation under a perfectly matched layer absorbing boundary condition. The analytical and simulation results indicate that the proposed structure can be utilized for filtering and splitting applications.
基金the National Natural Science Foundation of China (Grant No. 50608022)the Foundation of National Science and Technology(GrantNo.2006BAJ03B04)
文摘Aiming at the dynamic response of reticulated shell structures under wind load,systematic parameter analyses on wind-induced responses of Kiewitt6-6 type single-layer spherical reticulated shell structures and three-way grid single-layer cylindrical reticulated shell structures were performed with the random simulation method in time domain,including geometric parameters,structural parameters and aerodynamic parameters.Moreover,a wind-induced vibration coefficient was obtained,which can be a reference to the wind-resistance design of reticulated shell structures.The results indicate that the geometric parameters are the most important factor influencing wind-induced responses of the reticulated shell structures;the wind-induced vibration coeffi-cient is 3.0-3.2 for the spherical reticulated shell structures and that is 2.8-3.0 for the cylindrical reticula-ted shell structures,which shows that the wind-induced vibration coefficients of these two kinds of space frames are well-proportioned.
基金The authors gratefully acknowledge the supports from China Postdoctoral Science Foundation (Grant No. 2014M550327), the opening foundation of the State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, and the National Natural Science Foundation of China (Grant No. 11372275). The authors are also grateful to Chengbo Yu and Liang Yu for their introduction of the choanoid fluidized bed bioreactor and helpful discussions.
文摘Choanoid fluidized bed bioreactors (CFBBs) are newly developed core devices used in bioartificial liver- support systems to detoxify blood plasma of patients with microencapsulated liver cells. Direct numerical simulations (DNS) with a direct-forcing/fictitious domain (DF/FD) method were conducted to study the hydrodynamic performance of a CFBB. The effects of particle-fluid density ratio, particle number, and fil- ter screens preventing particles flowing out of the reactor were investigated. Depending on density ratio, two flow patterns are evident: the circulation mode in which the suspension rises along one sidewall and descends along the other sidewall, and the non-circulation mode in which the whole suspension roughly flows upward. The circulation mode takes place under non-neutral-buoyancy where the particle sedimentation dominates, whereas the non-circulation mode occurs under pure or near-neutral buoy- ancy with particle-fluid density ratios of unity or near unity. With particle-fluid density ratio of 1.01, the bioartificial liver reactor performs optimally as the significant particle accumulation existing in the non-circulation mode and the large shear forces on particles in the circulation mode are avoided. At higher particle volume fractions, more particles accumulate at the filter screens and a secondary counter circulation to the primary flow is observed at the top of the bed. Modelled as porous media, the filter screens play a negative role on particle fluidization velocities; without screens, particles are fluidized faster because of the higher fluid velocities in the jet center region. This work extends the DF/FD-based DNS to a fluidized bed and accounts for effects from inclined side walls and porous media, providing some hydrodynamics insight that is important for CFBB design and operation optimization.
基金supported by the National Natural Science Foundation of China (Grant No.50477035).
文摘The Taylor model arithmetic is introduced to deal with uncertainty.The uncertainty of model parameters is described by Taylor models and each variable in functions is replaced with the Taylor model(TM).Thus,time domain simulation under uncertainty is transformed to the integration of TM-based differential equations.In this paper,the Taylor series method is employed to compute differential equations;moreover,power system time domain simulation under uncertainty based on Taylor model method is presented.This method allows a rigorous estimation of the influence of either form of uncertainty and only needs one simulation.It is computationally fast compared with the Monte Carlo method,which is another technique for uncertainty analysis.The proposed method has been tested on the 39-bus New England system.The test results illustrate the effectiveness and practical value of the approach by comparing with the results of Monte Carlo simulation and traditional time domain simulation.
基金supported by the National Natural Science Foundation of China(Grant Nos.11572276&11502225)Hunan Provincial Natural Science Foundation of China(Grant No.14JJ6015)
文摘Based on characteristic functions of variants, we developed an unconventional phase field modeling for investigating domains formation and evolution in tetragonal ferroelectrics. In order to develop this computational approach, we constructed the anisotropy energy of tetragonal variants, which is used instead of Landau-Devonshire potential in the conventional phase field method, resulting in that much fewer parameters are needed for simulations. This approach is advantageous in simulations of emerging ferroelectric materials. We employ it to study the formation and evolution of domains in tetragonal barium titanate single crystal, as well as the nonlinear behaviors under cyclical stress and electric field loading. A multi-rank laminated ferroelectric domain pattern, 90° domain switching accompanied by polarization rotation, and 180° domain switching accompanied by move of domain wall are predicted. It is found that the speed of 90° domain switching is slower than that of 180° domain switching, due to both polarization and transformation strain changed in 90° domain switching. It also suggests that large strain actuation can be generated in single crystal ferroelectrics via combined electromechanical loading inducing 90° domain switching. The good agreement between simulation results and experimental measurements is observed.
基金This work was financially supported by the Natural Science Foundation of China(No.51902101)Natural Science Foundation of Jiangsu Province(No.BK20201381)+1 种基金Science Foundation of Nanjing University of Posts and Telecommunications(No.NY219144)the National College Student Innovation and Entrepreneurship Training Program(No.202210293017Z).
文摘The unique structural features represented by micro-nanoneedle tip structure reflect wonderful physical and chemical properties.The tip effect includes the concentration of energy such as electrons,photons and magnetism in the tip region,which has promising applications in the fields of energy conversion,water capture,environmental restoration and so on.In this review,a comprehensive and systematic summary of the latest advances in the application of the tip effect in different fields is provided.Utilizing advanced Finite Difference Time Domain simulation,we further propose our understanding of the fundamental mechanism of the tip effect induced by micro-nanostructure.However,we need to forge the present study to further reveal the essential law of the tip effect from the perspective of theoretical calculations.This review would provide a solid foundation for further development and application of the tip effect.
基金support from the National Natural Science Foundation of China(no.10872181)the National Basic Research Program of China(no.2006CB705400)+1 种基金Chinese Universities Scientific Fundthe Major Program of the National Natural Science Foundation of China(no.10632070)
文摘In this paper, we combine the direct-forcing fictitious domain (DF/FD) method and the sharp interface method to resolve the problem of particle dielectrophoresis in two dimensions. The flow field and the motion of particles are solved with the DF/FD method, the electric field is solved with the sharp inter- face method, and the electrostatic force on the particles is computed using the Maxwell stress tensor method. The proposed method is validated via three problems: effective conductivity of particle compos- ite between two planar plates, cell trapping in a channel, and motion of particles due to both conventional and traveling wave dielectrophoretic forces.