The objective of this paper is to present and to validate a new hybrid coupling （HC） algorithm for modeling of fluid-structure interaction （FSI） in incompressible, viscous flows. The HC algorithm is able to avoid ...The objective of this paper is to present and to validate a new hybrid coupling （HC） algorithm for modeling of fluid-structure interaction （FSI） in incompressible, viscous flows. The HC algorithm is able to avoid numerical instability issues associated with artificial added mass effects, which are often encountered by standard loosely coupled （LC） and tightly coupled （TC） algorithms, when modeling the FSI response of flexible structures in incompressible flow. The artificial added mass effect is caused by the lag in exchange of interfacial displacements and forces between the fluid and solid solvers in partitioned algorithms. The artificial added mass effect is much more prominent for light/flexible struc- tures moving in water, because the fluid forces are in the same order of magnitude as the solid forces, and because the speed at which numerical errors propagate in an incom- pressible fluid. The new HC algorithm avoids numerical instability issues associated with artificial added mass effects by embedding Theodorsen＇s analytical approximation of the hydroelastic forces in the solution process to obtain better initial estimates of the displacements. Details of the new HC algorithm are presented. Numerical validation studies are shown for the forced pitching response of a steel and a plastic hydrofoil. The results show that the HC algorithm is able to converge faster, and is able to avoid numerical insta- bility issues, compared to standard LC and TC algorithms, when modeling the transient FSI response of a plastic hydrofoil. Although the HC algorithm is only demonstrated for a NACA0009 hydrofoil subject to pure pitching motion, the method can be easily extended to model general 3-D FSI response and stability of complex, flexible structures in turbulent, incompressible, multiphase flows.展开更多
In this paper we deal with a nonlinear interaction problem between an incompressible viscous fluid and a nonlinear thermoelastic plate.The nonlinearity in the plate equation corresponds to nonlinear elastic force in v...In this paper we deal with a nonlinear interaction problem between an incompressible viscous fluid and a nonlinear thermoelastic plate.The nonlinearity in the plate equation corresponds to nonlinear elastic force in various physically relevant semilinear and quasilinear plate models.We prove the existence of a weak solution for this problem by constructing a hybrid approximation scheme that,via operator splitting,decouples the system into two sub-problems,one piece-wise stationary for the fluid and one time-continuous and in a finite basis for the structure.To prove the convergence of the approximate quasilinear elastic force,we develop a compensated compactness method that relies on the maximal monotonicity property of this nonlinear function.展开更多
Design for structural integrity requires an appreciation of where stress singularities can occur in structural configurations. While there is a rich literature devoted to the identification of such singular behavior i...Design for structural integrity requires an appreciation of where stress singularities can occur in structural configurations. While there is a rich literature devoted to the identification of such singular behavior in solid mechanics, to date there has been relatively little explicit identification of stress singularities caused by fluid flows. In this study, stress and pressure singularities induced by steady flows of viscous incompressible fluids are asymptotically identified. This is done by taking advantage of an earlier result that the Navier-Stokes equations are locally governed by Stokes flow in angular corners. Findings for power singularities are confirmed by developing and using an analogy with solid mechanics. This analogy also facilitates the identification of flow-induced log singularities. Both types of singularity are further confirmed for two global configurations by applying convergence-divergence checks to numerical results. Even though these flow-induced stress singularities are analogous to singularities in solid mechanics, they nonetheless render a number of structural configurations singular that were not previously appreciated as such from identifications within solid mechanics alone.展开更多
Centrifugal compressors with parallel-wall and contracting wall vaneless diffuser are designed by using centrifugal compressor computer-aided integrated design system. The internal flow fields of the compressor are ca...Centrifugal compressors with parallel-wall and contracting wall vaneless diffuser are designed by using centrifugal compressor computer-aided integrated design system. The internal flow fields of the compressor are calculated by solving three-dimensional Navier-Stokes equation. Four aspects are investigated and calculation results show that the total efficiencies and total pressure ratios of the compressor with contracting wall vandess diffuser is higher than that of the compressor with parallel-wall. The jet and wake don＇t mix rapidly inside vandess diffuser. The outlet blade lean angle doesn＇t affect the compressor performance. The greater the mass flow rate through impeller, the more uneven the velocity distribution at impeller outlet is.展开更多
Design for structural integrity requires an appreciation of where stress singularities can occur in structural configurations. While there is a rich literature devoted to the identification of such singular behavior i...Design for structural integrity requires an appreciation of where stress singularities can occur in structural configurations. While there is a rich literature devoted to the identification of such singular behavior in solid mechanics, only of late has there been much in the way of corresponding identifications of flow-induced stress singularities in fluid mechanics. These recent asymptotic identifications are for a single incompressible viscous fluid： Here the asymptotic approach is extended to apply to a configuration entailing two such fluids, For this configuration, various specifications leading to power or log singularities are determined. These results demonstrate that flow-induced stress singularities can occur in a structural container at a location where no singularities are identified within solid mechanics alone.展开更多
Study of the how held around the large scale offshore structures under the action of waves and viscous currents is of primary importance for the scouring estimation and protection in the vicinity of the structures. Bu...Study of the how held around the large scale offshore structures under the action of waves and viscous currents is of primary importance for the scouring estimation and protection in the vicinity of the structures. But very little has been known in its mechanism when the viscous effects is taken into consideration. As a part of the efforts to tackle the problem, a numerical model is presented for the simulation of the how held around a fixed vertical truncated circular cylinder subjected to waves and viscous currents based on the depth-averaged Reynolds equations and depth-averaged k-epsilon turbulence model. Finite difference method with a suitable iteration defect correct method and an artificial open boundary condition are adopted in the numerical process. Numerical results presented relate to the interactions of a pure incident viscous current with Reynolds number Re = 10(5), a pure incident regular sinusoidal wave, and the coexisting of viscous current and wave with a circular cylinder, respectively. Flow fields associated with the hydrodynamic coefficients of the fixed cylinder, as well as corresponding free surface profiles and wave amplitudes, are discussed. The present method is found to be relatively straightforward, computationally effective and numerically stable for treating the problem of interactions among waves, viscous currents and bodies.展开更多
Interactions between different scales in turbulence were studied starting from the incompressible Navier-Stokes equations. The integral and differential formulae of the short-range viscous stresses, which express the ...Interactions between different scales in turbulence were studied starting from the incompressible Navier-Stokes equations. The integral and differential formulae of the short-range viscous stresses, which express the short-range interactions between contiguous scales in turbulence,were given. A concept of the resonant-range interactions between extreme contiguous scales was introduced and the differential formula of the resonant-range viscous stresses was obtained. The short- and resonant-range viscous stresses were applied to deduce the large-eddy simulation (LES) equations as well as the multiscale equations, which are approximately closed and do not contain any empirical constants or relations. The properties and advantages of using the multiscale equations to compute turbulent flows were discussed. The short-range character of the interactions between the scales in turbulence means that the multiscale simulation is a very valuable technique for the calculation of turbulent flows. A few numerical examples were also given.展开更多
An earlier asymptotic analysis identified the flow-induced power and log singularities present when two nonmiscible,incompressible,viscous fluids contact a wall with no slip permitted.Here a parallel asymptotic analys...An earlier asymptotic analysis identified the flow-induced power and log singularities present when two nonmiscible,incompressible,viscous fluids contact a wall with no slip permitted.Here a parallel asymptotic analysis identifies such singularities in pressures and stresses when slip is admitted.While the singularity exponents for the two types of contact conditions are similar,the angles of interface impingement that have either power or log singularities are usually quite distinct.展开更多
A tilt-rotor unsteady flow analytical method has been developed based upon viscous vortex-particle method.In this method,the vorticity field is divided into small assembled vortex particles.Vortex motion and diffusion...A tilt-rotor unsteady flow analytical method has been developed based upon viscous vortex-particle method.In this method,the vorticity field is divided into small assembled vortex particles.Vortex motion and diffusion are obtained by solving the velocity-vorticity-formed incompressible Navier-Stokes equations using agrid-free Lagrangian simulation method.Generation of the newly vortex particles is calculated by using the Weissinger-L lifting surface model.Furthermore,in order to significantly improve computational efficiency,a fast multiple method(FMM)is introduced into the calculation of induced velocity and its gradient.Finally,the joint vertical experimental(JVX)tilt-rotor is taken as numerical examples to analyze.The wake geometry and downwash are investigated for both hover and airplane modes.The proposed method for tilt-rotor flow analysis is verified by comparing its results with those available measured data.Comparison indicates that the current method can accurately capture the complicated tilt-rotor wake variation and be suitable for aerodynamic interaction simulation in complex environments.Additionally,the aerodynamic interactional characteristics of dual-rotor wake are discussed in different rotor distance.Results show that there are significant differences on interactional characteristics between hover mode and airplane mode.展开更多
To compute transonic flows over a complex 3D aircraft configuration, a viscous/inviscid interaction method is developed by coupling an integral boundary-layer solver with an Eluer solver in a ＂semi-inverse＂ manner. ...To compute transonic flows over a complex 3D aircraft configuration, a viscous/inviscid interaction method is developed by coupling an integral boundary-layer solver with an Eluer solver in a ＂semi-inverse＂ manner. For the turbulent boundary-layer, an integral method using Green＇s lag equation is coupled with the outer inviscid flow. A blowing velocity approach is used to simulate the displacement effects of the boundary layer. To predict the aerodynamic drag, it is developed a numerical technique called far-field method that is based on the momentum theorem, in which the total drag is divided into three component drags, i.e. viscous, induced and wave-formed. Consequently, it can provide more physical insight into the drag sources than the often-used surface integral technique. The drag decomposition can be achieved with help of the second law of thermodynamics, which implies that entropy increases and total pressure decreases only across shock wave along a streamline of an inviscid non-isentropic flow. This method has been applied to the DLR-F4 wing/body configuration showing results in good agreement with the wind tunnel data.展开更多
Gao＇s viscous/in-viscid interacting shear flows （ISF） theory, proposed by professor Gao Zhi in Institute of Mechanics, China Academy of Science, and its inferences and their applications in computational fluid dyna...Gao＇s viscous/in-viscid interacting shear flows （ISF） theory, proposed by professor Gao Zhi in Institute of Mechanics, China Academy of Science, and its inferences and their applications in computational fluid dynamics （CFD） are reviewed and some subjects worthy to be studied are pro- posed in this paper. The flow-field and motion law of ISF, mathematics definition of strong viscous shear layer flow in ISF, ISF equations, wall-surface compatibility criteria （Gao＇s criteria ）, space scale variety law of strong viscous shear layer reveals flow mechanism and local space small scale triggered by strong interaction that cause some abnormal severe local pneumatic heating phenomenon in hypersonic flow. Gao＇s ISF theory was used in near wall flow, free ISF flow simulation and design of computing grids, Gao＇s wall-surface criteria were used to verify calculation reliability and accuracy of near wall flows, ISF theory approximate analytical result of shock waves-boundary layer interac- tion and ISF equations were used to obtain the numerical exact solution of local area flow （ such as stationary point flow）. Some new subjects, such as, improving near-wall turbulent models according to the turbulent flow simulation satisfying the wall-criteria and illustrating relation between grid-con- vergence based on the wall criteria and other convergence tactics, are suggested. The necessity of applying Gao＇s ISF theory and wall criteria is revealed. Difficulties and importance of hypersonic vis- cous/in-viscid interaction phenomenon were also emphasized.展开更多
In this paper, we study the inviscid limit problem for the scalar viscous conservation laws on half plane. We prove that if the solution of the corresponding inviscid equation on half plane is piecewise smooth with a ...In this paper, we study the inviscid limit problem for the scalar viscous conservation laws on half plane. We prove that if the solution of the corresponding inviscid equation on half plane is piecewise smooth with a single shock satisfying the entropy condition, then there exist solutions to the viscous conservation laws which converge to the inviscid solution away from the shock discontinuity and the boundary at a rate of ε1 as the viscosity ε tends to zero.展开更多
基金the financial support provided by the Office of Naval Research(ONR) through grant number N00014-09-1-1204 (managed by Dr. Ki-Han Kim)supported in part by the National Research Foundation of Korea (NRF)grant funded by the Korea government (MEST) through the GCRC-SOP Grant No. 2012-0004783
文摘The objective of this paper is to present and to validate a new hybrid coupling （HC） algorithm for modeling of fluid-structure interaction （FSI） in incompressible, viscous flows. The HC algorithm is able to avoid numerical instability issues associated with artificial added mass effects, which are often encountered by standard loosely coupled （LC） and tightly coupled （TC） algorithms, when modeling the FSI response of flexible structures in incompressible flow. The artificial added mass effect is caused by the lag in exchange of interfacial displacements and forces between the fluid and solid solvers in partitioned algorithms. The artificial added mass effect is much more prominent for light/flexible struc- tures moving in water, because the fluid forces are in the same order of magnitude as the solid forces, and because the speed at which numerical errors propagate in an incom- pressible fluid. The new HC algorithm avoids numerical instability issues associated with artificial added mass effects by embedding Theodorsen＇s analytical approximation of the hydroelastic forces in the solution process to obtain better initial estimates of the displacements. Details of the new HC algorithm are presented. Numerical validation studies are shown for the forced pitching response of a steel and a plastic hydrofoil. The results show that the HC algorithm is able to converge faster, and is able to avoid numerical insta- bility issues, compared to standard LC and TC algorithms, when modeling the transient FSI response of a plastic hydrofoil. Although the HC algorithm is only demonstrated for a NACA0009 hydrofoil subject to pure pitching motion, the method can be easily extended to model general 3-D FSI response and stability of complex, flexible structures in turbulent, incompressible, multiphase flows.
基金partially supported by National Natural Science Foundation of China(11631008)。
文摘In this paper we deal with a nonlinear interaction problem between an incompressible viscous fluid and a nonlinear thermoelastic plate.The nonlinearity in the plate equation corresponds to nonlinear elastic force in various physically relevant semilinear and quasilinear plate models.We prove the existence of a weak solution for this problem by constructing a hybrid approximation scheme that,via operator splitting,decouples the system into two sub-problems,one piece-wise stationary for the fluid and one time-continuous and in a finite basis for the structure.To prove the convergence of the approximate quasilinear elastic force,we develop a compensated compactness method that relies on the maximal monotonicity property of this nonlinear function.
文摘Design for structural integrity requires an appreciation of where stress singularities can occur in structural configurations. While there is a rich literature devoted to the identification of such singular behavior in solid mechanics, to date there has been relatively little explicit identification of stress singularities caused by fluid flows. In this study, stress and pressure singularities induced by steady flows of viscous incompressible fluids are asymptotically identified. This is done by taking advantage of an earlier result that the Navier-Stokes equations are locally governed by Stokes flow in angular corners. Findings for power singularities are confirmed by developing and using an analogy with solid mechanics. This analogy also facilitates the identification of flow-induced log singularities. Both types of singularity are further confirmed for two global configurations by applying convergence-divergence checks to numerical results. Even though these flow-induced stress singularities are analogous to singularities in solid mechanics, they nonetheless render a number of structural configurations singular that were not previously appreciated as such from identifications within solid mechanics alone.
文摘Centrifugal compressors with parallel-wall and contracting wall vaneless diffuser are designed by using centrifugal compressor computer-aided integrated design system. The internal flow fields of the compressor are calculated by solving three-dimensional Navier-Stokes equation. Four aspects are investigated and calculation results show that the total efficiencies and total pressure ratios of the compressor with contracting wall vandess diffuser is higher than that of the compressor with parallel-wall. The jet and wake don＇t mix rapidly inside vandess diffuser. The outlet blade lean angle doesn＇t affect the compressor performance. The greater the mass flow rate through impeller, the more uneven the velocity distribution at impeller outlet is.
文摘Design for structural integrity requires an appreciation of where stress singularities can occur in structural configurations. While there is a rich literature devoted to the identification of such singular behavior in solid mechanics, only of late has there been much in the way of corresponding identifications of flow-induced stress singularities in fluid mechanics. These recent asymptotic identifications are for a single incompressible viscous fluid： Here the asymptotic approach is extended to apply to a configuration entailing two such fluids, For this configuration, various specifications leading to power or log singularities are determined. These results demonstrate that flow-induced stress singularities can occur in a structural container at a location where no singularities are identified within solid mechanics alone.
基金The project supported by the National Natural Science Foundation of China Foundation of State Key Laboratory of Ocean Engineering at Shanghai Jiao Tong University
文摘Study of the how held around the large scale offshore structures under the action of waves and viscous currents is of primary importance for the scouring estimation and protection in the vicinity of the structures. But very little has been known in its mechanism when the viscous effects is taken into consideration. As a part of the efforts to tackle the problem, a numerical model is presented for the simulation of the how held around a fixed vertical truncated circular cylinder subjected to waves and viscous currents based on the depth-averaged Reynolds equations and depth-averaged k-epsilon turbulence model. Finite difference method with a suitable iteration defect correct method and an artificial open boundary condition are adopted in the numerical process. Numerical results presented relate to the interactions of a pure incident viscous current with Reynolds number Re = 10(5), a pure incident regular sinusoidal wave, and the coexisting of viscous current and wave with a circular cylinder, respectively. Flow fields associated with the hydrodynamic coefficients of the fixed cylinder, as well as corresponding free surface profiles and wave amplitudes, are discussed. The present method is found to be relatively straightforward, computationally effective and numerically stable for treating the problem of interactions among waves, viscous currents and bodies.
文摘Interactions between different scales in turbulence were studied starting from the incompressible Navier-Stokes equations. The integral and differential formulae of the short-range viscous stresses, which express the short-range interactions between contiguous scales in turbulence,were given. A concept of the resonant-range interactions between extreme contiguous scales was introduced and the differential formula of the resonant-range viscous stresses was obtained. The short- and resonant-range viscous stresses were applied to deduce the large-eddy simulation (LES) equations as well as the multiscale equations, which are approximately closed and do not contain any empirical constants or relations. The properties and advantages of using the multiscale equations to compute turbulent flows were discussed. The short-range character of the interactions between the scales in turbulence means that the multiscale simulation is a very valuable technique for the calculation of turbulent flows. A few numerical examples were also given.
文摘An earlier asymptotic analysis identified the flow-induced power and log singularities present when two nonmiscible,incompressible,viscous fluids contact a wall with no slip permitted.Here a parallel asymptotic analysis identifies such singularities in pressures and stresses when slip is admitted.While the singularity exponents for the two types of contact conditions are similar,the angles of interface impingement that have either power or log singularities are usually quite distinct.
基金Supported by the National Natural Science Foundation of China(11302103)
文摘A tilt-rotor unsteady flow analytical method has been developed based upon viscous vortex-particle method.In this method,the vorticity field is divided into small assembled vortex particles.Vortex motion and diffusion are obtained by solving the velocity-vorticity-formed incompressible Navier-Stokes equations using agrid-free Lagrangian simulation method.Generation of the newly vortex particles is calculated by using the Weissinger-L lifting surface model.Furthermore,in order to significantly improve computational efficiency,a fast multiple method(FMM)is introduced into the calculation of induced velocity and its gradient.Finally,the joint vertical experimental(JVX)tilt-rotor is taken as numerical examples to analyze.The wake geometry and downwash are investigated for both hover and airplane modes.The proposed method for tilt-rotor flow analysis is verified by comparing its results with those available measured data.Comparison indicates that the current method can accurately capture the complicated tilt-rotor wake variation and be suitable for aerodynamic interaction simulation in complex environments.Additionally,the aerodynamic interactional characteristics of dual-rotor wake are discussed in different rotor distance.Results show that there are significant differences on interactional characteristics between hover mode and airplane mode.
文摘To compute transonic flows over a complex 3D aircraft configuration, a viscous/inviscid interaction method is developed by coupling an integral boundary-layer solver with an Eluer solver in a ＂semi-inverse＂ manner. For the turbulent boundary-layer, an integral method using Green＇s lag equation is coupled with the outer inviscid flow. A blowing velocity approach is used to simulate the displacement effects of the boundary layer. To predict the aerodynamic drag, it is developed a numerical technique called far-field method that is based on the momentum theorem, in which the total drag is divided into three component drags, i.e. viscous, induced and wave-formed. Consequently, it can provide more physical insight into the drag sources than the often-used surface integral technique. The drag decomposition can be achieved with help of the second law of thermodynamics, which implies that entropy increases and total pressure decreases only across shock wave along a streamline of an inviscid non-isentropic flow. This method has been applied to the DLR-F4 wing/body configuration showing results in good agreement with the wind tunnel data.
基金Supported by the National Natural Science Foundation(10702009)
文摘Gao＇s viscous/in-viscid interacting shear flows （ISF） theory, proposed by professor Gao Zhi in Institute of Mechanics, China Academy of Science, and its inferences and their applications in computational fluid dynamics （CFD） are reviewed and some subjects worthy to be studied are pro- posed in this paper. The flow-field and motion law of ISF, mathematics definition of strong viscous shear layer flow in ISF, ISF equations, wall-surface compatibility criteria （Gao＇s criteria ）, space scale variety law of strong viscous shear layer reveals flow mechanism and local space small scale triggered by strong interaction that cause some abnormal severe local pneumatic heating phenomenon in hypersonic flow. Gao＇s ISF theory was used in near wall flow, free ISF flow simulation and design of computing grids, Gao＇s wall-surface criteria were used to verify calculation reliability and accuracy of near wall flows, ISF theory approximate analytical result of shock waves-boundary layer interac- tion and ISF equations were used to obtain the numerical exact solution of local area flow （ such as stationary point flow）. Some new subjects, such as, improving near-wall turbulent models according to the turbulent flow simulation satisfying the wall-criteria and illustrating relation between grid-con- vergence based on the wall criteria and other convergence tactics, are suggested. The necessity of applying Gao＇s ISF theory and wall criteria is revealed. Difficulties and importance of hypersonic vis- cous/in-viscid interaction phenomenon were also emphasized.
基金Acknowledgments The author is supported by Tianyuan Foundation （No. 11026093） and the National Natural Science Foundation of China （Nos. 11101162, 11071086）.
文摘In this paper, we study the inviscid limit problem for the scalar viscous conservation laws on half plane. We prove that if the solution of the corresponding inviscid equation on half plane is piecewise smooth with a single shock satisfying the entropy condition, then there exist solutions to the viscous conservation laws which converge to the inviscid solution away from the shock discontinuity and the boundary at a rate of ε1 as the viscosity ε tends to zero.
