Existing numerical methods for complex composites, such as multiscale simulation and neural network algorithms, face significant limitations. Multiscale techniques are often prohibitively expensive for large models, w...Existing numerical methods for complex composites, such as multiscale simulation and neural network algorithms, face significant limitations. Multiscale techniques are often prohibitively expensive for large models, while neural networks struggle to represent underlying microscopic material properties. To overcome these challenges, a meso-micro scale numerical method using a virtual node approach is developed in this study. A Wbraid/Al/Epoxy functional structural material is fabricated, and a representative periodic unit cell is identified based on its architecture. The complex structure is then discretized into nodes, and mechanical interactions are governed by pre-defined computation rules. This virtual node method is systematically compared against both multiscale simulation and a neural network algorithm, with validation provided through mechanical experiments. The results demonstrate that the nodal operation strategy significantly reduces computational resource requirements. By quantifying microscopic bonding with coefficients, explicit interface treatment is avoided, granting the method strong adaptability to lattice materials. The method can simulate extremely complex structures using parameters from simple tests and is suited for large systems. Compared to three-point bending experiments, errors for multiscale, virtual node, and neural network methods were 12.4%, 6.9%, and 34.5%, respectively. Under dynamic compression, the errors were 2.7%, 9.3%, and 15.43%. The virtual node method demonstrated superior accuracy under static conditions, enabling efficient prediction and auxiliary development of complex structural materials.展开更多
Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics h...Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.展开更多
Numerical simulation of meso-β-scale convective cloud systems associated with a PRE-STORM MCC case has been carried out using a 2-D version of the CSU Regional Atmospheric Modeling System (RAMS) nonhydrostatic model ...Numerical simulation of meso-β-scale convective cloud systems associated with a PRE-STORM MCC case has been carried out using a 2-D version of the CSU Regional Atmospheric Modeling System (RAMS) nonhydrostatic model with parameterized microphysics. It is found that the predicted meso-r-scale convective phenomena arc basically unsteady under the situation of strong shear at low-levels, while the meso-β-scale convective system is maintained up to 3 hours or more. The meso -β- scale cloud system exhibits characteristics of a multi-celled convective storm in which the meso-r-scale convective cells have lifetime of about 30 min. Pressure perturbation depicts a meso-low after a half hour in the low levels. As the cloud system evolves, the meso-low intensifies and extends to the upshear side and covers the entire domain in the mid-lower levels with the peak values of 5-8 hPa. Temperature perturbation depicts a warm region in the middle levels through the entire simulation period. The meso-r-scale warm cores with peak values of 4-8 ℃ are associated with strong convective cells. The cloud top evaporation causes a stronger cold layer around the cloud top levels.Simulation of microphysics exhibits that graupel is primarily concentrated in the strong convective cells forming the main source of convective rainfall after one hour of simulation time. Aggregates are mainly located in the stratiform region and decaying convective cells which produce the stratiform rainfall. Riming of the ice crystals is the predominant precipitation formation mechanism in the convection region, whereas aggregation of ice crystals is the predominant one in the stratiform region, which is consistent with observations. Sensitivity experiments of ice-phase mierophysical processes show that the microphysical structures of the convective cloud system can be simulated better with the diagnosed aggregation collection efficiencies.展开更多
The monopile is the most common foundation to support offshore wind turbines.In the marine environment,local scour due to combined currents and waves is a significant issue that must be considered in the design of win...The monopile is the most common foundation to support offshore wind turbines.In the marine environment,local scour due to combined currents and waves is a significant issue that must be considered in the design of wind turbine foundations.In this paper,a full-scale numerical model was developed and validated based on field data from Rudong,China.The scour development around monopiles was investigated,and the effects of waves and the Reynolds number Re were analyzed.Several formulas for predicting the scour depth in the literature have been evaluated.It is found that waves can accelerate scour development even if the KC number is small(0.78<KC<1.57).The formula obtained from small-scale model tests may be unsafe or wasteful when it is applied in practical design due to the scale effect.A new equation for predicting the scour depth based on the average pile Reynolds number(Rea)is proposed and validated with field data.The equilibrium scour depth predicted using the proposed equation is evaluated and compared with those from nine equations in the literature.It is demonstrated that the values predicted from the proposed equation and from the S/M(Sheppard/Melville)equation are closer to the field data.展开更多
Pore structure of porous media, including pore size and topology, is rather complex. In immiscible twophase displacement process, the capillary force affected by pore size dominates the two-phase flow in the porous me...Pore structure of porous media, including pore size and topology, is rather complex. In immiscible twophase displacement process, the capillary force affected by pore size dominates the two-phase flow in the porous media, affecting displacement results. Direct observation of the flow patterns in the porous media is difficult, and therefore knowledge about the two-phase displacement flow is insufficient. In this paper, a two-dimensional(2D) pore structure was extracted from a sandstone sample, and the flow process that CO_2 displaces resident brine in the extracted pore structure was simulated using the Navier eStokes equation combined with the conservative level set method. The simulation results reveal that the pore throat is a crucial factor for determining CO_2 displacement process in the porous media. The two-phase meniscuses in each pore throat were in a self-adjusting process. In the displacement process,CO_2 preferentially broke through the maximum pore throat. Before breaking through the maximum pore throat, the pressure of CO_2 continually increased, and the curvature and position of two-phase interfaces in the other pore throats adjusted accordingly. Once the maximum pore throat was broken through by the CO_2, the capillary force in the other pore throats released accordingly; subsequently, the interfaces withdrew under the effect of capillary fore, preparing for breaking through the next pore throat.Therefore, the two-phase displacement in CO_2 injection is accompanied by the breaking through and adjusting of the two-phase interfaces.展开更多
A solution scheme is proposed in this paper for an existing RTDHT system to simulate large-scale finite element （FE） numerical substructures. The analysis of the FE numerical substructure is split into response anal...A solution scheme is proposed in this paper for an existing RTDHT system to simulate large-scale finite element （FE） numerical substructures. The analysis of the FE numerical substructure is split into response analysis and signal generation tasks, and executed in two different target computers in real-time. One target computer implements the response analysis task, wherein a large time-step is used to solve the FE substructure, and another target computer implements the signal generation task, wherein an interpolation program is used to generate control signals in a small time-step to meet the input demand of the controller. By using this strategy, the scale of the FE numerical substructure simulation may be increased significantly. The proposed scheme is initially verified by two FE numerical substructure models with 98 and 1240 degrees of freedom （DOFs）. Thereafter, RTDHTs of a single frame-foundation structure are implemented where the foundation, considered as the numerical substructure, is simulated by the FE model with 1240 DOFs. Good agreements between the results of the RTDHT and those from the FE analysis in ABAQUS are obtained.展开更多
The Regional Atmospheric Modeling System （RAMS） and the computational fluid dynamics （CFD） codes known as FLUENT are combinatorially applied in a multi-scale numerical simulation of the urban surface layer （USL）...The Regional Atmospheric Modeling System （RAMS） and the computational fluid dynamics （CFD） codes known as FLUENT are combinatorially applied in a multi-scale numerical simulation of the urban surface layer （USL）. RAMS and FLUENT are combined as a multi-scale numerical modeling system, in which the RAMS simulated data are delivered to the computational model for FLUENT simulation in an offline way. Numerical simulations are performed to present and preliminarily validate the capability of the multi-scale modeling system, and the results show that the modeling system can reasonably provide information on the meteorological elements in an urban area from the urban scale to the city-block scale, especially the details of the turbulent flows within the USL.展开更多
Impacts of regional sea surface temperature(SST)anomalies on the interdecadal variation of the cross-equatorial flows(CEFs)in Eastern Hemisphere are studied using numerical simulations with a global atmospheric circul...Impacts of regional sea surface temperature(SST)anomalies on the interdecadal variation of the cross-equatorial flows(CEFs)in Eastern Hemisphere are studied using numerical simulations with a global atmospheric circulation model(NCAR CAM3)driven with 1950-2000 monthly SSTs in different marine areas(the globe,extratropics,tropics,tropical Indian Ocean-Pacific,and tropical Pacific)and ERA-40reanalysis data.Results show that all simulations,except the one driven with extratropical SSTs,can simulate the interdecadal strengthening of CEFs around Somali,120oE,and 150oE that occurred in the midand late-1970s.Among those simulated CEFs,the interdecadal variability in Somali and its interdecadal relationship with the East Asian summer monsoon are in better agreement with the observations,suggesting that changes in the SSTs of tropical oceans,especially the tropical Pacific,play a crucial role in the interdecadal variability of CEFs in Somali.The interdecadal change of CEFs in Somali is highly associated with the interdecadal variation of tropical Pacific SST.As the interdecadal warmer(colder)SST happens in the tropical Pacific,a"sandwich"pattern of SST anomalies,i.e."+,-,+"("-,+,-"),will occur in the eastern tropical Pacific from north to south with a pair of anomalous anticyclone(cyclone)at the lower troposphere;the pair links to another pair of anomalous cyclone(anticyclone)in the tropical Indian Ocean through an atmospheric bridge,and thus strengthens(weakens)the CEFs in Somali.展开更多
A high resolution, nonhydrostatic, three dimensional diagnostic PBL model over small scale concave terrain was established in this paper. A two dimensional prognostic model was developed based on the diagnostic mo...A high resolution, nonhydrostatic, three dimensional diagnostic PBL model over small scale concave terrain was established in this paper. A two dimensional prognostic model was developed based on the diagnostic model. The hydrostatic approximation was abandoned and the simple energy ( E e ) closure scheme was used in both models. Using the two models, characteristics of PBL structure and its evolution were fully studied. The main characteristic of the PBL is the circulation, and it fairly affects the distribution of the pollutant in the pit.展开更多
A spectrum method is used to simulate the time-developing free mixing layerwith cross shear which is introduced in different stages. The results show that the properties of flow are nearly the same for situations whet...A spectrum method is used to simulate the time-developing free mixing layerwith cross shear which is introduced in different stages. The results show that the properties of flow are nearly the same for situations whether the cross shear is introduced in theinitial time or in early stage. If cross shear is introduced in the stage that the roll-up ofmixing layer occurs, the turbulent intensities of now will increase and mixture of now willbe enhanced.展开更多
Thin film is a widely used structure in the present microelectromechanical systems (MEMS) and plays a vital role in many functional devices. However, the great size difference between the film's thickness and its ...Thin film is a widely used structure in the present microelectromechanical systems (MEMS) and plays a vital role in many functional devices. However, the great size difference between the film's thickness and its planar dimensions makes it difficult to study the thin film performance numerically. In this work, a scaling transformation was presented to make the different dimensional sizes equivalent, and thereby, to improve the grid quality considerably. Two numerical experiments were studied to validate the present scaling transformation method. The numerical results indicated that the largest grid size difference can be decreased to one to two orders of magnitude by using the present scaling transformation, and the memory required by the numerical simulation, i.e., the total grid number, could be reduced by about two to three orders of magnitude, while the numerical accuracies with and without this scaling transformation were nearly the same.展开更多
Cloud micro-physical structures in a precipitation system associated with the Meiyu front are observed using the balloon-borne Precipitation Particle Image Sensor at Baoshan observatory station, Shanghai during June a...Cloud micro-physical structures in a precipitation system associated with the Meiyu front are observed using the balloon-borne Precipitation Particle Image Sensor at Baoshan observatory station, Shanghai during June and July 1999. The vertical distributions of various cloud particle size, number density, and mass density are retrieved from the observations. Analyses of observations show that ice-phase particles (ice crystals, graupel, snowflakes, and frozen drops) often exist in the cloud of torrential rain associated with the Meiyu front. Among the various particles, ice crystals and graupel are the most numerous, but graupel and snow have the highest mass density. Ice-phase particles coexist with liquid water droplets near the 0°C level. The graupel is similarly distributed with height as the ice crystals. Raindrops below the 0°C level are mainly from melted grauple, snowflakes and frozen drops. They may further grow larger by coalescence with smaller ones as they fall from the cloud base. Numerical simulations using the non-hydrostatic meso-scale model MM5 with the Reisner graupel explicit moisture scheme confirm the main observational results. Rain water at the lower level is mainly generated from the melting of snow and graupel falling from the upper level where snow and graupel are generated and grown from collection with cloud and rain water. Thus the mixed-phase cloud process, in which ice phase coexists and interacts with liquid phase (cloud and rain drops), plays the most important role in the formation and development of heavy convective rainfall in the Meiyu frontal system.展开更多
Large eddy simulation (LES) explicitly calculates the large-scale vortex field and parameterizes the small-scale vortices.In this study,LES and κ-ε models were developed for a specific geometrical configuration of b...Large eddy simulation (LES) explicitly calculates the large-scale vortex field and parameterizes the small-scale vortices.In this study,LES and κ-ε models were developed for a specific geometrical configuration of backward-facing step (BFS).The simulation results were validated with particle image velocimetry (PIV) measurements and direct numerical simulation (DNS).This LES simulation was carried out with a Reynolds number of 9000 in a pressurized water tunnel with an expansion ratio of 2.00.The results indicate that the LES model can reveal largescale vortex motion although with a larger grid-cell size.However,the LES model tends to overestimate the top wall separation and the Reynolds stress components for the BFS flow simulation without a sufficiently fine grid.Overall,LES is a potential tool for simulating separated flow controlled by large-scale vortices.展开更多
A multi-scale continuous-discrete model based on the effects of the p27 gene control is built to simulate the avascular tumor growth. At the tissue level, the continuous Eulerian model is adopted to determine the dist...A multi-scale continuous-discrete model based on the effects of the p27 gene control is built to simulate the avascular tumor growth. At the tissue level, the continuous Eulerian model is adopted to determine the distribution of the concentration of oxygen, the extracellular matrix （ECM）, and the matrix-degradative enzyme （MDE）. At the cellular level, the discrete Lagrangien model is adopted to determine the movement, the proliferation, and the death of single tumor cells （TCs）. At the genetic level, whether a cell is committed to mitosis is determined by solving a set of equations modeling the effects of the p27 gene control. The avascular morphological evolution of the solid tumor growth is simulated, including the radius the oxygen distribution over time, and the expression. of the solid tumor, the number of the TCs, inhibiting effect＇ of the up-regulating p27 gene展开更多
In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommend...In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommended to determine the blast load on the structures via numerical simulation. Since the mesh size of the numerical model highly depends on the explosion scenario, there is no generally applicable approach for the mesh size selection. An efficient method to determine the mesh size of the numerical model of near-ground detonation based on explosion scenarios is proposed in this study. The effect of mesh size on the propagation of blast wave under different explosive weights was studied, and the correlations between the mesh size effect and the charge weight or the scaled distance was described. Based on the principle of the finite element method and Hopkinson-Cranz scaling law, a mesh size measurement unit related to the explosive weight was proposed as the criterion for determining the mesh size in the numerical simulation. Finally, the applicability of the method proposed in this paper was verified by comparing the results from numerical simulation and the explosion tests and was verified in AUTODYN.展开更多
Scale effect of ISWs loads on Floating Production Storage and Offloading(FPSO) is studied in this paper. The application conditions of KdV, eKdV and MCC ISWs theories are used in the numerical method. The depthaverage...Scale effect of ISWs loads on Floating Production Storage and Offloading(FPSO) is studied in this paper. The application conditions of KdV, eKdV and MCC ISWs theories are used in the numerical method. The depthaveraged velocities induced by ISWs are used for the velocity-inlet boundary. Three scale ratio numerical models λ=1, 20 and 300 were selected, which the scale ratio is the size ratio of numerical models to the experimental model.The comparisons between the numerical and former experimental results are performed to verify the feasibility of numerical method. The comparisons between the numerical and simplified theoretical results are performed to discuss the applicability of the simplified theoretical model summarized from the load experiments. Firstly, the numerical results of λ=1 numerical model showed a good agreement with former experimental and simplified theoretical results. It is feasible to simulate the ISWs loads on FPSO by the numerical method. Secondly, the comparisons between the results of three scale ratio numerical models and experimental results indicated that the scale ratios have more significant influence on the experimental horizontal forces than the vertical forces. The scale effect of horizontal forces mainly results from the different viscosity effects associated with the model’s dimension.Finally, through the comparisons between the numerical and simplified theoretical results for three scale ratio models, the simplified theoretical model of the pressure difference and friction forces exerted by ISWs on FPSO is applied for large-scale or full-scale FPSO.展开更多
A 3D numerical model considering the soil-structure interaction is presented in this paper to examine the ground movement and internal force during the construction of Qingdao North Metro Station, China with a special...A 3D numerical model considering the soil-structure interaction is presented in this paper to examine the ground movement and internal force during the construction of Qingdao North Metro Station, China with a special focus on the convex effect of the crossing excavation. The influence of intersection angles and soil resilience characteristics on deformation behavior is discussed, and the suitability of two alternative constitutive models applied in excavation simulation is also considered. The analysis results show that a notable convex effect appears to be associated with the crossing excavation, and the intersection is the key area requiring special attention. The displacements at the comer decrease with increasing crossing angles. The axial loads of struts along the retaining pile wall are unequal, and the values near the cross section are generally larger than the average loads of the left-sided ones. The modified Cam-Clay （MCC） model, which is capable of describing the loading-unloading criterion and identifying the stiffness difference of strain hardening between loading and unloading, can yield a relatively high accuracy of estimation for the behavior of excavations in comparison to the Mohr-Coulomb （MC） model. Furthermore, slight soil deformation resilience after unloading can reduce the ground surface settlement and enhance the ground stability.展开更多
Crystal plasticity theory was used to simulate upsetting tests of different dimensions and grain size micro copper cylinders in this study on the fluctuant flow stress scale effect. Results showed that with the decrea...Crystal plasticity theory was used to simulate upsetting tests of different dimensions and grain size micro copper cylinders in this study on the fluctuant flow stress scale effect. Results showed that with the decrease of billet grain quantity, flow stress fluctuation is not always increased, but there is a maximum. Through this study, the fluctuant flow stress scale effect can be understood deeper, and relevant necessary information was obtained for further prediction and control of this scale effect and to design the microforming process and die.展开更多
The effect of mandrel with the structure of ball and socket on the cross section quality of thin-walled tube numerical controlled(NC) bending was studied by numerical simulation method, combined with theoretical (anal...The effect of mandrel with the structure of ball and socket on the cross section quality of thin-walled tube numerical controlled(NC) bending was studied by numerical simulation method, combined with theoretical (analysis) and experiment. Influencing factors of the mandrel include the count of mandrel heads, the diameter of mandrel and its position. According to the principle of NC tube bending, quality defects possibly produced in thin-walled tube NC bending process were analyzed and two parameters were proposed in order to describe the cross section quality of thin-walled tube NC bending. According to the geometrical dimension of tube and dies, the range of mandrel protrusion was derived. The finite element model of thin-walled tube NC bending was established based on the DYNAFORM platform, and key technological problems were solved. The model was verified by experiment. The effect of the number of mandrel heads, the diameter of mandrel and the protrusion length of mandrel on the cross section quality of thin-walled tube NC bending was revealed and how to choose mandrel parameters was presented.展开更多
In modem manufacturing, a new type of sheet metal part with step cross-section in both inner hole and outer edge is proposed. The traditional stamping separating processes can only produce sheet metal part with vertic...In modem manufacturing, a new type of sheet metal part with step cross-section in both inner hole and outer edge is proposed. The traditional stamping separating processes can only produce sheet metal part with vertical cross-section. According to the latest developing theory and potential of cold pressure forming： combination of pressure and cold forging, a new flow control forming of sheet metal（FCF） is excogitated based on blanking process of general stamping and combined with cold forging processes such as extrusion and coining, etc, which is aiming at the above-mentioned new type of sheet metal part. With utilization of this new process, the new type of sheet metal parts can be manufactured. In order to shorten the testing period, the numerical simulation was carried out by using DEFORM-3D software, and both deformation and mechanics rules were analyzed. Based on the simulation, both punching part and blanked parts of this new type were successfully developed. Then a new conception of optimal distance between the step walls of inner hole and outside edge was proposed and the design principle for its numerical value was inferred. Furthermore, a mold set for combination of stamping ＆ cold forging was designed and manufactured, by which the technologic experiments were taken for validation with Aluminum plate of thickness 2.35 mm for power battery cover board, which verified the principle of the distance between the step walls. The research of cold pressure forming of thin sheet metal with step cross-section is significant, not only to the development of modem mechanical manufacture, but also to metal plastic forming science.展开更多
文摘Existing numerical methods for complex composites, such as multiscale simulation and neural network algorithms, face significant limitations. Multiscale techniques are often prohibitively expensive for large models, while neural networks struggle to represent underlying microscopic material properties. To overcome these challenges, a meso-micro scale numerical method using a virtual node approach is developed in this study. A Wbraid/Al/Epoxy functional structural material is fabricated, and a representative periodic unit cell is identified based on its architecture. The complex structure is then discretized into nodes, and mechanical interactions are governed by pre-defined computation rules. This virtual node method is systematically compared against both multiscale simulation and a neural network algorithm, with validation provided through mechanical experiments. The results demonstrate that the nodal operation strategy significantly reduces computational resource requirements. By quantifying microscopic bonding with coefficients, explicit interface treatment is avoided, granting the method strong adaptability to lattice materials. The method can simulate extremely complex structures using parameters from simple tests and is suited for large systems. Compared to three-point bending experiments, errors for multiscale, virtual node, and neural network methods were 12.4%, 6.9%, and 34.5%, respectively. Under dynamic compression, the errors were 2.7%, 9.3%, and 15.43%. The virtual node method demonstrated superior accuracy under static conditions, enabling efficient prediction and auxiliary development of complex structural materials.
文摘Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.
文摘Numerical simulation of meso-β-scale convective cloud systems associated with a PRE-STORM MCC case has been carried out using a 2-D version of the CSU Regional Atmospheric Modeling System (RAMS) nonhydrostatic model with parameterized microphysics. It is found that the predicted meso-r-scale convective phenomena arc basically unsteady under the situation of strong shear at low-levels, while the meso-β-scale convective system is maintained up to 3 hours or more. The meso -β- scale cloud system exhibits characteristics of a multi-celled convective storm in which the meso-r-scale convective cells have lifetime of about 30 min. Pressure perturbation depicts a meso-low after a half hour in the low levels. As the cloud system evolves, the meso-low intensifies and extends to the upshear side and covers the entire domain in the mid-lower levels with the peak values of 5-8 hPa. Temperature perturbation depicts a warm region in the middle levels through the entire simulation period. The meso-r-scale warm cores with peak values of 4-8 ℃ are associated with strong convective cells. The cloud top evaporation causes a stronger cold layer around the cloud top levels.Simulation of microphysics exhibits that graupel is primarily concentrated in the strong convective cells forming the main source of convective rainfall after one hour of simulation time. Aggregates are mainly located in the stratiform region and decaying convective cells which produce the stratiform rainfall. Riming of the ice crystals is the predominant precipitation formation mechanism in the convection region, whereas aggregation of ice crystals is the predominant one in the stratiform region, which is consistent with observations. Sensitivity experiments of ice-phase mierophysical processes show that the microphysical structures of the convective cloud system can be simulated better with the diagnosed aggregation collection efficiencies.
基金financially supported by the National Natural Science Foundation of China (Grant No.52378329)。
文摘The monopile is the most common foundation to support offshore wind turbines.In the marine environment,local scour due to combined currents and waves is a significant issue that must be considered in the design of wind turbine foundations.In this paper,a full-scale numerical model was developed and validated based on field data from Rudong,China.The scour development around monopiles was investigated,and the effects of waves and the Reynolds number Re were analyzed.Several formulas for predicting the scour depth in the literature have been evaluated.It is found that waves can accelerate scour development even if the KC number is small(0.78<KC<1.57).The formula obtained from small-scale model tests may be unsafe or wasteful when it is applied in practical design due to the scale effect.A new equation for predicting the scour depth based on the average pile Reynolds number(Rea)is proposed and validated with field data.The equilibrium scour depth predicted using the proposed equation is evaluated and compared with those from nine equations in the literature.It is demonstrated that the values predicted from the proposed equation and from the S/M(Sheppard/Melville)equation are closer to the field data.
基金funded by Key Laboratory of Coal-based CO_2 Capture and Geological Storage,Jiangsu Province,ChinaUS Advanced Coal Technology Consortium(No.2013 DFB60140-08)
文摘Pore structure of porous media, including pore size and topology, is rather complex. In immiscible twophase displacement process, the capillary force affected by pore size dominates the two-phase flow in the porous media, affecting displacement results. Direct observation of the flow patterns in the porous media is difficult, and therefore knowledge about the two-phase displacement flow is insufficient. In this paper, a two-dimensional(2D) pore structure was extracted from a sandstone sample, and the flow process that CO_2 displaces resident brine in the extracted pore structure was simulated using the Navier eStokes equation combined with the conservative level set method. The simulation results reveal that the pore throat is a crucial factor for determining CO_2 displacement process in the porous media. The two-phase meniscuses in each pore throat were in a self-adjusting process. In the displacement process,CO_2 preferentially broke through the maximum pore throat. Before breaking through the maximum pore throat, the pressure of CO_2 continually increased, and the curvature and position of two-phase interfaces in the other pore throats adjusted accordingly. Once the maximum pore throat was broken through by the CO_2, the capillary force in the other pore throats released accordingly; subsequently, the interfaces withdrew under the effect of capillary fore, preparing for breaking through the next pore throat.Therefore, the two-phase displacement in CO_2 injection is accompanied by the breaking through and adjusting of the two-phase interfaces.
基金National Natural Science Foundation under Grant Nos.51179093,91215301 and 41274106the Specialized Research Fund for the Doctoral Program of Higher Education under Grant No.20130002110032Tsinghua University Initiative Scientific Research Program under Grant No.20131089285
文摘A solution scheme is proposed in this paper for an existing RTDHT system to simulate large-scale finite element （FE） numerical substructures. The analysis of the FE numerical substructure is split into response analysis and signal generation tasks, and executed in two different target computers in real-time. One target computer implements the response analysis task, wherein a large time-step is used to solve the FE substructure, and another target computer implements the signal generation task, wherein an interpolation program is used to generate control signals in a small time-step to meet the input demand of the controller. By using this strategy, the scale of the FE numerical substructure simulation may be increased significantly. The proposed scheme is initially verified by two FE numerical substructure models with 98 and 1240 degrees of freedom （DOFs）. Thereafter, RTDHTs of a single frame-foundation structure are implemented where the foundation, considered as the numerical substructure, is simulated by the FE model with 1240 DOFs. Good agreements between the results of the RTDHT and those from the FE analysis in ABAQUS are obtained.
基金This study was supported by the National Natural Science Foundation of China （Grant Nos. 40233030, 40405004, 40405014）.
文摘The Regional Atmospheric Modeling System （RAMS） and the computational fluid dynamics （CFD） codes known as FLUENT are combinatorially applied in a multi-scale numerical simulation of the urban surface layer （USL）. RAMS and FLUENT are combined as a multi-scale numerical modeling system, in which the RAMS simulated data are delivered to the computational model for FLUENT simulation in an offline way. Numerical simulations are performed to present and preliminarily validate the capability of the multi-scale modeling system, and the results show that the modeling system can reasonably provide information on the meteorological elements in an urban area from the urban scale to the city-block scale, especially the details of the turbulent flows within the USL.
基金National Basic Research Program of China(2013CB4302022012CB955200)+2 种基金Special Program for China Meteorology Trade(GYHY201306020)Basic Natural Science Research of Jiangsu High Institutions(12KJB170007,10KJB170007)Priority Academic Program Development of Jiangsu Higher Education Institutions(PAPD)
文摘Impacts of regional sea surface temperature(SST)anomalies on the interdecadal variation of the cross-equatorial flows(CEFs)in Eastern Hemisphere are studied using numerical simulations with a global atmospheric circulation model(NCAR CAM3)driven with 1950-2000 monthly SSTs in different marine areas(the globe,extratropics,tropics,tropical Indian Ocean-Pacific,and tropical Pacific)and ERA-40reanalysis data.Results show that all simulations,except the one driven with extratropical SSTs,can simulate the interdecadal strengthening of CEFs around Somali,120oE,and 150oE that occurred in the midand late-1970s.Among those simulated CEFs,the interdecadal variability in Somali and its interdecadal relationship with the East Asian summer monsoon are in better agreement with the observations,suggesting that changes in the SSTs of tropical oceans,especially the tropical Pacific,play a crucial role in the interdecadal variability of CEFs in Somali.The interdecadal change of CEFs in Somali is highly associated with the interdecadal variation of tropical Pacific SST.As the interdecadal warmer(colder)SST happens in the tropical Pacific,a"sandwich"pattern of SST anomalies,i.e."+,-,+"("-,+,-"),will occur in the eastern tropical Pacific from north to south with a pair of anomalous anticyclone(cyclone)at the lower troposphere;the pair links to another pair of anomalous cyclone(anticyclone)in the tropical Indian Ocean through an atmospheric bridge,and thus strengthens(weakens)the CEFs in Somali.
文摘A high resolution, nonhydrostatic, three dimensional diagnostic PBL model over small scale concave terrain was established in this paper. A two dimensional prognostic model was developed based on the diagnostic model. The hydrostatic approximation was abandoned and the simple energy ( E e ) closure scheme was used in both models. Using the two models, characteristics of PBL structure and its evolution were fully studied. The main characteristic of the PBL is the circulation, and it fairly affects the distribution of the pollutant in the pit.
文摘A spectrum method is used to simulate the time-developing free mixing layerwith cross shear which is introduced in different stages. The results show that the properties of flow are nearly the same for situations whether the cross shear is introduced in theinitial time or in early stage. If cross shear is introduced in the stage that the roll-up ofmixing layer occurs, the turbulent intensities of now will increase and mixture of now willbe enhanced.
基金National Natural Science Foundation of China(No.60576020,No.60606014).
文摘Thin film is a widely used structure in the present microelectromechanical systems (MEMS) and plays a vital role in many functional devices. However, the great size difference between the film's thickness and its planar dimensions makes it difficult to study the thin film performance numerically. In this work, a scaling transformation was presented to make the different dimensional sizes equivalent, and thereby, to improve the grid quality considerably. Two numerical experiments were studied to validate the present scaling transformation method. The numerical results indicated that the largest grid size difference can be decreased to one to two orders of magnitude by using the present scaling transformation, and the memory required by the numerical simulation, i.e., the total grid number, could be reduced by about two to three orders of magnitude, while the numerical accuracies with and without this scaling transformation were nearly the same.
基金This study was supported by the State Key Basic Program:Research on the Formation Mechanism and Prediction Theory of Severe Synoptic Dis- asters in China No.G1998040907 and the National NaturalSciences Foundation of China under Grant No.49735180
文摘Cloud micro-physical structures in a precipitation system associated with the Meiyu front are observed using the balloon-borne Precipitation Particle Image Sensor at Baoshan observatory station, Shanghai during June and July 1999. The vertical distributions of various cloud particle size, number density, and mass density are retrieved from the observations. Analyses of observations show that ice-phase particles (ice crystals, graupel, snowflakes, and frozen drops) often exist in the cloud of torrential rain associated with the Meiyu front. Among the various particles, ice crystals and graupel are the most numerous, but graupel and snow have the highest mass density. Ice-phase particles coexist with liquid water droplets near the 0°C level. The graupel is similarly distributed with height as the ice crystals. Raindrops below the 0°C level are mainly from melted grauple, snowflakes and frozen drops. They may further grow larger by coalescence with smaller ones as they fall from the cloud base. Numerical simulations using the non-hydrostatic meso-scale model MM5 with the Reisner graupel explicit moisture scheme confirm the main observational results. Rain water at the lower level is mainly generated from the melting of snow and graupel falling from the upper level where snow and graupel are generated and grown from collection with cloud and rain water. Thus the mixed-phase cloud process, in which ice phase coexists and interacts with liquid phase (cloud and rain drops), plays the most important role in the formation and development of heavy convective rainfall in the Meiyu frontal system.
基金supported by the National Natural Science Foundation of China(Grant No.51379128)
文摘Large eddy simulation (LES) explicitly calculates the large-scale vortex field and parameterizes the small-scale vortices.In this study,LES and κ-ε models were developed for a specific geometrical configuration of backward-facing step (BFS).The simulation results were validated with particle image velocimetry (PIV) measurements and direct numerical simulation (DNS).This LES simulation was carried out with a Reynolds number of 9000 in a pressurized water tunnel with an expansion ratio of 2.00.The results indicate that the LES model can reveal largescale vortex motion although with a larger grid-cell size.However,the LES model tends to overestimate the top wall separation and the Reynolds stress components for the BFS flow simulation without a sufficiently fine grid.Overall,LES is a potential tool for simulating separated flow controlled by large-scale vortices.
基金Project supported by the National Natural Science Foundation of China （Nos. 10372026 and 10772751）
文摘A multi-scale continuous-discrete model based on the effects of the p27 gene control is built to simulate the avascular tumor growth. At the tissue level, the continuous Eulerian model is adopted to determine the distribution of the concentration of oxygen, the extracellular matrix （ECM）, and the matrix-degradative enzyme （MDE）. At the cellular level, the discrete Lagrangien model is adopted to determine the movement, the proliferation, and the death of single tumor cells （TCs）. At the genetic level, whether a cell is committed to mitosis is determined by solving a set of equations modeling the effects of the p27 gene control. The avascular morphological evolution of the solid tumor growth is simulated, including the radius the oxygen distribution over time, and the expression. of the solid tumor, the number of the TCs, inhibiting effect＇ of the up-regulating p27 gene
基金the funding supports of the National Key Research and Development Plan,China(Grant No.2022YFC3801800)National Natural Science Foundation of China(Grant Nos.52038010 and 52078368)。
文摘In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommended to determine the blast load on the structures via numerical simulation. Since the mesh size of the numerical model highly depends on the explosion scenario, there is no generally applicable approach for the mesh size selection. An efficient method to determine the mesh size of the numerical model of near-ground detonation based on explosion scenarios is proposed in this study. The effect of mesh size on the propagation of blast wave under different explosive weights was studied, and the correlations between the mesh size effect and the charge weight or the scaled distance was described. Based on the principle of the finite element method and Hopkinson-Cranz scaling law, a mesh size measurement unit related to the explosive weight was proposed as the criterion for determining the mesh size in the numerical simulation. Finally, the applicability of the method proposed in this paper was verified by comparing the results from numerical simulation and the explosion tests and was verified in AUTODYN.
基金financially supported by the National Natural Science Foundation of China(Grant No.11372184)the National Basic Research Program of China(973 Program,Grant Nos.2015CB251203-3 and 2013CB036103)
文摘Scale effect of ISWs loads on Floating Production Storage and Offloading(FPSO) is studied in this paper. The application conditions of KdV, eKdV and MCC ISWs theories are used in the numerical method. The depthaveraged velocities induced by ISWs are used for the velocity-inlet boundary. Three scale ratio numerical models λ=1, 20 and 300 were selected, which the scale ratio is the size ratio of numerical models to the experimental model.The comparisons between the numerical and former experimental results are performed to verify the feasibility of numerical method. The comparisons between the numerical and simplified theoretical results are performed to discuss the applicability of the simplified theoretical model summarized from the load experiments. Firstly, the numerical results of λ=1 numerical model showed a good agreement with former experimental and simplified theoretical results. It is feasible to simulate the ISWs loads on FPSO by the numerical method. Secondly, the comparisons between the results of three scale ratio numerical models and experimental results indicated that the scale ratios have more significant influence on the experimental horizontal forces than the vertical forces. The scale effect of horizontal forces mainly results from the different viscosity effects associated with the model’s dimension.Finally, through the comparisons between the numerical and simplified theoretical results for three scale ratio models, the simplified theoretical model of the pressure difference and friction forces exerted by ISWs on FPSO is applied for large-scale or full-scale FPSO.
基金Project (No. 2007AA11Z134) supported by the National High-Tech R&D Program (863) of China
文摘A 3D numerical model considering the soil-structure interaction is presented in this paper to examine the ground movement and internal force during the construction of Qingdao North Metro Station, China with a special focus on the convex effect of the crossing excavation. The influence of intersection angles and soil resilience characteristics on deformation behavior is discussed, and the suitability of two alternative constitutive models applied in excavation simulation is also considered. The analysis results show that a notable convex effect appears to be associated with the crossing excavation, and the intersection is the key area requiring special attention. The displacements at the comer decrease with increasing crossing angles. The axial loads of struts along the retaining pile wall are unequal, and the values near the cross section are generally larger than the average loads of the left-sided ones. The modified Cam-Clay （MCC） model, which is capable of describing the loading-unloading criterion and identifying the stiffness difference of strain hardening between loading and unloading, can yield a relatively high accuracy of estimation for the behavior of excavations in comparison to the Mohr-Coulomb （MC） model. Furthermore, slight soil deformation resilience after unloading can reduce the ground surface settlement and enhance the ground stability.
文摘Crystal plasticity theory was used to simulate upsetting tests of different dimensions and grain size micro copper cylinders in this study on the fluctuant flow stress scale effect. Results showed that with the decrease of billet grain quantity, flow stress fluctuation is not always increased, but there is a maximum. Through this study, the fluctuant flow stress scale effect can be understood deeper, and relevant necessary information was obtained for further prediction and control of this scale effect and to design the microforming process and die.
文摘The effect of mandrel with the structure of ball and socket on the cross section quality of thin-walled tube numerical controlled(NC) bending was studied by numerical simulation method, combined with theoretical (analysis) and experiment. Influencing factors of the mandrel include the count of mandrel heads, the diameter of mandrel and its position. According to the principle of NC tube bending, quality defects possibly produced in thin-walled tube NC bending process were analyzed and two parameters were proposed in order to describe the cross section quality of thin-walled tube NC bending. According to the geometrical dimension of tube and dies, the range of mandrel protrusion was derived. The finite element model of thin-walled tube NC bending was established based on the DYNAFORM platform, and key technological problems were solved. The model was verified by experiment. The effect of the number of mandrel heads, the diameter of mandrel and the protrusion length of mandrel on the cross section quality of thin-walled tube NC bending was revealed and how to choose mandrel parameters was presented.
文摘In modem manufacturing, a new type of sheet metal part with step cross-section in both inner hole and outer edge is proposed. The traditional stamping separating processes can only produce sheet metal part with vertical cross-section. According to the latest developing theory and potential of cold pressure forming： combination of pressure and cold forging, a new flow control forming of sheet metal（FCF） is excogitated based on blanking process of general stamping and combined with cold forging processes such as extrusion and coining, etc, which is aiming at the above-mentioned new type of sheet metal part. With utilization of this new process, the new type of sheet metal parts can be manufactured. In order to shorten the testing period, the numerical simulation was carried out by using DEFORM-3D software, and both deformation and mechanics rules were analyzed. Based on the simulation, both punching part and blanked parts of this new type were successfully developed. Then a new conception of optimal distance between the step walls of inner hole and outside edge was proposed and the design principle for its numerical value was inferred. Furthermore, a mold set for combination of stamping ＆ cold forging was designed and manufactured, by which the technologic experiments were taken for validation with Aluminum plate of thickness 2.35 mm for power battery cover board, which verified the principle of the distance between the step walls. The research of cold pressure forming of thin sheet metal with step cross-section is significant, not only to the development of modem mechanical manufacture, but also to metal plastic forming science.
