A novel gappy technology, gappy autoencoder with proper orthogonal decomposition(Gappy POD-AE), is proposed for reconstructing physical fields from sparse data. High-dimensional data are reduced via proper orthogonal ...A novel gappy technology, gappy autoencoder with proper orthogonal decomposition(Gappy POD-AE), is proposed for reconstructing physical fields from sparse data. High-dimensional data are reduced via proper orthogonal decomposition(POD),and low-dimensional data are used to train an autoencoder(AE). By integrating the POD operator with the decoder, a nonlinear solution form is established and incorporated into a new maximum-a-posteriori(MAP)-based objective for online reconstruction.The numerical results on the two-dimensional(2D) Bhatnagar-Gross-Krook-Boltzmann(BGK-Boltzmann) equation, wave equation, shallow-water equation, and satellite data show that Gappy POD-AE achieves higher accuracy than gappy proper orthogonal decomposition(Gappy POD), especially for the data with slowly decaying singular values,and is more efficient in training than gappy autoencoder(Gappy AE). The MAP-based formulation and new gappy procedure further enhance the reconstruction accuracy.展开更多
Proper Orthogonal Decomposition （POD） provides a powerful modal transformation tool for stochastic dynamics. In this paper, coherency matrix-based proper orthogonal decomposition （CPOD） is presented as an innovati...Proper Orthogonal Decomposition （POD） provides a powerful modal transformation tool for stochastic dynamics. In this paper, coherency matrix-based proper orthogonal decomposition （CPOD） is presented as an innovative form of the POD based on cross power spectral density matrices. By introducing a discretizing scheme, the CPOD-based spectral representation method is obtained for use in stochastic simulation. Moreover, some criteria are proposed that allow the truncation order of CPOD to be conveniently determined. A numerical example to illustrate the application of the proposed method for the simulation of a wind velocity field is provided.展开更多
The hydraulic fracturing is a nonlinear,fluid-solid coupling and transient problem,in most cases it is always time-consuming to simulate this process numerically.In recent years,although many numerical methods were pr...The hydraulic fracturing is a nonlinear,fluid-solid coupling and transient problem,in most cases it is always time-consuming to simulate this process numerically.In recent years,although many numerical methods were proposed to settle this problem,most of them still require a large amount of computer resources.Thus it is a high demand to develop more efficient numerical approaches to achieve the real-time monitoring of the fracture geometry during the hydraulic fracturing treatment.In this study,a reduced order modeling technique namely Proper Generalized Decomposition(PGD),is applied to accelerate the simulations of the transient,non-linear coupled system of hydraulic fracturing problem,to match this extremely tight response time constraint.The separability of the solution in space and time dimensions is studied for a simplified model problem.The solid and fluid equations are coupled explicitly by inverting the solid discrete problem,and a simple iterative procedure to handle the non-linear characteristic of the hydraulic fracturing problem is proposed in this work.Numeral validation illustrates that the results of PGD match well with these of standard finite element method in terms o f fracture opening and fluid pressure in the hydro-fracture.Moreover,after the off-line calculations,the numerical results can be obtained in real time.展开更多
Time-resolved particle image velocimetry（TRPIV） experiments are performed to investigate the coherent structure＇s performance of riblets in a turbulent boundary layer（TBL） at a friction Reynolds number of 185. To...Time-resolved particle image velocimetry（TRPIV） experiments are performed to investigate the coherent structure＇s performance of riblets in a turbulent boundary layer（TBL） at a friction Reynolds number of 185. To visualize the energetic large-scale coherent structures（CSs） over a smooth surface and riblets, the proper orthogonal decomposition（POD） and finite-time Lyapunov exponent（FTLE） are used to identify the CSs in the TBL. Spatial-temporal correlation is implemented to obtain the characters and transport properties of typical CSs in the FTLE fields. The results demonstrate that the generic flow structures, such as hairpin-like vortices, are also observed in the boundary layer flow over the riblets, consistent with its smooth counterpart. Low-order POD modes are more sensitive to the riblets in comparison with the high-order ones,and the wall-normal movement of the most energy-containing structures are suppressed over riblets. The spatial correlation analysis of the FTLE fields indicates that the evolution process of the hairpin vortex over riblets are inhibited. An apparent decrease of the convection velocity over riblets is noted, which is believed to reduce the ejection/sweep motions associated with high shear stress from the viscous sublayer. These reductions exhibit inhibition of momentum transfer among the structures near the wall in the TBL flows.展开更多
The proper orthogonal decomposition (POD) method for the instationary Navier-Stokes equations is considered. Several numerical approaches to evaluating the POD eigenfunctions are presented. The POD eigenfunctions are ...The proper orthogonal decomposition (POD) method for the instationary Navier-Stokes equations is considered. Several numerical approaches to evaluating the POD eigenfunctions are presented. The POD eigenfunctions are applied as a basis for a Galerkin projection of the instationary Navier-Stokes equations. And a low-dimensional ordinary differential models for fluid flows governed by the instationary Navier-Stokes equations are constructed. The numerical examples show that the method is feasible and efficient for optimal control of fluids.展开更多
The reduced-order model （ROM） for the two-dimensional supersonic cavity flow based on proper orthogonal decomposition （POD） and Galerkin projection is investigated. Presently, popular ROMs in cavity flows are base...The reduced-order model （ROM） for the two-dimensional supersonic cavity flow based on proper orthogonal decomposition （POD） and Galerkin projection is investigated. Presently, popular ROMs in cavity flows are based on an isentropic assumption, valid only for flows at low or moderate Mach numbers. A new ROM is constructed involving primitive variables of the fully compressible Navier-Stokes （N-S） equations, which is suitable for flows at high Mach numbers. Compared with the direct numerical simulation （DNS） results, the proposed model predicts flow dynamics （e.g., dominant frequency and amplitude） accurately for supersonic cavity flows, and is robust. The comparison between the present transient flow fields and those of the DNS shows that the proposed ROM can capture self-sustained oscillations of a shear layer. In addition, the present model reduction method can be easily extended to other supersonic flows.展开更多
A three-dimensional incompressible annular jet is simulated by the large eddy simulation(LES)method at a Reynolds number Re=8500.The time-averaged velocity field shows an asymmetric wake behind the central bluff-body ...A three-dimensional incompressible annular jet is simulated by the large eddy simulation(LES)method at a Reynolds number Re=8500.The time-averaged velocity field shows an asymmetric wake behind the central bluff-body although the flow geometry is symmetric.The proper orthogonal decomposition(POD)analysis of the velocity fluctuation vectors is conducted to study the flow dynamics of the wake flow.The distribution of turbulent kinetic energy across the three-dimensional POD modes shows that the first four eigenmodes each capture more than 1%of the turbulent kinetic energy,and hence their impact on the wake dynamics is studied.The results demonstrate that the asymmetric mean flow in the near-field of the annular jet is related to the first two POD modes which correspond to a radial shift of the stagnation point.The modes 3 and 4 involve the stretching or squeezing effects of the recirculation region in the radial direction.In addition,the spatial structure of these four POD eigenmodes also shows the counter-rotating vortices in the streamwise direction downstream of the flow reversal region.展开更多
Proper orthogonal decomposition （POD） is an effective statistical technique for data reduction and feature extraction of the random field including the wind field. This paper introduces the theory of the POD and ill...Proper orthogonal decomposition （POD） is an effective statistical technique for data reduction and feature extraction of the random field including the wind field. This paper introduces the theory of the POD and illustrates engineering of structures. Using the POD technique, it is shown that wind pressure data can be accurately reconstructed with a limited number of modes using the wind pressure data from wind tunnel test. Comparing the reconstructed values by POD with the original measured values from the wind tunnel test both in the time and frequency domains, it is concluded that the proper orthogonal decomposition（POD） is an efficient and practical technique for deriving the random wind pressure field from limited known data as shown in the pitched roof example in this paper.展开更多
The large storage requirement is a critical issue in cross-correlation imaging-condition based reverse time migration(RTM),because it requires the operation of the source and receiver wavefields at the same time.The b...The large storage requirement is a critical issue in cross-correlation imaging-condition based reverse time migration(RTM),because it requires the operation of the source and receiver wavefields at the same time.The boundary value method(BVM),based on the finite difference method(FDM),can be used to reconstruct the source wavefield in the reverse time propagation in the same way as the receiver wavefield,which can reduce the storage burden of the RTM data.Considering that the FDM cannot well handle models with discontinuous material properties and rough interfaces,we develop a source wavefield reconstruction strategy based on the finite element method(FEM),using proper orthogonal decomposition(POD)to enhance computational efficiency.In this method,we divide the whole time period into several segments,and construct the POD basis functions to get a reduced order model(ROM)for the source wavefield reconstruction in each segment.We show the corresponding quantitative analysis of the storage requirement of the POD-FEM.Numerical tests on the homogeneous model show the effectiveness of the proposed method,while the layered model and part of the Marmousi model tests indicate that the POD-FEM can keep an excellent balance between computational efficiency and memory usage compared with the full-stored method(FSM)and the BVM,and can be effectively applied in imaging.展开更多
Over the recent years there has been an increased trend in the use of Reduced Order Models （ROM） for modeling the coupled aeroelastic system. Of all the ROM models, the Proper Orthogonal Decomposition Method （POD）...Over the recent years there has been an increased trend in the use of Reduced Order Models （ROM） for modeling the coupled aeroelastic system. Of all the ROM models, the Proper Orthogonal Decomposition Method （POD） has been the most widely used, reason being the relative simplicity of implementation and the physical insight that it offers towards the physical problem. In this paper we begin by briefly recalling the recent work using POD for the computational aeroelasticity followed by the mathematical formulation. Mathematical formulation is important as it provides understanding of how POD method works. Implementation issues related to the POD method are presented next. Since POD is an empirical technique therefore, it is marred by the robustness issues as is the case with all the ROMs. In the end the variations of POD method, developed over the years are presented along with the most recent trend of using hybrid ROM.展开更多
The present work uses dynamic mode decomposition(DMD) to analyze wake flow of NACA0015 airfoil with Gurney flap.The physics of DMD is first introduced.Then the PIV-measured wake flow velocity field is decomposed into ...The present work uses dynamic mode decomposition(DMD) to analyze wake flow of NACA0015 airfoil with Gurney flap.The physics of DMD is first introduced.Then the PIV-measured wake flow velocity field is decomposed into dynamical modes.The vortex shedding pattern behind the trailing edge and its high-order harmonics have been captured with abundant information such as frequency,wavelength and convection speed.It is observed that high-order dynamic modes convect faster than low-order modes;moreover the wavelength of the dynamic modes scales with the corresponding frequency in power law.展开更多
Bends contribute to a flexible layout of pipeline system,but also lead to intensive energy costs due to the complex flow characteristic.This experimental study is conducted to investigate the impact of a single coarse...Bends contribute to a flexible layout of pipeline system,but also lead to intensive energy costs due to the complex flow characteristic.This experimental study is conducted to investigate the impact of a single coarse particle on the flow field in a bend.The velocity profiles of fluid on the axial symmetry plane of the bend are measured using time-resolved particle image velocimetry.The flow structures are extracted using the proper orthogonal decomposition method.The results reveal that there is a shear-layer flow in the bend during the transportation.With the increase in particle size,the particle has a dominant influence on the flow energy distribution of the overall flow.The impact of particles on the first few energetic flows is mainly in the latter part of the transportation,both temporally and spatially.As the particle size decreases,the most energetic unsteady flow within the bend changes from the convective flow to the flow of the shear layer.展开更多
This study experimentally investigates the hydrodynamic characteristics,geometric configurations,fluttering motions of the codend,and the instantaneous flow fields inside and around the codend,with and without a liner...This study experimentally investigates the hydrodynamic characteristics,geometric configurations,fluttering motions of the codend,and the instantaneous flow fields inside and around the codend,with and without a liner,under varying catch sizes and inflow velocities.A proper orthogonal decomposition method is employed to extract phase-averaged mean properties of unsteady turbulent flows from flow measurement data obtained using an electromagnetic current velocity meter inside and around the codend.The results reveal that as catch size increases,the drag force,codend motion,Reynolds number,and codend volume increase while the drag coefficient decreases.Owing to the codend shape and pronounced motion,a complex fluid–structure interaction occurs,demonstrating a strong correlation between drag force and codend volume.The oscillation amplitudes of the hydrodynamic forces and codend motions increase with increasing catch size,and their oscillations mainly involve low-frequency activity.A significant reduction in the flow field occurs inside and around the unlined codend without a catch.The flow field is 5.81%,14.39%,and 27.01%lower than the unlined codend with a catch,the codend with a liner but without a catch,and the codend with both a liner and a catch,respectively.Fourier analysis reveals that the codend motions and hydrodynamic forces are mainly characterized by low-frequency activity and are synchronized with the unsteady turbulent flow street.Furthermore,the proper orthogonal decomposition results reveal the development of unsteady turbulent flow inside and around the codend,driven by flow passage blockage caused by the presence of the liner,intense codend motions,and the catch.Understanding the hydrodynamic characteristics and flow instabilities inside and around the codend,particularly those associated with its fluttering motions,is crucial for optimizing trawl design and improving trawl selectivity.展开更多
Jacket platforms constitute the foundational infrastructure of offshore oil and gas field exploitation.How to efficiently and accurately monitor the mechanical properties of jacket structures is one of the key problem...Jacket platforms constitute the foundational infrastructure of offshore oil and gas field exploitation.How to efficiently and accurately monitor the mechanical properties of jacket structures is one of the key problems to be solved to ensure the safe operation of the platform.To address the practical engineering problem that it is difficult to monitor the stress response of the tubular joints of jacket platforms online,a digital twin reduced-order method for real-time prediction of the stress response of tubular joints is proposed.In the offline construction phase,multi-scale modeling and multi-parameter experimental design methods are used to obtain the stress response data set of the jacket structure.Proper orthogonal decomposition is employed to extract the main feature information from the snapshot matrix,resulting in a reduced-order basis.The leave-one-out cross-validation method is used to select the optimal modal order for constructing the reduced-order model(ROM).In the online prediction phase,a digital twin model of the tubular joint is established,and the prediction performance of the ROM is analyzed and verified through using random environmental load and field environmental monitoring data.The results indicate that,compared with traditional numerical simulations of tubular joints,the ROM based on the proposed reduced-order method is more efficient in predicting the stress response of tubular joints while ensuring accuracy and robustness.展开更多
Digital twin is considered the key technique for real-time monitoring and life-cycle management of electric equipment.To construct the digital twin model of electric equipment,a multi-parameter electromagnetic analysi...Digital twin is considered the key technique for real-time monitoring and life-cycle management of electric equipment.To construct the digital twin model of electric equipment,a multi-parameter electromagnetic analysis is needed to generate a large amount of high fidelity data under various working condition.However,repeated solving such multi-parameter electromagnetic problems based on full order finite element method may lead to extreme scale calculations.To address this issue,a hybrid approach that combines tensor decomposition and proper orthogonal decomposition(POD)is introduced,which can effectively establish a reduced order model for multi-parameter electromagnetic field problems.The performance of the proposed approach is illustrated through three numerical examples,namely an electrical motor,a transformer and a voice coil actuator including parameter variations of operating conditions,geometric parameters,and material parameters.The numerical results show that the proposed hybrid approach has significant advantage compared to conventional reduced order model in situation where the solution changes dramatically within the parameter variation range and even more apparent when the parameter dimension is high.展开更多
A proper orthogonal decomposition (POD) method is applied to a usual finite element (FE) formulation for parabolic equations so that it is reduced into a POD FE formulation with lower dimensions and enough high accura...A proper orthogonal decomposition (POD) method is applied to a usual finite element (FE) formulation for parabolic equations so that it is reduced into a POD FE formulation with lower dimensions and enough high accuracy. The errors between the reduced POD FE solution and the usual FE solution are analyzed. It is shown by numerical examples that the results of numerical computations are consistent with theoretical conclusions. Moreover, it is also shown that this validates the feasibility and efficiency of POD method.展开更多
The proper orthogonal decomposition(POD)and the singular value decomposition(SVD) are used to study the finite difference scheme(FDS)for the nonstationary Navier-Stokes equations. Ensembles of data are compiled from t...The proper orthogonal decomposition(POD)and the singular value decomposition(SVD) are used to study the finite difference scheme(FDS)for the nonstationary Navier-Stokes equations. Ensembles of data are compiled from the transient solutions computed from the discrete equation system derived from the FDS for the nonstationary Navier-Stokes equations.The optimal orthogonal bases are reconstructed by the elements of the ensemble with POD and SVD.Combining the above procedures with a Galerkin projection approach yields a new optimizing FDS model with lower dimensions and a high accuracy for the nonstationary Navier-Stokes equations.The errors between POD approximate solutions and FDS solutions are analyzed.It is shown by considering the results obtained for numerical simulations of cavity flows that the error between POD approximate solution and FDS solution is consistent with theoretical results.Moreover,it is also shown that this validates the feasibility and efficiency of POD method.展开更多
In aerodynamic optimization, global optimization methods such as genetic algorithms are preferred in many cases because of their advantage on reaching global optimum. However,for complex problems in which large number...In aerodynamic optimization, global optimization methods such as genetic algorithms are preferred in many cases because of their advantage on reaching global optimum. However,for complex problems in which large number of design variables are needed, the computational cost becomes prohibitive, and thus original global optimization strategies are required. To address this need, data dimensionality reduction method is combined with global optimization methods, thus forming a new global optimization system, aiming to improve the efficiency of conventional global optimization. The new optimization system involves applying Proper Orthogonal Decomposition（POD） in dimensionality reduction of design space while maintaining the generality of original design space. Besides, an acceleration approach for samples calculation in surrogate modeling is applied to reduce the computational time while providing sufficient accuracy. The optimizations of a transonic airfoil RAE2822 and the transonic wing ONERA M6 are performed to demonstrate the effectiveness of the proposed new optimization system. In both cases, we manage to reduce the number of design variables from 20 to 10 and from 42 to 20 respectively. The new design optimization system converges faster and it takes 1/3 of the total time of traditional optimization to converge to a better design, thus significantly reducing the overall optimization time and improving the efficiency of conventional global design optimization method.展开更多
Large eddy simulation(LES)is used to calculate the in-cylinder turbulent flow field in a direct injection spark ignition(DISI)engine.The computations are carried out for three different maximum valve lifts(MVL)and thr...Large eddy simulation(LES)is used to calculate the in-cylinder turbulent flow field in a direct injection spark ignition(DISI)engine.The computations are carried out for three different maximum valve lifts(MVL)and throughout 100 consecutive engine cycles.The simulated results as well as corresponding particle image velocimetry(PIV)measurement database are analyzed by the proper orthogonal decomposition(POD)method.Through a new developed POD quadruple decomposition the instantaneous in-cylinder flow fields are decomposed into four parts,named mean field,coherent field,transition field and turbulent field,respectively.Then the in-cylinder turbulent flow characteristics and cycle-to-cycle variations(CCV)are studied separately upon the four part flow fields.Results indicate that each part exhibits its specific characteristics and has close connection with others.The mean part contains more than 50%of the total kinetic energy and the energy cascade phenomenon occurs among the four part fields;the coherent field part possesses the highest CCV level which dominates CCV of the bulk flow.In addition,it is observed that a change in MVL affects significantly the in-cylinder flow behavior including CCV,especially for the coherent part.Furthermore,the POD analysis demonstrates that at least 25 sample cycles for the mean velocity and 50 sample cycles for the RMS velocity are necessary for obtaining converged and correct results in CCV.展开更多
In this article, a reduced mixed finite element （MFE） formulation based on proper orthogonal decomposition （POD） for the non-stationary conduction-convection problems is presented. Also the error estimates between...In this article, a reduced mixed finite element （MFE） formulation based on proper orthogonal decomposition （POD） for the non-stationary conduction-convection problems is presented. Also the error estimates between the reduced MFE solutions based on POD and usual MFE solutions are derived. It is shown by numerical examples that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced MFE formulation based on POD is feasible and efficient in finding numerical solutions for the non-stationary conduction-convection problems.展开更多
基金supported by the National Natural Science Foundation of China(No.12472197)。
文摘A novel gappy technology, gappy autoencoder with proper orthogonal decomposition(Gappy POD-AE), is proposed for reconstructing physical fields from sparse data. High-dimensional data are reduced via proper orthogonal decomposition(POD),and low-dimensional data are used to train an autoencoder(AE). By integrating the POD operator with the decoder, a nonlinear solution form is established and incorporated into a new maximum-a-posteriori(MAP)-based objective for online reconstruction.The numerical results on the two-dimensional(2D) Bhatnagar-Gross-Krook-Boltzmann(BGK-Boltzmann) equation, wave equation, shallow-water equation, and satellite data show that Gappy POD-AE achieves higher accuracy than gappy proper orthogonal decomposition(Gappy POD), especially for the data with slowly decaying singular values,and is more efficient in training than gappy autoencoder(Gappy AE). The MAP-based formulation and new gappy procedure further enhance the reconstruction accuracy.
基金Research Fund for Communications in Western China Under Grant No. 200431800028
文摘Proper Orthogonal Decomposition （POD） provides a powerful modal transformation tool for stochastic dynamics. In this paper, coherency matrix-based proper orthogonal decomposition （CPOD） is presented as an innovative form of the POD based on cross power spectral density matrices. By introducing a discretizing scheme, the CPOD-based spectral representation method is obtained for use in stochastic simulation. Moreover, some criteria are proposed that allow the truncation order of CPOD to be conveniently determined. A numerical example to illustrate the application of the proposed method for the simulation of a wind velocity field is provided.
基金the National Science Foundation of China(Grant Nos.51804033 and 51936001)China Postdoctoral Science and Foundation(Grant No.2018M641254)+3 种基金Beijing Postdoctoral Research Foundation(2018-ZZ-045)the Project of Construction of Innovative Teams and Teacher Career Development for Universities and Colleges Under Beijing Municipality(Grant No.IDHT20170507)Program of Great Wall Scholar(Grant No.CIT&TCD20180313)Jointly Projects of Beijing Natural Science Foundation and Beijing Municipal Education Commission(Grant No.KZ201810017023).
文摘The hydraulic fracturing is a nonlinear,fluid-solid coupling and transient problem,in most cases it is always time-consuming to simulate this process numerically.In recent years,although many numerical methods were proposed to settle this problem,most of them still require a large amount of computer resources.Thus it is a high demand to develop more efficient numerical approaches to achieve the real-time monitoring of the fracture geometry during the hydraulic fracturing treatment.In this study,a reduced order modeling technique namely Proper Generalized Decomposition(PGD),is applied to accelerate the simulations of the transient,non-linear coupled system of hydraulic fracturing problem,to match this extremely tight response time constraint.The separability of the solution in space and time dimensions is studied for a simplified model problem.The solid and fluid equations are coupled explicitly by inverting the solid discrete problem,and a simple iterative procedure to handle the non-linear characteristic of the hydraulic fracturing problem is proposed in this work.Numeral validation illustrates that the results of PGD match well with these of standard finite element method in terms o f fracture opening and fluid pressure in the hydro-fracture.Moreover,after the off-line calculations,the numerical results can be obtained in real time.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11332006,11732010,11572221,and 11502066)the Natural Science Foundation of Tianjin City(Grant No.18JCQNJC5100)
文摘Time-resolved particle image velocimetry（TRPIV） experiments are performed to investigate the coherent structure＇s performance of riblets in a turbulent boundary layer（TBL） at a friction Reynolds number of 185. To visualize the energetic large-scale coherent structures（CSs） over a smooth surface and riblets, the proper orthogonal decomposition（POD） and finite-time Lyapunov exponent（FTLE） are used to identify the CSs in the TBL. Spatial-temporal correlation is implemented to obtain the characters and transport properties of typical CSs in the FTLE fields. The results demonstrate that the generic flow structures, such as hairpin-like vortices, are also observed in the boundary layer flow over the riblets, consistent with its smooth counterpart. Low-order POD modes are more sensitive to the riblets in comparison with the high-order ones,and the wall-normal movement of the most energy-containing structures are suppressed over riblets. The spatial correlation analysis of the FTLE fields indicates that the evolution process of the hairpin vortex over riblets are inhibited. An apparent decrease of the convection velocity over riblets is noted, which is believed to reduce the ejection/sweep motions associated with high shear stress from the viscous sublayer. These reductions exhibit inhibition of momentum transfer among the structures near the wall in the TBL flows.
基金National Natural Science Foundation of China (No.10671153)
文摘The proper orthogonal decomposition (POD) method for the instationary Navier-Stokes equations is considered. Several numerical approaches to evaluating the POD eigenfunctions are presented. The POD eigenfunctions are applied as a basis for a Galerkin projection of the instationary Navier-Stokes equations. And a low-dimensional ordinary differential models for fluid flows governed by the instationary Navier-Stokes equations are constructed. The numerical examples show that the method is feasible and efficient for optimal control of fluids.
基金Project supported by the National Natural Science Foundation of China(Nos.11232011,11402262,11572314,and 11621202)
文摘The reduced-order model （ROM） for the two-dimensional supersonic cavity flow based on proper orthogonal decomposition （POD） and Galerkin projection is investigated. Presently, popular ROMs in cavity flows are based on an isentropic assumption, valid only for flows at low or moderate Mach numbers. A new ROM is constructed involving primitive variables of the fully compressible Navier-Stokes （N-S） equations, which is suitable for flows at high Mach numbers. Compared with the direct numerical simulation （DNS） results, the proposed model predicts flow dynamics （e.g., dominant frequency and amplitude） accurately for supersonic cavity flows, and is robust. The comparison between the present transient flow fields and those of the DNS shows that the proposed ROM can capture self-sustained oscillations of a shear layer. In addition, the present model reduction method can be easily extended to other supersonic flows.
文摘A three-dimensional incompressible annular jet is simulated by the large eddy simulation(LES)method at a Reynolds number Re=8500.The time-averaged velocity field shows an asymmetric wake behind the central bluff-body although the flow geometry is symmetric.The proper orthogonal decomposition(POD)analysis of the velocity fluctuation vectors is conducted to study the flow dynamics of the wake flow.The distribution of turbulent kinetic energy across the three-dimensional POD modes shows that the first four eigenmodes each capture more than 1%of the turbulent kinetic energy,and hence their impact on the wake dynamics is studied.The results demonstrate that the asymmetric mean flow in the near-field of the annular jet is related to the first two POD modes which correspond to a radial shift of the stagnation point.The modes 3 and 4 involve the stretching or squeezing effects of the recirculation region in the radial direction.In addition,the spatial structure of these four POD eigenmodes also shows the counter-rotating vortices in the streamwise direction downstream of the flow reversal region.
基金Acknowledgements The authors are grateful for the support of this research by the Committee of National Science Foundation of China （50908077） and Foundation of Heilongjiang Province Educational Committee （11551368）.
文摘Proper orthogonal decomposition （POD） is an effective statistical technique for data reduction and feature extraction of the random field including the wind field. This paper introduces the theory of the POD and illustrates engineering of structures. Using the POD technique, it is shown that wind pressure data can be accurately reconstructed with a limited number of modes using the wind pressure data from wind tunnel test. Comparing the reconstructed values by POD with the original measured values from the wind tunnel test both in the time and frequency domains, it is concluded that the proper orthogonal decomposition（POD） is an efficient and practical technique for deriving the random wind pressure field from limited known data as shown in the pitched roof example in this paper.
基金This work was supported by Natural Science Basic Research Program of Shaanxi(Program No.2023-JC-YB-269)the National Natural Science Foundation of China(Grant No.41974122).
文摘The large storage requirement is a critical issue in cross-correlation imaging-condition based reverse time migration(RTM),because it requires the operation of the source and receiver wavefields at the same time.The boundary value method(BVM),based on the finite difference method(FDM),can be used to reconstruct the source wavefield in the reverse time propagation in the same way as the receiver wavefield,which can reduce the storage burden of the RTM data.Considering that the FDM cannot well handle models with discontinuous material properties and rough interfaces,we develop a source wavefield reconstruction strategy based on the finite element method(FEM),using proper orthogonal decomposition(POD)to enhance computational efficiency.In this method,we divide the whole time period into several segments,and construct the POD basis functions to get a reduced order model(ROM)for the source wavefield reconstruction in each segment.We show the corresponding quantitative analysis of the storage requirement of the POD-FEM.Numerical tests on the homogeneous model show the effectiveness of the proposed method,while the layered model and part of the Marmousi model tests indicate that the POD-FEM can keep an excellent balance between computational efficiency and memory usage compared with the full-stored method(FSM)and the BVM,and can be effectively applied in imaging.
文摘Over the recent years there has been an increased trend in the use of Reduced Order Models （ROM） for modeling the coupled aeroelastic system. Of all the ROM models, the Proper Orthogonal Decomposition Method （POD） has been the most widely used, reason being the relative simplicity of implementation and the physical insight that it offers towards the physical problem. In this paper we begin by briefly recalling the recent work using POD for the computational aeroelasticity followed by the mathematical formulation. Mathematical formulation is important as it provides understanding of how POD method works. Implementation issues related to the POD method are presented next. Since POD is an empirical technique therefore, it is marred by the robustness issues as is the case with all the ROMs. In the end the variations of POD method, developed over the years are presented along with the most recent trend of using hybrid ROM.
基金supported by National Natural Science Foundation of China(Grant No.10832001)Vision Foundation of Beijing University of Aeronautics and Astronautics (Grant No.YWF-10-20-003)
文摘The present work uses dynamic mode decomposition(DMD) to analyze wake flow of NACA0015 airfoil with Gurney flap.The physics of DMD is first introduced.Then the PIV-measured wake flow velocity field is decomposed into dynamical modes.The vortex shedding pattern behind the trailing edge and its high-order harmonics have been captured with abundant information such as frequency,wavelength and convection speed.It is observed that high-order dynamic modes convect faster than low-order modes;moreover the wavelength of the dynamic modes scales with the corresponding frequency in power law.
基金supported by the Hainan Provincial Joint Project of Sanya Yazhou Bay Science and Technology City(Grant No.2021CXLH0003)the Innovational Fund for Scientific and Technological Personnel of Hainan Province(Grant No.KJRC2023D37).
文摘Bends contribute to a flexible layout of pipeline system,but also lead to intensive energy costs due to the complex flow characteristic.This experimental study is conducted to investigate the impact of a single coarse particle on the flow field in a bend.The velocity profiles of fluid on the axial symmetry plane of the bend are measured using time-resolved particle image velocimetry.The flow structures are extracted using the proper orthogonal decomposition method.The results reveal that there is a shear-layer flow in the bend during the transportation.With the increase in particle size,the particle has a dominant influence on the flow energy distribution of the overall flow.The impact of particles on the first few energetic flows is mainly in the latter part of the transportation,both temporally and spatially.As the particle size decreases,the most energetic unsteady flow within the bend changes from the convective flow to the flow of the shear layer.
基金sponsored by the National Natural Science Foundation of China(Grant No.32373187)the Research Fund for International Scientists of the National Natural Science Foundation of China(Grant No.32350410404)the Natural Science Foundation of Shanghai(Grant No.23ZR1427000).
文摘This study experimentally investigates the hydrodynamic characteristics,geometric configurations,fluttering motions of the codend,and the instantaneous flow fields inside and around the codend,with and without a liner,under varying catch sizes and inflow velocities.A proper orthogonal decomposition method is employed to extract phase-averaged mean properties of unsteady turbulent flows from flow measurement data obtained using an electromagnetic current velocity meter inside and around the codend.The results reveal that as catch size increases,the drag force,codend motion,Reynolds number,and codend volume increase while the drag coefficient decreases.Owing to the codend shape and pronounced motion,a complex fluid–structure interaction occurs,demonstrating a strong correlation between drag force and codend volume.The oscillation amplitudes of the hydrodynamic forces and codend motions increase with increasing catch size,and their oscillations mainly involve low-frequency activity.A significant reduction in the flow field occurs inside and around the unlined codend without a catch.The flow field is 5.81%,14.39%,and 27.01%lower than the unlined codend with a catch,the codend with a liner but without a catch,and the codend with both a liner and a catch,respectively.Fourier analysis reveals that the codend motions and hydrodynamic forces are mainly characterized by low-frequency activity and are synchronized with the unsteady turbulent flow street.Furthermore,the proper orthogonal decomposition results reveal the development of unsteady turbulent flow inside and around the codend,driven by flow passage blockage caused by the presence of the liner,intense codend motions,and the catch.Understanding the hydrodynamic characteristics and flow instabilities inside and around the codend,particularly those associated with its fluttering motions,is crucial for optimizing trawl design and improving trawl selectivity.
基金financially supported by the National Natural Science Foundation of China(Grant No.11472076).
文摘Jacket platforms constitute the foundational infrastructure of offshore oil and gas field exploitation.How to efficiently and accurately monitor the mechanical properties of jacket structures is one of the key problems to be solved to ensure the safe operation of the platform.To address the practical engineering problem that it is difficult to monitor the stress response of the tubular joints of jacket platforms online,a digital twin reduced-order method for real-time prediction of the stress response of tubular joints is proposed.In the offline construction phase,multi-scale modeling and multi-parameter experimental design methods are used to obtain the stress response data set of the jacket structure.Proper orthogonal decomposition is employed to extract the main feature information from the snapshot matrix,resulting in a reduced-order basis.The leave-one-out cross-validation method is used to select the optimal modal order for constructing the reduced-order model(ROM).In the online prediction phase,a digital twin model of the tubular joint is established,and the prediction performance of the ROM is analyzed and verified through using random environmental load and field environmental monitoring data.The results indicate that,compared with traditional numerical simulations of tubular joints,the ROM based on the proposed reduced-order method is more efficient in predicting the stress response of tubular joints while ensuring accuracy and robustness.
基金supported by the UCAS-ULille joint PhD training program and also partially supported within the frame of the EE4.0(Electrical Energy 4.0)project,which is co-financed by the European Union,with the financial support of the European Regional Development Fund(ERDF),the French state,and the French Region of Hautsde-France.
文摘Digital twin is considered the key technique for real-time monitoring and life-cycle management of electric equipment.To construct the digital twin model of electric equipment,a multi-parameter electromagnetic analysis is needed to generate a large amount of high fidelity data under various working condition.However,repeated solving such multi-parameter electromagnetic problems based on full order finite element method may lead to extreme scale calculations.To address this issue,a hybrid approach that combines tensor decomposition and proper orthogonal decomposition(POD)is introduced,which can effectively establish a reduced order model for multi-parameter electromagnetic field problems.The performance of the proposed approach is illustrated through three numerical examples,namely an electrical motor,a transformer and a voice coil actuator including parameter variations of operating conditions,geometric parameters,and material parameters.The numerical results show that the proposed hybrid approach has significant advantage compared to conventional reduced order model in situation where the solution changes dramatically within the parameter variation range and even more apparent when the parameter dimension is high.
基金supported by National Natural Science Foundation of China (Grant Nos. 10871022, 10771065,and 60573158)Natural Science Foundation of Hebei Province (Grant No. A2007001027)
文摘A proper orthogonal decomposition (POD) method is applied to a usual finite element (FE) formulation for parabolic equations so that it is reduced into a POD FE formulation with lower dimensions and enough high accuracy. The errors between the reduced POD FE solution and the usual FE solution are analyzed. It is shown by numerical examples that the results of numerical computations are consistent with theoretical conclusions. Moreover, it is also shown that this validates the feasibility and efficiency of POD method.
基金the National Natural Science Foundation of China(Grant Nos.10471100,40437017,and 60573158)Beijing Jiaotong University Science and Technology Foundation
文摘The proper orthogonal decomposition(POD)and the singular value decomposition(SVD) are used to study the finite difference scheme(FDS)for the nonstationary Navier-Stokes equations. Ensembles of data are compiled from the transient solutions computed from the discrete equation system derived from the FDS for the nonstationary Navier-Stokes equations.The optimal orthogonal bases are reconstructed by the elements of the ensemble with POD and SVD.Combining the above procedures with a Galerkin projection approach yields a new optimizing FDS model with lower dimensions and a high accuracy for the nonstationary Navier-Stokes equations.The errors between POD approximate solutions and FDS solutions are analyzed.It is shown by considering the results obtained for numerical simulations of cavity flows that the error between POD approximate solution and FDS solution is consistent with theoretical results.Moreover,it is also shown that this validates the feasibility and efficiency of POD method.
基金supported by the National Natural Science Foundation of China (No. 11502211)
文摘In aerodynamic optimization, global optimization methods such as genetic algorithms are preferred in many cases because of their advantage on reaching global optimum. However,for complex problems in which large number of design variables are needed, the computational cost becomes prohibitive, and thus original global optimization strategies are required. To address this need, data dimensionality reduction method is combined with global optimization methods, thus forming a new global optimization system, aiming to improve the efficiency of conventional global optimization. The new optimization system involves applying Proper Orthogonal Decomposition（POD） in dimensionality reduction of design space while maintaining the generality of original design space. Besides, an acceleration approach for samples calculation in surrogate modeling is applied to reduce the computational time while providing sufficient accuracy. The optimizations of a transonic airfoil RAE2822 and the transonic wing ONERA M6 are performed to demonstrate the effectiveness of the proposed new optimization system. In both cases, we manage to reduce the number of design variables from 20 to 10 and from 42 to 20 respectively. The new design optimization system converges faster and it takes 1/3 of the total time of traditional optimization to converge to a better design, thus significantly reducing the overall optimization time and improving the efficiency of conventional global design optimization method.
基金supported by the National Natural Science Foundation of China(Grant Nos.51176020 and 51376029)
文摘Large eddy simulation(LES)is used to calculate the in-cylinder turbulent flow field in a direct injection spark ignition(DISI)engine.The computations are carried out for three different maximum valve lifts(MVL)and throughout 100 consecutive engine cycles.The simulated results as well as corresponding particle image velocimetry(PIV)measurement database are analyzed by the proper orthogonal decomposition(POD)method.Through a new developed POD quadruple decomposition the instantaneous in-cylinder flow fields are decomposed into four parts,named mean field,coherent field,transition field and turbulent field,respectively.Then the in-cylinder turbulent flow characteristics and cycle-to-cycle variations(CCV)are studied separately upon the four part flow fields.Results indicate that each part exhibits its specific characteristics and has close connection with others.The mean part contains more than 50%of the total kinetic energy and the energy cascade phenomenon occurs among the four part fields;the coherent field part possesses the highest CCV level which dominates CCV of the bulk flow.In addition,it is observed that a change in MVL affects significantly the in-cylinder flow behavior including CCV,especially for the coherent part.Furthermore,the POD analysis demonstrates that at least 25 sample cycles for the mean velocity and 50 sample cycles for the RMS velocity are necessary for obtaining converged and correct results in CCV.
基金supported by the National Science Foundation of China (10871022 11061009+5 种基金 40821092)the National Basic Research Program (2010CB428403 2009CB421407 2010CB951001)Natural Science Foundation of Hebei Province (A2010001663)Chinese Universities Scientific Fund (2009-2-05)
文摘In this article, a reduced mixed finite element （MFE） formulation based on proper orthogonal decomposition （POD） for the non-stationary conduction-convection problems is presented. Also the error estimates between the reduced MFE solutions based on POD and usual MFE solutions are derived. It is shown by numerical examples that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced MFE formulation based on POD is feasible and efficient in finding numerical solutions for the non-stationary conduction-convection problems.
