Large language models(LLMs)and natural language processing(NLP)have significant promise to improve efficiency and refine healthcare decision-making and clinical results.Numerous domains,including healthcare,are rapidl...Large language models(LLMs)and natural language processing(NLP)have significant promise to improve efficiency and refine healthcare decision-making and clinical results.Numerous domains,including healthcare,are rapidly adopting LLMs for the classification of biomedical textual data in medical research.The LLM can derive insights from intricate,extensive,unstructured training data.Variants need to be accurately identified and classified to advance genetic research,provide individualized treatment,and assist physicians in making better choices.However,the sophisticated and perplexing language of medical reports is often beyond the capabilities of the devices we now utilize.Such an approach may result in incorrect diagnoses,which could affect a patient’s prognosis and course of therapy.This study evaluated the efficacy of the proposed model by looking at publicly accessible textual clinical data.We have cleaned the clinical textual data using various text preprocessing methods,including stemming,tokenization,and stop word removal.The important features are extracted using Bag of Words(BoW)and Term Frequency-Inverse Document Frequency(TFIDF)feature engineering methods.The important motive of this study is to predict the genetic variants based on the clinical evidence using a novel method with minimal error.According to the experimental results,the random forest model achieved 61%accuracy with 67%precision for class 9 using TFIDF features and 63%accuracy and a 73%F1 score for class 9 using Bag of Words features.The accuracy of the proposed BERT(Bidirectional Encoder Representations from Transformers)model was 70%with 5-fold cross-validation and 71%with 10-fold cross-validation.The research results provide a comprehensive overview of current LLM methods in healthcare,benefiting academics as well as professionals in the discipline.展开更多
With the advancement of deep learning in the automotive domain,more and more researchers are focusing on autonomous driving.Among these tasks,free space detection is particularly crucial.Currently,many model-based app...With the advancement of deep learning in the automotive domain,more and more researchers are focusing on autonomous driving.Among these tasks,free space detection is particularly crucial.Currently,many model-based approaches have achieved autonomous driving on well-structured urban roads,but these efforts primarily focus on urban road environments.In contrast,there are fewer deep learningmethods specifically designed for off-road traversable area detection,and their effectiveness is not yet satisfactory.This is because detecting traversable areas in complex outdoor environments poses significant challenges,and current methods often rely on single-image inputs,which do not align with contemporary multimodal approaches.Therefore,in this study,we propose a CFH-Net model for off-road traversable area detection.This model employs a Transformer architecture to enhance its capability of capturing global information.For multimodal feature extraction and fusion,we integrate the CM-FRM module for feature extraction and introduce the novel FFX module for feature fusion,thereby improving the perception capability of autonomous vehicles on unstructured roads.To address upsampling,we propose a new convolution precorrection method to reduce model parameters and computational complexity while enhancing the model’s ability to capture complex features.Finally,we conducted experiments on the ORFD off-road dataset and achieved outstanding results.展开更多
There is no unified planning standard for unstructured roads,and the morphological structures of these roads are complex and varied.It is important to maintain a balance between accuracy and speed for unstructured roa...There is no unified planning standard for unstructured roads,and the morphological structures of these roads are complex and varied.It is important to maintain a balance between accuracy and speed for unstructured road extraction models.Unstructured road extraction algorithms based on deep learning have problems such as high model complexity,high computational cost,and the inability to adapt to current edge computing devices.Therefore,it is best to use lightweight network models.Considering the need for lightweight models and the characteristics of unstructured roads with different pattern shapes,such as blocks and strips,a TMB(Triple Multi-Block)feature extraction module is proposed,and the overall structure of the TMBNet network is described.The TMB module was compared with SS-nbt,Non-bottleneck-1D,and other modules via experiments.The feasibility and effectiveness of the TMB module design were proven through experiments and visualizations.The comparison experiment,using multiple convolution kernel categories,proved that the TMB module can improve the segmentation accuracy of the network.The comparison with different semantic segmentation networks demonstrates that the TMBNet network has advantages in terms of unstructured road extraction.展开更多
Storm surge is often the marine disaster that poses the greatest threat to life and property in coastal areas.Accurate and timely issuance of storm surge warnings to take appropriate countermeasures is an important me...Storm surge is often the marine disaster that poses the greatest threat to life and property in coastal areas.Accurate and timely issuance of storm surge warnings to take appropriate countermeasures is an important means to reduce storm surge-related losses.Storm surge numerical models are important for storm surge forecasting.To further improve the performance of the storm surge forecast models,we developed a numerical storm surge forecast model based on an unstructured spherical centroidal Voronoi tessellation(SCVT)grid.The model is based on shallow water equations in vector-invariant form,and is discretized by Arakawa C grid.The SCVT grid can not only better describe the coastline information but also avoid rigid transitions,and it has a better global consistency by generating high-resolution grids in the key areas through transition refinement.In addition,the simulation speed of the model is accelerated by using the openACC-based GPU acceleration technology to meet the timeliness requirements of operational ensemble forecast.It only takes 37 s to simulate a day in the coastal waters of China.The newly developed storm surge model was applied to simulate typhoon-induced storm surges in the coastal waters of China.The hindcast experiments on the selected representative typhoon-induced storm surge processes indicate that the model can reasonably simulate the distribution characteristics of storm surges.The simulated maximum storm surges and their occurrence times are consistent with the observed data at the representative tide gauge stations,and the mean absolute errors are 3.5 cm and 0.6 h respectively,showing high accuracy and application prospects.展开更多
An upwind scheme based on the unstructured mesh is developed to solve ideal 2-D magnetohydrodynamics （MHD） equations. The inviscid fluxes are approximated by using the modified advection upstream splitting method （...An upwind scheme based on the unstructured mesh is developed to solve ideal 2-D magnetohydrodynamics （MHD） equations. The inviscid fluxes are approximated by using the modified advection upstream splitting method （AUSM） scheme, and a 5-stage explicit Runge-Kutta scheme is adopted in the time integration. To avoid the influence of the magnetic field divergence created during the simulation, the hyperbolic divergence cleaning method is introduced. The shock-capturing properties of the method are verified by solving the MHD shock-tube problem. Then the 2-D nozzle flow with the magnetic field is numerically simulated on the unstructured mesh. Computational results demonstrate the effects of the magnetic field and agree well with those from references.展开更多
Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficien...Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.展开更多
A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C-unstructured non-orthogonal grid in the hor- izontal d...A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C-unstructured non-orthogonal grid in the hor- izontal direction and z-level grid in the vertical direction. The three-dimensional model is discretized by the semi-implicit finite-volume method, in that the free-surface and the vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly by an integral method. The partial cell method is used for resolving topography, which enables the model to better represent irregular topography. The model has been tested against analytical cases for wind and tidal oscillation circulation, and is applied to simulating the tidal flow in the Bohal Sea. The results are in good agreement both with the analytical solutions and measurement results.展开更多
The lower-upper symmetric Gauss-Seidel （LU-SGS） implicit relaxation has been widely used because it has the merits of less dependency on grid topology, low numerical complexity and modest memory requirements. In ori...The lower-upper symmetric Gauss-Seidel （LU-SGS） implicit relaxation has been widely used because it has the merits of less dependency on grid topology, low numerical complexity and modest memory requirements. In original LU-SGS scheme, the implicit system matrix is constructed based on the splitting of convective flux Jacobian according to its spectral radius. Although this treatment has the merit of reducing computational complexity and helps to ensure the diagonally dominant property of the implicit system matrix, it can also cause serious distortions on the implicit system matrix because too many approximations are introduced by this splitting method if the contravariant velocity is small or close to sonic speed. To overcome this shortcoming, an improved LU-SGS scheme with a hybrid construction method for the implicit system matrix is developed in this paper. The hybrid way is that： on the cell faces having small contravariant velocity or transonic contravariant velocity, the accurate derivative of the convective flux term is used to construct more accurate implicit system matrix, while the original Jacobian splitting method is adopted on the other cell faces to reduce computational complexity and ensure the diagonally dominant property of the implicit system matrix. To investigate the convergence performance of the improved LU-SGS scheme, 2D and 3D turbulent flows around the NACA0012 airfoil, RAE2822 airfoil and LANN wing are simulated on hybrid unstructured meshes. The nu- merical results show that the improved LU-SGS scheme is significantly more efficient than the original LU-SGS scheme.展开更多
Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured...Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured-grid-based methods can model complex underground structures with high accuracy and overcome the defects of traditional methods, such as the high computational cost for improving model accuracy and the difficulty of inverting with topography. In this paper, we used the limited-memory quasi-Newton(L-BFGS) method with an unstructured finite-element grid to perform 3D MT inversions. This method avoids explicitly calculating Hessian matrices, which greatly reduces the memory requirements. After the first iteration, the approximate inverse Hessian matrix well approximates the true one, and the Newton step(set to 1) can meet the sufficient descent condition. Only one calculation of the objective function and its gradient are needed for each iteration, which greatly improves its computational efficiency. This approach is well-suited for large-scale 3D MT inversions. We have tested our algorithm on data with and without topography, and the results matched the real models well. We can recommend performing inversions based on an unstructured finite-element method and the L-BFGS method for situations with topography and complex underground structures.展开更多
We propose a new 3D inversion scheme to invert the near-and transition-zone data of CSAMT with topography accurately.In this new method,the earth was discretized into unstructured tetrahedra to fit the ragged topograp...We propose a new 3D inversion scheme to invert the near-and transition-zone data of CSAMT with topography accurately.In this new method,the earth was discretized into unstructured tetrahedra to fit the ragged topography and the vector fi nite-element method was adopted to obtain precise responses and good sensitivity.To simulate the attitude and shape of the transmitter,we divided a long-grounded transmitter into dipoles and integrated these dipoles to obtain good responses in the near-and transition-fi eld zones.Next,we designed an L2 norm-based objective functional and applied a standard quasi-Newton method as the optimization method to solve the inverse problem and guarantee steady convergence.We tested our 3D inversion method first on synthetic data and then on a field dataset acquired from select sites near Changbai Mountain,China.In both tests,the new inversion algorithm achieved excellent fitting between the predicted and observed data,even in near-and transition-fi eld zones,and the inversion results agreed well with the true model.These fi ndings reveal that the proposed algorithm is eff ective for 3D inversion of CSAMT data.展开更多
Passive worms can passively propagate through embedding themselves into some sharing files, which can result in significant damage to unstructured P2P networks. To study the passive worm behaviors, this paper firstly ...Passive worms can passively propagate through embedding themselves into some sharing files, which can result in significant damage to unstructured P2P networks. To study the passive worm behaviors, this paper firstly analyzes and obtains the average delay for all peers in the whole transmitting process due to the limitation of network throughput, and then proposes a mathematical model for the propagation of passive worms over the unstructured P2P networks. The model mainly takes the effect of the network throughput into account, and applies a new healthy files dissemination-based defense strategy according to the file popularity which follows the Zipf distribution. The simulation results show that the propagation of passive worms is mainly governed by the number of hops, initially infected files and uninfected files. The larger the number of hops, the more rapidly the passive worms propagate. If the number of the initially infected files is increased by the attackers, the propagation speed of passive worms increases obviously. A larger size of the uninfected file results in a better attack performance. However, the number of files generated by passive worms is not an important factor governing the propagation of passive worms. The effectiveness of healthy files dissemination strategy is verified. This model can provide a guideline in the control of unstructured P2P networks as well as passive worm defense.展开更多
In this study, porosity was introduced into two-dimensional shallow water equations to reflect the effects of obstructions, leading to the modification of the expressions for the flux and source terms. An extra porosi...In this study, porosity was introduced into two-dimensional shallow water equations to reflect the effects of obstructions, leading to the modification of the expressions for the flux and source terms. An extra porosity source term appears in the momentum equation. The numerical model of the shallow water equations with porosity is presented with the finite volume method on unstructured grids and the modified Roe-type approximate Riemann solver. The source terms of the bed slope and porosity are both decomposed in the characteristic direction so that the numerical scheme can exactly satisfy the conservative property. The present model was tested with a dam break with discontinuous porosity and a flash flood in the Toce River Valley. The results show that the model can simulate the influence of obstructions, and the numerical scheme can maintain the flux balance at the interface with high efficiency and resolution.展开更多
A novel VOF-type volume-tracking method for two-dimensional free-surface flows based on the unstructured triangular mesh is presented. Owing to the inherent merit of the unstructured triangular mesh in fitting curved ...A novel VOF-type volume-tracking method for two-dimensional free-surface flows based on the unstructured triangular mesh is presented. Owing to the inherent merit of the unstructured triangular mesh in fitting curved boundaries, this method can handle the free-surface problems with complex geometries accurately and directly, without introducing any complicated boundary treatment or artificial diffusion. The method solves the volume transport equation geometrically through the Modified Lagrangian-Eulerian Re-map （MLER） method, which is applied to advective fluid volumes. Moreover, the PLIC method is adopted to give a second-order reconstructed interface approximation. To validate this method, two advection tests were performed for the establishment of the accuracy and convergence rate of the solutions. Numerical results for these complex tests provide convincing evidence for the excellent solution quality and fidelity of the method.展开更多
This paper presents a method to generate unstructured adaptive meshes with moving boundaries and its application to CFD. Delaunay triangulation criterion in conjunction with the automatic point creation is used to gen...This paper presents a method to generate unstructured adaptive meshes with moving boundaries and its application to CFD. Delaunay triangulation criterion in conjunction with the automatic point creation is used to generate 2 D and 3 D unstructured grids. A local grid regeneration method is proposed to cope with moving boundaries. Numerical examples include the interactions of shock waves with movable bodies and the movement of a projectile within a ram accelerator, illustrating an efficient and robust mesh generation method developed.展开更多
A three-dimensional(3D)step-by-step inversion strategy for audio magnetotellurics(AMT)is investigated in this study.The objective function is minimized by iteratively solving the Gauss-Newton normal equation,and the i...A three-dimensional(3D)step-by-step inversion strategy for audio magnetotellurics(AMT)is investigated in this study.The objective function is minimized by iteratively solving the Gauss-Newton normal equation,and the inversion region is discretized with unstructured tetrahedral elements.The inversion proceeds step-by-step from a coarse mesh to a fine mesh.In the inversion iteration process,a mesh is adaptively optimized according to the spatial gradient information about the model resistivity to fine delineate the boundaries of abnormal bodies.In the early stage of inversion execution,a coarse mesh is used for inversion,and the inversion stability is improved by reducing the number of inversion elements.In addition,mesh refinement is performed in the iterative inversion process.The inversion results obtained from the previous mesh are used as the reference and initial models for the next mesh iterative inversion.The step-by-step inversion strategy can ensure that the inversion is performed in the correct direction,improving the inversion stability and results gradually.Synthetic results show that the step-by-step inversion strategy with a Gauss-Newton method for 3D AMT inversion is stable and reliable,which lays a foundation for further practical 3D AMT data inversion.展开更多
A two-dimensional coastal ocean model based on unstructured C-grid is built, in which the momentum equation is discretized on the faces of each cell, and the continuity equation is discretized on the cell. The model i...A two-dimensional coastal ocean model based on unstructured C-grid is built, in which the momentum equation is discretized on the faces of each cell, and the continuity equation is discretized on the cell. The model is discretized by semi-implicit finite volume method, in that the free surface is semi-implicit and the bottom friction is implicit, thereby removing stability limitations associated with the surface gravity wave and friction. The remaining terms in the momentum equations are discretized explicitly by integral finite volume method and second-order Adams-Bashforth method. Tidal flow in the polar quadrant with known analytic solution is employed to test the proposed model. Finally, the performance of the present model to simulate tidal flow in a geometrically complex domain is examined by simulation of tidal currents in the Pearl River Estuary.展开更多
In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstructio...In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstruction technique in terms of compact stencil is improved to reduce local condition numbers. To further improve the efficiency of computation, the adaptive mesh refinement technique is implemented in the framework of high-order finite volume methods. Mesh refinement and coarsening criteria are chosen to be the indicators for certain flow structures. One important challenge of the adaptive mesh refinement technique on unstructured grids is the dynamic load balancing in parallel computation. To solve this problem, the open-source library p4 est based on the forest of octrees is adopted. Several two-and three-dimensional test cases are computed to verify the accuracy and robustness of the proposed numerical schemes.展开更多
Unstructured P2P has power-law link distribution, and the random walk in power-law networks is analyzed. The analysis results show that the probability that a random walker walks through the high degree nodes is high ...Unstructured P2P has power-law link distribution, and the random walk in power-law networks is analyzed. The analysis results show that the probability that a random walker walks through the high degree nodes is high in the power-law network, and the information on the high degree nodes can be easily found through random walk. Random walk spread and random walk search method （RWSS） is proposed based on the analysis result. Simulation results show that RWSS achieves high success rates at low cost and is robust to high degree node failure.展开更多
A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method o...A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method on triangular grids. The developed HWENO methodology utilizes high-order derivative information to keep WENO re- construction stencils in the von Neumann neighborhood. A simple and efficient technique is also proposed to enhance the smoothness of the existing stencils, making higher-order scheme stable and simplifying the reconstruction process at the same time. The resulting HWENO-based limiters are as compact as the underlying DG schemes and therefore easy to implement. Numerical results for a wide range of flow conditions demonstrate that for DG schemes of up to fourth order of accuracy, the designed HWENO limiters can simul- taneously obtain uniform high order accuracy and sharp, es- sentially non-oscillatory shock transition.展开更多
An efficient high-order numerical method for supersonic reactive flows is proposed in this article.The reactive source term and convection term are solved separately by splitting scheme.In the reaction step,an adaptiv...An efficient high-order numerical method for supersonic reactive flows is proposed in this article.The reactive source term and convection term are solved separately by splitting scheme.In the reaction step,an adaptive time-step method is presented,which can improve the efficiency greatly.In the convection step,a third-order accurate weighted essentially non-oscillatory（WENO）method is adopted to reconstruct the solution in the unstructured grids.Numerical results show that our new method can capture the correct propagation speed of the detonation wave exactly even in coarse grids,while high order accuracy can be achieved in the smooth region.In addition,the proposed adaptive splitting method can reduce the computational cost greatly compared with the traditional splitting method.展开更多
基金funded by Princess Nourah bint Abdulrahman University and Researchers Supporting Project number(PNURSP2025R346),Princess Nourah bint Abdulrahman University,Riyadh,Saudi Arabia.
文摘Large language models(LLMs)and natural language processing(NLP)have significant promise to improve efficiency and refine healthcare decision-making and clinical results.Numerous domains,including healthcare,are rapidly adopting LLMs for the classification of biomedical textual data in medical research.The LLM can derive insights from intricate,extensive,unstructured training data.Variants need to be accurately identified and classified to advance genetic research,provide individualized treatment,and assist physicians in making better choices.However,the sophisticated and perplexing language of medical reports is often beyond the capabilities of the devices we now utilize.Such an approach may result in incorrect diagnoses,which could affect a patient’s prognosis and course of therapy.This study evaluated the efficacy of the proposed model by looking at publicly accessible textual clinical data.We have cleaned the clinical textual data using various text preprocessing methods,including stemming,tokenization,and stop word removal.The important features are extracted using Bag of Words(BoW)and Term Frequency-Inverse Document Frequency(TFIDF)feature engineering methods.The important motive of this study is to predict the genetic variants based on the clinical evidence using a novel method with minimal error.According to the experimental results,the random forest model achieved 61%accuracy with 67%precision for class 9 using TFIDF features and 63%accuracy and a 73%F1 score for class 9 using Bag of Words features.The accuracy of the proposed BERT(Bidirectional Encoder Representations from Transformers)model was 70%with 5-fold cross-validation and 71%with 10-fold cross-validation.The research results provide a comprehensive overview of current LLM methods in healthcare,benefiting academics as well as professionals in the discipline.
文摘With the advancement of deep learning in the automotive domain,more and more researchers are focusing on autonomous driving.Among these tasks,free space detection is particularly crucial.Currently,many model-based approaches have achieved autonomous driving on well-structured urban roads,but these efforts primarily focus on urban road environments.In contrast,there are fewer deep learningmethods specifically designed for off-road traversable area detection,and their effectiveness is not yet satisfactory.This is because detecting traversable areas in complex outdoor environments poses significant challenges,and current methods often rely on single-image inputs,which do not align with contemporary multimodal approaches.Therefore,in this study,we propose a CFH-Net model for off-road traversable area detection.This model employs a Transformer architecture to enhance its capability of capturing global information.For multimodal feature extraction and fusion,we integrate the CM-FRM module for feature extraction and introduce the novel FFX module for feature fusion,thereby improving the perception capability of autonomous vehicles on unstructured roads.To address upsampling,we propose a new convolution precorrection method to reduce model parameters and computational complexity while enhancing the model’s ability to capture complex features.Finally,we conducted experiments on the ORFD off-road dataset and achieved outstanding results.
基金Supported by National Natural Science Foundation of China(Grant Nos.62261160575,61991414,61973036)Technical Field Foundation of the National Defense Science and Technology 173 Program of China(Grant Nos.20220601053,20220601030)。
文摘There is no unified planning standard for unstructured roads,and the morphological structures of these roads are complex and varied.It is important to maintain a balance between accuracy and speed for unstructured road extraction models.Unstructured road extraction algorithms based on deep learning have problems such as high model complexity,high computational cost,and the inability to adapt to current edge computing devices.Therefore,it is best to use lightweight network models.Considering the need for lightweight models and the characteristics of unstructured roads with different pattern shapes,such as blocks and strips,a TMB(Triple Multi-Block)feature extraction module is proposed,and the overall structure of the TMBNet network is described.The TMB module was compared with SS-nbt,Non-bottleneck-1D,and other modules via experiments.The feasibility and effectiveness of the TMB module design were proven through experiments and visualizations.The comparison experiment,using multiple convolution kernel categories,proved that the TMB module can improve the segmentation accuracy of the network.The comparison with different semantic segmentation networks demonstrates that the TMBNet network has advantages in terms of unstructured road extraction.
基金The National Natural Science Foundation of China under contract No.42076214.
文摘Storm surge is often the marine disaster that poses the greatest threat to life and property in coastal areas.Accurate and timely issuance of storm surge warnings to take appropriate countermeasures is an important means to reduce storm surge-related losses.Storm surge numerical models are important for storm surge forecasting.To further improve the performance of the storm surge forecast models,we developed a numerical storm surge forecast model based on an unstructured spherical centroidal Voronoi tessellation(SCVT)grid.The model is based on shallow water equations in vector-invariant form,and is discretized by Arakawa C grid.The SCVT grid can not only better describe the coastline information but also avoid rigid transitions,and it has a better global consistency by generating high-resolution grids in the key areas through transition refinement.In addition,the simulation speed of the model is accelerated by using the openACC-based GPU acceleration technology to meet the timeliness requirements of operational ensemble forecast.It only takes 37 s to simulate a day in the coastal waters of China.The newly developed storm surge model was applied to simulate typhoon-induced storm surges in the coastal waters of China.The hindcast experiments on the selected representative typhoon-induced storm surge processes indicate that the model can reasonably simulate the distribution characteristics of storm surges.The simulated maximum storm surges and their occurrence times are consistent with the observed data at the representative tide gauge stations,and the mean absolute errors are 3.5 cm and 0.6 h respectively,showing high accuracy and application prospects.
文摘An upwind scheme based on the unstructured mesh is developed to solve ideal 2-D magnetohydrodynamics （MHD） equations. The inviscid fluxes are approximated by using the modified advection upstream splitting method （AUSM） scheme, and a 5-stage explicit Runge-Kutta scheme is adopted in the time integration. To avoid the influence of the magnetic field divergence created during the simulation, the hyperbolic divergence cleaning method is introduced. The shock-capturing properties of the method are verified by solving the MHD shock-tube problem. Then the 2-D nozzle flow with the magnetic field is numerically simulated on the unstructured mesh. Computational results demonstrate the effects of the magnetic field and agree well with those from references.
基金National High Technology Research and Development Program(863 Program)(No.2006AA06Z105,2007AA06Z134)
文摘Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.
基金The Major State Basic Research Program of China under contract No. 2012CB417002the National Natural Science Foundation of China under contract Nos 50909065 and 51109039
文摘A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C-unstructured non-orthogonal grid in the hor- izontal direction and z-level grid in the vertical direction. The three-dimensional model is discretized by the semi-implicit finite-volume method, in that the free-surface and the vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly by an integral method. The partial cell method is used for resolving topography, which enables the model to better represent irregular topography. The model has been tested against analytical cases for wind and tidal oscillation circulation, and is applied to simulating the tidal flow in the Bohal Sea. The results are in good agreement both with the analytical solutions and measurement results.
基金Foundation item： National Natural Science Foundation of China （10802067）
文摘The lower-upper symmetric Gauss-Seidel （LU-SGS） implicit relaxation has been widely used because it has the merits of less dependency on grid topology, low numerical complexity and modest memory requirements. In original LU-SGS scheme, the implicit system matrix is constructed based on the splitting of convective flux Jacobian according to its spectral radius. Although this treatment has the merit of reducing computational complexity and helps to ensure the diagonally dominant property of the implicit system matrix, it can also cause serious distortions on the implicit system matrix because too many approximations are introduced by this splitting method if the contravariant velocity is small or close to sonic speed. To overcome this shortcoming, an improved LU-SGS scheme with a hybrid construction method for the implicit system matrix is developed in this paper. The hybrid way is that： on the cell faces having small contravariant velocity or transonic contravariant velocity, the accurate derivative of the convective flux term is used to construct more accurate implicit system matrix, while the original Jacobian splitting method is adopted on the other cell faces to reduce computational complexity and ensure the diagonally dominant property of the implicit system matrix. To investigate the convergence performance of the improved LU-SGS scheme, 2D and 3D turbulent flows around the NACA0012 airfoil, RAE2822 airfoil and LANN wing are simulated on hybrid unstructured meshes. The nu- merical results show that the improved LU-SGS scheme is significantly more efficient than the original LU-SGS scheme.
基金financially supported by the National Natural Science Foundation of China(No.41774125)Key Program of National Natural Science Foundation of China(No.41530320)+1 种基金the Key National Research Project of China(Nos.2016YFC0303100 and 2017YFC0601900)the Strategic Priority Research Program of Chinese Academy of Sciences Pilot Special(No.XDA 14020102)
文摘Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured-grid-based methods can model complex underground structures with high accuracy and overcome the defects of traditional methods, such as the high computational cost for improving model accuracy and the difficulty of inverting with topography. In this paper, we used the limited-memory quasi-Newton(L-BFGS) method with an unstructured finite-element grid to perform 3D MT inversions. This method avoids explicitly calculating Hessian matrices, which greatly reduces the memory requirements. After the first iteration, the approximate inverse Hessian matrix well approximates the true one, and the Newton step(set to 1) can meet the sufficient descent condition. Only one calculation of the objective function and its gradient are needed for each iteration, which greatly improves its computational efficiency. This approach is well-suited for large-scale 3D MT inversions. We have tested our algorithm on data with and without topography, and the results matched the real models well. We can recommend performing inversions based on an unstructured finite-element method and the L-BFGS method for situations with topography and complex underground structures.
基金the Strategic Priority Research Program of the Chinese Academy of Sciences(No.XDA14020102)the National Natural Science Foundation of China(Nos.41774125,41530320,41904104)+1 种基金the Key National Research Project of China(No.2018YFC0603300)the S&T Program of Beijing(No.Z181100005718001).
文摘We propose a new 3D inversion scheme to invert the near-and transition-zone data of CSAMT with topography accurately.In this new method,the earth was discretized into unstructured tetrahedra to fit the ragged topography and the vector fi nite-element method was adopted to obtain precise responses and good sensitivity.To simulate the attitude and shape of the transmitter,we divided a long-grounded transmitter into dipoles and integrated these dipoles to obtain good responses in the near-and transition-fi eld zones.Next,we designed an L2 norm-based objective functional and applied a standard quasi-Newton method as the optimization method to solve the inverse problem and guarantee steady convergence.We tested our 3D inversion method first on synthetic data and then on a field dataset acquired from select sites near Changbai Mountain,China.In both tests,the new inversion algorithm achieved excellent fitting between the predicted and observed data,even in near-and transition-fi eld zones,and the inversion results agreed well with the true model.These fi ndings reveal that the proposed algorithm is eff ective for 3D inversion of CSAMT data.
基金National Natural Science Foundation of China (No.60633020 and No. 90204012)Natural Science Foundation of Hebei Province (No. F2006000177)
文摘Passive worms can passively propagate through embedding themselves into some sharing files, which can result in significant damage to unstructured P2P networks. To study the passive worm behaviors, this paper firstly analyzes and obtains the average delay for all peers in the whole transmitting process due to the limitation of network throughput, and then proposes a mathematical model for the propagation of passive worms over the unstructured P2P networks. The model mainly takes the effect of the network throughput into account, and applies a new healthy files dissemination-based defense strategy according to the file popularity which follows the Zipf distribution. The simulation results show that the propagation of passive worms is mainly governed by the number of hops, initially infected files and uninfected files. The larger the number of hops, the more rapidly the passive worms propagate. If the number of the initially infected files is increased by the attackers, the propagation speed of passive worms increases obviously. A larger size of the uninfected file results in a better attack performance. However, the number of files generated by passive worms is not an important factor governing the propagation of passive worms. The effectiveness of healthy files dissemination strategy is verified. This model can provide a guideline in the control of unstructured P2P networks as well as passive worm defense.
基金supported by the National Natural Science Foundation of China (Grants No. 50909065 and 51109039)the National Basic Research Program of China (973 Program, Grant No. 2012CB417002)
文摘In this study, porosity was introduced into two-dimensional shallow water equations to reflect the effects of obstructions, leading to the modification of the expressions for the flux and source terms. An extra porosity source term appears in the momentum equation. The numerical model of the shallow water equations with porosity is presented with the finite volume method on unstructured grids and the modified Roe-type approximate Riemann solver. The source terms of the bed slope and porosity are both decomposed in the characteristic direction so that the numerical scheme can exactly satisfy the conservative property. The present model was tested with a dam break with discontinuous porosity and a flash flood in the Toce River Valley. The results show that the model can simulate the influence of obstructions, and the numerical scheme can maintain the flux balance at the interface with high efficiency and resolution.
文摘A novel VOF-type volume-tracking method for two-dimensional free-surface flows based on the unstructured triangular mesh is presented. Owing to the inherent merit of the unstructured triangular mesh in fitting curved boundaries, this method can handle the free-surface problems with complex geometries accurately and directly, without introducing any complicated boundary treatment or artificial diffusion. The method solves the volume transport equation geometrically through the Modified Lagrangian-Eulerian Re-map （MLER） method, which is applied to advective fluid volumes. Moreover, the PLIC method is adopted to give a second-order reconstructed interface approximation. To validate this method, two advection tests were performed for the establishment of the accuracy and convergence rate of the solutions. Numerical results for these complex tests provide convincing evidence for the excellent solution quality and fidelity of the method.
基金the Natural Science Foundation of China(No.5 99760 13 and1983 2 0 3 0 )
文摘This paper presents a method to generate unstructured adaptive meshes with moving boundaries and its application to CFD. Delaunay triangulation criterion in conjunction with the automatic point creation is used to generate 2 D and 3 D unstructured grids. A local grid regeneration method is proposed to cope with moving boundaries. Numerical examples include the interactions of shock waves with movable bodies and the movement of a projectile within a ram accelerator, illustrating an efficient and robust mesh generation method developed.
文摘A three-dimensional(3D)step-by-step inversion strategy for audio magnetotellurics(AMT)is investigated in this study.The objective function is minimized by iteratively solving the Gauss-Newton normal equation,and the inversion region is discretized with unstructured tetrahedral elements.The inversion proceeds step-by-step from a coarse mesh to a fine mesh.In the inversion iteration process,a mesh is adaptively optimized according to the spatial gradient information about the model resistivity to fine delineate the boundaries of abnormal bodies.In the early stage of inversion execution,a coarse mesh is used for inversion,and the inversion stability is improved by reducing the number of inversion elements.In addition,mesh refinement is performed in the iterative inversion process.The inversion results obtained from the previous mesh are used as the reference and initial models for the next mesh iterative inversion.The step-by-step inversion strategy can ensure that the inversion is performed in the correct direction,improving the inversion stability and results gradually.Synthetic results show that the step-by-step inversion strategy with a Gauss-Newton method for 3D AMT inversion is stable and reliable,which lays a foundation for further practical 3D AMT data inversion.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.50909065 and 51109039)the Major State Basic Research Program of China(973 Program,Grant No.2012CB417002)
文摘A two-dimensional coastal ocean model based on unstructured C-grid is built, in which the momentum equation is discretized on the faces of each cell, and the continuity equation is discretized on the cell. The model is discretized by semi-implicit finite volume method, in that the free surface is semi-implicit and the bottom friction is implicit, thereby removing stability limitations associated with the surface gravity wave and friction. The remaining terms in the momentum equations are discretized explicitly by integral finite volume method and second-order Adams-Bashforth method. Tidal flow in the polar quadrant with known analytic solution is employed to test the proposed model. Finally, the performance of the present model to simulate tidal flow in a geometrically complex domain is examined by simulation of tidal currents in the Pearl River Estuary.
基金supported by the National Natural Science Foundation of China(Nos.91752114 and 11672160)
文摘In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstruction technique in terms of compact stencil is improved to reduce local condition numbers. To further improve the efficiency of computation, the adaptive mesh refinement technique is implemented in the framework of high-order finite volume methods. Mesh refinement and coarsening criteria are chosen to be the indicators for certain flow structures. One important challenge of the adaptive mesh refinement technique on unstructured grids is the dynamic load balancing in parallel computation. To solve this problem, the open-source library p4 est based on the forest of octrees is adopted. Several two-and three-dimensional test cases are computed to verify the accuracy and robustness of the proposed numerical schemes.
文摘Unstructured P2P has power-law link distribution, and the random walk in power-law networks is analyzed. The analysis results show that the probability that a random walker walks through the high degree nodes is high in the power-law network, and the information on the high degree nodes can be easily found through random walk. Random walk spread and random walk search method （RWSS） is proposed based on the analysis result. Simulation results show that RWSS achieves high success rates at low cost and is robust to high degree node failure.
基金supported by the National Basic Research Program of China (2009CB724104)the National Natural Science Foundation of China (90716010)
文摘A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method on triangular grids. The developed HWENO methodology utilizes high-order derivative information to keep WENO re- construction stencils in the von Neumann neighborhood. A simple and efficient technique is also proposed to enhance the smoothness of the existing stencils, making higher-order scheme stable and simplifying the reconstruction process at the same time. The resulting HWENO-based limiters are as compact as the underlying DG schemes and therefore easy to implement. Numerical results for a wide range of flow conditions demonstrate that for DG schemes of up to fourth order of accuracy, the designed HWENO limiters can simul- taneously obtain uniform high order accuracy and sharp, es- sentially non-oscillatory shock transition.
基金supported by the National Natural Science Foundation of China(Grants 51476152,11302213,and 11572336)
文摘An efficient high-order numerical method for supersonic reactive flows is proposed in this article.The reactive source term and convection term are solved separately by splitting scheme.In the reaction step,an adaptive time-step method is presented,which can improve the efficiency greatly.In the convection step,a third-order accurate weighted essentially non-oscillatory（WENO）method is adopted to reconstruct the solution in the unstructured grids.Numerical results show that our new method can capture the correct propagation speed of the detonation wave exactly even in coarse grids,while high order accuracy can be achieved in the smooth region.In addition,the proposed adaptive splitting method can reduce the computational cost greatly compared with the traditional splitting method.
