Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise re...Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.展开更多
Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of un...Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of unsteady flows.To enhance computational efficiency,we propose the Implicit-Explicit Two-Step Runge-Kutta(IMEX-TSRK)time-stepping discretization methods for unsteady flows,and develop a novel adaptive algorithm that correctly partitions spatial regions to apply implicit or explicit methods.The novel adaptive IMEX-TSRK schemes effectively handle the numerical stiffness of the small grid size and improve computational efficiency.Compared to implicit and explicit Runge-Kutta(RK)schemes,the IMEX-TSRK methods achieve the same order of accuracy with fewer first derivative calculations.Numerical case tests demonstrate that the IMEX-TSRK methods maintain numerical stability while enhancing computational efficiency.Specifically,in high Reynolds number flows,the computational efficiency of the IMEX-TSRK methods surpasses that of explicit RK schemes by more than one order of magnitude,and that of implicit RK schemes several times over.展开更多
Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation o...Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation of new holes.However,most pertinent investigations in the field rely predominantly on fixed background mesh,which is never remeshed.Consequently,the mesh element partitioned by material interface during the optimization process necessitates approximation by using artificial interpolation models to obtain its element stiffness or other properties.This paper introduces a novel approach to topology op-timization by integrating the PLSM with body-fitted adaptive mesh and Helmholtz-type filter.Primarily,combining the PLSM with body-fitted adaptive mesh enables the regeneration of mesh based on the zero level-set interface.This not only precludes the direct traversal of the material interface through the mesh element during the topology optimization process,but also improves the accuracy of calculation.Additionally,the incorporation of a Helmholtz-type partial differential equation filter,relying solely on mesh information essential for finite element discretization,serves to regulate the topological complexity and the minimum feature size of the optimized structure.Leveraging these advantages,the topology optimization program demonstrates its versa-tility by successfully addressing various design problems,encompassing the minimum mean compliance problem and minimum energy dissipation problem.Ultimately,the result of numerical example indicates that the optimized structure exhibits a dis-tinct and smooth boundary,affirming the effective control over both topological complexity and the minimum feature size of the optimized structure.展开更多
We present a gain adaptive tuning method for fiber Raman amplifier(FRA) using two-stage neural networks(NNs) and double weights updates. After training the connection weights of two-stage NNs separately in training ph...We present a gain adaptive tuning method for fiber Raman amplifier(FRA) using two-stage neural networks(NNs) and double weights updates. After training the connection weights of two-stage NNs separately in training phase, the connection weights of the unified NN are updated again in verification phase according to error between the predicted and target gains to eliminate the inherent error of the NNs. The simulation results show that the mean of root mean square error(RMSE) and maximum error of gains are 0.131 d B and 0.281 d B, respectively. It shows that the method can realize adaptive adjustment function of FRA gain with high accuracy.展开更多
To accomplish the reliability analyses of the correlation of multi-analytical objectives,an innovative framework of Dimensional Synchronous Modeling(DSM)and correlation analysis is developed based on the stepwise mode...To accomplish the reliability analyses of the correlation of multi-analytical objectives,an innovative framework of Dimensional Synchronous Modeling(DSM)and correlation analysis is developed based on the stepwise modeling strategy,cell array operation principle,and Copula theory.Under this framework,we propose a DSM-based Enhanced Kriging(DSMEK)algorithm to synchronously derive the modeling of multi-objective,and explore an adaptive Copula function approach to analyze the correlation among multiple objectives and to assess the synthetical reliability level.In the proposed DSMEK and adaptive Copula methods,the Kriging model is treated as the basis function of DSMEK model,the Multi-Objective Snake Optimizer(MOSO)algorithm is used to search the optimal values of hyperparameters of basis functions,the cell array operation principle is adopted to establish a whole model of multiple objectives,the goodness of fit is utilized to determine the forms of Copula functions,and the determined Copula functions are employed to perform the reliability analyses of the correlation of multi-analytical objectives.Furthermore,three examples,including multi-objective complex function approximation,aeroengine turbine bladeddisc multi-failure mode reliability analyses and aircraft landing gear system brake temperature reliability analyses,are performed to verify the effectiveness of the proposed methods,from the viewpoints of mathematics and engineering.The results show that the DSMEK and adaptive Copula approaches hold obvious advantages in terms of modeling features and simulation performance.The efforts of this work provide a useful way for the modeling of multi-analytical objectives and synthetical reliability analyses of complex structure/system with multi-output responses.展开更多
To address the issues of single warning indicators,fixed thresholds,and insufficient adaptability in coal and gas outburst early warning models,this study proposes a dynamic early warning model for gas outbursts based...To address the issues of single warning indicators,fixed thresholds,and insufficient adaptability in coal and gas outburst early warning models,this study proposes a dynamic early warning model for gas outbursts based on adaptive fractal dimension characterization.By analyzing the nonlinear characteristics of gas concentration data,an adaptive window fractal analysis method is introduced.Combined with boxcounting dimension and variation of box dimension metrics,a cross-scale dynamic warning model for disaster prevention is established.The implementation involves three key phases:First,wavelet denoising and interpolation methods are employed for raw data preprocessing,followed by validation of fractal characteristics.Second,an adaptive window cross-scale fractal dimension method is proposed to calculate the box-counting dimension of gas concentration,enabling effective capture of multi-scale complex features.Finally,dynamic threshold partitioning is achieved through membership functions and the 3σprinciple,establishing a graded classification standard for the mine gas disaster(MGD)index.Validated through engineering applications at Shoushan#1 Coal Mine in Henan Province,the results demonstrate that the adaptive window fractal dimension curve exhibits significantly enhanced fluctuation characteristics compared to fixed window methods,with local feature detection capability improved and warning accuracy reaching 86.9%.The research reveals that this model effectively resolves the limitations of traditional methods in capturing local features and dependency on subjective thresholds through multiindicator fusion and threshold optimization,providing both theoretical foundation and practical tool for coal mine gas outburst early warning.展开更多
This study explored the observation strategy and effectiveness of synoptic-scale adaptive observations for improving sea fog prediction in coastal regions around the Bohai Sea based on a poorly predicted fog event wit...This study explored the observation strategy and effectiveness of synoptic-scale adaptive observations for improving sea fog prediction in coastal regions around the Bohai Sea based on a poorly predicted fog event with cold-front synoptic pattern(CFSP).An ensemble Kalman filter data assimilation system for the Weather Research and Forecasting model was adopted with ensemble sensitivity analysis(ESA).By comparing observation impacts(estimated from a 40-member ensemble with ESA)among different meteorological observation variables and pressure levels,the temperature at 850 hPa and surface layer(850 hPa-and-surface temperature)was selected as the target observation type.Additionally,the area with large observation impacts for this observation type was predicted in the transition region of the surface low–high system.This area developed southward with the low and moved eastward with the low–high system,which could be explained by the main features of CFSP.Moreover,both experiments assimilating synthetic and real observations showed that assimilating 850 hPa-and-surface temperature observations generally yielded better fog coverage forecasts in areas with greater observation impacts than areas with smaller impacts.However,the effectiveness of adaptive observations was reduced when real observations rather than synthetic observations were assimilated,which is possibly due to factors such as observation and model errors.The main conclusions above were verified by another typical fog event with CFSP characteristics.Results of this study highlight the importance of improved initial conditions in the transition region of the low–high system for improving fog prediction and provide scientific guidance for implementing an observation network for fog forecasting over the Bohai Sea.展开更多
3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic m...3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.展开更多
A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic gr...A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.展开更多
High fidelity analysis models,which are beneficial to improving the design quality,have been more and more widely utilized in the modern engineering design optimization problems.However,the high fidelity analysis mode...High fidelity analysis models,which are beneficial to improving the design quality,have been more and more widely utilized in the modern engineering design optimization problems.However,the high fidelity analysis models are so computationally expensive that the time required in design optimization is usually unacceptable.In order to improve the efficiency of optimization involving high fidelity analysis models,the optimization efficiency can be upgraded through applying surrogates to approximate the computationally expensive models,which can greately reduce the computation time.An efficient heuristic global optimization method using adaptive radial basis function（RBF） based on fuzzy clustering（ARFC） is proposed.In this method,a novel algorithm of maximin Latin hypercube design using successive local enumeration（SLE） is employed to obtain sample points with good performance in both space-filling and projective uniformity properties,which does a great deal of good to metamodels accuracy.RBF method is adopted for constructing the metamodels,and with the increasing the number of sample points the approximation accuracy of RBF is gradually enhanced.The fuzzy c-means clustering method is applied to identify the reduced attractive regions in the original design space.The numerical benchmark examples are used for validating the performance of ARFC.The results demonstrates that for most application examples the global optima are effectively obtained and comparison with adaptive response surface method（ARSM） proves that the proposed method can intuitively capture promising design regions and can efficiently identify the global or near-global design optimum.This method improves the efficiency and global convergence of the optimization problems,and gives a new optimization strategy for engineering design optimization problems involving computationally expensive models.展开更多
The paper presents the utilization of the adaptive Delaunay triangulation in the finite element modeling of two dimensional crack propagation problems, including detailed description of the proposed procedure which co...The paper presents the utilization of the adaptive Delaunay triangulation in the finite element modeling of two dimensional crack propagation problems, including detailed description of the proposed procedure which consists of the Delaunay triangulation algorithm and an adaptive remeshing technique. The adaptive remeshing technique generates small elements around crack tips and large elements in the other regions. The resulting stress intensity factors and simulated crack propagation behavior are used to evaluate the effectiveness of the procedure. Three sample problems of a center cracked plate, a single edge cracked plate and a compact tension specimen, are simulated and their results assessed.展开更多
An automatically adaptive element free method is presented to analyze the seismic response of liquefiable soils. The method is based on the element free Galerkin method （EFGM） and the fission procedure that is part ...An automatically adaptive element free method is presented to analyze the seismic response of liquefiable soils. The method is based on the element free Galerkin method （EFGM） and the fission procedure that is part of h-refinement, indicated by error estimation. In the proposed method, a posteriori error estimate procedure that depends on the energy norm of stress and the T-Belytschko （TB） stress recovery scheme is incorporated. The effective cyclic elasto-plastic constitutive model is used to describe the nonlinear behavior of the saturated soil. The governing equations are established by u-p formulation. The proposed method can effectively avoid the volumetric locking due to large deformation that usually occurs in numerical computations using the finite element method （FEM）. The efficiency of the proposed method is demonstrated by evaluating the seismic response of an embankment and comparing it to results obtained through FEM. It is shown that the proposed method provides an accurate seismic analysis of saturated soil that includes the effects of liquefaction .展开更多
In this paper,a robust adaptive controller is designed for a guided spinning rocket,whose dynamics presents the characteristics of pitch-yaw cross coupling,fast time-varying aerodynamics parameters and wide flight env...In this paper,a robust adaptive controller is designed for a guided spinning rocket,whose dynamics presents the characteristics of pitch-yaw cross coupling,fast time-varying aerodynamics parameters and wide flight envelop.First,a coupled nonlinear six-degree-of-freedom equation of motion for a guided spinning rocket is developed,and the lateral acceleration motion is modeled as a control plant with time-varying matched uncertainties and unmodeled dynamics.Then,a robust adaptive control method is proposed by combining Bregman divergence and variational method to achieve fast adaption and maintain bounded tracking.The stability of the resulting closed-loop system is proved,and the ultimate bound and convergence rate are also analyzed.Finally,numerical simulations are performed for a single operating point and the whole flight trajectory to show the robustness and adaptability of the proposed method with respect to timevarying uncertainties and unmodeled dynamics.展开更多
In order to establish an adaptive turbo-shaft engine model with high accuracy, a new modeling method based on parameter selection （PS） algorithm and multi-input multi-output recursive reduced least square support ve...In order to establish an adaptive turbo-shaft engine model with high accuracy, a new modeling method based on parameter selection （PS） algorithm and multi-input multi-output recursive reduced least square support vector regression （MRR-LSSVR） machine is proposed. Firstly, the PS algorithm is designed to choose the most reasonable inputs of the adaptive module. During this process, a wrapper criterion based on least square support vector regression （LSSVR） machine is adopted, which can not only reduce computational complexity but also enhance generalization performance. Secondly, with the input variables determined by the PS algorithm, a mapping model of engine parameter estimation is trained off-line using MRR-LSSVR, which has a satisfying accuracy within 5＆. Finally, based on a numerical simulation platform of an integrated helicopter/ turbo-shaft engine system, an adaptive turbo-shaft engine model is developed and tested in a certain flight envelope. Under the condition of single or multiple engine components being degraded, many simulation experiments are carried out, and the simulation results show the effectiveness and validity of the proposed adaptive modeling method.展开更多
Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid...Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.展开更多
An adaptive finite element-element-free Galerkin (FE-EFG) coupling method is proposed and developed for the numerical simulation of bulk metal forming processes. This approach is able to adaptively convert distorted F...An adaptive finite element-element-free Galerkin (FE-EFG) coupling method is proposed and developed for the numerical simulation of bulk metal forming processes. This approach is able to adaptively convert distorted FE elements to EFG domain in analysis. A new scheme to implement adaptive conversion and coupling is presented. The coupling method takes both advantages of finite element method (FEM) and meshless methods. It is capable of handling large deformations with no need of remeshing procedures, while it is computationally more efficient than those full meshless methods. The effectiveness of the proposed method is demonstrated with the numerical simulations of the bulk metal forming processes including forging and extrusion.展开更多
The trajectory tracking control problem is addressed for autonomous underwater vehicle(AUV) in marine environ?ment, with presence of the influence of the uncertain factors including ocean current disturbance, dynamic ...The trajectory tracking control problem is addressed for autonomous underwater vehicle(AUV) in marine environ?ment, with presence of the influence of the uncertain factors including ocean current disturbance, dynamic modeling uncertainty, and thrust model errors. To improve the trajectory tracking accuracy of AUV, an adaptive backstepping terminal sliding mode control based on recurrent neural networks(RNN) is proposed. Firstly, considering the inaccu?rate of thrust model of thruster, a Taylor’s polynomial is used to obtain the thrust model errors. And then, the dynamic modeling uncertainty and thrust model errors are combined into the system model uncertainty(SMU) of AUV; through the RNN, the SMU and ocean current disturbance are classified, approximated online. Finally, the weights of RNN and other control parameters are adjusted online based on the backstepping terminal sliding mode controller. In addition, a chattering?reduction method is proposed based on sigmoid function. In chattering?reduction method, the sigmoid function is used to realize the continuity of the sliding mode switching function, and the sliding mode switching gain is adjusted online based on the exponential form of the sliding mode function. Based on the Lyapu?nov theory and Barbalat’s lemma, it is theoretically proved that the AUV trajectory tracking error can quickly converge to zero in the finite time. This research proposes a trajectory tracking control method of AUV, which can e ectively achieve high?precision trajectory tracking control of AUV under the influence of the uncertain factors. The feasibility and e ectiveness of the proposed method is demonstrated with trajectory tracking simulations and pool?experi?ments of AUV.展开更多
There are two separate traditional approaches to model contact problems: continuum and atomistic theory. Continuum theory is successfully used in many domains, but when the scale of the model comes to nanometer, conti...There are two separate traditional approaches to model contact problems: continuum and atomistic theory. Continuum theory is successfully used in many domains, but when the scale of the model comes to nanometer, continuum approximation meets challenges. Atomistic theory can catch the detailed behaviors of an individual atom by using molecular dynamics (MD) or quantum mechanics, although accurately, it is usually time-consuming. A multiscale method coupled MD and finite element (FE) is presented. To mesh the FE region automatically, an adaptive method based on the strain energy gradient is introduced to the multiscale method to constitute an adaptive multiscale method. Utilizing the proposed method, adhesive contacts between a rigid cylinder and an elastic substrate are studied, and the results are compared with full MD simulations. The process of FE meshes refinement shows that adaptive multiscale method can make FE mesh generation more flexible. Comparison of the displacements of boundary atoms in the overlap region with the results from full MD simulations indicates that adaptive multiscale method can transfer displacements effectively. Displacements of atoms and FE nodes on the center line of the multiscale model agree well with that of atoms in full MD simulations, which shows the continuity in the overlap region. Furthermore, the Von Mises stress contours and contact force distributions in the contact region are almost same as full MD simulations. The method presented combines multiscale method and adaptive technique, and can provide a more effective way to multiscale method and to the investigation on nanoscale contact problems.展开更多
In this paper, the generalized nonlinear Schrodinger equation （GNLSE） is solved by an adaptive split-step Fourier method （ASSFM）. It is found that ASSFM must be used to solve GNLSE to ensure precision when the sol...In this paper, the generalized nonlinear Schrodinger equation （GNLSE） is solved by an adaptive split-step Fourier method （ASSFM）. It is found that ASSFM must be used to solve GNLSE to ensure precision when the soliton selffrequency shift is remarkable and the photonic crystal fibre （PCF） parameters vary with the frequency considerably. The precision of numerical simulation by using ASSFM is higher than that by using split-step Fourier method in the process of laser pulse propagation in PCFs due to the fact that the variation of fibre parameters with the peak frequency in the pulse spectrum can be taken into account fully.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.12302435 and 12221002)。
文摘Shock wave caused by a sudden release of high-energy,such as explosion and blast,usually affects a significant range of areas.The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource.This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size.In this work,a variable-domain-size adaptive mesh enlargement(vAME)method is developed based on the proposed adaptive mesh enlargement(AME)method for modeling multi-explosives explosion problems.The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions,rather than in all directions as in AME method.A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy.An estimate of speedup ratio is analyzed for further efficiency comparison.Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave,reflected wave,and Mach wave.Additionally,the vAME method is employed to validate the accuracy,as well as to investigate the performance of the fluid field and shock wave propagation,considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass.The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations.The vAME method yields a competitive efficiency,increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method,respectively.It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.
基金supported by the National Natural Science Foundation of China(No.92252201)the Fundamental Research Funds for the Central Universitiesthe Academic Excellence Foundation of Beihang University(BUAA)for PhD Students。
文摘Efficient and accurate simulation of unsteady flow presents a significant challenge that needs to be overcome in computational fluid dynamics.Temporal discretization method plays a crucial role in the simulation of unsteady flows.To enhance computational efficiency,we propose the Implicit-Explicit Two-Step Runge-Kutta(IMEX-TSRK)time-stepping discretization methods for unsteady flows,and develop a novel adaptive algorithm that correctly partitions spatial regions to apply implicit or explicit methods.The novel adaptive IMEX-TSRK schemes effectively handle the numerical stiffness of the small grid size and improve computational efficiency.Compared to implicit and explicit Runge-Kutta(RK)schemes,the IMEX-TSRK methods achieve the same order of accuracy with fewer first derivative calculations.Numerical case tests demonstrate that the IMEX-TSRK methods maintain numerical stability while enhancing computational efficiency.Specifically,in high Reynolds number flows,the computational efficiency of the IMEX-TSRK methods surpasses that of explicit RK schemes by more than one order of magnitude,and that of implicit RK schemes several times over.
基金supported by the National Natural Science Foundation of China(Grant Nos.12372200 and 12072242).
文摘Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation of new holes.However,most pertinent investigations in the field rely predominantly on fixed background mesh,which is never remeshed.Consequently,the mesh element partitioned by material interface during the optimization process necessitates approximation by using artificial interpolation models to obtain its element stiffness or other properties.This paper introduces a novel approach to topology op-timization by integrating the PLSM with body-fitted adaptive mesh and Helmholtz-type filter.Primarily,combining the PLSM with body-fitted adaptive mesh enables the regeneration of mesh based on the zero level-set interface.This not only precludes the direct traversal of the material interface through the mesh element during the topology optimization process,but also improves the accuracy of calculation.Additionally,the incorporation of a Helmholtz-type partial differential equation filter,relying solely on mesh information essential for finite element discretization,serves to regulate the topological complexity and the minimum feature size of the optimized structure.Leveraging these advantages,the topology optimization program demonstrates its versa-tility by successfully addressing various design problems,encompassing the minimum mean compliance problem and minimum energy dissipation problem.Ultimately,the result of numerical example indicates that the optimized structure exhibits a dis-tinct and smooth boundary,affirming the effective control over both topological complexity and the minimum feature size of the optimized structure.
基金supported by the Natural Science Research Project of Colleges and Universities in Anhui Province (No.KJ2021A0479)the Science Research Program of Anhui University of Finance and Economics (No.ACKYC22082)。
文摘We present a gain adaptive tuning method for fiber Raman amplifier(FRA) using two-stage neural networks(NNs) and double weights updates. After training the connection weights of two-stage NNs separately in training phase, the connection weights of the unified NN are updated again in verification phase according to error between the predicted and target gains to eliminate the inherent error of the NNs. The simulation results show that the mean of root mean square error(RMSE) and maximum error of gains are 0.131 d B and 0.281 d B, respectively. It shows that the method can realize adaptive adjustment function of FRA gain with high accuracy.
基金co-supported by the National Natural Science Foundation of China(Nos.52405293,52375237)China Postdoctoral Science Foundation(No.2024M754219)Shaanxi Province Postdoctoral Research Project Funding,China。
文摘To accomplish the reliability analyses of the correlation of multi-analytical objectives,an innovative framework of Dimensional Synchronous Modeling(DSM)and correlation analysis is developed based on the stepwise modeling strategy,cell array operation principle,and Copula theory.Under this framework,we propose a DSM-based Enhanced Kriging(DSMEK)algorithm to synchronously derive the modeling of multi-objective,and explore an adaptive Copula function approach to analyze the correlation among multiple objectives and to assess the synthetical reliability level.In the proposed DSMEK and adaptive Copula methods,the Kriging model is treated as the basis function of DSMEK model,the Multi-Objective Snake Optimizer(MOSO)algorithm is used to search the optimal values of hyperparameters of basis functions,the cell array operation principle is adopted to establish a whole model of multiple objectives,the goodness of fit is utilized to determine the forms of Copula functions,and the determined Copula functions are employed to perform the reliability analyses of the correlation of multi-analytical objectives.Furthermore,three examples,including multi-objective complex function approximation,aeroengine turbine bladeddisc multi-failure mode reliability analyses and aircraft landing gear system brake temperature reliability analyses,are performed to verify the effectiveness of the proposed methods,from the viewpoints of mathematics and engineering.The results show that the DSMEK and adaptive Copula approaches hold obvious advantages in terms of modeling features and simulation performance.The efforts of this work provide a useful way for the modeling of multi-analytical objectives and synthetical reliability analyses of complex structure/system with multi-output responses.
基金funded by the National Key Research and Development ProgramFund for Young Scientists(No.2021YFC2900400)+5 种基金the National Natural Science Foundation of China(No.52304123)Fundamental Research Funds for the Central Universities(No.2024CDJXY025)Sichuan-Chongqing Science and Technology Innovation Cooperation Program Project(No.CSTB2024TIAD-CYKJCXX0016)Postdoctoral Research Foundation of China(No.2023M730412)Postdoctoral Fellowship Program of China Postdoctoral Science Foundation(No.GZB20230914)Chongqing Outstanding Youth Science Foundation Program(No.CSTB2023NSCQ-JQX0027)。
文摘To address the issues of single warning indicators,fixed thresholds,and insufficient adaptability in coal and gas outburst early warning models,this study proposes a dynamic early warning model for gas outbursts based on adaptive fractal dimension characterization.By analyzing the nonlinear characteristics of gas concentration data,an adaptive window fractal analysis method is introduced.Combined with boxcounting dimension and variation of box dimension metrics,a cross-scale dynamic warning model for disaster prevention is established.The implementation involves three key phases:First,wavelet denoising and interpolation methods are employed for raw data preprocessing,followed by validation of fractal characteristics.Second,an adaptive window cross-scale fractal dimension method is proposed to calculate the box-counting dimension of gas concentration,enabling effective capture of multi-scale complex features.Finally,dynamic threshold partitioning is achieved through membership functions and the 3σprinciple,establishing a graded classification standard for the mine gas disaster(MGD)index.Validated through engineering applications at Shoushan#1 Coal Mine in Henan Province,the results demonstrate that the adaptive window fractal dimension curve exhibits significantly enhanced fluctuation characteristics compared to fixed window methods,with local feature detection capability improved and warning accuracy reaching 86.9%.The research reveals that this model effectively resolves the limitations of traditional methods in capturing local features and dependency on subjective thresholds through multiindicator fusion and threshold optimization,providing both theoretical foundation and practical tool for coal mine gas outburst early warning.
基金supported by the National Natural Science Foundation of China(Grant No.41705081)the Shandong Natural Science Foundation Project(Grant No.ZR2019ZD12)the Laoshan Laboratory(Grant No.LSKJ202202203).
文摘This study explored the observation strategy and effectiveness of synoptic-scale adaptive observations for improving sea fog prediction in coastal regions around the Bohai Sea based on a poorly predicted fog event with cold-front synoptic pattern(CFSP).An ensemble Kalman filter data assimilation system for the Weather Research and Forecasting model was adopted with ensemble sensitivity analysis(ESA).By comparing observation impacts(estimated from a 40-member ensemble with ESA)among different meteorological observation variables and pressure levels,the temperature at 850 hPa and surface layer(850 hPa-and-surface temperature)was selected as the target observation type.Additionally,the area with large observation impacts for this observation type was predicted in the transition region of the surface low–high system.This area developed southward with the low and moved eastward with the low–high system,which could be explained by the main features of CFSP.Moreover,both experiments assimilating synthetic and real observations showed that assimilating 850 hPa-and-surface temperature observations generally yielded better fog coverage forecasts in areas with greater observation impacts than areas with smaller impacts.However,the effectiveness of adaptive observations was reduced when real observations rather than synthetic observations were assimilated,which is possibly due to factors such as observation and model errors.The main conclusions above were verified by another typical fog event with CFSP characteristics.Results of this study highlight the importance of improved initial conditions in the transition region of the low–high system for improving fog prediction and provide scientific guidance for implementing an observation network for fog forecasting over the Bohai Sea.
基金The authors thank the funds supported by the China National Nuclear Corporation under Grants Nos.WUQNYC2101 and WUHTLM2101-04National Natural Science Foundation of China(42074132,42274154).
文摘3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.
基金Projects(51161011,11364024)supported by the National Natural Science Foundation of ChinaProject(1204GKCA065)supported by the Key Technology R&D Program of Gansu Province,China+1 种基金Project(201210)supported by the Fundamental Research Funds for the Universities of Gansu Province,ChinaProject(J201304)supported by the Funds for Distinguished Young Scientists of Lanzhou University of Technology,China
文摘A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.
基金supported by National Natural Science Foundation of China (Grant Nos. 50875024,51105040)Excellent Young Scholars Research Fund of Beijing Institute of Technology,China (Grant No.2010Y0102)Defense Creative Research Group Foundation of China(Grant No. GFTD0803)
文摘High fidelity analysis models,which are beneficial to improving the design quality,have been more and more widely utilized in the modern engineering design optimization problems.However,the high fidelity analysis models are so computationally expensive that the time required in design optimization is usually unacceptable.In order to improve the efficiency of optimization involving high fidelity analysis models,the optimization efficiency can be upgraded through applying surrogates to approximate the computationally expensive models,which can greately reduce the computation time.An efficient heuristic global optimization method using adaptive radial basis function（RBF） based on fuzzy clustering（ARFC） is proposed.In this method,a novel algorithm of maximin Latin hypercube design using successive local enumeration（SLE） is employed to obtain sample points with good performance in both space-filling and projective uniformity properties,which does a great deal of good to metamodels accuracy.RBF method is adopted for constructing the metamodels,and with the increasing the number of sample points the approximation accuracy of RBF is gradually enhanced.The fuzzy c-means clustering method is applied to identify the reduced attractive regions in the original design space.The numerical benchmark examples are used for validating the performance of ARFC.The results demonstrates that for most application examples the global optima are effectively obtained and comparison with adaptive response surface method（ARSM） proves that the proposed method can intuitively capture promising design regions and can efficiently identify the global or near-global design optimum.This method improves the efficiency and global convergence of the optimization problems,and gives a new optimization strategy for engineering design optimization problems involving computationally expensive models.
文摘The paper presents the utilization of the adaptive Delaunay triangulation in the finite element modeling of two dimensional crack propagation problems, including detailed description of the proposed procedure which consists of the Delaunay triangulation algorithm and an adaptive remeshing technique. The adaptive remeshing technique generates small elements around crack tips and large elements in the other regions. The resulting stress intensity factors and simulated crack propagation behavior are used to evaluate the effectiveness of the procedure. Three sample problems of a center cracked plate, a single edge cracked plate and a compact tension specimen, are simulated and their results assessed.
基金the National Natural Science Foundation of China Under Grant No. 50779003and50639010Natural Science Foundation of Liaoning Proviance Under Grant No. 20071085
文摘An automatically adaptive element free method is presented to analyze the seismic response of liquefiable soils. The method is based on the element free Galerkin method （EFGM） and the fission procedure that is part of h-refinement, indicated by error estimation. In the proposed method, a posteriori error estimate procedure that depends on the energy norm of stress and the T-Belytschko （TB） stress recovery scheme is incorporated. The effective cyclic elasto-plastic constitutive model is used to describe the nonlinear behavior of the saturated soil. The governing equations are established by u-p formulation. The proposed method can effectively avoid the volumetric locking due to large deformation that usually occurs in numerical computations using the finite element method （FEM）. The efficiency of the proposed method is demonstrated by evaluating the seismic response of an embankment and comparing it to results obtained through FEM. It is shown that the proposed method provides an accurate seismic analysis of saturated soil that includes the effects of liquefaction .
基金supported by the National Natural Science Foundation of China (No. 11532002)。
文摘In this paper,a robust adaptive controller is designed for a guided spinning rocket,whose dynamics presents the characteristics of pitch-yaw cross coupling,fast time-varying aerodynamics parameters and wide flight envelop.First,a coupled nonlinear six-degree-of-freedom equation of motion for a guided spinning rocket is developed,and the lateral acceleration motion is modeled as a control plant with time-varying matched uncertainties and unmodeled dynamics.Then,a robust adaptive control method is proposed by combining Bregman divergence and variational method to achieve fast adaption and maintain bounded tracking.The stability of the resulting closed-loop system is proved,and the ultimate bound and convergence rate are also analyzed.Finally,numerical simulations are performed for a single operating point and the whole flight trajectory to show the robustness and adaptability of the proposed method with respect to timevarying uncertainties and unmodeled dynamics.
基金co-supported by Aeronautical Science Foundation of China (No. 2010ZB52011)Funding of Jiangsu Innovation Program for Graduate Education (No.CXLX11_0213)
文摘In order to establish an adaptive turbo-shaft engine model with high accuracy, a new modeling method based on parameter selection （PS） algorithm and multi-input multi-output recursive reduced least square support vector regression （MRR-LSSVR） machine is proposed. Firstly, the PS algorithm is designed to choose the most reasonable inputs of the adaptive module. During this process, a wrapper criterion based on least square support vector regression （LSSVR） machine is adopted, which can not only reduce computational complexity but also enhance generalization performance. Secondly, with the input variables determined by the PS algorithm, a mapping model of engine parameter estimation is trained off-line using MRR-LSSVR, which has a satisfying accuracy within 5＆. Finally, based on a numerical simulation platform of an integrated helicopter/ turbo-shaft engine system, an adaptive turbo-shaft engine model is developed and tested in a certain flight envelope. Under the condition of single or multiple engine components being degraded, many simulation experiments are carried out, and the simulation results show the effectiveness and validity of the proposed adaptive modeling method.
基金Projects(2006AA06Z105, 2007AA06Z134) supported by the National High-Tech Research and Development Program of ChinaProjects(2007, 2008) supported by China Scholarship Council (CSC)
文摘Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.
基金Project (No. 50575143) supported by the National Natural Science Foundation of China
文摘An adaptive finite element-element-free Galerkin (FE-EFG) coupling method is proposed and developed for the numerical simulation of bulk metal forming processes. This approach is able to adaptively convert distorted FE elements to EFG domain in analysis. A new scheme to implement adaptive conversion and coupling is presented. The coupling method takes both advantages of finite element method (FEM) and meshless methods. It is capable of handling large deformations with no need of remeshing procedures, while it is computationally more efficient than those full meshless methods. The effectiveness of the proposed method is demonstrated with the numerical simulations of the bulk metal forming processes including forging and extrusion.
基金Basic Research Program of Ministry of Industry and Information Technology of China(Grant No.B2420133003)National Natural Science Foundation of China(Grant Nos.51779060,51679054)
文摘The trajectory tracking control problem is addressed for autonomous underwater vehicle(AUV) in marine environ?ment, with presence of the influence of the uncertain factors including ocean current disturbance, dynamic modeling uncertainty, and thrust model errors. To improve the trajectory tracking accuracy of AUV, an adaptive backstepping terminal sliding mode control based on recurrent neural networks(RNN) is proposed. Firstly, considering the inaccu?rate of thrust model of thruster, a Taylor’s polynomial is used to obtain the thrust model errors. And then, the dynamic modeling uncertainty and thrust model errors are combined into the system model uncertainty(SMU) of AUV; through the RNN, the SMU and ocean current disturbance are classified, approximated online. Finally, the weights of RNN and other control parameters are adjusted online based on the backstepping terminal sliding mode controller. In addition, a chattering?reduction method is proposed based on sigmoid function. In chattering?reduction method, the sigmoid function is used to realize the continuity of the sliding mode switching function, and the sliding mode switching gain is adjusted online based on the exponential form of the sliding mode function. Based on the Lyapu?nov theory and Barbalat’s lemma, it is theoretically proved that the AUV trajectory tracking error can quickly converge to zero in the finite time. This research proposes a trajectory tracking control method of AUV, which can e ectively achieve high?precision trajectory tracking control of AUV under the influence of the uncertain factors. The feasibility and e ectiveness of the proposed method is demonstrated with trajectory tracking simulations and pool?experi?ments of AUV.
基金supported by National Natural Science Foundation of China (Grant Nos. 51205313, 50975232)Northwestern Polytechnical University Foundation for Fundamental Research of China (Grant No.JC20110249)
文摘There are two separate traditional approaches to model contact problems: continuum and atomistic theory. Continuum theory is successfully used in many domains, but when the scale of the model comes to nanometer, continuum approximation meets challenges. Atomistic theory can catch the detailed behaviors of an individual atom by using molecular dynamics (MD) or quantum mechanics, although accurately, it is usually time-consuming. A multiscale method coupled MD and finite element (FE) is presented. To mesh the FE region automatically, an adaptive method based on the strain energy gradient is introduced to the multiscale method to constitute an adaptive multiscale method. Utilizing the proposed method, adhesive contacts between a rigid cylinder and an elastic substrate are studied, and the results are compared with full MD simulations. The process of FE meshes refinement shows that adaptive multiscale method can make FE mesh generation more flexible. Comparison of the displacements of boundary atoms in the overlap region with the results from full MD simulations indicates that adaptive multiscale method can transfer displacements effectively. Displacements of atoms and FE nodes on the center line of the multiscale model agree well with that of atoms in full MD simulations, which shows the continuity in the overlap region. Furthermore, the Von Mises stress contours and contact force distributions in the contact region are almost same as full MD simulations. The method presented combines multiscale method and adaptive technique, and can provide a more effective way to multiscale method and to the investigation on nanoscale contact problems.
文摘In this paper, the generalized nonlinear Schrodinger equation （GNLSE） is solved by an adaptive split-step Fourier method （ASSFM）. It is found that ASSFM must be used to solve GNLSE to ensure precision when the soliton selffrequency shift is remarkable and the photonic crystal fibre （PCF） parameters vary with the frequency considerably. The precision of numerical simulation by using ASSFM is higher than that by using split-step Fourier method in the process of laser pulse propagation in PCFs due to the fact that the variation of fibre parameters with the peak frequency in the pulse spectrum can be taken into account fully.
