In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over ...In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over large-scale and complex terrain requires additional methods. To this end, we have proposed a new topography-capable By performing a two-dimensional Fourier transform in the horizontal directions, threedimensional partial differential equations in the spatial domain were transformed into a group of independent, one-dimensional differential equations engaged with different wave numbers. These independent differential equations are highly parallel across different wave numbers. differential equations with different wave numbers, and the efficiency of solving fixedbandwidth linear equations was further improved by a chasing method. In a synthetic test, a prism model was used to verify the accuracy and reliability of the proposed algorithm by comparing the numerical solution with the analytical solution. We studied the computational precision and efficiency with and without topography using different Fourier transform methods. The results showed that the Guass-FFT method has higher numerical precision, while the standard FFT method is superior, in terms of computation time, for inversion and quantitative interpretation under complicated terrain.展开更多
Applying the methods of on-site observation and dynamic model, the research on the fronts at the Jiulong Estuary has been carried out, during which spatial and temporal distribution, dynamic characteristics and format...Applying the methods of on-site observation and dynamic model, the research on the fronts at the Jiulong Estuary has been carried out, during which spatial and temporal distribution, dynamic characteristics and formation mechanism of salinity fronts are analyzed and discussed. The research shows that the estuarine fronts mainly lie in the area from the Jiyu Islet to the Haimen Island, outside of Yuweizai to Hulishan cross-section, the near coast of Yuweizai and the south of the Songyu-Gulangyu Channel. The fronts in the former two regions are formed directly by plume, while the one near the coast of Yuweizai is a tidal intrusion front caused by flood current and the one at the south of the Songyu-Gulangyu Channel is the result of current shear transformation. Under normal circumstances, fresh water of the Jiulong River mainly influences the inside of the Xiamen Bay, and when it is in typhoon seasons, plume front can affect the Taiwan Strait and has an effect on the biogeochemical Drocesses in the strait.展开更多
The characteristics of three-dimensional （3-D） tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, ma...The characteristics of three-dimensional （3-D） tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, σ 3-D numerical model oil an unstructured triangular grid has been degeloped. The σ coordinate transforination, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simulations in the Oujing Estuary with submerged dyke and moveable boundary problems. The model is evaluated with the in situ data, and the results show that the calculated water elevations at 19 stations and currents at 19 profiler stations are in good agreement with measured data both in magnitude and phase. This numerical model is applied to the 3-D tidal circulation simulations of experiments in stopping flow transport through the South Branch of the Oujiang Estuary, and the feasibility to cutoff the flow in the South Branch of the Oujiang Estuary is demonstrated by numerical simulation experiments. The developed numerical model simulated the 3-D tidal current circulations in complicated coastal and estuarine waters very well.展开更多
Saline intrusion into marine sewage ouffalls will greatly decrease the efficiency of sewage disposal. In order to investigate the mechanisms of this flow, in this paper, a three-dimensional numerical model based on F...Saline intrusion into marine sewage ouffalls will greatly decrease the efficiency of sewage disposal. In order to investigate the mechanisms of this flow, in this paper, a three-dimensional numerical model based on FVM （Finite Volume Method） is established, The RNG κ-ε model is selected for turbulence modeling. The time-averaged vohtme fraction equations are introduced to simulate the stratification and inteffaeial exchange of sewage and seawater in outfalls. Validity of the established three-dimensional numerical model is evaluated by comparisons of numerical results with experimental data. With this three-dimensional numerical model, the internal flow characteristics in ouffalls for different sewage discharges are simulated. The results indicate that for a low sewage discharge, saline circulates in the outfall due to intrusion and both the inflowing momentum and the inteffaeial turbulent mixing are important mechanisms to extrude the saline. For a high sewage discharge, saline intrusion could be avoided. The inflow momentum is the main mechanism to extrude the saline and the inteffacial turbulent mixing is nut important relatively. Even at a high sewage discharge, the saline wedge would be retained in the main ouffall pipe after the risers are purged. It takes a long time for this saline wedge to be extruded by interracial turbulent mixing.展开更多
In order to make the numerical calculation of viscous flows more convenient for the flows in channel with complicated profile governing equations expressed in the arbitrary curvilinear coordinates were derived by mean...In order to make the numerical calculation of viscous flows more convenient for the flows in channel with complicated profile governing equations expressed in the arbitrary curvilinear coordinates were derived by means of Favre density-weighted averaged method, and a turbulent model with effect of curvature modification was also derived. The numerical calculation of laminar and turbulent flown in divergent curved channels was carried out by means of parabolizeil computation method. The calculating results were used to analyze and investigate the aerodynamic performance of talor cascades in compressors preliminarily.展开更多
The motion of gas bubbles beneath a free surface will lead to a spike of fluid on the free surface. The distance of the bubbles to the free surface is the key factor to different phenomena. When the inception distance...The motion of gas bubbles beneath a free surface will lead to a spike of fluid on the free surface. The distance of the bubbles to the free surface is the key factor to different phenomena. When the inception distance varies in some range, crown phenomenon would happen after the impact of weak buoyancy bubbles, so this kind of spike is defined as crown spike in the present paper. Based on potential flow theory, a three-dimensional numerical model is established to simulate the motion of the free-surface spike generated by one bubble or a horizontal line of two in-phase bubbles. After the downward jet formed near the end of the collapse phase, the simulation of the free surface is performed to study the crown spike without regard to the toroidal bubble＇s effect. Calculations about the interaction between one bubble and free surface agree well with the experimental results conducted with a high-speed camera, and relative error is within 15%. Crown spike in both single- and two-bubble cases are simulated numerically. Different features and laws of the motion of crown spike, depending on the bubble-boundary distances and the inter-bubble distances, have been investigated.展开更多
Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wa...Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.展开更多
Rogue waves pose a significant threat to the safety of ships and offshore structures,making it crucial to understand their physical mechanisms,such as spatial-temporal focusing,which can lead to their formation.This s...Rogue waves pose a significant threat to the safety of ships and offshore structures,making it crucial to understand their physical mechanisms,such as spatial-temporal focusing,which can lead to their formation.This study investigates three-dimensional focused waves using the newly developed deep-water high-level Green-Naghdi(HLGN)model.Through numerical simulations,we evaluate the selection of the involved wave numbers within the HLGN model and present the algorithm for the three-dimensional implementation.Validation of the model is conducted through numerical reproduction of the three-dimensional focused waves considered in other’s laboratory measurements.The simulated wave profiles and velocity fields are compared with experimental data,demonstrating strong agreement.Discussion is provided about the robustness and accuracy of the HLGN model in simulating three-dimensional focused waves under deep-water conditions.展开更多
A three-dimensional coupled physical and water quality model was developed and applied to the Jiaozhou Bay to study water quality involving nutrients, biochemical oxygen demand, dissolved oxygen, and phytoplankton tha...A three-dimensional coupled physical and water quality model was developed and applied to the Jiaozhou Bay to study water quality involving nutrients, biochemical oxygen demand, dissolved oxygen, and phytoplankton that are closely related to eutrophication process. The physical model is a modified ECOM-si version with inclusion of flooding/draining processes over the intertidal zone. The water quality model is based on WASP5 which quantifies processes governing internal nutrients cycling, dissolved oxygen balance and phytoplankton growth. The model was used to simulate the spatial distribution and the temporal variation of water quality in the Jiaozhou Bay for the period of May 2005 to May 2006. In addition, the effect of reduction of riverine nutrients load was simulated and evaluated. The simulated results show that under the influence of nutrients discharged from river, the concentrations of nutrients and phytoplankton were higher in the northwest and northeast of the bay, and decreased from the inner bay to the outer. Affected by strong tidal mixing, the concentrations of all state variables were vertically homogeneous except in the deeper regions where a small gradient was found. Obvious seasonal variation of phytoplankton biomass was found, which exhibited two peaks in March and July, respectively. The variation of riverine waste loads had remarkable impact on nutrients concentration in coastal areas, but slightly altered the distribution in the center of the bay.展开更多
The Global and Regional Assimilation and Prediction System (GRAPES), a limited-area regional model, was used to simulate the onset of South China Sea summer monsoon. In view of the relatively insufficient information ...The Global and Regional Assimilation and Prediction System (GRAPES), a limited-area regional model, was used to simulate the onset of South China Sea summer monsoon. In view of the relatively insufficient information about the initial field in simulation predictions, the Advanced Microwave Sounding Unit-B (AMSU-B) data from a NOAA satellite were introduced to improve the initial values. By directly using the 3-dimensional variational data assimilation system of GRAPES, two schemes for assimilation tests were designed. In the design, Test 1 (T1) assimilates both sounding and AMSU-B data, and Test 2 (T2) assimilates only the conventional sounding data, before applying the model in simulation forecasts. Comparative experiments showed that the model was very sensitive to initial fields and successful in reproducing the monsoon onset, allocation of high- and low-level wind fields during the pentad of onset, and the northward advancement of the monsoon and monsoon rain bands. The scheme, however, simulated rainfall and the location of the subtropical high with deviations from observations. The simulated location of the subtropical high was more westward and northward and the simulated rainfall for the South China Sea was larger and covered a broader area.展开更多
The optimal evacuation scheme is studied based on the dam-break flood numerical simulation. A three- dimensional dam-break mathematical model combined with the volume of fluid （VOF） method is adopted. According to t...The optimal evacuation scheme is studied based on the dam-break flood numerical simulation. A three- dimensional dam-break mathematical model combined with the volume of fluid （VOF） method is adopted. According to the hydraulic information obtained from numerical simulation and selecting principles of evacuation emergency scheme, evacuation route analysis model is proposed, which consists of the road right model and random degree model. The road right model is used to calculate the consumption time in roads, and the random degree model is used to judge whether the roads are blocked. Then the shortest evacuation route is obtained based on Dijstra algorithm. Gongming Reservoir located in Shenzhen is taken as a case to study. The results show that industrial area I is flooded at 2 500 s, and after 5 500 s, most of industrial area II is submerged. The Hushan, Loucun Forest and Chaishan are not flooded around industrial area I and II. Based on the above analysis, the optimal evacuation scheme is determined.展开更多
By using the lattice model combined with finite element methods andstatistical techniques, a numerical approach is developed to establish mechanical models ofthree-dimensional heterogeneous brittle materials. A specia...By using the lattice model combined with finite element methods andstatistical techniques, a numerical approach is developed to establish mechanical models ofthree-dimensional heterogeneous brittle materials. A special numerical code is introduced, in whicha lattice model and statistical approaches are used to simulate the initial heterogeneity ofmaterial properties. The size of displacement-load step is adap-tively determined so that only fewelements would fail in each load step. When the tensile principal strain in an element exceeds theultimate strain of this element, the element is considered broken and its Young's modulus is set tobe very low. Some important behaviors of heterogeneous brittle materials are indicated using thiscode. Load-displacement curves and figures of three-dimensional fracture patterns are alsonumerically obtained, which are similar to those observed in laboratory tests.展开更多
The use of the mathematical models so far for three-dimensional flow has some limitations because of their simplifications. Many characteristics of the flow field can not be predicted by these models. In this paper th...The use of the mathematical models so far for three-dimensional flow has some limitations because of their simplifications. Many characteristics of the flow field can not be predicted by these models. In this paper the three dimensional elliptic governing equations are solved by finite-volume methods; the buoyancy extensions of the widely tested k-?. model is adapted. The method is first applied to calculate the field of side discharge into open channel flow. The results are in good agreement with those of ref. [7]. Then it is further used to the intake discharge problem which is of a typical layout in cooling-water projects, and the calculated results, which predict in detail the charactreistics of flow field, are reasonable.展开更多
The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially i...The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.展开更多
Submerged dikes were frequently built in curved estuarine channels to improve navigational conditions. A three-dimensional numerical model of tidal motions was developed to simulate flows through the submerged dike al...Submerged dikes were frequently built in curved estuarine channels to improve navigational conditions. A three-dimensional numerical model of tidal motions was developed to simulate flows through the submerged dike alignment in curved estuarine reach of complex bathymetry on the basis of satisfactory calibrations. Through the analysis of the morphological characteristics of the curved channel of the Huangpu River mouth in the Yangtze River Delta, a submerged dike alignment project was examined. The effect of the navigation channel improvement project and its influence on the flow regimes, the navigational conditions, the flood control and the operating conditions of harbours within the adjacent water area were evaluated, and the technical feasibility of the navigation channel improvement scheme was demonstrated.展开更多
Saline intrusion into sewage outfalls will greatly decrease the efficiency of the structures. The numerical model for this flow has been limited to one- and two-dimensional ones. In this article, a three-dimensional n...Saline intrusion into sewage outfalls will greatly decrease the efficiency of the structures. The numerical model for this flow has been limited to one- and two-dimensional ones. In this article, a three-dimensional numerical model for saline intrusion and purging in sewage outfalls was developed. The flow was modeled in three dimensions under turbulent conditions with the RNG k-ε turbulence model. The numerical results provided quantitative evidence of the fundamental flow mechanisms that took place during saline intrusion and purging. The comparisons of the results with that of two-dimensional model and that of experiments indicate that the three-dimensional numerical model developed in this article is more effective in predicting the internal flow in outfalls.展开更多
A three-dimensional meso-γ,scale elastic atmospheric numerical model has been established and used to simulate a severe storm case.The important characteristics of the modeled storm are close to those of the observed...A three-dimensional meso-γ,scale elastic atmospheric numerical model has been established and used to simulate a severe storm case.The important characteristics of the modeled storm are close to those of the observed case.展开更多
基金supported by the Natural Science Foundation of China(No.41574127)the China Postdoctoral Science Foundation(No.2017M622608)the project for the independent exploration of graduate students at Central South University(No.2017zzts008)
文摘In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over large-scale and complex terrain requires additional methods. To this end, we have proposed a new topography-capable By performing a two-dimensional Fourier transform in the horizontal directions, threedimensional partial differential equations in the spatial domain were transformed into a group of independent, one-dimensional differential equations engaged with different wave numbers. These independent differential equations are highly parallel across different wave numbers. differential equations with different wave numbers, and the efficiency of solving fixedbandwidth linear equations was further improved by a chasing method. In a synthetic test, a prism model was used to verify the accuracy and reliability of the proposed algorithm by comparing the numerical solution with the analytical solution. We studied the computational precision and efficiency with and without topography using different Fourier transform methods. The results showed that the Guass-FFT method has higher numerical precision, while the standard FFT method is superior, in terms of computation time, for inversion and quantitative interpretation under complicated terrain.
文摘Applying the methods of on-site observation and dynamic model, the research on the fronts at the Jiulong Estuary has been carried out, during which spatial and temporal distribution, dynamic characteristics and formation mechanism of salinity fronts are analyzed and discussed. The research shows that the estuarine fronts mainly lie in the area from the Jiyu Islet to the Haimen Island, outside of Yuweizai to Hulishan cross-section, the near coast of Yuweizai and the south of the Songyu-Gulangyu Channel. The fronts in the former two regions are formed directly by plume, while the one near the coast of Yuweizai is a tidal intrusion front caused by flood current and the one at the south of the Songyu-Gulangyu Channel is the result of current shear transformation. Under normal circumstances, fresh water of the Jiulong River mainly influences the inside of the Xiamen Bay, and when it is in typhoon seasons, plume front can affect the Taiwan Strait and has an effect on the biogeochemical Drocesses in the strait.
基金The Natural Science Foundation of Tianjin, China under contract No.08JCZDZT00200
文摘The characteristics of three-dimensional （3-D） tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, σ 3-D numerical model oil an unstructured triangular grid has been degeloped. The σ coordinate transforination, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simulations in the Oujing Estuary with submerged dyke and moveable boundary problems. The model is evaluated with the in situ data, and the results show that the calculated water elevations at 19 stations and currents at 19 profiler stations are in good agreement with measured data both in magnitude and phase. This numerical model is applied to the 3-D tidal circulation simulations of experiments in stopping flow transport through the South Branch of the Oujiang Estuary, and the feasibility to cutoff the flow in the South Branch of the Oujiang Estuary is demonstrated by numerical simulation experiments. The developed numerical model simulated the 3-D tidal current circulations in complicated coastal and estuarine waters very well.
文摘Saline intrusion into marine sewage ouffalls will greatly decrease the efficiency of sewage disposal. In order to investigate the mechanisms of this flow, in this paper, a three-dimensional numerical model based on FVM （Finite Volume Method） is established, The RNG κ-ε model is selected for turbulence modeling. The time-averaged vohtme fraction equations are introduced to simulate the stratification and inteffaeial exchange of sewage and seawater in outfalls. Validity of the established three-dimensional numerical model is evaluated by comparisons of numerical results with experimental data. With this three-dimensional numerical model, the internal flow characteristics in ouffalls for different sewage discharges are simulated. The results indicate that for a low sewage discharge, saline circulates in the outfall due to intrusion and both the inflowing momentum and the inteffaeial turbulent mixing are important mechanisms to extrude the saline. For a high sewage discharge, saline intrusion could be avoided. The inflow momentum is the main mechanism to extrude the saline and the inteffacial turbulent mixing is nut important relatively. Even at a high sewage discharge, the saline wedge would be retained in the main ouffall pipe after the risers are purged. It takes a long time for this saline wedge to be extruded by interracial turbulent mixing.
文摘In order to make the numerical calculation of viscous flows more convenient for the flows in channel with complicated profile governing equations expressed in the arbitrary curvilinear coordinates were derived by means of Favre density-weighted averaged method, and a turbulent model with effect of curvature modification was also derived. The numerical calculation of laminar and turbulent flown in divergent curved channels was carried out by means of parabolizeil computation method. The calculating results were used to analyze and investigate the aerodynamic performance of talor cascades in compressors preliminarily.
基金Project supported by the Major Basic Research Project of National Security of China(Grant No.613157)the Excellent Young Scientists Fund of China(Grant No.51222904)
文摘The motion of gas bubbles beneath a free surface will lead to a spike of fluid on the free surface. The distance of the bubbles to the free surface is the key factor to different phenomena. When the inception distance varies in some range, crown phenomenon would happen after the impact of weak buoyancy bubbles, so this kind of spike is defined as crown spike in the present paper. Based on potential flow theory, a three-dimensional numerical model is established to simulate the motion of the free-surface spike generated by one bubble or a horizontal line of two in-phase bubbles. After the downward jet formed near the end of the collapse phase, the simulation of the free surface is performed to study the crown spike without regard to the toroidal bubble＇s effect. Calculations about the interaction between one bubble and free surface agree well with the experimental results conducted with a high-speed camera, and relative error is within 15%. Crown spike in both single- and two-bubble cases are simulated numerically. Different features and laws of the motion of crown spike, depending on the bubble-boundary distances and the inter-bubble distances, have been investigated.
基金National Natural Science Foundation (Project number 40604013).
文摘Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.
基金Project supported by the National Natural Science Foundation of China(Grant No.12172099).
文摘Rogue waves pose a significant threat to the safety of ships and offshore structures,making it crucial to understand their physical mechanisms,such as spatial-temporal focusing,which can lead to their formation.This study investigates three-dimensional focused waves using the newly developed deep-water high-level Green-Naghdi(HLGN)model.Through numerical simulations,we evaluate the selection of the involved wave numbers within the HLGN model and present the algorithm for the three-dimensional implementation.Validation of the model is conducted through numerical reproduction of the three-dimensional focused waves considered in other’s laboratory measurements.The simulated wave profiles and velocity fields are compared with experimental data,demonstrating strong agreement.Discussion is provided about the robustness and accuracy of the HLGN model in simulating three-dimensional focused waves under deep-water conditions.
文摘A three-dimensional coupled physical and water quality model was developed and applied to the Jiaozhou Bay to study water quality involving nutrients, biochemical oxygen demand, dissolved oxygen, and phytoplankton that are closely related to eutrophication process. The physical model is a modified ECOM-si version with inclusion of flooding/draining processes over the intertidal zone. The water quality model is based on WASP5 which quantifies processes governing internal nutrients cycling, dissolved oxygen balance and phytoplankton growth. The model was used to simulate the spatial distribution and the temporal variation of water quality in the Jiaozhou Bay for the period of May 2005 to May 2006. In addition, the effect of reduction of riverine nutrients load was simulated and evaluated. The simulated results show that under the influence of nutrients discharged from river, the concentrations of nutrients and phytoplankton were higher in the northwest and northeast of the bay, and decreased from the inner bay to the outer. Affected by strong tidal mixing, the concentrations of all state variables were vertically homogeneous except in the deeper regions where a small gradient was found. Obvious seasonal variation of phytoplankton biomass was found, which exhibited two peaks in March and July, respectively. The variation of riverine waste loads had remarkable impact on nutrients concentration in coastal areas, but slightly altered the distribution in the center of the bay.
基金National 863 Special Project (2006AA01A123)Research on Interpretation Techniques for High-Resolution Numerical Prediction of Hunan ProvinceResearch on Rainstorm Forecast System with GRAPES and Its Application and Accessment in Hunan and Key Project of Science in Hunan (2008FJ1006)
文摘The Global and Regional Assimilation and Prediction System (GRAPES), a limited-area regional model, was used to simulate the onset of South China Sea summer monsoon. In view of the relatively insufficient information about the initial field in simulation predictions, the Advanced Microwave Sounding Unit-B (AMSU-B) data from a NOAA satellite were introduced to improve the initial values. By directly using the 3-dimensional variational data assimilation system of GRAPES, two schemes for assimilation tests were designed. In the design, Test 1 (T1) assimilates both sounding and AMSU-B data, and Test 2 (T2) assimilates only the conventional sounding data, before applying the model in simulation forecasts. Comparative experiments showed that the model was very sensitive to initial fields and successful in reproducing the monsoon onset, allocation of high- and low-level wind fields during the pentad of onset, and the northward advancement of the monsoon and monsoon rain bands. The scheme, however, simulated rainfall and the location of the subtropical high with deviations from observations. The simulated location of the subtropical high was more westward and northward and the simulated rainfall for the South China Sea was larger and covered a broader area.
基金Supported by Natural Science Foundation of Tianjin (No.09JCYBJC08700)the Foundation for Innovative Research Groups of National Natural Science Foundation of China (No.51021004)National Natural Science Foundation of China (No.90815019)
文摘The optimal evacuation scheme is studied based on the dam-break flood numerical simulation. A three- dimensional dam-break mathematical model combined with the volume of fluid （VOF） method is adopted. According to the hydraulic information obtained from numerical simulation and selecting principles of evacuation emergency scheme, evacuation route analysis model is proposed, which consists of the road right model and random degree model. The road right model is used to calculate the consumption time in roads, and the random degree model is used to judge whether the roads are blocked. Then the shortest evacuation route is obtained based on Dijstra algorithm. Gongming Reservoir located in Shenzhen is taken as a case to study. The results show that industrial area I is flooded at 2 500 s, and after 5 500 s, most of industrial area II is submerged. The Hushan, Loucun Forest and Chaishan are not flooded around industrial area I and II. Based on the above analysis, the optimal evacuation scheme is determined.
文摘By using the lattice model combined with finite element methods andstatistical techniques, a numerical approach is developed to establish mechanical models ofthree-dimensional heterogeneous brittle materials. A special numerical code is introduced, in whicha lattice model and statistical approaches are used to simulate the initial heterogeneity ofmaterial properties. The size of displacement-load step is adap-tively determined so that only fewelements would fail in each load step. When the tensile principal strain in an element exceeds theultimate strain of this element, the element is considered broken and its Young's modulus is set tobe very low. Some important behaviors of heterogeneous brittle materials are indicated using thiscode. Load-displacement curves and figures of three-dimensional fracture patterns are alsonumerically obtained, which are similar to those observed in laboratory tests.
基金The Project supported by the Foundation of Science Research of the State Education Commission of the People's Republic of China
文摘The use of the mathematical models so far for three-dimensional flow has some limitations because of their simplifications. Many characteristics of the flow field can not be predicted by these models. In this paper the three dimensional elliptic governing equations are solved by finite-volume methods; the buoyancy extensions of the widely tested k-?. model is adapted. The method is first applied to calculate the field of side discharge into open channel flow. The results are in good agreement with those of ref. [7]. Then it is further used to the intake discharge problem which is of a typical layout in cooling-water projects, and the calculated results, which predict in detail the charactreistics of flow field, are reasonable.
基金the National Supercomputer Center in Tianjin for their patient assistance in providing the compilation environment.We thank the editor,Huajian Yao,for handling the manuscript and Mingming Li and another anonymous reviewer for their constructive comments.The research leading to these results has received funding from National Natural Science Foundation of China projects(Grant Nos.92355302 and 42121005)Taishan Scholar projects(Grant No.tspd20210305)others(Grant Nos.XDB0710000,L2324203,XK2023DXC001,LSKJ202204400,and ZR2021ZD09).
文摘The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.
文摘Submerged dikes were frequently built in curved estuarine channels to improve navigational conditions. A three-dimensional numerical model of tidal motions was developed to simulate flows through the submerged dike alignment in curved estuarine reach of complex bathymetry on the basis of satisfactory calibrations. Through the analysis of the morphological characteristics of the curved channel of the Huangpu River mouth in the Yangtze River Delta, a submerged dike alignment project was examined. The effect of the navigation channel improvement project and its influence on the flow regimes, the navigational conditions, the flood control and the operating conditions of harbours within the adjacent water area were evaluated, and the technical feasibility of the navigation channel improvement scheme was demonstrated.
文摘Saline intrusion into sewage outfalls will greatly decrease the efficiency of the structures. The numerical model for this flow has been limited to one- and two-dimensional ones. In this article, a three-dimensional numerical model for saline intrusion and purging in sewage outfalls was developed. The flow was modeled in three dimensions under turbulent conditions with the RNG k-ε turbulence model. The numerical results provided quantitative evidence of the fundamental flow mechanisms that took place during saline intrusion and purging. The comparisons of the results with that of two-dimensional model and that of experiments indicate that the three-dimensional numerical model developed in this article is more effective in predicting the internal flow in outfalls.
文摘A three-dimensional meso-γ,scale elastic atmospheric numerical model has been established and used to simulate a severe storm case.The important characteristics of the modeled storm are close to those of the observed case.
