By means of difference calculus, a multifold analog of Gould-Hsu inversion is established. Asconcrete examples, some multifold inverse relations are presented. By developing a standard prooftechnique for combinatorial...By means of difference calculus, a multifold analog of Gould-Hsu inversion is established. Asconcrete examples, some multifold inverse relations are presented. By developing a standard prooftechnique for combinatorial identities, we demonstrate Mohanty-Hanna's convolution formulas almostmechanically. In addition, an application of reciprocal formulas to interpolation process is alsosketched.展开更多
In this study,an inverse design framework was established to find lightweight honeycomb structures(HCSs)with high impact resistance.The hybrid HCS,composed of re-entrant(RE)and elliptical annular re-entrant(EARE)honey...In this study,an inverse design framework was established to find lightweight honeycomb structures(HCSs)with high impact resistance.The hybrid HCS,composed of re-entrant(RE)and elliptical annular re-entrant(EARE)honeycomb cells,was created by constructing arrangement matrices to achieve structural lightweight.The machine learning(ML)framework consisted of a neural network(NN)forward regression model for predicting impact resistance and a multi-objective optimization algorithm for generating high-performance designs.The surrogate of the local design space was initially realized by establishing the NN in the small sample dataset,and the active learning strategy was used to continuously extended the local optimal design until the model converged in the global space.The results indicated that the active learning strategy significantly improved the inference capability of the NN model in unknown design domains.By guiding the iteration direction of the optimization algorithm,lightweight designs with high impact resistance were identified.The energy absorption capacity of the optimal design reached 94.98%of the EARE honeycomb,while the initial peak stress and mass decreased by 28.85%and 19.91%,respectively.Furthermore,Shapley Additive Explanations(SHAP)for global explanation of the NN indicated a strong correlation between the arrangement mode of HCS and its impact resistance.By reducing the stiffness of the cells at the top boundary of the structure,the initial impact damage sustained by the structure can be significantly improved.Overall,this study proposed a general lightweight design method for array structures under impact loads,which is beneficial for the widespread application of honeycomb-based protective structures.展开更多
Co-assembling chiral molecules with achiral compounds via non-covalent interactions like areneperfluoroarene(AP) interactions offers an effective approach for fabricating chiral functional materials.Herein,chiral mole...Co-assembling chiral molecules with achiral compounds via non-covalent interactions like areneperfluoroarene(AP) interactions offers an effective approach for fabricating chiral functional materials.Herein,chiral molecules L/D-PF1 and L/D-PF2 with pyrene groups were synthesized and its chiroptical properties upon co-assembly with achiral compound octafluoronaphthalene(OFN) through AP interaction were systemically studied.The co-assembly of L/D-PF1/OFN and L/D-PF2/OFN exhibited distinct chiroptical properties such as circular dichroism(CD) and circularly polarized luminescence(CPL) signals.Chirality transfer from the chirality center of L/D-PF1 and L/D-PF2 to the achiral OFN and chiral amplification were successfully achieved.Besides,no significant CPL signal was observed in the self-assembly of L/DPF1 or L/D-PF2 while co-assembly with OFN exhibited obvious CPL amplification induced by AP interaction.Notably,a reversal CD signal and CPL signal could be observed in L/D-PF2/OFN when the molar ratio changed from 1:1 to 1:2 while not found in L/D-PF1/OFN,indicating that that minor structural changes of molecules could cause large changes in assembly.In addition,a series of computational calculations were conducted to verify the AP interaction between L-PF1/L-PF2 and OFN.This work demonstrated that arene-perfluoroarene interaction could drive chiral transfer,chiral amplification and chiral inversion and provided a new method for the preparation of chiroptical materials.展开更多
Patient-specific finite element analysis(FEA)is a promising tool for noninvasive quantification of cardiac and vascular structural mechanics in vivo.However,inverse material property identification using FEA,which req...Patient-specific finite element analysis(FEA)is a promising tool for noninvasive quantification of cardiac and vascular structural mechanics in vivo.However,inverse material property identification using FEA,which requires iteratively solving nonlinear hyperelasticity problems,is computationally expensive which limits the ability to provide timely patient-specific insights to clinicians.In this study,we present an inverse material parameter identification strategy that integrates deep neural networks(DNNs)with FEA,namely inverse DNN-FEA.In this framework,a DNN encodes the spatial distribution of material parameters and effectively regularizes the inverse solution,which aims to reduce susceptibility to local optima that often arise in heterogeneous nonlinear hyperelastic problems.Consequently,inverse DNN-FEA enables identification of material parameters at the element level.For validation,we applied DNN-FEA to identify four spatially varying passive Holzapfel-Ogden material parameters of the left ventricular myocardium in synthetic benchmark cases with a clinically-derived geometry.To evaluate the benefit of DNN integration,a baseline FEA-only solver implemented in PyTorch was used for comparison.Results demonstrated that DNN-FEA achieved substantially lower average errors in parameter identification compared to FEA(case 1,DNN-FEA:0.37%~2.15%vs.FEA:2.64%~12.91%).The results also demonstrate that the same DNN architecture is capable of identifying a different spatial material property distribution(case 2,DNN-FEA:0.03%~0.60%vs.FEA:0.93%~16.25%).These findings suggest that DNN-FEA provides an accurate framework for inverse identification of heterogeneous myocardial material properties.This approach may facilitate future applications in patient-specific modeling based on in vivo clinical imaging and could be extended to other biomechanical simulation problems.展开更多
Knowledge of the statistical characteristics of inversions and their effects on aerosols under different large-scale synoptic circulations is important for studying and modeling the diffusion of pollutants in the boun...Knowledge of the statistical characteristics of inversions and their effects on aerosols under different large-scale synoptic circulations is important for studying and modeling the diffusion of pollutants in the boundary layer. Based on results gen- erated using the self-organizing map （SOM） weather classification method, this study compares the statistical characteristics of surface-based inversions （SBIs） and elevated inversions （EIs）, and quantitatively evaluates the effect of SBIs on aerosol condensation nuclei （CN） concentrations and the relationship between temperature gradients and aerosols for six prevailing synoptic patterns over the the Southern Great Plains （SGP） site during 2001-10. Large-scale synoptic patterns strongly influ- ence the statistical characteristics of inversions and the accumulation of aerosols in the low-level atmosphere. The activity, frequency, intensity, and vertical distribution of inversions are significantly different among these synoptic patterns. The verti- cal distribution of inversions varies diurnally and is significantly different among the different synoptic patterns. Anticyclonic patterns affect the accumulation of aerosols near the ground more strongly than cyclonic patterns. Mean aerosol CN con- centrations increase during SBIs compared to no inversion cases by 16.1%, 22.6%, 24.5%, 58.7%, 29.8% and 23.7% for the six synoptic patterns. This study confirms that there is a positive correlation between temperature gradients and aerosol CN concentrations near the ground at night under similar large-scale synoptic patterns. The relationship is different for different synoptic patterns and can be described by linear functions. These findings suggest that large-scale synoptic patterns change the static stability of the atmosphere and inversions in the lower atmosphere, thereby influencing the diffusion of aerosols near the ground.展开更多
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.展开更多
High-quality and continuous radiosonde, aerosol and surface meteorology datasets are used to investigate the statistical characteristics of meteorological parameters and their effects on aerosols. The data were collec...High-quality and continuous radiosonde, aerosol and surface meteorology datasets are used to investigate the statistical characteristics of meteorological parameters and their effects on aerosols. The data were collected at the Atmospheric Radiation Measurement Southern Great Plains climate research facility during 2000–15. The parameters and vertical distribution of temperature inversion layers were found to have strong diurnal and seasonal changes. For surface-based temperature inversion (SBI), the mean frequency and depth of temperature inversion layers were 39.4% and 198 m, respectively. The temperature difference between the top and bottom of SBI was 4.8℃, and so the temperature gradient was 2.4℃(100 m)^-1. The detailed vertical distributions of temperature inversion had been determined, and only the temperature inversion layers below 1000 m showed diurnal and seasonal variations. Mean surface aerosol number concentrations increased by 43.0%, 21.9% and 49.2% when SBIs were present at 0530, 1730 and 2330 LST, respectively. The effect of SBI on surface aerosol concentration was weakest in summer (18.1%) and strongest in winter (58.4%). During elevated temperature inversion events, there was no noticeable difference in surface aerosol number concentrations. Temperature differences and temperature gradients across SBIs correlated fairly well with aerosol number concentrations, especially for temperature gradients. The vertical distribution of aerosol optical properties with and without temperature inversions was different. Surface aerosol measurements were representative of the air within (below), but not above, SBIs and EIs. These results provide a basis for developing a boundary layer aerosol accumulation model and for improving radiative transfer models in the lower atmosphere.展开更多
We conducted rapid inversions of rupture process for the 2023 earthquake doublet occurred in SE Türkiye,the first with a magnitude of M_(W)7.8 and the second with a magnitude of M_(W)7.6,using teleseismic and str...We conducted rapid inversions of rupture process for the 2023 earthquake doublet occurred in SE Türkiye,the first with a magnitude of M_(W)7.8 and the second with a magnitude of M_(W)7.6,using teleseismic and strong-motion data.The teleseismic rupture models of the both events were obtained approximately 88 and 55 minutes after their occurrences,respectively.The rupture models indicated that the first event was an asymmetric bilateral event with ruptures mainly propagating to the northeast,while the second one was a unilateral event with ruptures propagating to the west.This information could be useful in locating the meizoseismal areas.Compared with teleseismic models,the strong-motion models showed relatively higher resolution.A noticeable difference was found for the M_(W)7.6 earthquake,for which the strong-motion models shows a bilateral event,rather than a unilateral event,but the dominant rupture direction is still westward.Nevertheless,all strong-motion models are consistent with the teleseismic models in terms of magnitudes,durations,and dominant rupture directions.This suggests that both teleseismic and strong-motion data can be used for fast determination of major source characteristics.In contrast,the strong-motion data would be preferable in future emergency responses since they are recorded earlier and have a better resolution ability on the source ruptures.展开更多
In this paper, by means of Gould-Hsu inverse series relations, we establish several Gould-Hsu inversion chains. As consequence, some new transformation formulae as well as some famous hypergeometric series identities ...In this paper, by means of Gould-Hsu inverse series relations, we establish several Gould-Hsu inversion chains. As consequence, some new transformation formulae as well as some famous hypergeometric series identities are derived.展开更多
The duplicate form of the generalized Gould-Hsu inversions has been obtained by Shi and Zhang. In this paper, we present a simple proof of this duplicate form. With the same method, we construct the duplicate form of ...The duplicate form of the generalized Gould-Hsu inversions has been obtained by Shi and Zhang. In this paper, we present a simple proof of this duplicate form. With the same method, we construct the duplicate form of the generalized Carlitz inversions. Using this duplicate form, we obtain several terminating basic hypergeometric identities and some limiting cases.展开更多
It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regio...It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regional earthquake tomography, but seldom are all three classes of model parameters updated simultaneously. This is mainly due to the trade-off between the different types of model parameters and the lack of different seismic phases to constrain the model parameters. Using a spherical-coordinate ray tracing algorithm for first and later（primary reflected） arrival tracing algorithm in combination with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model parameters in regional or global tomographic studies. In this paper we incorporate the multistage irregular shortest-path ray tracing algorithm（in a spherical coordinate system） with a subspace inversion solver to formulate a simultaneous inversion algorithm for triple model parameters updating using direct and later arrival time information.Comparison tests for two sets of data（noise free and added noise） indicate that the new triple-class parameter inversion algorithm is capable of obtaining nearly the same results as the double-class parameter inversion scheme. Furthermore,the proposed multi-parameter type inversion method is not sensitive to a modest level of picking error in the traveltime data, and also performs well with a relatively large uncertainty in earthquake hypocentral locations. This shows it to be a feasible and promising approach in regional or global tomographic applications.展开更多
In the paper, a 2D symmetrical anisotropic medium whose strike agrees with one of the horizontal principal axes is considered to develop a corresponding inversion technique. In the specified conditions, if we assume a...In the paper, a 2D symmetrical anisotropic medium whose strike agrees with one of the horizontal principal axes is considered to develop a corresponding inversion technique. In the specified conditions, if we assume an equivalent conductivity anisotropy in both the vertical and dipping directions, i.e., σzz=σyy, the differential equations obtained are formally the same as that for TE and TM modes in the 2D isotropic geoelectrical media. The same inversion technique as that in the 2D isotropic media can be employed to obtain the anisotropic conductivities. It means that the TE and TM inversion results in the isotropic media can be respectively thought as the resistivities in the two principal directions of the symmetrically anisotropic media, which has offered a new approach and a theoretical guidance for interpreting magnetotelluric data. And the inversion technique developed here is used to test the magnetotelluric data in the area of Tianzhu and Yongdeng in Gansu Province, so that the crust anisotropic geoelectrical structures in this region can be obtained.展开更多
The Bayesian inversion method is a stochastic approach based on the Bayesian theory.With the development of sampling algorithms and computer technologies,the Bayesian inversion method has been widely used in geophysic...The Bayesian inversion method is a stochastic approach based on the Bayesian theory.With the development of sampling algorithms and computer technologies,the Bayesian inversion method has been widely used in geophysical inversion problems.In this study,we conduct inversion experiments using crosshole seismic travel-time data to examine the characteristics and performance of the stochastic Bayesian inversion based on the Markov chain Monte Carlo sampling scheme and the traditional deterministic inversion with Tikhonov regularization.Velocity structures with two different spatial variations are considered,one with a chessboard pattern and the other with an interface mimicking the Mohorovicicdiscontinuity(Moho).Inversions are carried out with different scenarios of model discretization and source–receiver configurations.Results show that the Bayesian method yields more robust single-model estimations than the deterministic method,with smaller model errors.In addition,the Bayesian method provides the posterior probabilistic distribution function of the model space,which can help us evaluate the quality of the inversion result.展开更多
HF (high frequency) radar sounder technology has been developed for several missions of Mars surface/subsurface exploration. This paper presents a model of rough surface and stratified sub-surfaces to describe the mul...HF (high frequency) radar sounder technology has been developed for several missions of Mars surface/subsurface exploration. This paper presents a model of rough surface and stratified sub-surfaces to describe the multi-layer structure of Mars polar deposits. Based on numerical simulation of radar echoes from rough surface/stratified interfaces, an inversion approach is developed to obtain the parameters of Polar Layered Deposits, i.e. layers thickness and dielectric constants. As a validation example, the SHARAD radar sounder data of the Promethei Lingula of Mars South Polar region is adopted for parameters inversion. The result of stratification is also analyzed and compared with the optical photo of the deep cliff of Chasma Australe canyon. Dielectric inversions show that the deposit media are not uniform, and the dielectric constants of the Promethei Lingula surfaces are large, and become reduced around the depth of 20 m - 30 m, below where most of the deposits are nearly pure ice, except a few thin layers with a lot of dust.展开更多
As is well known, Greece has a significant number of earthquakes each year. Ιn recent years, several earthquakes have occurred in Greece. For this scope, a methodology was used to determine the source parameters. Thi...As is well known, Greece has a significant number of earthquakes each year. Ιn recent years, several earthquakes have occurred in Greece. For this scope, a methodology was used to determine the source parameters. This methodology is based on minimizing the difference between the observed and the synthetic waveforms, using the method Source Parameters Calculation—SPCa <a href="#ref1" target="_blank">[1]</a>. The source parameters, using the proposed methodology, are calculated by comparing observed seismograms and synthetic by inverting data. The synthetics are calculated using the reflectivity method (Kennett, 1983) as implemented by Randall et al. (1994) for a given earth structure. This study includes inversion results for the strongest events that occurred in Greece from 2008 to 2014. For the same events calculated the main fault plane, using the method of Hypocenter Centroid-plot (HC-plot) <a href="#ref2" target="_blank">[2]</a> <a href="#ref3" target="_blank">[3]</a>. This methodology is a simple geometrical method based on the combination between the hypocentral position and the two possible fault planes.展开更多
Four exotic chiral organocatalysts, 9-amino-（9-deoxy） cinchona alkaloids with （8S, 9R） and （8R, 9S）- configurations, were conveniently synthesized for the first time in 27-72% total yields through two conversion...Four exotic chiral organocatalysts, 9-amino-（9-deoxy） cinchona alkaloids with （8S, 9R） and （8R, 9S）- configurations, were conveniently synthesized for the first time in 27-72% total yields through two conversions of configuration at the 9-stereogenic centers of commercially available cinchona alkaloids.展开更多
Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppr...Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppritz equations to estimate Young's modulus,which can introduce cumulative errors and reduce the accuracy of inversion results.To address these issues,this paper introduces the analytical solution of the Zoeppritz equation into the inversion process.The equation is re-derived and expressed in terms of Young's modulus,Poisson's ratio,and density.Within the Bayesian framework,we construct an objective function for the joint inversion of PP and PS waves.Traditional gradient-based algorithms often suffer from low precision and the computational complexity.In this study,we address limitations of conventional approaches related to low precision and complicated code by using Circle chaotic mapping,Levy flights,and Gaussian mutation to optimize the quantum particle swarm optimization(QPSO),named improved quantum particle swarm optimization(IQPSO).The IQPSO demonstrates superior global optimization capabilities.We test the proposed inversion method with both synthetic and field data.The test results demonstrate the proposed method's feasibility and effectiveness,indicating an improvement in inversion accuracy over traditional methods.展开更多
Variations in incoming shortwave radiation influence the net surface heat flux,contributing to the formation of a temperature inversion.The effects of shortwave radiation on the temperature inversions in the Bay of Be...Variations in incoming shortwave radiation influence the net surface heat flux,contributing to the formation of a temperature inversion.The effects of shortwave radiation on the temperature inversions in the Bay of Bengal and eastern equatorial Indian Ocean have never been investigated.Thus,a high-resolution(horizontal resolution of 0.07°×0.07° with 50 vertical layers) Regional Ocean Modeling System(ROMS) model is utilized to quantify the contributions of shortwave radiation to the temperature inversions in the study domain.Analyses of the mixed layer heat and salt budgets are performed,and different model simulations are compared.The model results suggest that a 30% change in shortwave radiation can change approximately 3% of the temperature inversion area in the Bay of Bengal.Low shortwave radiation reduces the net surface heat flux and cools the mixed layer substantially;it also reduces the evaporation rate,causing less evaporative water vapor losses from the ocean than the typical situation,and ultimately enhances haline stratification.Thus,the rudimentary outcome of this research is that a decrease in shortwave radiation produces more temperature inversion in the study region,which is primarily driven by the net surface cooling and supported by the intensive haline stratification.Moreover,low shortwave radiation eventually intensifies the temperature inversion layer by thickening the barrier layer.This study could be an important reference for predicting how the Indian Ocean climate will respond to future changes in shortwave radiation.展开更多
We shall give natural generalized solutions of Hadamard and tensor products equations for matrices by the concept of the Tikhonov regularization combined with the theory of reproducing kernels.
文摘By means of difference calculus, a multifold analog of Gould-Hsu inversion is established. Asconcrete examples, some multifold inverse relations are presented. By developing a standard prooftechnique for combinatorial identities, we demonstrate Mohanty-Hanna's convolution formulas almostmechanically. In addition, an application of reciprocal formulas to interpolation process is alsosketched.
基金the financial supports from National Key R&D Program for Young Scientists of China(Grant No.2022YFC3080900)National Natural Science Foundation of China(Grant No.52374181)+1 种基金BIT Research and Innovation Promoting Project(Grant No.2024YCXZ017)supported by Science and Technology Innovation Program of Beijing institute of technology under Grant No.2022CX01025。
文摘In this study,an inverse design framework was established to find lightweight honeycomb structures(HCSs)with high impact resistance.The hybrid HCS,composed of re-entrant(RE)and elliptical annular re-entrant(EARE)honeycomb cells,was created by constructing arrangement matrices to achieve structural lightweight.The machine learning(ML)framework consisted of a neural network(NN)forward regression model for predicting impact resistance and a multi-objective optimization algorithm for generating high-performance designs.The surrogate of the local design space was initially realized by establishing the NN in the small sample dataset,and the active learning strategy was used to continuously extended the local optimal design until the model converged in the global space.The results indicated that the active learning strategy significantly improved the inference capability of the NN model in unknown design domains.By guiding the iteration direction of the optimization algorithm,lightweight designs with high impact resistance were identified.The energy absorption capacity of the optimal design reached 94.98%of the EARE honeycomb,while the initial peak stress and mass decreased by 28.85%and 19.91%,respectively.Furthermore,Shapley Additive Explanations(SHAP)for global explanation of the NN indicated a strong correlation between the arrangement mode of HCS and its impact resistance.By reducing the stiffness of the cells at the top boundary of the structure,the initial impact damage sustained by the structure can be significantly improved.Overall,this study proposed a general lightweight design method for array structures under impact loads,which is beneficial for the widespread application of honeycomb-based protective structures.
基金financially supported by the National Natural Science Foundation of China (Nos.22171165 and 22371170)Natural Science Foundation of Shandong Province (No.ZR2022MB080)Scientific and Technological Frontiers in Project of Henan Province(No.242102110192)。
文摘Co-assembling chiral molecules with achiral compounds via non-covalent interactions like areneperfluoroarene(AP) interactions offers an effective approach for fabricating chiral functional materials.Herein,chiral molecules L/D-PF1 and L/D-PF2 with pyrene groups were synthesized and its chiroptical properties upon co-assembly with achiral compound octafluoronaphthalene(OFN) through AP interaction were systemically studied.The co-assembly of L/D-PF1/OFN and L/D-PF2/OFN exhibited distinct chiroptical properties such as circular dichroism(CD) and circularly polarized luminescence(CPL) signals.Chirality transfer from the chirality center of L/D-PF1 and L/D-PF2 to the achiral OFN and chiral amplification were successfully achieved.Besides,no significant CPL signal was observed in the self-assembly of L/DPF1 or L/D-PF2 while co-assembly with OFN exhibited obvious CPL amplification induced by AP interaction.Notably,a reversal CD signal and CPL signal could be observed in L/D-PF2/OFN when the molar ratio changed from 1:1 to 1:2 while not found in L/D-PF1/OFN,indicating that that minor structural changes of molecules could cause large changes in assembly.In addition,a series of computational calculations were conducted to verify the AP interaction between L-PF1/L-PF2 and OFN.This work demonstrated that arene-perfluoroarene interaction could drive chiral transfer,chiral amplification and chiral inversion and provided a new method for the preparation of chiroptical materials.
基金supported in part by the National Science Foundation under GrantsDMS 2436630 and 2436629.
文摘Patient-specific finite element analysis(FEA)is a promising tool for noninvasive quantification of cardiac and vascular structural mechanics in vivo.However,inverse material property identification using FEA,which requires iteratively solving nonlinear hyperelasticity problems,is computationally expensive which limits the ability to provide timely patient-specific insights to clinicians.In this study,we present an inverse material parameter identification strategy that integrates deep neural networks(DNNs)with FEA,namely inverse DNN-FEA.In this framework,a DNN encodes the spatial distribution of material parameters and effectively regularizes the inverse solution,which aims to reduce susceptibility to local optima that often arise in heterogeneous nonlinear hyperelastic problems.Consequently,inverse DNN-FEA enables identification of material parameters at the element level.For validation,we applied DNN-FEA to identify four spatially varying passive Holzapfel-Ogden material parameters of the left ventricular myocardium in synthetic benchmark cases with a clinically-derived geometry.To evaluate the benefit of DNN integration,a baseline FEA-only solver implemented in PyTorch was used for comparison.Results demonstrated that DNN-FEA achieved substantially lower average errors in parameter identification compared to FEA(case 1,DNN-FEA:0.37%~2.15%vs.FEA:2.64%~12.91%).The results also demonstrate that the same DNN architecture is capable of identifying a different spatial material property distribution(case 2,DNN-FEA:0.03%~0.60%vs.FEA:0.93%~16.25%).These findings suggest that DNN-FEA provides an accurate framework for inverse identification of heterogeneous myocardial material properties.This approach may facilitate future applications in patient-specific modeling based on in vivo clinical imaging and could be extended to other biomechanical simulation problems.
基金sponsored by the U.S. Department of Energy (DOE)supported by the Ministry of Science and Technology of China (Grant Nos. 2010CB950804 and 2013CB955801)+1 种基金the "Strategic Priority Research Program" of the Chinese Academy of Sciences (Grant No. XDA05100300)the National Natural Science Foundation of China (Grant No. 41305011)
文摘Knowledge of the statistical characteristics of inversions and their effects on aerosols under different large-scale synoptic circulations is important for studying and modeling the diffusion of pollutants in the boundary layer. Based on results gen- erated using the self-organizing map （SOM） weather classification method, this study compares the statistical characteristics of surface-based inversions （SBIs） and elevated inversions （EIs）, and quantitatively evaluates the effect of SBIs on aerosol condensation nuclei （CN） concentrations and the relationship between temperature gradients and aerosols for six prevailing synoptic patterns over the the Southern Great Plains （SGP） site during 2001-10. Large-scale synoptic patterns strongly influ- ence the statistical characteristics of inversions and the accumulation of aerosols in the low-level atmosphere. The activity, frequency, intensity, and vertical distribution of inversions are significantly different among these synoptic patterns. The verti- cal distribution of inversions varies diurnally and is significantly different among the different synoptic patterns. Anticyclonic patterns affect the accumulation of aerosols near the ground more strongly than cyclonic patterns. Mean aerosol CN con- centrations increase during SBIs compared to no inversion cases by 16.1%, 22.6%, 24.5%, 58.7%, 29.8% and 23.7% for the six synoptic patterns. This study confirms that there is a positive correlation between temperature gradients and aerosol CN concentrations near the ground at night under similar large-scale synoptic patterns. The relationship is different for different synoptic patterns and can be described by linear functions. These findings suggest that large-scale synoptic patterns change the static stability of the atmosphere and inversions in the lower atmosphere, thereby influencing the diffusion of aerosols near the ground.
基金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.
基金supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDA17010101)the National Natural Science Foundation of China (Grant Nos. 41305011, 41775033, 41575033 and 41675034)+1 种基金the China Postdoctoral Science Foundation (Grant No.2014M550797)the National Key R&D Program of China (Grant No. 2017YFA0603504)
文摘High-quality and continuous radiosonde, aerosol and surface meteorology datasets are used to investigate the statistical characteristics of meteorological parameters and their effects on aerosols. The data were collected at the Atmospheric Radiation Measurement Southern Great Plains climate research facility during 2000–15. The parameters and vertical distribution of temperature inversion layers were found to have strong diurnal and seasonal changes. For surface-based temperature inversion (SBI), the mean frequency and depth of temperature inversion layers were 39.4% and 198 m, respectively. The temperature difference between the top and bottom of SBI was 4.8℃, and so the temperature gradient was 2.4℃(100 m)^-1. The detailed vertical distributions of temperature inversion had been determined, and only the temperature inversion layers below 1000 m showed diurnal and seasonal variations. Mean surface aerosol number concentrations increased by 43.0%, 21.9% and 49.2% when SBIs were present at 0530, 1730 and 2330 LST, respectively. The effect of SBI on surface aerosol concentration was weakest in summer (18.1%) and strongest in winter (58.4%). During elevated temperature inversion events, there was no noticeable difference in surface aerosol number concentrations. Temperature differences and temperature gradients across SBIs correlated fairly well with aerosol number concentrations, especially for temperature gradients. The vertical distribution of aerosol optical properties with and without temperature inversions was different. Surface aerosol measurements were representative of the air within (below), but not above, SBIs and EIs. These results provide a basis for developing a boundary layer aerosol accumulation model and for improving radiative transfer models in the lower atmosphere.
基金supported by the National Key Research and Development Program of China(2022YFF0800603).
文摘We conducted rapid inversions of rupture process for the 2023 earthquake doublet occurred in SE Türkiye,the first with a magnitude of M_(W)7.8 and the second with a magnitude of M_(W)7.6,using teleseismic and strong-motion data.The teleseismic rupture models of the both events were obtained approximately 88 and 55 minutes after their occurrences,respectively.The rupture models indicated that the first event was an asymmetric bilateral event with ruptures mainly propagating to the northeast,while the second one was a unilateral event with ruptures propagating to the west.This information could be useful in locating the meizoseismal areas.Compared with teleseismic models,the strong-motion models showed relatively higher resolution.A noticeable difference was found for the M_(W)7.6 earthquake,for which the strong-motion models shows a bilateral event,rather than a unilateral event,but the dominant rupture direction is still westward.Nevertheless,all strong-motion models are consistent with the teleseismic models in terms of magnitudes,durations,and dominant rupture directions.This suggests that both teleseismic and strong-motion data can be used for fast determination of major source characteristics.In contrast,the strong-motion data would be preferable in future emergency responses since they are recorded earlier and have a better resolution ability on the source ruptures.
基金Supported by the Science Foundation of Dalian Nationalities University (Grant No. DC110305)
文摘In this paper, by means of Gould-Hsu inverse series relations, we establish several Gould-Hsu inversion chains. As consequence, some new transformation formulae as well as some famous hypergeometric series identities are derived.
文摘The duplicate form of the generalized Gould-Hsu inversions has been obtained by Shi and Zhang. In this paper, we present a simple proof of this duplicate form. With the same method, we construct the duplicate form of the generalized Carlitz inversions. Using this duplicate form, we obtain several terminating basic hypergeometric identities and some limiting cases.
基金partially supported by the Doctoral Programming Research Fund of Higher Education, Chinese Ministry of Education (No. 20110205110010)
文摘It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regional earthquake tomography, but seldom are all three classes of model parameters updated simultaneously. This is mainly due to the trade-off between the different types of model parameters and the lack of different seismic phases to constrain the model parameters. Using a spherical-coordinate ray tracing algorithm for first and later（primary reflected） arrival tracing algorithm in combination with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model parameters in regional or global tomographic studies. In this paper we incorporate the multistage irregular shortest-path ray tracing algorithm（in a spherical coordinate system） with a subspace inversion solver to formulate a simultaneous inversion algorithm for triple model parameters updating using direct and later arrival time information.Comparison tests for two sets of data（noise free and added noise） indicate that the new triple-class parameter inversion algorithm is capable of obtaining nearly the same results as the double-class parameter inversion scheme. Furthermore,the proposed multi-parameter type inversion method is not sensitive to a modest level of picking error in the traveltime data, and also performs well with a relatively large uncertainty in earthquake hypocentral locations. This shows it to be a feasible and promising approach in regional or global tomographic applications.
基金National Natural Science Foundation of China (40074010).
文摘In the paper, a 2D symmetrical anisotropic medium whose strike agrees with one of the horizontal principal axes is considered to develop a corresponding inversion technique. In the specified conditions, if we assume an equivalent conductivity anisotropy in both the vertical and dipping directions, i.e., σzz=σyy, the differential equations obtained are formally the same as that for TE and TM modes in the 2D isotropic geoelectrical media. The same inversion technique as that in the 2D isotropic media can be employed to obtain the anisotropic conductivities. It means that the TE and TM inversion results in the isotropic media can be respectively thought as the resistivities in the two principal directions of the symmetrically anisotropic media, which has offered a new approach and a theoretical guidance for interpreting magnetotelluric data. And the inversion technique developed here is used to test the magnetotelluric data in the area of Tianzhu and Yongdeng in Gansu Province, so that the crust anisotropic geoelectrical structures in this region can be obtained.
基金supported by the National Natural Science Foundation of China (grant nos. 41930103 and 41674052)
文摘The Bayesian inversion method is a stochastic approach based on the Bayesian theory.With the development of sampling algorithms and computer technologies,the Bayesian inversion method has been widely used in geophysical inversion problems.In this study,we conduct inversion experiments using crosshole seismic travel-time data to examine the characteristics and performance of the stochastic Bayesian inversion based on the Markov chain Monte Carlo sampling scheme and the traditional deterministic inversion with Tikhonov regularization.Velocity structures with two different spatial variations are considered,one with a chessboard pattern and the other with an interface mimicking the Mohorovicicdiscontinuity(Moho).Inversions are carried out with different scenarios of model discretization and source–receiver configurations.Results show that the Bayesian method yields more robust single-model estimations than the deterministic method,with smaller model errors.In addition,the Bayesian method provides the posterior probabilistic distribution function of the model space,which can help us evaluate the quality of the inversion result.
文摘HF (high frequency) radar sounder technology has been developed for several missions of Mars surface/subsurface exploration. This paper presents a model of rough surface and stratified sub-surfaces to describe the multi-layer structure of Mars polar deposits. Based on numerical simulation of radar echoes from rough surface/stratified interfaces, an inversion approach is developed to obtain the parameters of Polar Layered Deposits, i.e. layers thickness and dielectric constants. As a validation example, the SHARAD radar sounder data of the Promethei Lingula of Mars South Polar region is adopted for parameters inversion. The result of stratification is also analyzed and compared with the optical photo of the deep cliff of Chasma Australe canyon. Dielectric inversions show that the deposit media are not uniform, and the dielectric constants of the Promethei Lingula surfaces are large, and become reduced around the depth of 20 m - 30 m, below where most of the deposits are nearly pure ice, except a few thin layers with a lot of dust.
文摘As is well known, Greece has a significant number of earthquakes each year. Ιn recent years, several earthquakes have occurred in Greece. For this scope, a methodology was used to determine the source parameters. This methodology is based on minimizing the difference between the observed and the synthetic waveforms, using the method Source Parameters Calculation—SPCa <a href="#ref1" target="_blank">[1]</a>. The source parameters, using the proposed methodology, are calculated by comparing observed seismograms and synthetic by inverting data. The synthetics are calculated using the reflectivity method (Kennett, 1983) as implemented by Randall et al. (1994) for a given earth structure. This study includes inversion results for the strongest events that occurred in Greece from 2008 to 2014. For the same events calculated the main fault plane, using the method of Hypocenter Centroid-plot (HC-plot) <a href="#ref2" target="_blank">[2]</a> <a href="#ref3" target="_blank">[3]</a>. This methodology is a simple geometrical method based on the combination between the hypocentral position and the two possible fault planes.
基金Research support from the National Science Foundation ofChina(No.21071116)Chongqing Scientific Foundation,China(No.2010BB4126)
文摘Four exotic chiral organocatalysts, 9-amino-（9-deoxy） cinchona alkaloids with （8S, 9R） and （8R, 9S）- configurations, were conveniently synthesized for the first time in 27-72% total yields through two conversions of configuration at the 9-stereogenic centers of commercially available cinchona alkaloids.
基金supported by Fundamental Research Funds for the Central Universities,CHD300102264715National Key Research and Development Program of China under Grant 2021YFA0716902Natural Science Basic Research Program of Shaanxi 2024JCYBMS-199。
文摘Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppritz equations to estimate Young's modulus,which can introduce cumulative errors and reduce the accuracy of inversion results.To address these issues,this paper introduces the analytical solution of the Zoeppritz equation into the inversion process.The equation is re-derived and expressed in terms of Young's modulus,Poisson's ratio,and density.Within the Bayesian framework,we construct an objective function for the joint inversion of PP and PS waves.Traditional gradient-based algorithms often suffer from low precision and the computational complexity.In this study,we address limitations of conventional approaches related to low precision and complicated code by using Circle chaotic mapping,Levy flights,and Gaussian mutation to optimize the quantum particle swarm optimization(QPSO),named improved quantum particle swarm optimization(IQPSO).The IQPSO demonstrates superior global optimization capabilities.We test the proposed inversion method with both synthetic and field data.The test results demonstrate the proposed method's feasibility and effectiveness,indicating an improvement in inversion accuracy over traditional methods.
基金The Marine Scholarship of ChinaChina Scholarship Council for International Doctoral Students under contract No.2017SOA016552the National Natural Science Foundation of China under contract Nos U2106204 and 41676003。
文摘Variations in incoming shortwave radiation influence the net surface heat flux,contributing to the formation of a temperature inversion.The effects of shortwave radiation on the temperature inversions in the Bay of Bengal and eastern equatorial Indian Ocean have never been investigated.Thus,a high-resolution(horizontal resolution of 0.07°×0.07° with 50 vertical layers) Regional Ocean Modeling System(ROMS) model is utilized to quantify the contributions of shortwave radiation to the temperature inversions in the study domain.Analyses of the mixed layer heat and salt budgets are performed,and different model simulations are compared.The model results suggest that a 30% change in shortwave radiation can change approximately 3% of the temperature inversion area in the Bay of Bengal.Low shortwave radiation reduces the net surface heat flux and cools the mixed layer substantially;it also reduces the evaporation rate,causing less evaporative water vapor losses from the ocean than the typical situation,and ultimately enhances haline stratification.Thus,the rudimentary outcome of this research is that a decrease in shortwave radiation produces more temperature inversion in the study region,which is primarily driven by the net surface cooling and supported by the intensive haline stratification.Moreover,low shortwave radiation eventually intensifies the temperature inversion layer by thickening the barrier layer.This study could be an important reference for predicting how the Indian Ocean climate will respond to future changes in shortwave radiation.
文摘We shall give natural generalized solutions of Hadamard and tensor products equations for matrices by the concept of the Tikhonov regularization combined with the theory of reproducing kernels.
