Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an...Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an approach for the long-term extreme-response analysis of floating structures. A modified gradient-based retrieval algorithm in conjunction with the inverse first-order reliability method(IFORM) is proposed to enable the use of convolution models in long-term extreme analysis of structures with an analytical formula of response amplitude operator(RAO). The proposed algorithm ensures convergence stability and iteration accuracy and exhibits a higher computational efficiency than the traditional backtracking method. However, when the RAO of general offshore structures cannot be analytically expressed, the convolutional integration method fails to function properly. A numerical discretization approach is further proposed for offshore structures in the case when the analytical expression of the RAO is not feasible. Through iterative discretization of environmental contours(ECs) and RAOs, a detailed procedure is proposed to calculate the long-term response extremes of offshore structures. The validity and accuracy of the proposed approach are tested using a floating offshore wind turbine as a numerical example. The long-term extreme heave responses of various return periods are calculated via the IFORM in conjunction with a numerical discretization approach. The environmental data corresponding to N-year structural responses are located inside the ECs, which indicates that the selection of design points directly along the ECs yields conservative design results.展开更多
Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this cha...Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.展开更多
Serial-parallel manipulators are of great interest to academic community in recent years,especially those composed of classical parallel mechanisms.There have been many studies around 2(3RPS)and 2(3SPR)S-PMs,but unfor...Serial-parallel manipulators are of great interest to academic community in recent years,especially those composed of classical parallel mechanisms.There have been many studies around 2(3RPS)and 2(3SPR)S-PMs,but unfortunately their inverse kinematics have not yet been resolved.This paper discovers that the unknown kinematic parameters of middle platform are responsible for the unresolvable of inverse kinematics,meanwhile the unknown kinematic parameters of middle platform also have huge coupling relationships.Therefore,to break through this challenges,the huge coupling relationships are decoupled layer by layer,the kinematic parameters of middle platform are solved by combining Sylvester's elimination method,and the inverse displacements of 2(3RPS)and 2(3SPR)S-PMs are obtained subsequently.This paper not only solves the inverse kinematics of classical 2(3RPS)and 2(3SPR)S-PMs,but also reveals the essence of the inverse kinematics of general(3-DOF)+(3-DOF)6-DOF S-PMs and proposes a corresponding solution.展开更多
How to utilize existing flow control mechanisms to make profiled end wall design more flexible,efficient,and physical is a meaningful challenge.This study presents a three-dimensional inverse method for profiled end w...How to utilize existing flow control mechanisms to make profiled end wall design more flexible,efficient,and physical is a meaningful challenge.This study presents a three-dimensional inverse method for profiled end wall design to achieve the application of flow control mechanisms.The predetermined pressure distribution on the end wall is reached by modifying the end wall geometry during flow field calculation.A motion velocity model is derived from the normal momentum equation of the moving no-slip boundary to modify the end wall geometry.A Reynolds-Averaged Navier-Stokes(RANS)solver based on the Semi-Implicit Method for Pressure Linked Equations(SIMPLE)algorithm is adopted to simulate the flow field.Based on the mechanism understanding obtained through numerical optimization results,this study adopts the inverse method to redesign an optimized end wall in a compressor cascade.The results indicate that the redesigned end wall exhibits better loss reduction,reducing the overall total pressure loss by 5.5%,whereas the optimized end wall reduces it by 3%.The inverse method allows the imposition of desired influences on the end wall flow without constructing a database,making it highly flexible,efficient,and physical.展开更多
The ^(12)C+^(12)C fusion is one of the most important reactions in modern nuclear astrophysics.The trend and magnitude of the reaction rate within the Gamow window strongly influence various astrophysical processes.Ho...The ^(12)C+^(12)C fusion is one of the most important reactions in modern nuclear astrophysics.The trend and magnitude of the reaction rate within the Gamow window strongly influence various astrophysical processes.However,direct measurement of this reaction is extremely difficult,which makes it necessary to develop indirect methods.In this study,the ^(23)Na+p reaction system was used to study the compound nucleus ^(24)Mg.We employed a thick-target inverse kinematics method combined with theγ-charged-particle coincidence technique to measure the proton andα exit channels of ^(24)Mg.Technical details of the ^(23)Na+p thick-target inverse kinematics experiment and analysis are presented herein.展开更多
This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeabi...This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.展开更多
Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods rema...Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods remains challenging,as it requires the estimation of more than eight parameters.Assuming the reservoir to be a weakly anisotropic ORT medium with small contrasts in the background elastic parameters,a new azimuthal elastic impedance equation was first derived using parameter combinations and mathematical approximations.This equation exhibited almost the same accuracy as the original equation and contained only six model parameters:the compression modulus,anisotropic shear modulus,anisotropic compression modulus,density,normal fracture weakness,and tangential fracture weakness.Subsequently,a stepwise inversion method using second-order derivatives of the elastic impedance was developed to estimate these parameters.Moreover,the Thomsen anisotropy parameter,epsilon,was estimated from the inversion results using the ratio of the anisotropic compression modulus to the compression modulus.Synthetic examples with moderate noise and field data examples confirm the feasibility and effectiveness of the inversion method.The proposed method exhibited accuracy similar to that of previous inversion strategies and could predict richer vertical fracture information.Ultimately,the method was applied to a three-dimensional work area,and the predictions were consistent with logging and geological a priori information,confirming the effectiveness of this method.Summarily,the proposed stepwise inversion method can alleviate the uncertainty of multi-parameter inversion in ORT medium,thereby improving the reliability of fracture detection.展开更多
An inverse analysis is presented to estimate line heat source in two-dimensional steady-state and transient heat transfer problems.A constant heat source is considered in the steady-state heat transfer problem(a param...An inverse analysis is presented to estimate line heat source in two-dimensional steady-state and transient heat transfer problems.A constant heat source is considered in the steady-state heat transfer problem(a parameter estimation problem)and a time-varying heat source is considered in the transient heat transfer problem(a function estimation problem).Since a general irregular 2D heat conducting body is considered,a body-fitted grid generation is used to mesh the domain.Then governing equations and associated boundary and initial conditions are transformed from the physical domain to the computational domain and finite difference method is used to solve the governing equations to obtain the temperature distribution in the body.Using an efficient,accurate,and very easy to implement sensitivity analysis incorporated in a gradient based minimization method(here,steepest descentmethod),the unknown heat source is estimated accurately.In the function estimation part,it is assumed that there is no prior information on the functional form of the heat source and the estimation process can be performed with a reasonable initial guess for the heat source.The main advantage of the proposed inverse analysis is that the sensitivity matrix(and hence,the objective function gradient with respect to the unknown variables)can be computed during the direct heat transfer solution through newyet simple explicit expressions with no need to solve extra equations such as the sensitivity and adjoint problems and impose additional computational costs comparable to the direct problem solution ones.Some test cases are presented to investigate the accuracy,efficiency,and effect of measurement error on the estimated parameter and function for the line heat source.展开更多
The reliable estimation of the wavenumber space(k-space)of the plates remains a longterm concern for acoustic modeling and structural dynamic behavior characterization.Most current analyses of wavenumber identificatio...The reliable estimation of the wavenumber space(k-space)of the plates remains a longterm concern for acoustic modeling and structural dynamic behavior characterization.Most current analyses of wavenumber identification methods are based on the deterministic hypothesis.To this end,an inverse method is proposed for identifying wave propagation characteristics of twodimensional structures under stochastic conditions,such as wavenumber space,dispersion curves,and band gaps.The proposed method is developed based on an algebraic identification scheme in the polar coordinate system framework,thus named Algebraic K-Space Identification(AKSI)technique.Additionally,a model order estimation strategy and a wavenumber filter are proposed to ensure that AKSI is successfully applied.The main benefit of AKSI is that it is a reliable and fast method under four stochastic conditions:(A)High level of signal noise;(B)Small perturbation caused by uncertainties in measurement points’coordinates;(C)Non-periodic sampling;(D)Unknown structural periodicity.To validate the proposed method,we numerically benchmark AKSI and three other inverse methods to extract dispersion curves on three plates under stochastic conditions.One experiment is then performed on an isotropic steel plate.These investigations demonstrate that AKSI is a good in-situ k-space estimator under stochastic conditions.展开更多
This study presents a novel two-step approach to assess plate-like structural laminar damages,particularly for delamination damage detection of composite structures.Firstly,a 2-D continuous wavelet transform is employ...This study presents a novel two-step approach to assess plate-like structural laminar damages,particularly for delamination damage detection of composite structures.Firstly,a 2-D continuous wavelet transform is employed to identify the damage location and sizes from vibration curvature data.An inverse method is subsequently then used to determine the bending stiffness reduction ratio along a specified direction,enabling the quantification of the delamination severity.The method employed in this study is an extension of the one-dimensional inverse method developed in a previous work of the authors.The applicability of the two-step inverse approach is demonstrated in a simulation analysis and by an experimental study on a cantilever composite plate containing a single delamination.The inverse method is shown to have the capacity to reveal the detailed damage information of delamination within a constrained searching space and can be used to determine the effective flexural stiffness of composite plate structures,even in cases of complex delamination damage.展开更多
[Objective] The research aimed to isolate flanking sequences adjacent to the transgenic T-DNA in Brassica napus by an improved inverse PCR method.[Method] Using single clone of transgenic FS4 in Brassica napus as the ...[Objective] The research aimed to isolate flanking sequences adjacent to the transgenic T-DNA in Brassica napus by an improved inverse PCR method.[Method] Using single clone of transgenic FS4 in Brassica napus as the research materials,total DNA was extracted from transgenic Brassica napus by using modified CTAB method.After enzyme digestion and purification,self-joining was made.Two circles of nested PCR and the sequence alignment were carried out.[Result] A fragement with the size of 4.0 kb was amplified ...展开更多
The deflection of the vertical(DOV)is the key information in the study of ocean gravity field.However,in most areas,the precision of the prime component of DoV is significantly lower than that of the meridian componen...The deflection of the vertical(DOV)is the key information in the study of ocean gravity field.However,in most areas,the precision of the prime component of DoV is significantly lower than that of the meridian component.To obtain higher accuracy and resolution of ocean gravity information,researchers have proposed a novel altimeter called the wide-swath altimeter.This altimeter allows for the simultaneous acquisition of high-precision and high-resolution two-dimensional measurements of sea surface height(SSH).In this paper,the Surface Water and Ocean Topography(SWOT)mission with a wide-swath altimeter on board is selected for research.One cycle of swoT sea surface height data is simulated to inverse the DOV in the Arabian Sea(45°E—80°E,0°-30°N),and the inversion results are compared with those of conventional altimeter data.The results demonstrate that the difference between the meridian and prime components derived from the inversion of swoT wide-swath data is minimal,significantly outperforming the inversion results of conventional nadir altimeter data.The advantage of swoT wide-swath altimeter lies in its ability to use the multi-directional geoid slope at any sea surface measurement point to invert the components in the meridian and prime directions.To investigate the impact of this advantage on inversion precision,this paper employs a method to calculate the gradient of the geoid in multiple directions to invert DoV components.The improvement effect of calculating the gradient of the geoid in multiple directions on the precision of DoV component is analyzed.It is found that the accuracy of DoV inversion has significantly improved with the increase of geodetic gradient calculation direction.In addition,the effects of various errors and grid spacing in SwoT wide sea surface height data on the precision of Dov inversion are also analyzed.展开更多
To improve the accuracy and reduce the calcu- lation cost for the inverse problem of centrifugal pump impeller, the new inverse method based on proper orthog- onal decomposition (POD) is proposed. The pump blade sha...To improve the accuracy and reduce the calcu- lation cost for the inverse problem of centrifugal pump impeller, the new inverse method based on proper orthog- onal decomposition (POD) is proposed. The pump blade shape is parameterized by quartic Bezier curve, and the initial snapshots is generated by introducing the perturbation of the blade shape control parameters. The internal flow field and its hydraulic performance is predicted by CFD method. The snapshots vector includes the blade shape parameter and the distribution of blade load. The POD basis for the snap- shots set are deduced by proper orthogonal decomposition. The sample vector set is expressed in terms of the linear combination of the orthogonal basis. The objective blade shape corresponding to the objective distribution of blade load is obtained by least square fit. The Iterative correction algorithm for the centrifugal pump blade inverse method based on POD is proposed. The objective blade load dis- tributions are corrected according to the difference of the CFD result and the POD result. The two dimensional and three dimensional blade calculation cases show that the proposed centrifugal pump blade inverse method based on POD have good convergence and high accuracy, and thecalculation cost is greatly reduced. After two iterations, the deviation of the blade load and the pump hydraulic perfor- mance are limited within 4.0% and 6.0% individually for most of the flow rate range. This paper provides a promising inverse method for centrifugal pump impeller, which will benefit the hydraulic optimization of centrifugal pump.展开更多
Inverse method was used in single crystal superalloy DD6 processing simulation during solidification. Numerical modeling coupled with experiments has been used to estimate the interface heat transfer coefficient (IHT...Inverse method was used in single crystal superalloy DD6 processing simulation during solidification. Numerical modeling coupled with experiments has been used to estimate the interface heat transfer coefficient (IHTC) between the surface of slab casting and inner mold. Calculated temperature dependent values of IHTC were obtained from a numerical solution. The calculated temperatures agreed well with the measurement of cooling profile.展开更多
Autonomous underwater gliders are highly effcient,buoyancy-driven,winged autonomous underwater vehicles. Their dynamics are multivariable nonlinear systems. In addition,the gliders are underactuated and diffcult to ma...Autonomous underwater gliders are highly effcient,buoyancy-driven,winged autonomous underwater vehicles. Their dynamics are multivariable nonlinear systems. In addition,the gliders are underactuated and diffcult to maneuver,and also dependent on their operational environment. To confront these problems and to design an effective controller,the inverse system method was used to decouple the original system into two independent single variable linear subsystems. The stability of the zero dynamics was analyzed,and an additional closed-loop controller for each linear subsystem was designed by sliding mode control method to form a type of composite controller. Simulation results demonstrate that the derived nonlinear controller is able to cope with the aforementioned problems simultaneously and satisfactorily.展开更多
The heat transfer coefficient in a multidimensional heat conduction problem is obtained from the solution of the inverse heat conduction problem based on the thermographic temperature measurement. The modified one-dim...The heat transfer coefficient in a multidimensional heat conduction problem is obtained from the solution of the inverse heat conduction problem based on the thermographic temperature measurement. The modified one-dimensional correction method (MODCM), along with the finite volume method, is employed for both two- and three-dimensional inverse problems. A series of numerical experiments are conducted in order to verify the effectiveness of the method. In addition, the effect of the temperature measurement error, the ending criterion of the iteration, etc. on the result of the inverse problem is investigated. It is proved that the method is a simple, stable and accurate one that can solve successfully the inverse heat conduction problem.展开更多
On the basis of hydrographic data obtained in August 2000 cruise, the circulation in the South China Sea (SCS) is computed by the modified inverse method in combination with SSH data from TOPEX/ERS-2 analysis. For stu...On the basis of hydrographic data obtained in August 2000 cruise, the circulation in the South China Sea (SCS) is computed by the modified inverse method in combination with SSH data from TOPEX/ERS-2 analysis. For study of the dynamical mechanism, which causes the pattern of summer circulation in the SCS, the diagnostic model (Yuan et al. 1982. Acta Oceanologica Sinica,4(1):1-11; Yuan and Su. 1992. Numerical Computation of Physical Oceanography.474-542) is used to simulate numerically the summer circulation in the SCS. The following results have been obtained. (1) The central and southwestern SCSs are dominated mainly by anticy-clonic circulation systems. They are mainly as follows. 1) There is strong anticyclonic eddy southeast of Vietnam (W1). Its horizontal scale is about 300 km, and it extends vertically from the surface to the about 1 000 m level. 2) There are a warm eddy W2 southeast of Zhongsha Islands and the anticyclonic circulation system W3 west off the Luzon Island. 3) There is a stronger cyclonic eddy C1 between the anticyclonic eddies W1 and W2.4) A strong northward coastal jet is present near the coast of Vietnam, and separates from the coast of Vietnam at about 12° N to the northeast.(2)The northern SCS is dominated mainly by a cyclonic circulation system. There is a cyclonic circulation system near and north of Section N2. (3) The southeastern SCS is dominated mainly by the cyclonic circulation system. (4) Comparing the results of circulation in the SCS during the summer of 2000 with those during the summer of 1998, it is found that they agree qualitatively, but there is the some difference between them in quantity.This shows that the circulation in the SCS has obviously seasonal feature. (5) The dynamical mechanism which products the basic pattern of summer circulation is because the following two reasons: 1) the joint effect of the baroclinity and relief (JEBAR) is essential dynamical cause; and 2) it is next important dynamical cause that the interaction between the wind stress and bottom topography under the southerly monsoon. (6) Comparing the hydrographic structure and distribution of stream functions with the SSH data from TOPEX/ERS-2 analysis in the SCS during August of 2000, they agree qualitatively.展开更多
Inverse heat conduction method (IHCM) is one of the most effective approaches to obtaining the boiling heat transfer coefficient from measured results. This paper focuses on its application in cryogenic boiling heat t...Inverse heat conduction method (IHCM) is one of the most effective approaches to obtaining the boiling heat transfer coefficient from measured results. This paper focuses on its application in cryogenic boiling heat transfer. Experiments were conducted on the heat transfer of a stainless steel block in a liquid nitrogen bath, with the assumption of a 1D conduction condition to realize fast acquisition of the temperature of the test points inside the block. With the inverse-heat conduction theory and the explicit finite difference model, a solving program was developed to calculate the heat flux and the boiling heat transfer coefficient of a stainless steel block in liquid nitrogen bath based on the temperature acquisition data. Considering the oscillating data and some unsmooth transition points in the inverse-heat-conduction calculation result of the heat-transfer coefficient, a two-step data-fitting procedure was proposed to obtain the expression for the boiling heat transfer coefficients. The coefficient was then verified for accuracy by a comparison between the simulation results using this expression and the verifying experimental results of a stainless steel block. The maximum error with a revised segment fitting is around 6%, which verifies the feasibility of using IHCM to measure the boiling heat transfer coefficient in liquid nitrogen bath.展开更多
The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise w...The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.展开更多
This paper is concerned with the inverse scattering problems for Schrdinger equations with compactly supported potentials.For purpose of reconstructing the support of the potential,we derive a factorization of the sca...This paper is concerned with the inverse scattering problems for Schrdinger equations with compactly supported potentials.For purpose of reconstructing the support of the potential,we derive a factorization of the scattering amplitude operator A and prove that the ranges of (A* A) ^1/4 and G which maps more general incident fields than plane waves into the scattering amplitude coincide.As an application we characterize the support of the potential using only the spectral data of the operator A.展开更多
基金Supported by the National Natural Science Foundation of China (Grant Nos.52088102 and 51879287)National Key Research and Development Program of China (Grant No.2022YFB2602301)。
文摘Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an approach for the long-term extreme-response analysis of floating structures. A modified gradient-based retrieval algorithm in conjunction with the inverse first-order reliability method(IFORM) is proposed to enable the use of convolution models in long-term extreme analysis of structures with an analytical formula of response amplitude operator(RAO). The proposed algorithm ensures convergence stability and iteration accuracy and exhibits a higher computational efficiency than the traditional backtracking method. However, when the RAO of general offshore structures cannot be analytically expressed, the convolutional integration method fails to function properly. A numerical discretization approach is further proposed for offshore structures in the case when the analytical expression of the RAO is not feasible. Through iterative discretization of environmental contours(ECs) and RAOs, a detailed procedure is proposed to calculate the long-term response extremes of offshore structures. The validity and accuracy of the proposed approach are tested using a floating offshore wind turbine as a numerical example. The long-term extreme heave responses of various return periods are calculated via the IFORM in conjunction with a numerical discretization approach. The environmental data corresponding to N-year structural responses are located inside the ECs, which indicates that the selection of design points directly along the ECs yields conservative design results.
基金funded by the National Key R&D Program of China(Grant No.2022YFC2903904)the National Natural Science Foundation of China(Grant Nos.51904057 and U1906208).
文摘Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.
基金Supported by National Natural Science Foundation of China(Grant No.52275033)National Natural Science Youth Foundation of China(Grant No.52205033)Hebei Provincial Natural Science Foundation of China(Grant No.E2021203019)。
文摘Serial-parallel manipulators are of great interest to academic community in recent years,especially those composed of classical parallel mechanisms.There have been many studies around 2(3RPS)and 2(3SPR)S-PMs,but unfortunately their inverse kinematics have not yet been resolved.This paper discovers that the unknown kinematic parameters of middle platform are responsible for the unresolvable of inverse kinematics,meanwhile the unknown kinematic parameters of middle platform also have huge coupling relationships.Therefore,to break through this challenges,the huge coupling relationships are decoupled layer by layer,the kinematic parameters of middle platform are solved by combining Sylvester's elimination method,and the inverse displacements of 2(3RPS)and 2(3SPR)S-PMs are obtained subsequently.This paper not only solves the inverse kinematics of classical 2(3RPS)and 2(3SPR)S-PMs,but also reveals the essence of the inverse kinematics of general(3-DOF)+(3-DOF)6-DOF S-PMs and proposes a corresponding solution.
基金Supported by the National Natural Science Foundation of China(No.52376021).
文摘How to utilize existing flow control mechanisms to make profiled end wall design more flexible,efficient,and physical is a meaningful challenge.This study presents a three-dimensional inverse method for profiled end wall design to achieve the application of flow control mechanisms.The predetermined pressure distribution on the end wall is reached by modifying the end wall geometry during flow field calculation.A motion velocity model is derived from the normal momentum equation of the moving no-slip boundary to modify the end wall geometry.A Reynolds-Averaged Navier-Stokes(RANS)solver based on the Semi-Implicit Method for Pressure Linked Equations(SIMPLE)algorithm is adopted to simulate the flow field.Based on the mechanism understanding obtained through numerical optimization results,this study adopts the inverse method to redesign an optimized end wall in a compressor cascade.The results indicate that the redesigned end wall exhibits better loss reduction,reducing the overall total pressure loss by 5.5%,whereas the optimized end wall reduces it by 3%.The inverse method allows the imposition of desired influences on the end wall flow without constructing a database,making it highly flexible,efficient,and physical.
基金supported by the National Key Research and Development Project (No. 2022YFA1602301)the National Natural Science Foundation of China (Nos. U2267205, 12275361, 12125509, 12222514, 11961141003, and 12005304)+2 种基金the CAST Young Talent Support Planthe CNNC Science Fund for Talented Young Scholarsthe Continuous-Support Basic Scientific Research Project
文摘The ^(12)C+^(12)C fusion is one of the most important reactions in modern nuclear astrophysics.The trend and magnitude of the reaction rate within the Gamow window strongly influence various astrophysical processes.However,direct measurement of this reaction is extremely difficult,which makes it necessary to develop indirect methods.In this study,the ^(23)Na+p reaction system was used to study the compound nucleus ^(24)Mg.We employed a thick-target inverse kinematics method combined with theγ-charged-particle coincidence technique to measure the proton andα exit channels of ^(24)Mg.Technical details of the ^(23)Na+p thick-target inverse kinematics experiment and analysis are presented herein.
文摘This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.
基金sponsorship of the National Natural Science Foundation of China(42430809,42274157,42030103,42404132)the Fund of State Key Laboratory of Deep Oil and Gas,China University of Petroleum(East China)(SKLDOG2024-ZYTS-02)+5 种基金the Postdoctoral Fellowship Program of CPSF(GZB20240850)the Postdoctoral Project of Qingdao(QDBSH20240102082)the Fundamental Research Funds for the Central Universities(24CX07004A,24CX06036A)the CNPC Innovation Fund(2024DQ02-0505,2024DQ02-0136)the Innovation fund project for graduate student of China University of Petroleum(East China)the Fundamental Research Funds for the Central Universities(24CX04002A).
文摘Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods remains challenging,as it requires the estimation of more than eight parameters.Assuming the reservoir to be a weakly anisotropic ORT medium with small contrasts in the background elastic parameters,a new azimuthal elastic impedance equation was first derived using parameter combinations and mathematical approximations.This equation exhibited almost the same accuracy as the original equation and contained only six model parameters:the compression modulus,anisotropic shear modulus,anisotropic compression modulus,density,normal fracture weakness,and tangential fracture weakness.Subsequently,a stepwise inversion method using second-order derivatives of the elastic impedance was developed to estimate these parameters.Moreover,the Thomsen anisotropy parameter,epsilon,was estimated from the inversion results using the ratio of the anisotropic compression modulus to the compression modulus.Synthetic examples with moderate noise and field data examples confirm the feasibility and effectiveness of the inversion method.The proposed method exhibited accuracy similar to that of previous inversion strategies and could predict richer vertical fracture information.Ultimately,the method was applied to a three-dimensional work area,and the predictions were consistent with logging and geological a priori information,confirming the effectiveness of this method.Summarily,the proposed stepwise inversion method can alleviate the uncertainty of multi-parameter inversion in ORT medium,thereby improving the reliability of fracture detection.
文摘An inverse analysis is presented to estimate line heat source in two-dimensional steady-state and transient heat transfer problems.A constant heat source is considered in the steady-state heat transfer problem(a parameter estimation problem)and a time-varying heat source is considered in the transient heat transfer problem(a function estimation problem).Since a general irregular 2D heat conducting body is considered,a body-fitted grid generation is used to mesh the domain.Then governing equations and associated boundary and initial conditions are transformed from the physical domain to the computational domain and finite difference method is used to solve the governing equations to obtain the temperature distribution in the body.Using an efficient,accurate,and very easy to implement sensitivity analysis incorporated in a gradient based minimization method(here,steepest descentmethod),the unknown heat source is estimated accurately.In the function estimation part,it is assumed that there is no prior information on the functional form of the heat source and the estimation process can be performed with a reasonable initial guess for the heat source.The main advantage of the proposed inverse analysis is that the sensitivity matrix(and hence,the objective function gradient with respect to the unknown variables)can be computed during the direct heat transfer solution through newyet simple explicit expressions with no need to solve extra equations such as the sensitivity and adjoint problems and impose additional computational costs comparable to the direct problem solution ones.Some test cases are presented to investigate the accuracy,efficiency,and effect of measurement error on the estimated parameter and function for the line heat source.
基金supported by the Lyon Acoustics Center of Lyon University,Francefunded by the China Scholarship Council(CSC)。
文摘The reliable estimation of the wavenumber space(k-space)of the plates remains a longterm concern for acoustic modeling and structural dynamic behavior characterization.Most current analyses of wavenumber identification methods are based on the deterministic hypothesis.To this end,an inverse method is proposed for identifying wave propagation characteristics of twodimensional structures under stochastic conditions,such as wavenumber space,dispersion curves,and band gaps.The proposed method is developed based on an algebraic identification scheme in the polar coordinate system framework,thus named Algebraic K-Space Identification(AKSI)technique.Additionally,a model order estimation strategy and a wavenumber filter are proposed to ensure that AKSI is successfully applied.The main benefit of AKSI is that it is a reliable and fast method under four stochastic conditions:(A)High level of signal noise;(B)Small perturbation caused by uncertainties in measurement points’coordinates;(C)Non-periodic sampling;(D)Unknown structural periodicity.To validate the proposed method,we numerically benchmark AKSI and three other inverse methods to extract dispersion curves on three plates under stochastic conditions.One experiment is then performed on an isotropic steel plate.These investigations demonstrate that AKSI is a good in-situ k-space estimator under stochastic conditions.
文摘This study presents a novel two-step approach to assess plate-like structural laminar damages,particularly for delamination damage detection of composite structures.Firstly,a 2-D continuous wavelet transform is employed to identify the damage location and sizes from vibration curvature data.An inverse method is subsequently then used to determine the bending stiffness reduction ratio along a specified direction,enabling the quantification of the delamination severity.The method employed in this study is an extension of the one-dimensional inverse method developed in a previous work of the authors.The applicability of the two-step inverse approach is demonstrated in a simulation analysis and by an experimental study on a cantilever composite plate containing a single delamination.The inverse method is shown to have the capacity to reveal the detailed damage information of delamination within a constrained searching space and can be used to determine the effective flexural stiffness of composite plate structures,even in cases of complex delamination damage.
基金Supported by National 863 Program of China(2006AA10A113)Natural Science foundation of Zhejiang Province(Y306097)~~
文摘[Objective] The research aimed to isolate flanking sequences adjacent to the transgenic T-DNA in Brassica napus by an improved inverse PCR method.[Method] Using single clone of transgenic FS4 in Brassica napus as the research materials,total DNA was extracted from transgenic Brassica napus by using modified CTAB method.After enzyme digestion and purification,self-joining was made.Two circles of nested PCR and the sequence alignment were carried out.[Result] A fragement with the size of 4.0 kb was amplified ...
基金support from the National Natural Science Foundation of China(No.42274006,42192535,42242015).
文摘The deflection of the vertical(DOV)is the key information in the study of ocean gravity field.However,in most areas,the precision of the prime component of DoV is significantly lower than that of the meridian component.To obtain higher accuracy and resolution of ocean gravity information,researchers have proposed a novel altimeter called the wide-swath altimeter.This altimeter allows for the simultaneous acquisition of high-precision and high-resolution two-dimensional measurements of sea surface height(SSH).In this paper,the Surface Water and Ocean Topography(SWOT)mission with a wide-swath altimeter on board is selected for research.One cycle of swoT sea surface height data is simulated to inverse the DOV in the Arabian Sea(45°E—80°E,0°-30°N),and the inversion results are compared with those of conventional altimeter data.The results demonstrate that the difference between the meridian and prime components derived from the inversion of swoT wide-swath data is minimal,significantly outperforming the inversion results of conventional nadir altimeter data.The advantage of swoT wide-swath altimeter lies in its ability to use the multi-directional geoid slope at any sea surface measurement point to invert the components in the meridian and prime directions.To investigate the impact of this advantage on inversion precision,this paper employs a method to calculate the gradient of the geoid in multiple directions to invert DoV components.The improvement effect of calculating the gradient of the geoid in multiple directions on the precision of DoV component is analyzed.It is found that the accuracy of DoV inversion has significantly improved with the increase of geodetic gradient calculation direction.In addition,the effects of various errors and grid spacing in SwoT wide sea surface height data on the precision of Dov inversion are also analyzed.
基金Supported by National Natural Science Foundation of China(Grant Nos.51469014,51676003)National Key Research and Development Program of China(Grant No.20016YFB0200901)
文摘To improve the accuracy and reduce the calcu- lation cost for the inverse problem of centrifugal pump impeller, the new inverse method based on proper orthog- onal decomposition (POD) is proposed. The pump blade shape is parameterized by quartic Bezier curve, and the initial snapshots is generated by introducing the perturbation of the blade shape control parameters. The internal flow field and its hydraulic performance is predicted by CFD method. The snapshots vector includes the blade shape parameter and the distribution of blade load. The POD basis for the snap- shots set are deduced by proper orthogonal decomposition. The sample vector set is expressed in terms of the linear combination of the orthogonal basis. The objective blade shape corresponding to the objective distribution of blade load is obtained by least square fit. The Iterative correction algorithm for the centrifugal pump blade inverse method based on POD is proposed. The objective blade load dis- tributions are corrected according to the difference of the CFD result and the POD result. The two dimensional and three dimensional blade calculation cases show that the proposed centrifugal pump blade inverse method based on POD have good convergence and high accuracy, and thecalculation cost is greatly reduced. After two iterations, the deviation of the blade load and the pump hydraulic perfor- mance are limited within 4.0% and 6.0% individually for most of the flow rate range. This paper provides a promising inverse method for centrifugal pump impeller, which will benefit the hydraulic optimization of centrifugal pump.
基金supported by National Basic Research Program of China(No.2005CB724105)National Natural Science Foundation of China (No.10477010)National High Technical Research and Development Program of China(No.2007AA04Z141)
文摘Inverse method was used in single crystal superalloy DD6 processing simulation during solidification. Numerical modeling coupled with experiments has been used to estimate the interface heat transfer coefficient (IHTC) between the surface of slab casting and inner mold. Calculated temperature dependent values of IHTC were obtained from a numerical solution. The calculated temperatures agreed well with the measurement of cooling profile.
基金the National Natural Science Foundation of China (No. 50979058)the Special Research Fund for the Doctoral Program of Higher Education (No. 20090073110012)
文摘Autonomous underwater gliders are highly effcient,buoyancy-driven,winged autonomous underwater vehicles. Their dynamics are multivariable nonlinear systems. In addition,the gliders are underactuated and diffcult to maneuver,and also dependent on their operational environment. To confront these problems and to design an effective controller,the inverse system method was used to decouple the original system into two independent single variable linear subsystems. The stability of the zero dynamics was analyzed,and an additional closed-loop controller for each linear subsystem was designed by sliding mode control method to form a type of composite controller. Simulation results demonstrate that the derived nonlinear controller is able to cope with the aforementioned problems simultaneously and satisfactorily.
文摘The heat transfer coefficient in a multidimensional heat conduction problem is obtained from the solution of the inverse heat conduction problem based on the thermographic temperature measurement. The modified one-dimensional correction method (MODCM), along with the finite volume method, is employed for both two- and three-dimensional inverse problems. A series of numerical experiments are conducted in order to verify the effectiveness of the method. In addition, the effect of the temperature measurement error, the ending criterion of the iteration, etc. on the result of the inverse problem is investigated. It is proved that the method is a simple, stable and accurate one that can solve successfully the inverse heat conduction problem.
基金the Major State Basic Research Program of China un der contract No.G1999043805.
文摘On the basis of hydrographic data obtained in August 2000 cruise, the circulation in the South China Sea (SCS) is computed by the modified inverse method in combination with SSH data from TOPEX/ERS-2 analysis. For study of the dynamical mechanism, which causes the pattern of summer circulation in the SCS, the diagnostic model (Yuan et al. 1982. Acta Oceanologica Sinica,4(1):1-11; Yuan and Su. 1992. Numerical Computation of Physical Oceanography.474-542) is used to simulate numerically the summer circulation in the SCS. The following results have been obtained. (1) The central and southwestern SCSs are dominated mainly by anticy-clonic circulation systems. They are mainly as follows. 1) There is strong anticyclonic eddy southeast of Vietnam (W1). Its horizontal scale is about 300 km, and it extends vertically from the surface to the about 1 000 m level. 2) There are a warm eddy W2 southeast of Zhongsha Islands and the anticyclonic circulation system W3 west off the Luzon Island. 3) There is a stronger cyclonic eddy C1 between the anticyclonic eddies W1 and W2.4) A strong northward coastal jet is present near the coast of Vietnam, and separates from the coast of Vietnam at about 12° N to the northeast.(2)The northern SCS is dominated mainly by a cyclonic circulation system. There is a cyclonic circulation system near and north of Section N2. (3) The southeastern SCS is dominated mainly by the cyclonic circulation system. (4) Comparing the results of circulation in the SCS during the summer of 2000 with those during the summer of 1998, it is found that they agree qualitatively, but there is the some difference between them in quantity.This shows that the circulation in the SCS has obviously seasonal feature. (5) The dynamical mechanism which products the basic pattern of summer circulation is because the following two reasons: 1) the joint effect of the baroclinity and relief (JEBAR) is essential dynamical cause; and 2) it is next important dynamical cause that the interaction between the wind stress and bottom topography under the southerly monsoon. (6) Comparing the hydrographic structure and distribution of stream functions with the SSH data from TOPEX/ERS-2 analysis in the SCS during August of 2000, they agree qualitatively.
基金supported by the National Natural Sciences Foundation of China (No. 50776075)
文摘Inverse heat conduction method (IHCM) is one of the most effective approaches to obtaining the boiling heat transfer coefficient from measured results. This paper focuses on its application in cryogenic boiling heat transfer. Experiments were conducted on the heat transfer of a stainless steel block in a liquid nitrogen bath, with the assumption of a 1D conduction condition to realize fast acquisition of the temperature of the test points inside the block. With the inverse-heat conduction theory and the explicit finite difference model, a solving program was developed to calculate the heat flux and the boiling heat transfer coefficient of a stainless steel block in liquid nitrogen bath based on the temperature acquisition data. Considering the oscillating data and some unsmooth transition points in the inverse-heat-conduction calculation result of the heat-transfer coefficient, a two-step data-fitting procedure was proposed to obtain the expression for the boiling heat transfer coefficients. The coefficient was then verified for accuracy by a comparison between the simulation results using this expression and the verifying experimental results of a stainless steel block. The maximum error with a revised segment fitting is around 6%, which verifies the feasibility of using IHCM to measure the boiling heat transfer coefficient in liquid nitrogen bath.
文摘The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.
基金The Major State Basic Research Development Program Grant (2005CB321701)the Heilongjiang Education Committee Grant (11551364) of China
文摘This paper is concerned with the inverse scattering problems for Schrdinger equations with compactly supported potentials.For purpose of reconstructing the support of the potential,we derive a factorization of the scattering amplitude operator A and prove that the ranges of (A* A) ^1/4 and G which maps more general incident fields than plane waves into the scattering amplitude coincide.As an application we characterize the support of the potential using only the spectral data of the operator A.