A finite element reconstruction algorithm for ultrasound tomography based on the Helmholtz equation in frequency domain is presented to monitor the grouting defects in reinforced concrete structures.In this algorithm,...A finite element reconstruction algorithm for ultrasound tomography based on the Helmholtz equation in frequency domain is presented to monitor the grouting defects in reinforced concrete structures.In this algorithm,a hybrid regularizations-based iterative Newton method is implemented to provide stable inverse solutions.Furthermore,a dual mesh scheme and an adjoint method are adopted to reduce the computation cost and improve the efficiency of reconstruction.Simultaneous reconstruction of both acoustic velocity and attenuation coefficient for a reinforced concrete model is achieved with multiple frequency data.The algorithm is evaluated with numerical simulation under various practical scenarios including varied transmission/receiving modes,different noise levels,different source/detector numbers,and different contrast levels between the heterogeneity and background region.Results obtained suggest that the algorithm is insensitive to noise,and the reconstructions are quantitatively accurate in terms of the location,size and acoustic properties of the target over a range of contrast levels.展开更多
Computational models provide additional tools for studying the brain,however,many techniques are currently disconnected from each other.There is a need for new computational approaches that span the range of physics o...Computational models provide additional tools for studying the brain,however,many techniques are currently disconnected from each other.There is a need for new computational approaches that span the range of physics operating in the brain.In this review paper,we offer some new perspectives on how the embedded element method can fill this gap and has the potential to connect a myriad of modeling genre.The embedded element method is a mesh superposition technique used within finite element analysis.This method allows for the incorporation of axonal fiber tracts to be explicitly represented.Here,we explore the use of the approach beyond its original goal of predicting axonal strain in brain injury.We explore the potential application of the embedded element method in areas of electrophysiology,neurodegeneration,neuropharmacology and mechanobiology.We conclude that this method has the potential to provide us with an integrated computational framework that can assist in developing improved diagnostic tools and regeneration technologies.展开更多
In a recent paper, an efficient semi-implicit finite element scheme for 2-dimensional tidal flow computations is proposed. In that scheme, each term of the governing equations, rather than each dependent variable, is ...In a recent paper, an efficient semi-implicit finite element scheme for 2-dimensional tidal flow computations is proposed. In that scheme, each term of the governing equations, rather than each dependent variable, is ex- panded in terms of the unknown nodal values. Simpson's rule ix used for numerical integration to make the mass matrix diagonal. The friction terms are represented semi-implicitly to improve stability, but no additional compu- tational effort is required. The shortcomings of this scheme are that the time-stepping scheme is only first-order ae- curate and artificial smoothing is required to control the numerical noise. In this paper, the previous scheme is im- proved by including the eddy viscosity terms in the governing equations to replace artificial smoothing in noise con- trol and the time-stepping scheme is modified to make it second-order accurate. These improvements can be achieved with only a slight increase in computational effort. The test cases used previously to validate the former scheme are again employed to test the present scheme.展开更多
The wave-induced seepage force is investigated on marine structures resting on or buried in the seabed.The bed is modelled as a poroelastic medium containing a nearly saturated water.The governing equations are solved...The wave-induced seepage force is investigated on marine structures resting on or buried in the seabed.The bed is modelled as a poroelastic medium containing a nearly saturated water.The governing equations are solved with Finite Element Method.For a pipeline buried in the seabed,agreement between the present numerical results and that of Cheng H.D.(1986)is quite satisfactory.展开更多
Region partition(RP) is the key technique to the finite element parallel computing(FEPC),and its performance has a decisive influence on the entire process of analysis and computation.The performance evaluation index ...Region partition(RP) is the key technique to the finite element parallel computing(FEPC),and its performance has a decisive influence on the entire process of analysis and computation.The performance evaluation index of RP method for the three-dimensional finite element model(FEM) has been given.By taking the electric field of aluminum reduction cell(ARC) as the research object,the performance of two classical RP methods,which are Al-NASRA and NGUYEN partition(ANP) algorithm and the multi-level partition(MLP) method,has been analyzed and compared.The comparison results indicate a sound performance of ANP algorithm,but to large-scale models,the computing time of ANP algorithm increases notably.This is because the ANP algorithm determines only one node based on the minimum weight and just adds the elements connected to the node into the sub-region during each iteration.To obtain the satisfied speed and the precision,an improved dynamic self-adaptive ANP(DSA-ANP) algorithm has been proposed.With consideration of model scale,complexity and sub-RP stage,the improved algorithm adaptively determines the number of nodes and selects those nodes with small enough weight,and then dynamically adds these connected elements.The proposed algorithm has been applied to the finite element analysis(FEA) of the electric field simulation of ARC.Compared with the traditional ANP algorithm,the computational efficiency of the proposed algorithm has been shortened approximately from 260 s to 13 s.This proves the superiority of the improved algorithm on computing time performance.展开更多
In this paper,a bioheat model of temperature distribution in the human eye is studied,the mathematical formulation of this model is described using adequate mathematical tools.The existence and the uniqueness of the s...In this paper,a bioheat model of temperature distribution in the human eye is studied,the mathematical formulation of this model is described using adequate mathematical tools.The existence and the uniqueness of the solution of this problem is proven and four algorithms based on finite element method approximation and domain decomposition methods are presented in details.The validation of all algorithm is done using a numerical application for an example where the analytical solution is known.The properties and parameters reported in the open literature for the human eye are used to approximate numerically the temperature for bioheat model by finite element approximation and nonoverlapping domain decomposition method.The obtained results that are verified using the experimental results recorded in the literature revealed a better accuracy by the use of algorithm proposed.展开更多
基金Project(31200748)supported by the National Natural Science Foundation of China
文摘A finite element reconstruction algorithm for ultrasound tomography based on the Helmholtz equation in frequency domain is presented to monitor the grouting defects in reinforced concrete structures.In this algorithm,a hybrid regularizations-based iterative Newton method is implemented to provide stable inverse solutions.Furthermore,a dual mesh scheme and an adjoint method are adopted to reduce the computation cost and improve the efficiency of reconstruction.Simultaneous reconstruction of both acoustic velocity and attenuation coefficient for a reinforced concrete model is achieved with multiple frequency data.The algorithm is evaluated with numerical simulation under various practical scenarios including varied transmission/receiving modes,different noise levels,different source/detector numbers,and different contrast levels between the heterogeneity and background region.Results obtained suggest that the algorithm is insensitive to noise,and the reconstructions are quantitatively accurate in terms of the location,size and acoustic properties of the target over a range of contrast levels.
基金support provided by Computational Fluid Dynamics Research Corporation(CFDRC)under a sub-contract funded by the Department of Defense,Department of Health Program through contract W81XWH-14-C-0045
文摘Computational models provide additional tools for studying the brain,however,many techniques are currently disconnected from each other.There is a need for new computational approaches that span the range of physics operating in the brain.In this review paper,we offer some new perspectives on how the embedded element method can fill this gap and has the potential to connect a myriad of modeling genre.The embedded element method is a mesh superposition technique used within finite element analysis.This method allows for the incorporation of axonal fiber tracts to be explicitly represented.Here,we explore the use of the approach beyond its original goal of predicting axonal strain in brain injury.We explore the potential application of the embedded element method in areas of electrophysiology,neurodegeneration,neuropharmacology and mechanobiology.We conclude that this method has the potential to provide us with an integrated computational framework that can assist in developing improved diagnostic tools and regeneration technologies.
文摘In a recent paper, an efficient semi-implicit finite element scheme for 2-dimensional tidal flow computations is proposed. In that scheme, each term of the governing equations, rather than each dependent variable, is ex- panded in terms of the unknown nodal values. Simpson's rule ix used for numerical integration to make the mass matrix diagonal. The friction terms are represented semi-implicitly to improve stability, but no additional compu- tational effort is required. The shortcomings of this scheme are that the time-stepping scheme is only first-order ae- curate and artificial smoothing is required to control the numerical noise. In this paper, the previous scheme is im- proved by including the eddy viscosity terms in the governing equations to replace artificial smoothing in noise con- trol and the time-stepping scheme is modified to make it second-order accurate. These improvements can be achieved with only a slight increase in computational effort. The test cases used previously to validate the former scheme are again employed to test the present scheme.
文摘The wave-induced seepage force is investigated on marine structures resting on or buried in the seabed.The bed is modelled as a poroelastic medium containing a nearly saturated water.The governing equations are solved with Finite Element Method.For a pipeline buried in the seabed,agreement between the present numerical results and that of Cheng H.D.(1986)is quite satisfactory.
基金Project(61273187)supported by the National Natural Science Foundation of ChinaProject(61321003)supported by the Foundation for Innovative Research Groups of the National Natural Science Foundation of China
文摘Region partition(RP) is the key technique to the finite element parallel computing(FEPC),and its performance has a decisive influence on the entire process of analysis and computation.The performance evaluation index of RP method for the three-dimensional finite element model(FEM) has been given.By taking the electric field of aluminum reduction cell(ARC) as the research object,the performance of two classical RP methods,which are Al-NASRA and NGUYEN partition(ANP) algorithm and the multi-level partition(MLP) method,has been analyzed and compared.The comparison results indicate a sound performance of ANP algorithm,but to large-scale models,the computing time of ANP algorithm increases notably.This is because the ANP algorithm determines only one node based on the minimum weight and just adds the elements connected to the node into the sub-region during each iteration.To obtain the satisfied speed and the precision,an improved dynamic self-adaptive ANP(DSA-ANP) algorithm has been proposed.With consideration of model scale,complexity and sub-RP stage,the improved algorithm adaptively determines the number of nodes and selects those nodes with small enough weight,and then dynamically adds these connected elements.The proposed algorithm has been applied to the finite element analysis(FEA) of the electric field simulation of ARC.Compared with the traditional ANP algorithm,the computational efficiency of the proposed algorithm has been shortened approximately from 260 s to 13 s.This proves the superiority of the improved algorithm on computing time performance.
文摘In this paper,a bioheat model of temperature distribution in the human eye is studied,the mathematical formulation of this model is described using adequate mathematical tools.The existence and the uniqueness of the solution of this problem is proven and four algorithms based on finite element method approximation and domain decomposition methods are presented in details.The validation of all algorithm is done using a numerical application for an example where the analytical solution is known.The properties and parameters reported in the open literature for the human eye are used to approximate numerically the temperature for bioheat model by finite element approximation and nonoverlapping domain decomposition method.The obtained results that are verified using the experimental results recorded in the literature revealed a better accuracy by the use of algorithm proposed.