For the second-order finite volume method,implicit schemes and reconstruction methods are two main algorithms which influence the robustness and efficiency of the numerical simulations of compressible turbulent flows....For the second-order finite volume method,implicit schemes and reconstruction methods are two main algorithms which influence the robustness and efficiency of the numerical simulations of compressible turbulent flows.In this paper,a compact least-squares reconstruction method is proposed to calculate the gradients for the distribution of flow field variables approximation.The compactness of the new reconstruction method is reflected in the gradient calculation process.The geometries of the face-neighboring elements are no longer utilized,and the weighted average values at the centroid of the interfaces are used to calculate the gradients instead of the values at the centroid of the face-neighboring elements.Meanwhile,an exact Jacobian solving strategy is developed for implicit temporal discretization.The accurate processing of Jacobian matrix can extensively improve the invertibility of the Jacobian matrix and avoid introducing extra numerical errors.In addition,a modified Venkatakrishnan limiter is applied to deal with the shock which may appear in transonic flows and the applicability of the mentioned methods is enhanced further.The combination of the proposed methods makes the numerical simulations of turbulent flow converge rapidly and steadily with an adaptive increasing CFL number.The numerical results of several benchmarks indicate that the proposed methods perform well in terms of robustness,efficiency and accuracy,and have good application potential in turbulent flow simulations of complex configurations.展开更多
The main purpose of reverse engineering is to convert discrete data pointsinto piecewise smooth, continuous surface models. Before carrying out model reconstruction it issignificant to extract geometric features becau...The main purpose of reverse engineering is to convert discrete data pointsinto piecewise smooth, continuous surface models. Before carrying out model reconstruction it issignificant to extract geometric features because the quality of modeling greatly depends on therepresentation of features. Some fitting techniques of natural quadric surfaces with least-squaresmethod are described. And these techniques can be directly used to extract quadric surfaces featuresduring the process of segmentation for point cloud.展开更多
The Galerkin and least-squares methods are two classes of the most popular Krylov subspace methOds for solving large linear systems of equations. Unfortunately, both the methods may suffer from serious breakdowns of t...The Galerkin and least-squares methods are two classes of the most popular Krylov subspace methOds for solving large linear systems of equations. Unfortunately, both the methods may suffer from serious breakdowns of the same type: In a breakdown situation the Galerkin method is unable to calculate an approximate solution, while the least-squares method, although does not really break down, is unsucessful in reducing the norm of its residual. In this paper we first establish a unified theorem which gives a relationship between breakdowns in the two methods. We further illustrate theoretically and experimentally that if the coefficient matrix of a lienar system is of high defectiveness with the associated eigenvalues less than 1, then the restarted Galerkin and least-squares methods will be in great risks of complete breakdowns. It appears that our findings may help to understand phenomena observed practically and to derive treatments for breakdowns of this type.展开更多
A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercriti...A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.展开更多
Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of s...Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.展开更多
The purpose of this article is to develop and analyze least-squares approximations for the incompressible magnetohydrodynamic equations. The major advantage of the least-squares finite element method is that it is not...The purpose of this article is to develop and analyze least-squares approximations for the incompressible magnetohydrodynamic equations. The major advantage of the least-squares finite element method is that it is not subjected to the so-called Ladyzhenskaya-Babuska-Brezzi （LBB） condition. The authors employ least-squares functionals which involve a discrete inner product which is related to the inner product in H^-1（Ω）.展开更多
With the development of computational power, there has been an increased focus on data-fitting related seismic inversion techniques for high fidelity seismic velocity model and image, such as full-waveform inversion a...With the development of computational power, there has been an increased focus on data-fitting related seismic inversion techniques for high fidelity seismic velocity model and image, such as full-waveform inversion and least squares migration. However, though more advanced than conventional methods, these data fitting methods can be very expensive in terms of computational cost. Recently, various techniques to optimize these data-fitting seismic inversion problems have been implemented to cater for the industrial need for much improved efficiency. In this study, we propose a general stochastic conjugate gradient method for these data-fitting related inverse problems. We first prescribe the basic theory of our method and then give synthetic examples. Our numerical experiments illustrate the potential of this method for large-size seismic inversion application.展开更多
Elevation change monitoring of the Antarctic ice sheet has been a key issue in global change research.Satellite altimetry has been proven to be effective in detecting ice sheet variations. With the development of ICES...Elevation change monitoring of the Antarctic ice sheet has been a key issue in global change research.Satellite altimetry has been proven to be effective in detecting ice sheet variations. With the development of ICESat-2, many elevation observations can be used to derive elevation changes. However, the large amount of multitemporal data may include anomalous data points, increasing the uncertainty of the results. In this work, we improved the traditional repeat track method by introducing the Institute of Geodesy and Geophysics Ⅲ(IGGⅢ) method to obtain high-accuracy estimates of elevation change. The improved method was applied to analyze elevation changes along the transect from Zhongshan Station to Dome A in East Antarctica via ICESat-2 satellite altimetry data. The results show that the improved and traditional methods yield consistent numerical and spatial elevation change distributions. The elevation change calculated via the traditional method is 0.033 ± 0.131 m/yr, whereas the elevation change estimated via the IGGⅢ robust estimation method is 0.033 ± 0.109 m/yr from March 2019 to December 2021.In terms of spatial distribution, elevation changes in inland areas remain close to equilibrium, whereas regions with steeper ice sheet margins exhibit positive accumulation trends in elevation changes. The improved method reduces the standard error of the adjustment function from 0.975 to 0.691 m/yr. The improvement is particularly remarkable in the area between 72°S and 77°S. The results demonstrate that the IGGⅢ method effectively reduces errors caused by the inclusion of anomalous data and maintains the high data utilization rate of repeat-orbit methods.展开更多
Physics-informed neural networks(PINNs)have shown considerable promise for performing numerical simulations in fluid mechanics.They provide mesh-free,end-to-end approaches by embedding physical laws into their loss fu...Physics-informed neural networks(PINNs)have shown considerable promise for performing numerical simulations in fluid mechanics.They provide mesh-free,end-to-end approaches by embedding physical laws into their loss functions.However,when addressing complex flow problems,PINNs still face some challenges such as activation saturation and vanishing gradients in deep network training,leading to slow convergence and insufficient prediction accuracy.We present physics-informed neural networks incorporating lattice Boltzmann method optimized by tanh robust weight initialization(T-PINN-LBM)to address these challenges.This approach fuses the mesoscopic lattice Boltzmann model with the automatic differentiation framework of PINNs.It also implements a tanh robust weight initialization method derived from fixed point analysis.This model effectively mitigates activation and gradient decay in deep networks,improving convergence speed and data efficiency in multiscale flow simulations.We validate the effectiveness of the model on the classical arithmetic example of lid-driven cavity flow.Compared to the traditional Xavier initialized PINN and PINN-LBM,T-PINNLBM reduces the mean absolute error(MAE)by one order of magnitude at the same network depth and maintains stable convergence in deeper networks.The results demonstrate that this model can accurately capture complex flow structures without prior data,providing a new feasible pathway for data-free driven fluid simulation.展开更多
A least-squares mixed finite element (LSMFE) method for the numerical solution of fourth order parabolic problems analyzed and developed in this paper. The Ciarlet-Raviart mixed finite element space is used to approxi...A least-squares mixed finite element (LSMFE) method for the numerical solution of fourth order parabolic problems analyzed and developed in this paper. The Ciarlet-Raviart mixed finite element space is used to approximate. The a posteriori error estimator which is needed in the adaptive refinement algorithm is proposed. The local evaluation of the least-squares functional serves as a posteriori error estimator. The posteriori errors are effectively estimated. The convergence of the adaptive least-squares mixed finite element method is proved.展开更多
In this paper,we present a novel initial costates solver for initializing time-optimal trajectory problems in relative motion with continuous low thrust.The proposed solver consists of two primary components:training ...In this paper,we present a novel initial costates solver for initializing time-optimal trajectory problems in relative motion with continuous low thrust.The proposed solver consists of two primary components:training a Multilayer Perceptron(MLP)for generating reference sequence and Time of Flight(TOF)to the target,and deriving a system of linear algebraic equations for obtaining the initial costates.To overcome the challenge of generating training samples for the MLP,the backward generation method is proposed to obtain five different training databases.The training database and sample form are determined by analyzing the input and output correlation using the Pearson correlation coefficient.The best-performing MLP is obtained by analyzing the training results with various hyper-parameter combinations.A reference sequence starting from the initial states is obtained by integrating forward with the near-optimal control vector from the output of MLP.Finally,a system of linear algebraic equations for estimating the initial costates is derived using the reference sequence and the necessary conditions for optimality.Simulation results demonstrate that the proposed initial costates solver improves the convergence ratio and reduce the function calls of the shooting function.Furthermore,Monte-Carlo simulation illustrates that the initial costates solver is applicable to different initial velocities,demonstrating excellent generalization ability.展开更多
A least-squares mixed finite element method was formulated for a class of Stokes equations in two dimensional domains. The steady state and the time-dependent Stokes' equations were considered. For the stationary ...A least-squares mixed finite element method was formulated for a class of Stokes equations in two dimensional domains. The steady state and the time-dependent Stokes' equations were considered. For the stationary equation, optimal H-t and L-2-error estimates are derived under the standard regularity assumption on the finite element partition ( the LBB-condition is not required). Far the evolutionary equation, optimal L-2 estimates are derived under the conventional Raviart-Thomas spaces.展开更多
In this paper, a least-squares finite element method for the upper-convected Maxell (UCM) fluid is proposed. We first linearize the constitutive and momentum equations and then apply a least-squares method to the line...In this paper, a least-squares finite element method for the upper-convected Maxell (UCM) fluid is proposed. We first linearize the constitutive and momentum equations and then apply a least-squares method to the linearized version of the viscoelastic UCM model. The L2 least-squares functional involves the residuals of each equation multiplied by proper weights. The corresponding homogeneous functional is equivalent to a natural norm. The error estimates of the finite element solution are analyzed when the conforming piecewise polynomial elements are used for the unknowns.展开更多
Background:This article investigates the Least-Squares Monte Carlo Method by using different polynomial basis in American Asian Options pricing.The standard approach in the option pricing literature is to choose the b...Background:This article investigates the Least-Squares Monte Carlo Method by using different polynomial basis in American Asian Options pricing.The standard approach in the option pricing literature is to choose the basis arbitrarily.By comparing four different polynomial basis we show that the choice of basis interferes in the option's price.Methods:We assess Least-Squares Method performance in pricing four different American Asian Options by using four polynomial basis:Power,Laguerre,Legendre and Hermite A.To every American Asian Option priced,three sets of parameters are used in order to evaluate it properly.Results:We show that the choice of the basis interferes in the option's price by showing that one of them converges to the option's value faster than any other by using fewer simulated paths.In the case of an Amerasian call option,for example,we find that the preferable polynomial basis is Hermite A.For an Amerasian put option,the Power polynomial basis is recommended.Such empirical outcome is theoretically unpredictable,since in principle all basis can be indistinctly used when pricing the derivative.Conclusion:In this article The Least-Squares Monte Carlo Method performance is assessed in pricing four different types of American Asian Options by using four different polynomial basis through three different sets of parameters.Our results suggest that one polynomial basis is best suited to perform the method when pricing an American Asian option.Theoretically all basis can be indistinctly used when pricing the derivative.However,our results does not confirm these.We find that when pricing an American Asian put option,Power A is better than the other basis we have studied here whereas when pricing an American Asian call,Hermite A is better.展开更多
The Nelder-Mead simplex method is a well-known algorithm enabling the minimization of functions that are not available in closed-form and that need not be differentiable or convex.Furthermore,it is particularly parsim...The Nelder-Mead simplex method is a well-known algorithm enabling the minimization of functions that are not available in closed-form and that need not be differentiable or convex.Furthermore,it is particularly parsimonious on the number of function evaluations,thus making it preferable to convex optimization paradigms in the case,common when dealing with control design problems,that the objective function of the optimization problem is non-differentiable,non-convex,and its closed-form is not available or difficult to be computed analytically.The main goal of this paper is to show how the joint use of the Nelder-Mead simplex method and the Morrison algorithm can be successfully used to solve relevant and challenging control problems that cannot be easily solved using analytic methods.In particular,it is shown how the problems of strong stabilization,static output feedback stabilization,and design of robust controllers having fixed structure can be framed as optimization problems,which,in turn,can be efficiently solved by coupling the two above mentioned algorithms.The performance of this procedure is compared with state-of-the-art techniques on dozens of static output feedback benchmark case studies,and its effectiveness is demonstrated by several examples.展开更多
The full-spectrum least-squares(FSLS) method is introduced to perform quantitative energy-dispersive X-ray fluorescence analysis for unknown solid samples.Based on the conventional least-squares principle, this spectr...The full-spectrum least-squares(FSLS) method is introduced to perform quantitative energy-dispersive X-ray fluorescence analysis for unknown solid samples.Based on the conventional least-squares principle, this spectrum evaluation method is able to obtain the background-corrected and interference-free net peaks, which is significant for quantization analyses. A variety of analytical parameters and functions to describe the features of the fluorescence spectra of pure elements are used and established, such as the mass absorption coefficient, the Gi factor, and fundamental fluorescence formulas. The FSLS iterative program was compiled in the C language. The content of each component should reach the convergence criterion at the end of the calculations. After a basic theory analysis and experimental preparation, 13 national standard soil samples were detected using a spectrometer to test the feasibility of using the algorithm. The results show that the calculated contents of Ti, Fe, Ni, Cu, and Zn have the same changing tendency as the corresponding standard content in the 13 reference samples. Accuracies of 0.35% and 14.03% are obtained, respectively, for Fe and Ti, whose standard concentrations are 8.82% and 0.578%, respectively. However, the calculated results of trace elements (only tens of lg/g) deviate from the standard values. This may be because of measurement accuracy and mutual effects between the elements.展开更多
A linear-correction least-squares（LCLS） estimation procedure is proposed for geolocation using frequency difference of arrival （FDOA） measurements only. We first analyze the measurements of FDOA, and further deriv...A linear-correction least-squares（LCLS） estimation procedure is proposed for geolocation using frequency difference of arrival （FDOA） measurements only. We first analyze the measurements of FDOA, and further derive the Cramer-Rao lower bound （CRLB） of geoloeation using FDOA measurements. For the localization model is a nonlinear least squares（LS） estimator with a nonlinear constrained, a linearizing method is used to convert the model to a linear least squares estimator with a nonlinear con- strained. The Gauss-Newton iteration method is developed to conquer the source localization problem. From the analysis of solving Lagrange multiplier, the algorithm is a generalization of linear-correction least squares estimation procedure under the condition of geolocation using FDOA measurements only. The algorithm is compared with common least squares estimation. Comparisons of their estimation accuracy and the CRLB are made, and the proposed method attains the CRLB. Simulation re- sults are included to corroborate the theoretical development.展开更多
The robust magnesium surfaces with multi-functions are highly desirable,and the simple and scalable methods to construct such surfaces are urgently indispensable.Herein,we conducted a one-step spraying method to facil...The robust magnesium surfaces with multi-functions are highly desirable,and the simple and scalable methods to construct such surfaces are urgently indispensable.Herein,we conducted a one-step spraying method to facilely fabricate the robust coating with multi-functions on magnesium alloys.The as-sprayed magnesium alloys surface is superhydrophobic with a static water contact angle(WCA)of 157.0°and a roll-off angle of 6.0°.Such surface has excellent mechanical,chemical and thermal stabilities,even undergoing various physical and chemical damages,including sand impact(10 gmin^(-1),≥20 min),water impact(2 impacts s^(-1),≥180 min),abrasion(1.00 kPa,≥25 cycles),peeling(≥2.15 kPa),high temperature(200°C,≥24 h),strong acidic/salty/basic media(p H=113)and organic-solvent immersion(ethanol and n-hexane,≥24 h),demonstrating brilliant robustness.Notably,the surface displays multi-functions of corrosion protection,anti-fouling and heat insulation,which will undoubtedly promote the much wider applications of magnesium alloys.展开更多
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.展开更多
The design optimization taking into account the impact of uncertainties favors improving the robustness of the design.A Surrogate-Assisted Gradient-Based(SAGB)method for the robust aerodynamic design optimization of t...The design optimization taking into account the impact of uncertainties favors improving the robustness of the design.A Surrogate-Assisted Gradient-Based(SAGB)method for the robust aerodynamic design optimization of turbomachinery blades considering large-scale uncertainty is introduced,verified and validated in the study.The gradient-based method is employed due to its high optimization efficiency and any one surrogate model with sufficient response accuracy can be employed to quantify the nonlinear performance changes.The gradients of objective performance function to the design parameters are calculated first for all the training samples,from which the gradients of cost function can be fast determined.To reveal the high efficiency and high accuracy of SAGB on gradient calculation,the number of flow computations needed is evaluated and compared with three other methods.Through the aerodynamic design optimization of a transonic turbine cascade minimizing total pressure loss at the outlet,the SAGB-based gradients of the base and optimized blades are compared with those obtained by the Monte Carlo-assisted finite difference method.Moreover,the results of both the robust and deterministic aerodynamic design optimizations are presented and compared to demonstrate the practicability of SAGB on improving the aerodynamic robustness of turbomachinery blades.展开更多
基金supported by the National Natural Science Foundation of China(Nos.11702329,12102247)the Shanghai Key Lab of Vehicle Aerodynamics and Vehicle Thermal Management Systems,China(No.VATLAB-2021-01)。
文摘For the second-order finite volume method,implicit schemes and reconstruction methods are two main algorithms which influence the robustness and efficiency of the numerical simulations of compressible turbulent flows.In this paper,a compact least-squares reconstruction method is proposed to calculate the gradients for the distribution of flow field variables approximation.The compactness of the new reconstruction method is reflected in the gradient calculation process.The geometries of the face-neighboring elements are no longer utilized,and the weighted average values at the centroid of the interfaces are used to calculate the gradients instead of the values at the centroid of the face-neighboring elements.Meanwhile,an exact Jacobian solving strategy is developed for implicit temporal discretization.The accurate processing of Jacobian matrix can extensively improve the invertibility of the Jacobian matrix and avoid introducing extra numerical errors.In addition,a modified Venkatakrishnan limiter is applied to deal with the shock which may appear in transonic flows and the applicability of the mentioned methods is enhanced further.The combination of the proposed methods makes the numerical simulations of turbulent flow converge rapidly and steadily with an adaptive increasing CFL number.The numerical results of several benchmarks indicate that the proposed methods perform well in terms of robustness,efficiency and accuracy,and have good application potential in turbulent flow simulations of complex configurations.
基金This project is supported by Research Foundation for Doctoral Program of Higher Education, China (No.98033532)
文摘The main purpose of reverse engineering is to convert discrete data pointsinto piecewise smooth, continuous surface models. Before carrying out model reconstruction it issignificant to extract geometric features because the quality of modeling greatly depends on therepresentation of features. Some fitting techniques of natural quadric surfaces with least-squaresmethod are described. And these techniques can be directly used to extract quadric surfaces featuresduring the process of segmentation for point cloud.
文摘The Galerkin and least-squares methods are two classes of the most popular Krylov subspace methOds for solving large linear systems of equations. Unfortunately, both the methods may suffer from serious breakdowns of the same type: In a breakdown situation the Galerkin method is unable to calculate an approximate solution, while the least-squares method, although does not really break down, is unsucessful in reducing the norm of its residual. In this paper we first establish a unified theorem which gives a relationship between breakdowns in the two methods. We further illustrate theoretically and experimentally that if the coefficient matrix of a lienar system is of high defectiveness with the associated eigenvalues less than 1, then the restarted Galerkin and least-squares methods will be in great risks of complete breakdowns. It appears that our findings may help to understand phenomena observed practically and to derive treatments for breakdowns of this type.
基金the National Science Council ot Taiwan,China for funding this research(Project no.:NSC 94-2218-E-035-011)
文摘A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.
基金the National Science Council of Taiwan for funding this research (NSC 96-2221-E-019-061).
文摘Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.
基金supported by the National Basic Research Program of China (2005CB321701)NSF of mathematics research special fund of Hebei Province(08M005)
文摘The purpose of this article is to develop and analyze least-squares approximations for the incompressible magnetohydrodynamic equations. The major advantage of the least-squares finite element method is that it is not subjected to the so-called Ladyzhenskaya-Babuska-Brezzi （LBB） condition. The authors employ least-squares functionals which involve a discrete inner product which is related to the inner product in H^-1（Ω）.
基金partially supported by the National Natural Science Foundation of China (No.41230318)
文摘With the development of computational power, there has been an increased focus on data-fitting related seismic inversion techniques for high fidelity seismic velocity model and image, such as full-waveform inversion and least squares migration. However, though more advanced than conventional methods, these data fitting methods can be very expensive in terms of computational cost. Recently, various techniques to optimize these data-fitting seismic inversion problems have been implemented to cater for the industrial need for much improved efficiency. In this study, we propose a general stochastic conjugate gradient method for these data-fitting related inverse problems. We first prescribe the basic theory of our method and then give synthetic examples. Our numerical experiments illustrate the potential of this method for large-size seismic inversion application.
基金supported by the National Key Research and Development Program of China under grant number 2023YFC2809103the Fundamental Research Funds for the Central Universities under grant numbers 2042022kf1204, 2042022kf1069, 2042023gf0012, 2042022dx0001+1 种基金the Hubei Provincial Natural Science Foundation of China under grant number 2022CFB081the State Key Laboratory of Geodesy and Earth's Dynamics, Innovation Academy for Precision Measurement Science and Technology under grant number SKLGED2023-2-6
文摘Elevation change monitoring of the Antarctic ice sheet has been a key issue in global change research.Satellite altimetry has been proven to be effective in detecting ice sheet variations. With the development of ICESat-2, many elevation observations can be used to derive elevation changes. However, the large amount of multitemporal data may include anomalous data points, increasing the uncertainty of the results. In this work, we improved the traditional repeat track method by introducing the Institute of Geodesy and Geophysics Ⅲ(IGGⅢ) method to obtain high-accuracy estimates of elevation change. The improved method was applied to analyze elevation changes along the transect from Zhongshan Station to Dome A in East Antarctica via ICESat-2 satellite altimetry data. The results show that the improved and traditional methods yield consistent numerical and spatial elevation change distributions. The elevation change calculated via the traditional method is 0.033 ± 0.131 m/yr, whereas the elevation change estimated via the IGGⅢ robust estimation method is 0.033 ± 0.109 m/yr from March 2019 to December 2021.In terms of spatial distribution, elevation changes in inland areas remain close to equilibrium, whereas regions with steeper ice sheet margins exhibit positive accumulation trends in elevation changes. The improved method reduces the standard error of the adjustment function from 0.975 to 0.691 m/yr. The improvement is particularly remarkable in the area between 72°S and 77°S. The results demonstrate that the IGGⅢ method effectively reduces errors caused by the inclusion of anomalous data and maintains the high data utilization rate of repeat-orbit methods.
基金supported by the National Natural Science Foundation of China(Grant Nos.12474453,12174085,and 12404530).
文摘Physics-informed neural networks(PINNs)have shown considerable promise for performing numerical simulations in fluid mechanics.They provide mesh-free,end-to-end approaches by embedding physical laws into their loss functions.However,when addressing complex flow problems,PINNs still face some challenges such as activation saturation and vanishing gradients in deep network training,leading to slow convergence and insufficient prediction accuracy.We present physics-informed neural networks incorporating lattice Boltzmann method optimized by tanh robust weight initialization(T-PINN-LBM)to address these challenges.This approach fuses the mesoscopic lattice Boltzmann model with the automatic differentiation framework of PINNs.It also implements a tanh robust weight initialization method derived from fixed point analysis.This model effectively mitigates activation and gradient decay in deep networks,improving convergence speed and data efficiency in multiscale flow simulations.We validate the effectiveness of the model on the classical arithmetic example of lid-driven cavity flow.Compared to the traditional Xavier initialized PINN and PINN-LBM,T-PINNLBM reduces the mean absolute error(MAE)by one order of magnitude at the same network depth and maintains stable convergence in deeper networks.The results demonstrate that this model can accurately capture complex flow structures without prior data,providing a new feasible pathway for data-free driven fluid simulation.
文摘A least-squares mixed finite element (LSMFE) method for the numerical solution of fourth order parabolic problems analyzed and developed in this paper. The Ciarlet-Raviart mixed finite element space is used to approximate. The a posteriori error estimator which is needed in the adaptive refinement algorithm is proposed. The local evaluation of the least-squares functional serves as a posteriori error estimator. The posteriori errors are effectively estimated. The convergence of the adaptive least-squares mixed finite element method is proved.
基金This study was funded by the National Natural Science Foundation of China(Nos.11972077 and 12272039).
文摘In this paper,we present a novel initial costates solver for initializing time-optimal trajectory problems in relative motion with continuous low thrust.The proposed solver consists of two primary components:training a Multilayer Perceptron(MLP)for generating reference sequence and Time of Flight(TOF)to the target,and deriving a system of linear algebraic equations for obtaining the initial costates.To overcome the challenge of generating training samples for the MLP,the backward generation method is proposed to obtain five different training databases.The training database and sample form are determined by analyzing the input and output correlation using the Pearson correlation coefficient.The best-performing MLP is obtained by analyzing the training results with various hyper-parameter combinations.A reference sequence starting from the initial states is obtained by integrating forward with the near-optimal control vector from the output of MLP.Finally,a system of linear algebraic equations for estimating the initial costates is derived using the reference sequence and the necessary conditions for optimality.Simulation results demonstrate that the proposed initial costates solver improves the convergence ratio and reduce the function calls of the shooting function.Furthermore,Monte-Carlo simulation illustrates that the initial costates solver is applicable to different initial velocities,demonstrating excellent generalization ability.
文摘A least-squares mixed finite element method was formulated for a class of Stokes equations in two dimensional domains. The steady state and the time-dependent Stokes' equations were considered. For the stationary equation, optimal H-t and L-2-error estimates are derived under the standard regularity assumption on the finite element partition ( the LBB-condition is not required). Far the evolutionary equation, optimal L-2 estimates are derived under the conventional Raviart-Thomas spaces.
文摘In this paper, a least-squares finite element method for the upper-convected Maxell (UCM) fluid is proposed. We first linearize the constitutive and momentum equations and then apply a least-squares method to the linearized version of the viscoelastic UCM model. The L2 least-squares functional involves the residuals of each equation multiplied by proper weights. The corresponding homogeneous functional is equivalent to a natural norm. The error estimates of the finite element solution are analyzed when the conforming piecewise polynomial elements are used for the unknowns.
文摘Background:This article investigates the Least-Squares Monte Carlo Method by using different polynomial basis in American Asian Options pricing.The standard approach in the option pricing literature is to choose the basis arbitrarily.By comparing four different polynomial basis we show that the choice of basis interferes in the option's price.Methods:We assess Least-Squares Method performance in pricing four different American Asian Options by using four polynomial basis:Power,Laguerre,Legendre and Hermite A.To every American Asian Option priced,three sets of parameters are used in order to evaluate it properly.Results:We show that the choice of the basis interferes in the option's price by showing that one of them converges to the option's value faster than any other by using fewer simulated paths.In the case of an Amerasian call option,for example,we find that the preferable polynomial basis is Hermite A.For an Amerasian put option,the Power polynomial basis is recommended.Such empirical outcome is theoretically unpredictable,since in principle all basis can be indistinctly used when pricing the derivative.Conclusion:In this article The Least-Squares Monte Carlo Method performance is assessed in pricing four different types of American Asian Options by using four different polynomial basis through three different sets of parameters.Our results suggest that one polynomial basis is best suited to perform the method when pricing an American Asian option.Theoretically all basis can be indistinctly used when pricing the derivative.However,our results does not confirm these.We find that when pricing an American Asian put option,Power A is better than the other basis we have studied here whereas when pricing an American Asian call,Hermite A is better.
基金partially supported by the Italian Ministry for Research in the framework of the 2020 Program for Research Projects of National Interest(2020RTWES4)。
文摘The Nelder-Mead simplex method is a well-known algorithm enabling the minimization of functions that are not available in closed-form and that need not be differentiable or convex.Furthermore,it is particularly parsimonious on the number of function evaluations,thus making it preferable to convex optimization paradigms in the case,common when dealing with control design problems,that the objective function of the optimization problem is non-differentiable,non-convex,and its closed-form is not available or difficult to be computed analytically.The main goal of this paper is to show how the joint use of the Nelder-Mead simplex method and the Morrison algorithm can be successfully used to solve relevant and challenging control problems that cannot be easily solved using analytic methods.In particular,it is shown how the problems of strong stabilization,static output feedback stabilization,and design of robust controllers having fixed structure can be framed as optimization problems,which,in turn,can be efficiently solved by coupling the two above mentioned algorithms.The performance of this procedure is compared with state-of-the-art techniques on dozens of static output feedback benchmark case studies,and its effectiveness is demonstrated by several examples.
基金supported by the National Key R&D Project of China(No.2017YFC0602100)the National Natural Science Foundation of China(No.41774147)Sichuan Science and Technology Support Program(No.2015GZ0272)
文摘The full-spectrum least-squares(FSLS) method is introduced to perform quantitative energy-dispersive X-ray fluorescence analysis for unknown solid samples.Based on the conventional least-squares principle, this spectrum evaluation method is able to obtain the background-corrected and interference-free net peaks, which is significant for quantization analyses. A variety of analytical parameters and functions to describe the features of the fluorescence spectra of pure elements are used and established, such as the mass absorption coefficient, the Gi factor, and fundamental fluorescence formulas. The FSLS iterative program was compiled in the C language. The content of each component should reach the convergence criterion at the end of the calculations. After a basic theory analysis and experimental preparation, 13 national standard soil samples were detected using a spectrometer to test the feasibility of using the algorithm. The results show that the calculated contents of Ti, Fe, Ni, Cu, and Zn have the same changing tendency as the corresponding standard content in the 13 reference samples. Accuracies of 0.35% and 14.03% are obtained, respectively, for Fe and Ti, whose standard concentrations are 8.82% and 0.578%, respectively. However, the calculated results of trace elements (only tens of lg/g) deviate from the standard values. This may be because of measurement accuracy and mutual effects between the elements.
基金National High-tech Research and Development Program of China (2011AA7072043)National Defense Key Laboratory Foundation of China (9140C860304)Innovation Fund of Graduate School of NUDT (B120406)
文摘A linear-correction least-squares（LCLS） estimation procedure is proposed for geolocation using frequency difference of arrival （FDOA） measurements only. We first analyze the measurements of FDOA, and further derive the Cramer-Rao lower bound （CRLB） of geoloeation using FDOA measurements. For the localization model is a nonlinear least squares（LS） estimator with a nonlinear constrained, a linearizing method is used to convert the model to a linear least squares estimator with a nonlinear con- strained. The Gauss-Newton iteration method is developed to conquer the source localization problem. From the analysis of solving Lagrange multiplier, the algorithm is a generalization of linear-correction least squares estimation procedure under the condition of geolocation using FDOA measurements only. The algorithm is compared with common least squares estimation. Comparisons of their estimation accuracy and the CRLB are made, and the proposed method attains the CRLB. Simulation re- sults are included to corroborate the theoretical development.
基金supported by the National Natural Science Foundation of China(21773019,21972012)the Graduate Research and Innovation Foundation of Chongqing(CYB18044)the sharing fund of Chongqing University s Large-scale Equipment
文摘The robust magnesium surfaces with multi-functions are highly desirable,and the simple and scalable methods to construct such surfaces are urgently indispensable.Herein,we conducted a one-step spraying method to facilely fabricate the robust coating with multi-functions on magnesium alloys.The as-sprayed magnesium alloys surface is superhydrophobic with a static water contact angle(WCA)of 157.0°and a roll-off angle of 6.0°.Such surface has excellent mechanical,chemical and thermal stabilities,even undergoing various physical and chemical damages,including sand impact(10 gmin^(-1),≥20 min),water impact(2 impacts s^(-1),≥180 min),abrasion(1.00 kPa,≥25 cycles),peeling(≥2.15 kPa),high temperature(200°C,≥24 h),strong acidic/salty/basic media(p H=113)and organic-solvent immersion(ethanol and n-hexane,≥24 h),demonstrating brilliant robustness.Notably,the surface displays multi-functions of corrosion protection,anti-fouling and heat insulation,which will undoubtedly promote the much wider applications of magnesium alloys.
基金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.
基金National Natural Science Foundation of China(Nos.51676003,51976183)National Science and Technology Major Project of China(No.J2019II-0012-0032)。
文摘The design optimization taking into account the impact of uncertainties favors improving the robustness of the design.A Surrogate-Assisted Gradient-Based(SAGB)method for the robust aerodynamic design optimization of turbomachinery blades considering large-scale uncertainty is introduced,verified and validated in the study.The gradient-based method is employed due to its high optimization efficiency and any one surrogate model with sufficient response accuracy can be employed to quantify the nonlinear performance changes.The gradients of objective performance function to the design parameters are calculated first for all the training samples,from which the gradients of cost function can be fast determined.To reveal the high efficiency and high accuracy of SAGB on gradient calculation,the number of flow computations needed is evaluated and compared with three other methods.Through the aerodynamic design optimization of a transonic turbine cascade minimizing total pressure loss at the outlet,the SAGB-based gradients of the base and optimized blades are compared with those obtained by the Monte Carlo-assisted finite difference method.Moreover,the results of both the robust and deterministic aerodynamic design optimizations are presented and compared to demonstrate the practicability of SAGB on improving the aerodynamic robustness of turbomachinery blades.
