How to get the rapid and stable inversion results and reconstruct the clear subsurface resistivity structures is a focus problem in current magnetotelluric inversion. A stable solution of an ill-posed inverse problem ...How to get the rapid and stable inversion results and reconstruct the clear subsurface resistivity structures is a focus problem in current magnetotelluric inversion. A stable solution of an ill-posed inverse problem was obtained by the regularization methods in which some desired structures were imposed to stabilize the inverse problem. By the smoothness-constrained model and approximate sensitivity method, the stable subsurface resistivity structures were reconstructed. The synthetic examples show that the smoothness-constrained regularized inversion method is effective and can be reasonable to reconstruct three-dimensional subsurface resistivity structures.展开更多
To comprehensively utilize the valuable geological map,exploration profile,borehole,and geochemical logging data and the knowledge on the formation of the Jinshan Ag-Au deposit for forecasting the exploration targets ...To comprehensively utilize the valuable geological map,exploration profile,borehole,and geochemical logging data and the knowledge on the formation of the Jinshan Ag-Au deposit for forecasting the exploration targets of concealed ore bodies,three-dimensional Mineral Prospectivity Modeling(MPM)of the deposit has been conducted using the weights-of-evidence(WofE)method.Conditional independence between evidence layers was tested,and the outline results using the prediction-volume(P-V)and Student's t-statistic methods for delineating favorable mineralization areas from continuous posterior probability map were critically compared.Four exploration targets delineated ultimately by the Student's t-statistic method for the discovery of minable ore bodies in each of the target areas were discussed in detail.The main conclusions include:(1)three-dimensional modeling of a deposit using multi-source reconnaissance data is useful for MPM in interpreting their relationships with known ore bodies;(2)WofE modeling can be used as a straightforward tool for integrating deposit model and reconnaissance data in MPM;(3)the Student's t-statistic method is more applicable in binarizing the continuous prospectivity map for exploration targeting than the PV approach;and(4)two target areas within high potential to find undiscovered ore bodies were diagnosed to guide future near-mine exploration activities of the Jinshan deposit.展开更多
To capture the nonlinear dynamics and gain evolution in chirped pulse amplification(CPA)systems,the split-step Fourier method and the fourth-order Runge–Kutta method are integrated to iteratively address the generali...To capture the nonlinear dynamics and gain evolution in chirped pulse amplification(CPA)systems,the split-step Fourier method and the fourth-order Runge–Kutta method are integrated to iteratively address the generalized nonlinear Schrödinger equation and the rate equations.However,this approach is burdened by substantial computational demands,resulting in significant time expenditures.In the context of intelligent laser optimization and inverse design,the necessity for numerous simulations further exacerbates this issue,highlighting the need for fast and accurate simulation methodologies.Here,we introduce an end-to-end model augmented with active learning(E2E-AL)with decent generalization through different dedicated embedding methods over various parameters.On an identical computational platform,the artificial intelligence–driven model is 2000 times faster than the conventional simulation method.Benefiting from the active learning strategy,the E2E-AL model achieves decent precision with only two-thirds of the training samples compared with the case without such a strategy.Furthermore,we demonstrate a multi-objective inverse design of the CPA systems enabled by the E2E-AL model.The E2E-AL framework manifests the potential of becoming a standard approach for the rapid and accurate modeling of ultrafast lasers and is readily extended to simulate other complex systems.展开更多
Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to en...Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to engineering requirements, aiming to optimize satellite heat dissipation while considering constraints on static stability, 3D geometric relationships between components, and special component positions. The 3D-SCALO problem is a challenging bilevel combinatorial optimization task, involving the optimization of discrete component assignment variables in the outer layer and continuous component position variables in the inner layer,with both influencing each other. To address this issue, first, a Mixed Integer Programming(MIP) model is proposed, which reformulates the original bilevel problem into a single-level optimization problem, enabling the exploration of a more comprehensive optimization space while avoiding iterative nested optimization. Then, to model the 3D geometric relationships between components within the MIP framework, a linearized 3D Phi-function method is proposed, which handles non-overlapping and safety distance constraints between cuboid components in an explicit and effective way. Subsequently, the Finite-Rectangle Method(FRM) is proposed to manage 3D geometric constraints for complex-shaped components by approximating them with a finite set of cuboids, extending the applicability of the geometric modeling approach. Finally, the feasibility and effectiveness of the proposed MIP model are demonstrated through two numerical examples"and a real-world engineering case, which confirms its suitability for complex-shaped components and real engineering applications.展开更多
Background:Traditional imaging approaches to keratoconus(KCN)have thus far failed to produce a standardized approach for diagnosis.While many diagnostic modalities and metrics exist,none have proven robust enough to b...Background:Traditional imaging approaches to keratoconus(KCN)have thus far failed to produce a standardized approach for diagnosis.While many diagnostic modalities and metrics exist,none have proven robust enough to be considered a gold standard.This study aims to introduce novel metrics to differentiate between KCN and healthy corneas using three-dimensional(3D)measurements of surface area and volume.Methods:This retrospective observational study examined KCN patients along with healthy control patients between the ages of 20 and 79 years old at the University of Maryland,Baltimore.The selected patients underwent a nine-line raster scan anterior segment optical coherence tomography(AS-OCT).ImageJ was used to determine the central 6 mm of each image and each corneal image was then divided into six 1 mm segments.Free-D software was then used to render the nine different images into a 3D model to calculate corneal surface area and volume.A two-tailed Mann-Whitney test was used to assess statistical significance when comparing these subsets.Results:Thirty-three eyes with KCN,along with 33 healthy control,were enrolled.There were statistically significant differences between the healthy and KCN groups in the metric of anterior corneal surface area(13.927 vs.13.991 mm^(2),P=0.046),posterior corneal surface area(14.045 vs.14.173 mm^(2),P<0.001),and volume(8.430 vs.7.773 mm3,P<0.001)within the central 6 mm.Conclusions:3D corneal models derived from AS-OCT can be used to measure anterior corneal surface area,posterior corneal surface area,and corneal volume.All three parameters are statistically different between corneas with KCN and healthy corneas.Further study and application of these parameters may yield new methodologies for the detection of KCN.展开更多
Currently,there are a limited number of dynamic models available for braided composite plates with large overall motions,despite the incorporation of three-dimensional(3D)braided composites into rotating blade compone...Currently,there are a limited number of dynamic models available for braided composite plates with large overall motions,despite the incorporation of three-dimensional(3D)braided composites into rotating blade components.In this paper,a dynamic model of 3D 4-directional braided composite thin plates considering braiding directions is established.Based on Kirchhoff's plate assumptions,the displacement variables of the plate are expressed.By incorporating the braiding directions into the constitutive equation of the braided composites,the dynamic model of the plate considering braiding directions is obtained.The effects of the speeds,braiding directions,and braided angles on the responses of the plate with fixed-axis rotation and translational motion,respectively,are investigated.This paper presents a dynamic theory for calculating the deformation of 3D braided composite structures undergoing both translational and rotational motions.It also provides a simulation method for investigating the dynamic behavior of non-isotropic material plates in various applications.展开更多
Borehole-to-surface electrical imaging （BSEI） uses a line source and a point source to generate a stable electric field in the ground. In order to study the surface potential of anomalies, three-dimensional forward ...Borehole-to-surface electrical imaging （BSEI） uses a line source and a point source to generate a stable electric field in the ground. In order to study the surface potential of anomalies, three-dimensional forward modeling of point and line sources was conducted by using the finite-difference method and the incomplete Cholesky conjugate gradient （ICCG） method. Then, the damping least square method was used in the 3D inversion of the formation resistivity data. Several geological models were considered in the forward modeling and inversion. The forward modeling results suggest that the potentials generated by the two sources have different surface signatures. The inversion data suggest that the low- resistivity anomaly is outlined better than the high-resistivity anomaly. Moreover, when the point source is under the anomaly, the resistivity anomaly boundaries are better outlined than when using a line source.展开更多
The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying ge...The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying geotechnical responses(e.g.consolidation settlement)in a 3D spatial domain.However,traditional 3D numerical model updating approaches are computationally prohibitive and therefore difficult to update the 3D responses in real time.To address these challenges,this study proposes a novel machine learning framework called sparse dictionary learning(T-3D-SDL)for real-time updating of time-varying 3D geotechnical responses.In T-3D-SDL,a concerned dataset(e.g.time-varying 3D settlement)is approximated as a linear superposition of dictionary atoms generated from 3D random FEM analyses.Field monitoring data are then used to identify non-trivial atoms and estimate their weights within a Bayesian framework for model updating and prediction.The proposed approach enables the real-time update of temporally varying settlements with a high 3D spatial resolution and quantified uncertainty as field monitoring data evolve.The proposed approach is illustrated using an embankment construction project.The results show that the proposed approach effectively improves settlement predictions along temporal and 3D spatial dimensions,with minimal latency(e.g.within minutes),as monitoring data appear.In addition,the proposed approach requires only a reasonably small number of 3D FEM model evaluations,avoids the use of widely adopted yet often criticized surrogate models,and effectively addresses the limitations(e.g.computational inefficiency)of existing 3D model updating approaches.展开更多
This paper presents a hybrid model for three-dimensional Geographical Information Systems which is an integration of surface- and volume-based models. The Triangulated Irregular Network (TIN) and octree models are int...This paper presents a hybrid model for three-dimensional Geographical Information Systems which is an integration of surface- and volume-based models. The Triangulated Irregular Network (TIN) and octree models are integrated in this hybrid models. The TIN model works as a surface-based model which mainly serves for surface presentation and visualization. On the other hand, the octree encoding supports volumetric analysis. The designed data structure brings a major advantage in the three-dimensional selective retrieval. This technique increases the efficiency of three-dimensional data operation.展开更多
This paper presents a unique and formal method of quantifying the similarity or distance between sedimentary facies successions from measured sections in outcrop or drilled wells and demonstrates its first application...This paper presents a unique and formal method of quantifying the similarity or distance between sedimentary facies successions from measured sections in outcrop or drilled wells and demonstrates its first application in inverse stratigraphic modeling. A sedimentary facies succession is represented with a string of symbols, or facies codes in its natural vertical order, in which each symbol brings with it one attribute such as thickness for the facies. These strings are called attributed strings. A similarity measure is defined between the attributed strings based on a syntactic pattern-recognition technique. A dynamic programming algorithm is used to calculate the similarity. Inverse stratigraphic modeling aims to generate quantitative 3D facies models based on forward stratigraphic modeling that honors observed datasets. One of the key techniques in inverse stratigraphic modeling is how to quantify the similarity or distance between simulated and observed sedimentary facies successions at data locations in order for the forward model to condition the simulation results to the observed dataset such as measured sections or drilled wells. This quantification technique comparing sedimentary successions is demonstrated in the form of a cost function based on the defined distance in our inverse stratigraphic modeling implemented with forward modeling optimization.展开更多
By using aluminum alloys,the properties of the material in sheet hydroforming were obtained based on the identification of parameters for constitutive models by inverse modeling in which the friction coefficients were...By using aluminum alloys,the properties of the material in sheet hydroforming were obtained based on the identification of parameters for constitutive models by inverse modeling in which the friction coefficients were also considered in 2D and 3D simulations.With consideration of identified simulation parameters by inverse modeling,some key process parameters including tool dimensions and pre-bulging on the forming processes in sheet hydroforming were investigated and optimized.Based on the optimized parameters,the sheet hydroforming process can be analyzed more accurately to improve the robust design.It proves that the results from simulation based on the identified parameters are in good agreement with those from experiments.展开更多
Identifying the unknown geometric and material information of a multi-shield object by analyzing the radiation signature measurements is always an important problem in national and global security. In order to identif...Identifying the unknown geometric and material information of a multi-shield object by analyzing the radiation signature measurements is always an important problem in national and global security. In order to identify the unknown shielding layer thicknesses of a source/shield system with gamma-ray spectra, we have developed a derivative-free inverse radiation transport model based on a differential evolution algorithm with global and local neighbourhoods(IRT-DEGL). In the present paper, the IRT-DEGL model is further extended for estimating the unknown thicknesses with random initial guesses and material mass densities of multi-shielding layers as well as their combinations. Using the detected gamma-ray spectra,the illustration of inverse studies is implemented and the main influence factors for inverse results are also analyzed.展开更多
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.展开更多
Errors inevitably exist in numerical weather prediction （NWP） due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the pred...Errors inevitably exist in numerical weather prediction （NWP） due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the prediction equations can be estimated inversely by using the past data, which are presumed to represent the imperfection of the NWP model （model error, denoted as ME）. In this first paper of a two-part series, an iteration method for obtaining the MEs in past intervals is presented, and the results from testing its convergence in idealized experiments are reported. Moreover, two batches of iteration tests were applied in the global forecast system of the Global and Regional Assimilation and Prediction System （GRAPES-GFS） for July-August 2009 and January-February 2010. The datasets associated with the initial conditions and sea surface temperature （SST） were both based on NCEP （National Centers for Environmental Prediction） FNL （final） data. The results showed that 6th h forecast errors were reduced to 10% of their original value after a 20-step iteration. Then, off-line forecast error corrections were estimated linearly based on the 2-month mean MEs and compared with forecast errors. The estimated error corrections agreed well with the forecast errors, but the linear growth rate of the estimation was steeper than the forecast error. The advantage of this iteration method is that the MEs can provide the foundation for online correction. A larger proportion of the forecast errors can be expected to be canceled out by properly introducing the model error correction into GRAPES-GFS.展开更多
This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zer...This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zero norm solution. The inversion approach mainly employs forward modeling; a depth weight function is introduced into the objective function of the zero norms. Sparse inversion results are obtained by the corresponding optimal mathematical method. To achieve the practical geophysical and geological significance of the results, penalty function is applied to constrain the density values. Results obtained by proposed provide clear boundary depth and density contrast distribution information. The method's accuracy, validity, and reliability are verified by comparing its results with those of synthetic models. To further explain its reliability, a practical gravity data is obtained for a region in Texas, USA is applied. Inversion results for this region are compared with those of previous studies, including a research of logging data in the same area. The depth of salt dome obtained by the inversion method is 4.2 km, which is in good agreement with the 4.4 km value from the logging data. From this, the practicality of the inversion method is also validated.展开更多
Model error is one of the key factors restricting the accuracy of numerical weather prediction （NWP）. Considering the continuous evolution of the atmosphere, the observed data （ignoring the measurement error） can ...Model error is one of the key factors restricting the accuracy of numerical weather prediction （NWP）. Considering the continuous evolution of the atmosphere, the observed data （ignoring the measurement error） can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers＇ equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.展开更多
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.展开更多
The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the tim...The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.展开更多
Two cardiac functional models are constructed in this paper. One is a single current model and the other is a current multipole model. Parameters denoting the properties of these two models are calculated by a least-s...Two cardiac functional models are constructed in this paper. One is a single current model and the other is a current multipole model. Parameters denoting the properties of these two models are calculated by a least-square fit to the measurements using a simulated annealing algorithm. The measured signals are detected at 36 observation nodes by a superconducting quantum interference device （SQUID）. By studying the trends of position, orientation and magnitude of the single current dipole model and the current multipole model in the QRS complex during one time span and comparing the reconstructed magnetocardiography （MCG） of these two cardiac models, we find that the current multipole model is a more appropriate model to represent cardiac electrophysiological activity.展开更多
Three-dimensional(3D)culture systems are becoming increasingly popular due to their ability to mimic tissue-like structures more effectively than the monolayer cultures.In cancer and stem cell research,the natural cel...Three-dimensional(3D)culture systems are becoming increasingly popular due to their ability to mimic tissue-like structures more effectively than the monolayer cultures.In cancer and stem cell research,the natural cell characteristics and architectures are closely mimicked by the 3D cell models.Thus,the 3D cell cultures are promising and suitable systems for various proposes,ranging from disease modeling to drug target identification as well as potential therapeutic substances that may transform our lives.This review provides a comprehensive compendium of recent advancements in culturing cells,in particular cancer and stem cells,using 3D culture techniques.The major approaches highlighted here include cell spheroids,hydrogel embedding,bioreactors,scaffolds,and bioprinting.In addition,the progress of employing 3D cell culture systems as a platform for cancer and stem cell research was addressed,and the prominent studies of 3D cell culture systems were discussed.展开更多
基金Project(20110162120064)supported by Higher School Doctor Subject Special Scientific Research Foundation of ChinaProject(10JJ6059)supported by the Natural Science Foundation of Hunan Province,China
文摘How to get the rapid and stable inversion results and reconstruct the clear subsurface resistivity structures is a focus problem in current magnetotelluric inversion. A stable solution of an ill-posed inverse problem was obtained by the regularization methods in which some desired structures were imposed to stabilize the inverse problem. By the smoothness-constrained model and approximate sensitivity method, the stable subsurface resistivity structures were reconstructed. The synthetic examples show that the smoothness-constrained regularized inversion method is effective and can be reasonable to reconstruct three-dimensional subsurface resistivity structures.
基金financially supported by the Ministry of Science and Technology of China(Nos.2022YFF0801201,2021YFC2900300)the National Natural Science Foundation of China(Nos.41872245,U1911202)the Guangdong Basic and Applied Basic Research Foundation(No.2020A1515010666)。
文摘To comprehensively utilize the valuable geological map,exploration profile,borehole,and geochemical logging data and the knowledge on the formation of the Jinshan Ag-Au deposit for forecasting the exploration targets of concealed ore bodies,three-dimensional Mineral Prospectivity Modeling(MPM)of the deposit has been conducted using the weights-of-evidence(WofE)method.Conditional independence between evidence layers was tested,and the outline results using the prediction-volume(P-V)and Student's t-statistic methods for delineating favorable mineralization areas from continuous posterior probability map were critically compared.Four exploration targets delineated ultimately by the Student's t-statistic method for the discovery of minable ore bodies in each of the target areas were discussed in detail.The main conclusions include:(1)three-dimensional modeling of a deposit using multi-source reconnaissance data is useful for MPM in interpreting their relationships with known ore bodies;(2)WofE modeling can be used as a straightforward tool for integrating deposit model and reconnaissance data in MPM;(3)the Student's t-statistic method is more applicable in binarizing the continuous prospectivity map for exploration targeting than the PV approach;and(4)two target areas within high potential to find undiscovered ore bodies were diagnosed to guide future near-mine exploration activities of the Jinshan deposit.
基金supported by the National Natural Science Foundation of China(Grant Nos.62227821,62025503,and 62205199).
文摘To capture the nonlinear dynamics and gain evolution in chirped pulse amplification(CPA)systems,the split-step Fourier method and the fourth-order Runge–Kutta method are integrated to iteratively address the generalized nonlinear Schrödinger equation and the rate equations.However,this approach is burdened by substantial computational demands,resulting in significant time expenditures.In the context of intelligent laser optimization and inverse design,the necessity for numerous simulations further exacerbates this issue,highlighting the need for fast and accurate simulation methodologies.Here,we introduce an end-to-end model augmented with active learning(E2E-AL)with decent generalization through different dedicated embedding methods over various parameters.On an identical computational platform,the artificial intelligence–driven model is 2000 times faster than the conventional simulation method.Benefiting from the active learning strategy,the E2E-AL model achieves decent precision with only two-thirds of the training samples compared with the case without such a strategy.Furthermore,we demonstrate a multi-objective inverse design of the CPA systems enabled by the E2E-AL model.The E2E-AL framework manifests the potential of becoming a standard approach for the rapid and accurate modeling of ultrafast lasers and is readily extended to simulate other complex systems.
基金supported by the National Natural Science Foundation of China(No.92371206)the Postgraduate Scientific Research Innovation Project of Hunan Province,China(No.CX2023063).
文摘Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to engineering requirements, aiming to optimize satellite heat dissipation while considering constraints on static stability, 3D geometric relationships between components, and special component positions. The 3D-SCALO problem is a challenging bilevel combinatorial optimization task, involving the optimization of discrete component assignment variables in the outer layer and continuous component position variables in the inner layer,with both influencing each other. To address this issue, first, a Mixed Integer Programming(MIP) model is proposed, which reformulates the original bilevel problem into a single-level optimization problem, enabling the exploration of a more comprehensive optimization space while avoiding iterative nested optimization. Then, to model the 3D geometric relationships between components within the MIP framework, a linearized 3D Phi-function method is proposed, which handles non-overlapping and safety distance constraints between cuboid components in an explicit and effective way. Subsequently, the Finite-Rectangle Method(FRM) is proposed to manage 3D geometric constraints for complex-shaped components by approximating them with a finite set of cuboids, extending the applicability of the geometric modeling approach. Finally, the feasibility and effectiveness of the proposed MIP model are demonstrated through two numerical examples"and a real-world engineering case, which confirms its suitability for complex-shaped components and real engineering applications.
文摘Background:Traditional imaging approaches to keratoconus(KCN)have thus far failed to produce a standardized approach for diagnosis.While many diagnostic modalities and metrics exist,none have proven robust enough to be considered a gold standard.This study aims to introduce novel metrics to differentiate between KCN and healthy corneas using three-dimensional(3D)measurements of surface area and volume.Methods:This retrospective observational study examined KCN patients along with healthy control patients between the ages of 20 and 79 years old at the University of Maryland,Baltimore.The selected patients underwent a nine-line raster scan anterior segment optical coherence tomography(AS-OCT).ImageJ was used to determine the central 6 mm of each image and each corneal image was then divided into six 1 mm segments.Free-D software was then used to render the nine different images into a 3D model to calculate corneal surface area and volume.A two-tailed Mann-Whitney test was used to assess statistical significance when comparing these subsets.Results:Thirty-three eyes with KCN,along with 33 healthy control,were enrolled.There were statistically significant differences between the healthy and KCN groups in the metric of anterior corneal surface area(13.927 vs.13.991 mm^(2),P=0.046),posterior corneal surface area(14.045 vs.14.173 mm^(2),P<0.001),and volume(8.430 vs.7.773 mm3,P<0.001)within the central 6 mm.Conclusions:3D corneal models derived from AS-OCT can be used to measure anterior corneal surface area,posterior corneal surface area,and corneal volume.All three parameters are statistically different between corneas with KCN and healthy corneas.Further study and application of these parameters may yield new methodologies for the detection of KCN.
基金Project supported by the National Natural Science Foundation of China(Nos.12372071 and 12372070)the Aeronautical Science Fund of China(No.2022Z055052001)the Foundation of China Scholarship Council(No.202306830079)。
文摘Currently,there are a limited number of dynamic models available for braided composite plates with large overall motions,despite the incorporation of three-dimensional(3D)braided composites into rotating blade components.In this paper,a dynamic model of 3D 4-directional braided composite thin plates considering braiding directions is established.Based on Kirchhoff's plate assumptions,the displacement variables of the plate are expressed.By incorporating the braiding directions into the constitutive equation of the braided composites,the dynamic model of the plate considering braiding directions is obtained.The effects of the speeds,braiding directions,and braided angles on the responses of the plate with fixed-axis rotation and translational motion,respectively,are investigated.This paper presents a dynamic theory for calculating the deformation of 3D braided composite structures undergoing both translational and rotational motions.It also provides a simulation method for investigating the dynamic behavior of non-isotropic material plates in various applications.
基金sponsored by the National Major Project(No.2016ZX05014-001)the National Natural Science Foundation of China(No.41172130 and U1403191)the Fundamental Research Funds for the Central Universities(No.2-9-2015-209)
文摘Borehole-to-surface electrical imaging （BSEI） uses a line source and a point source to generate a stable electric field in the ground. In order to study the surface potential of anomalies, three-dimensional forward modeling of point and line sources was conducted by using the finite-difference method and the incomplete Cholesky conjugate gradient （ICCG） method. Then, the damping least square method was used in the 3D inversion of the formation resistivity data. Several geological models were considered in the forward modeling and inversion. The forward modeling results suggest that the potentials generated by the two sources have different surface signatures. The inversion data suggest that the low- resistivity anomaly is outlined better than the high-resistivity anomaly. Moreover, when the point source is under the anomaly, the resistivity anomaly boundaries are better outlined than when using a line source.
基金supported by a grant from the Research Grant Council of Hong Kong Special Administrative Region(Project No.11207724).
文摘The development of digital twins for geotechnical structures necessitates the real-time updates of threedimensional(3D)virtual models(e.g.numerical finite element method(FEM)model)to accurately predict time-varying geotechnical responses(e.g.consolidation settlement)in a 3D spatial domain.However,traditional 3D numerical model updating approaches are computationally prohibitive and therefore difficult to update the 3D responses in real time.To address these challenges,this study proposes a novel machine learning framework called sparse dictionary learning(T-3D-SDL)for real-time updating of time-varying 3D geotechnical responses.In T-3D-SDL,a concerned dataset(e.g.time-varying 3D settlement)is approximated as a linear superposition of dictionary atoms generated from 3D random FEM analyses.Field monitoring data are then used to identify non-trivial atoms and estimate their weights within a Bayesian framework for model updating and prediction.The proposed approach enables the real-time update of temporally varying settlements with a high 3D spatial resolution and quantified uncertainty as field monitoring data evolve.The proposed approach is illustrated using an embankment construction project.The results show that the proposed approach effectively improves settlement predictions along temporal and 3D spatial dimensions,with minimal latency(e.g.within minutes),as monitoring data appear.In addition,the proposed approach requires only a reasonably small number of 3D FEM model evaluations,avoids the use of widely adopted yet often criticized surrogate models,and effectively addresses the limitations(e.g.computational inefficiency)of existing 3D model updating approaches.
文摘This paper presents a hybrid model for three-dimensional Geographical Information Systems which is an integration of surface- and volume-based models. The Triangulated Irregular Network (TIN) and octree models are integrated in this hybrid models. The TIN model works as a surface-based model which mainly serves for surface presentation and visualization. On the other hand, the octree encoding supports volumetric analysis. The designed data structure brings a major advantage in the three-dimensional selective retrieval. This technique increases the efficiency of three-dimensional data operation.
基金financially was supported by Colorado School of Minessupported by the Science and Technology Ministry of China (2016ZX05033003)+1 种基金China Academy of Sciences (XDA14010204)Sinopec (G5800-15-ZS-KJB016)
文摘This paper presents a unique and formal method of quantifying the similarity or distance between sedimentary facies successions from measured sections in outcrop or drilled wells and demonstrates its first application in inverse stratigraphic modeling. A sedimentary facies succession is represented with a string of symbols, or facies codes in its natural vertical order, in which each symbol brings with it one attribute such as thickness for the facies. These strings are called attributed strings. A similarity measure is defined between the attributed strings based on a syntactic pattern-recognition technique. A dynamic programming algorithm is used to calculate the similarity. Inverse stratigraphic modeling aims to generate quantitative 3D facies models based on forward stratigraphic modeling that honors observed datasets. One of the key techniques in inverse stratigraphic modeling is how to quantify the similarity or distance between simulated and observed sedimentary facies successions at data locations in order for the forward model to condition the simulation results to the observed dataset such as measured sections or drilled wells. This quantification technique comparing sedimentary successions is demonstrated in the form of a cost function based on the defined distance in our inverse stratigraphic modeling implemented with forward modeling optimization.
基金Project(9901351) supported by the Hydromechanical Deep Drawing without a Draw DieProject(1057001) supported by the National Natural Science Foundation of China
文摘By using aluminum alloys,the properties of the material in sheet hydroforming were obtained based on the identification of parameters for constitutive models by inverse modeling in which the friction coefficients were also considered in 2D and 3D simulations.With consideration of identified simulation parameters by inverse modeling,some key process parameters including tool dimensions and pre-bulging on the forming processes in sheet hydroforming were investigated and optimized.Based on the optimized parameters,the sheet hydroforming process can be analyzed more accurately to improve the robust design.It proves that the results from simulation based on the identified parameters are in good agreement with those from experiments.
基金supported by the National Natural Science Foundation of China(Nos.11605163 and 21504085)the China Academy of Engineering Physics Foundation for Development of Science and Technology(No.201580103014 and No.2015B0301063)+1 种基金the Foundation for Special Talents in China Academy of Engineering Physics(No.TP201502-3)the Sichuan Science and Technology Development Foundation for Young Scientists(No.2017Q0050)
文摘Identifying the unknown geometric and material information of a multi-shield object by analyzing the radiation signature measurements is always an important problem in national and global security. In order to identify the unknown shielding layer thicknesses of a source/shield system with gamma-ray spectra, we have developed a derivative-free inverse radiation transport model based on a differential evolution algorithm with global and local neighbourhoods(IRT-DEGL). In the present paper, the IRT-DEGL model is further extended for estimating the unknown thicknesses with random initial guesses and material mass densities of multi-shielding layers as well as their combinations. Using the detected gamma-ray spectra,the illustration of inverse studies is implemented and the main influence factors for inverse results are also analyzed.
基金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.
基金funded by the National Natural Science Foundation Science Fund for Youth (Grant No.41405095)the Key Projects in the National Science and Technology Pillar Program during the Twelfth Fiveyear Plan Period (Grant No.2012BAC22B02)the National Natural Science Foundation Science Fund for Creative Research Groups (Grant No.41221064)
文摘Errors inevitably exist in numerical weather prediction （NWP） due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the prediction equations can be estimated inversely by using the past data, which are presumed to represent the imperfection of the NWP model （model error, denoted as ME）. In this first paper of a two-part series, an iteration method for obtaining the MEs in past intervals is presented, and the results from testing its convergence in idealized experiments are reported. Moreover, two batches of iteration tests were applied in the global forecast system of the Global and Regional Assimilation and Prediction System （GRAPES-GFS） for July-August 2009 and January-February 2010. The datasets associated with the initial conditions and sea surface temperature （SST） were both based on NCEP （National Centers for Environmental Prediction） FNL （final） data. The results showed that 6th h forecast errors were reduced to 10% of their original value after a 20-step iteration. Then, off-line forecast error corrections were estimated linearly based on the 2-month mean MEs and compared with forecast errors. The estimated error corrections agreed well with the forecast errors, but the linear growth rate of the estimation was steeper than the forecast error. The advantage of this iteration method is that the MEs can provide the foundation for online correction. A larger proportion of the forecast errors can be expected to be canceled out by properly introducing the model error correction into GRAPES-GFS.
基金supported by the Development of airborne gravity gradiometer(No.2017YFC0601601)open subject of Key Laboratory of Petroleum Resources Research,Institute of Geology and Geophysics,Chinese Academy of Sciences(No.KLOR2018-8)
文摘This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zero norm solution. The inversion approach mainly employs forward modeling; a depth weight function is introduced into the objective function of the zero norms. Sparse inversion results are obtained by the corresponding optimal mathematical method. To achieve the practical geophysical and geological significance of the results, penalty function is applied to constrain the density values. Results obtained by proposed provide clear boundary depth and density contrast distribution information. The method's accuracy, validity, and reliability are verified by comparing its results with those of synthetic models. To further explain its reliability, a practical gravity data is obtained for a region in Texas, USA is applied. Inversion results for this region are compared with those of previous studies, including a research of logging data in the same area. The depth of salt dome obtained by the inversion method is 4.2 km, which is in good agreement with the 4.4 km value from the logging data. From this, the practicality of the inversion method is also validated.
基金Project supported by the Special Scientific Research Project for Public Interest(Grant No.GYHY201206009)the Fundamental Research Funds for the Central Universities,China(Grant Nos.lzujbky-2012-13 and lzujbky-2013-11)the National Basic Research Program of China(Grant Nos.2012CB955902 and 2013CB430204)
文摘Model error is one of the key factors restricting the accuracy of numerical weather prediction （NWP）. Considering the continuous evolution of the atmosphere, the observed data （ignoring the measurement error） can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers＇ equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.
基金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.
基金supported by the Joint Fund of Seismological Science(Grant No.U1839206)the National R&D Program on Monitoring,Early Warning and Prevention of Major Natural Disaster(Grant No.2017YFC1500301)+2 种基金supported by IGGCAS Research Start-up Funds(Grant No.E0515402)National Natural Science Foundation of China(Grant No.E1115401)supported by National Natural Science Foundation of China(Grant No.11971258).
文摘The nearly analytic discrete(NAD)method is a kind of finite difference method with advantages of high accuracy and stability.Previous studies have investigated the NAD method for simulating wave propagation in the time-domain.This study applies the NAD method to solving three-dimensional(3D)acoustic wave equations in the frequency-domain.This forward modeling approach is then used as the“engine”for implementing 3D frequency-domain full waveform inversion(FWI).In the numerical modeling experiments,synthetic examples are first given to show the superiority of the NAD method in forward modeling compared with traditional finite difference methods.Synthetic 3D frequency-domain FWI experiments are then carried out to examine the effectiveness of the proposed methods.The inversion results show that the NAD method is more suitable than traditional methods,in terms of computational cost and stability,for 3D frequency-domain FWI,and represents an effective approach for inversion of subsurface model structures.
基金Project supported by the State Key Development Program for Basic Research of China (Grant No 2006CB601007)the National Natural Science Foundation of China (Grant No 10674006)the National High Technology Research and Development Program of China (Grant No 2007AA03Z238)
文摘Two cardiac functional models are constructed in this paper. One is a single current model and the other is a current multipole model. Parameters denoting the properties of these two models are calculated by a least-square fit to the measurements using a simulated annealing algorithm. The measured signals are detected at 36 observation nodes by a superconducting quantum interference device （SQUID）. By studying the trends of position, orientation and magnitude of the single current dipole model and the current multipole model in the QRS complex during one time span and comparing the reconstructed magnetocardiography （MCG） of these two cardiac models, we find that the current multipole model is a more appropriate model to represent cardiac electrophysiological activity.
文摘Three-dimensional(3D)culture systems are becoming increasingly popular due to their ability to mimic tissue-like structures more effectively than the monolayer cultures.In cancer and stem cell research,the natural cell characteristics and architectures are closely mimicked by the 3D cell models.Thus,the 3D cell cultures are promising and suitable systems for various proposes,ranging from disease modeling to drug target identification as well as potential therapeutic substances that may transform our lives.This review provides a comprehensive compendium of recent advancements in culturing cells,in particular cancer and stem cells,using 3D culture techniques.The major approaches highlighted here include cell spheroids,hydrogel embedding,bioreactors,scaffolds,and bioprinting.In addition,the progress of employing 3D cell culture systems as a platform for cancer and stem cell research was addressed,and the prominent studies of 3D cell culture systems were discussed.
