The Lagrange multiplier method plays an important role in establishing generalized variational principles notonly in tluid mechallics. but also in elasticity. Sometimes, however, one may come across variational crisi...The Lagrange multiplier method plays an important role in establishing generalized variational principles notonly in tluid mechallics. but also in elasticity. Sometimes, however, one may come across variational crisis(somemultipliers vanish identically). failing to achieve his aim. The crisis is caused by the fact that the Inultipliers are treatedas independent variables in the process of variatioll. but after identification they become functions of the originalindependent variables. To overcome it, a Inodified Lagrange multiplier method or semi-inverse method has beenproposed to deduce generalized varistional principles. Some e-camples are given to illustrate its convenience andeffectiveness of the novel method.展开更多
This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the...This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as "fictitious" fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure' s rigid body shape and behavior, a rigid body constraint is enforced on the "fictitious" fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain (DLM/ FD) method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.展开更多
Using the concept of the base forces, a new finite element method (base force element method, BFEM) based on the complementary energy principle is presented for accurate modeling of structures with large displacemen...Using the concept of the base forces, a new finite element method (base force element method, BFEM) based on the complementary energy principle is presented for accurate modeling of structures with large displacements and large rotations. First, the complementary energy of an element is described by taking the base forces as state variables, and is then separated into deformation and rotation parts for the case of large deformation. Second, the control equations of the BFEM based on the complementary energy principle are derived using the Lagrange multiplier method. Nonlinear procedure of the BFEM is then developed. Finally, several examples are analyzed to illustrate the reliability and accuracy of the BFEM.展开更多
By redefining the multiplier associated with inequality constraint as a positive definite function of the originally-defined multiplier, say, u2_i, i=1, 2, ..., m, nonnegative constraints imposed on inequality constra...By redefining the multiplier associated with inequality constraint as a positive definite function of the originally-defined multiplier, say, u2_i, i=1, 2, ..., m, nonnegative constraints imposed on inequality constraints in Karush-Kuhn-Tucker necessary conditions are removed. For constructing the Lagrange neural network and Lagrange multiplier method, it is no longer necessary to convert inequality constraints into equality constraints by slack variables in order to reuse those results dedicated to equality constraints, and they can be similarly proved with minor modification. Utilizing this technique, a new type of Lagrange neural network and a new type of Lagrange multiplier method are devised, which both handle inequality constraints directly. Also, their stability and convergence are analyzed rigorously.展开更多
Roadways excavated in soft rocks at great depth are difficult to be maintained due to large deformation of surrounding rocks, which greatly influences the safety and efficiency of deep resources exploitation. During t...Roadways excavated in soft rocks at great depth are difficult to be maintained due to large deformation of surrounding rocks, which greatly influences the safety and efficiency of deep resources exploitation. During the excavation process of a deep soft rock tunnel, the rock wall may be compacted due to large deformation. In this paper, the technique to address this problem by a two-dimensional (2D) finite element software, large deformation engineering analyses software (LDEAS 1.0), is provided. By using the Lagrange multiplier method, the kinematic constraint of non-penetrating condition and static constraint of Coulomb friction are introduced to the governing equations in the form of incremental displacement. The numerical example demonstrates the efficiency of this technology. Deformations of a transportation tunnel in inclined soft rock strata at the depth of 1 000 m in Qishan coal mine and a tunnel excavated to three different depths are analyzed by two models, i.e. the additive decomposition model and polar decomposition model. It can be found that the deformation of the transportation tunnel is asymmetrical due to the inclination of rock strata. For extremely soft rock, large deformation can converge only for the additive decomposition model. The deformation of surrounding rocks increases with the increase in the tunnel depth for both models. At the same depth, the deformation calculated by the additive decomposition model is smaller than that by the polar decomposition model.展开更多
A chaotic system is bounded, and its trajectory is confined to a certain region which is called the chaotic attractor. No matter how unstable the interior of the system is, the trajectory never exceeds the chaotic att...A chaotic system is bounded, and its trajectory is confined to a certain region which is called the chaotic attractor. No matter how unstable the interior of the system is, the trajectory never exceeds the chaotic attractor. In the present paper, the sphere bound of the generalized Lorenz system is given, based on the Lyapunov function and the Lagrange multiplier method. Furthermore, we show the actual parameters and perform numerical simulations.展开更多
In this paper, by applying Lagrange, multiplier method and high order Lagrange multiplier method [1], we systematically derive coupled potential energy principle.coupled complementary energy principle,and generalized...In this paper, by applying Lagrange, multiplier method and high order Lagrange multiplier method [1], we systematically derive coupled potential energy principle.coupled complementary energy principle,and generalized coupled potential energy principles and generalized coupled complementary energy principles with two and three kinds of variables in photoelasticity.展开更多
The mantle unsteady flows, which are in an incompressible and isoviscous spherical shell, are investigated by using algorithms of the parallel Lagrange multiplier dissonant decomposition method (LMDDM) and the paralle...The mantle unsteady flows, which are in an incompressible and isoviscous spherical shell, are investigated by using algorithms of the parallel Lagrange multiplier dissonant decomposition method (LMDDM) and the parallel Lagrange multiplier discontinuous deformation analyses (LMDDA) in this paper. Some physical fields about mantle flows such as velocity, pressure, temperature, stress and the force to the crust of the Asian continent are calculated on a parallel computer.展开更多
An improved implementation of Distributed Lagrange multiplier/fictitious domain method was presented and used to simulate the interactions between two circular particles sedimenting in a two_dimensional channel. The s...An improved implementation of Distributed Lagrange multiplier/fictitious domain method was presented and used to simulate the interactions between two circular particles sedimenting in a two_dimensional channel. The simulation results were verified by comparison with experiments. The results show that the interactions between two particles with different sizes can be described as drafting, kissing, tumbling and separating. Only for small diameter ratio, the two particles will interact undergoing repeated DKT (Drafting, Kissing and Tumbling) process. Otherwise, the two particles will separate after their tumbling. The results also show that, during the interaction process, the motion of the small particle is strongly affected while the large particle is affected slightly.展开更多
An improved implementation of distributed multiplier/fictitious domain method is presented for the direct numerical simulation of particulate flow. The key improvement is to replace a finite element triangulation for...An improved implementation of distributed multiplier/fictitious domain method is presented for the direct numerical simulation of particulate flow. The key improvement is to replace a finite element triangulation for the velocity and a “twice coarser' triangulation for the pressure with a rectangular discretization for the velocity and pressure. For code validation, the sedimentation of a single particle in a two dimensional channel was simulated. The results showed that the simulation is independent of the mesh size as well as the time step. The comparison between experimental data and this simulation showed that our code can give a more accurate simulation on the motion of particles than previous DLM code. The code was then applied to simulate the sedimentation of 600 particles in a rectangular box. The falling course is presented and discussed. At the same time, this simulation also demonstrates that the method presented in this paper can be used for solving the initial problems involving a lager number of particles exactly with computing durations kept at acceptable levels.展开更多
.In this paper,an augmented Lagrangian Uzawa iterative method is developed and analyzed for solving a class of double saddle-point systems with semidefinite(2,2)block.Convergence of the iterativemethod is proved under....In this paper,an augmented Lagrangian Uzawa iterative method is developed and analyzed for solving a class of double saddle-point systems with semidefinite(2,2)block.Convergence of the iterativemethod is proved under the assumption that the double saddle-point problem exists a unique solution.An application of the iterative method to the double saddle-point systems arising from the distributed Lagrange multiplier/fictitious domain(DLM/FD)finite element method for solving elliptic interface problems is also presented,in which the existence and uniqueness of the double saddle-point system is guaranteed by the analysis of the DLM/FD finite element method.Numerical experiments are conducted to validate the theoretical results and to study the performance of the proposed iterative method.展开更多
In this paper, a new method-Distributed La-grange Multiplier/fictitiousdomain (DLM) method for partic-ulate flows was improved. A rectangular mesh was introduced todiscretize the domain, and the buoyant force was cons...In this paper, a new method-Distributed La-grange Multiplier/fictitiousdomain (DLM) method for partic-ulate flows was improved. A rectangular mesh was introduced todiscretize the domain, and the buoyant force was considered for predicting the positions ofparticles. In order to validate the presented algorithm, the sedimentation of a single circularparticle was simulated using different mesh sizes and time steps firstly. The results show that thesimulation is independent of the mesh size as well as the time step. And then, the interactionsbetween two falling particles, including drafting, kissing and tumbling, and the sedimentation of 18particles also have been simulated with the code.展开更多
The sedimentation of circular particles in a vertical channel filled withOldroyd ― B fluid was studied by an improved Distributed Lagrange Multiplier/fictitious domain(DLM) method. The sedimenting behaviors of two pa...The sedimentation of circular particles in a vertical channel filled withOldroyd ― B fluid was studied by an improved Distributed Lagrange Multiplier/fictitious domain(DLM) method. The sedimenting behaviors of two particles are presented firstly, which shows that,when the particles are dropped in a viscoealstic fluid, the stable configuration is the one wherethe particles are aligned parallel to the flow direction when the Mach number Mis less than 1 andthe elasticity number E is greater than 1. This agrees well with the known experimental in Ref. [1]and simulation results in Ref. [2]. Our simulations also show that, as in Newtonian fluid, thesedimentation of the particles will be accelerated due to the .interaction between particles in aviscoealstic fluid.展开更多
Distributed Lagrange Multiplier/Fictitious Domain (DLM/FD) method was employed to simulate the channel flow with polymer suspension. The polymer molecules were modeled as Finitely Extensible Nonlinear Elastic (FENE...Distributed Lagrange Multiplier/Fictitious Domain (DLM/FD) method was employed to simulate the channel flow with polymer suspension. The polymer molecules were modeled as Finitely Extensible Nonlinear Elastic (FENE) chains and FENE chain cluster. The coiled stretched transition of FENE chains was examined and the change in configuration of FENE chains was presented. The average velocity profile of the fully developed channel flow with 64 FENE chains was given and fitted well with the power-law curve. The change of chain cluster configuration was also simulated. These simulations show that DLM/FD method is capable of simulating the motion of not only FENE chain, but also FENE chain cluster.展开更多
文摘The Lagrange multiplier method plays an important role in establishing generalized variational principles notonly in tluid mechallics. but also in elasticity. Sometimes, however, one may come across variational crisis(somemultipliers vanish identically). failing to achieve his aim. The crisis is caused by the fact that the Inultipliers are treatedas independent variables in the process of variatioll. but after identification they become functions of the originalindependent variables. To overcome it, a Inodified Lagrange multiplier method or semi-inverse method has beenproposed to deduce generalized varistional principles. Some e-camples are given to illustrate its convenience andeffectiveness of the novel method.
基金This study is supported by the National Natural Science Foundation of China (Grant No50579046) the Science Foundation of Tianjin Municipal Commission of Science and Technology (Grant No043114711)
文摘This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as "fictitious" fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure' s rigid body shape and behavior, a rigid body constraint is enforced on the "fictitious" fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain (DLM/ FD) method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.
基金supported by the China Postdoctoral Science Foundation Funded Project (20080430038) the Funding Project for Academic Human Resources Development in Institutions of Higher Learning Under the Jurisdiction of Beijing Municipality (05004999200602)
文摘Using the concept of the base forces, a new finite element method (base force element method, BFEM) based on the complementary energy principle is presented for accurate modeling of structures with large displacements and large rotations. First, the complementary energy of an element is described by taking the base forces as state variables, and is then separated into deformation and rotation parts for the case of large deformation. Second, the control equations of the BFEM based on the complementary energy principle are derived using the Lagrange multiplier method. Nonlinear procedure of the BFEM is then developed. Finally, several examples are analyzed to illustrate the reliability and accuracy of the BFEM.
文摘By redefining the multiplier associated with inequality constraint as a positive definite function of the originally-defined multiplier, say, u2_i, i=1, 2, ..., m, nonnegative constraints imposed on inequality constraints in Karush-Kuhn-Tucker necessary conditions are removed. For constructing the Lagrange neural network and Lagrange multiplier method, it is no longer necessary to convert inequality constraints into equality constraints by slack variables in order to reuse those results dedicated to equality constraints, and they can be similarly proved with minor modification. Utilizing this technique, a new type of Lagrange neural network and a new type of Lagrange multiplier method are devised, which both handle inequality constraints directly. Also, their stability and convergence are analyzed rigorously.
基金Supported by the Fundamental Research Funds for the Central Universities of China (2009QL05)
文摘Roadways excavated in soft rocks at great depth are difficult to be maintained due to large deformation of surrounding rocks, which greatly influences the safety and efficiency of deep resources exploitation. During the excavation process of a deep soft rock tunnel, the rock wall may be compacted due to large deformation. In this paper, the technique to address this problem by a two-dimensional (2D) finite element software, large deformation engineering analyses software (LDEAS 1.0), is provided. By using the Lagrange multiplier method, the kinematic constraint of non-penetrating condition and static constraint of Coulomb friction are introduced to the governing equations in the form of incremental displacement. The numerical example demonstrates the efficiency of this technology. Deformations of a transportation tunnel in inclined soft rock strata at the depth of 1 000 m in Qishan coal mine and a tunnel excavated to three different depths are analyzed by two models, i.e. the additive decomposition model and polar decomposition model. It can be found that the deformation of the transportation tunnel is asymmetrical due to the inclination of rock strata. For extremely soft rock, large deformation can converge only for the additive decomposition model. The deformation of surrounding rocks increases with the increase in the tunnel depth for both models. At the same depth, the deformation calculated by the additive decomposition model is smaller than that by the polar decomposition model.
基金Project supported by the National Natural Science Foundation of China (Grant No. 60674059)Research Fund of University of Science and Technology Beijing, China (Grant No. 00009010)
文摘A chaotic system is bounded, and its trajectory is confined to a certain region which is called the chaotic attractor. No matter how unstable the interior of the system is, the trajectory never exceeds the chaotic attractor. In the present paper, the sphere bound of the generalized Lorenz system is given, based on the Lyapunov function and the Lagrange multiplier method. Furthermore, we show the actual parameters and perform numerical simulations.
文摘In this paper, by applying Lagrange, multiplier method and high order Lagrange multiplier method [1], we systematically derive coupled potential energy principle.coupled complementary energy principle,and generalized coupled potential energy principles and generalized coupled complementary energy principles with two and three kinds of variables in photoelasticity.
基金State Climbing Project (95-S-05-02) and State Natural Science Foundation of China (49724232).
文摘The mantle unsteady flows, which are in an incompressible and isoviscous spherical shell, are investigated by using algorithms of the parallel Lagrange multiplier dissonant decomposition method (LMDDM) and the parallel Lagrange multiplier discontinuous deformation analyses (LMDDA) in this paper. Some physical fields about mantle flows such as velocity, pressure, temperature, stress and the force to the crust of the Asian continent are calculated on a parallel computer.
文摘An improved implementation of Distributed Lagrange multiplier/fictitious domain method was presented and used to simulate the interactions between two circular particles sedimenting in a two_dimensional channel. The simulation results were verified by comparison with experiments. The results show that the interactions between two particles with different sizes can be described as drafting, kissing, tumbling and separating. Only for small diameter ratio, the two particles will interact undergoing repeated DKT (Drafting, Kissing and Tumbling) process. Otherwise, the two particles will separate after their tumbling. The results also show that, during the interaction process, the motion of the small particle is strongly affected while the large particle is affected slightly.
基金TheNationalNaturalSciencesFoundationforOutstandingYouthofChina (No .19925210)andZhejiangProvincialNaturalScienceFoundationofChina(No .10 10 4 7)
文摘An improved implementation of distributed multiplier/fictitious domain method is presented for the direct numerical simulation of particulate flow. The key improvement is to replace a finite element triangulation for the velocity and a “twice coarser' triangulation for the pressure with a rectangular discretization for the velocity and pressure. For code validation, the sedimentation of a single particle in a two dimensional channel was simulated. The results showed that the simulation is independent of the mesh size as well as the time step. The comparison between experimental data and this simulation showed that our code can give a more accurate simulation on the motion of particles than previous DLM code. The code was then applied to simulate the sedimentation of 600 particles in a rectangular box. The falling course is presented and discussed. At the same time, this simulation also demonstrates that the method presented in this paper can be used for solving the initial problems involving a lager number of particles exactly with computing durations kept at acceptable levels.
基金supported by the 10 plus 10 project of Tongji University(No.4260141304/004/010).
文摘.In this paper,an augmented Lagrangian Uzawa iterative method is developed and analyzed for solving a class of double saddle-point systems with semidefinite(2,2)block.Convergence of the iterativemethod is proved under the assumption that the double saddle-point problem exists a unique solution.An application of the iterative method to the double saddle-point systems arising from the distributed Lagrange multiplier/fictitious domain(DLM/FD)finite element method for solving elliptic interface problems is also presented,in which the existence and uniqueness of the double saddle-point system is guaranteed by the analysis of the DLM/FD finite element method.Numerical experiments are conducted to validate the theoretical results and to study the performance of the proposed iterative method.
文摘In this paper, a new method-Distributed La-grange Multiplier/fictitiousdomain (DLM) method for partic-ulate flows was improved. A rectangular mesh was introduced todiscretize the domain, and the buoyant force was considered for predicting the positions ofparticles. In order to validate the presented algorithm, the sedimentation of a single circularparticle was simulated using different mesh sizes and time steps firstly. The results show that thesimulation is independent of the mesh size as well as the time step. And then, the interactionsbetween two falling particles, including drafting, kissing and tumbling, and the sedimentation of 18particles also have been simulated with the code.
文摘The sedimentation of circular particles in a vertical channel filled withOldroyd ― B fluid was studied by an improved Distributed Lagrange Multiplier/fictitious domain(DLM) method. The sedimenting behaviors of two particles are presented firstly, which shows that,when the particles are dropped in a viscoealstic fluid, the stable configuration is the one wherethe particles are aligned parallel to the flow direction when the Mach number Mis less than 1 andthe elasticity number E is greater than 1. This agrees well with the known experimental in Ref. [1]and simulation results in Ref. [2]. Our simulations also show that, as in Newtonian fluid, thesedimentation of the particles will be accelerated due to the .interaction between particles in aviscoealstic fluid.
文摘Distributed Lagrange Multiplier/Fictitious Domain (DLM/FD) method was employed to simulate the channel flow with polymer suspension. The polymer molecules were modeled as Finitely Extensible Nonlinear Elastic (FENE) chains and FENE chain cluster. The coiled stretched transition of FENE chains was examined and the change in configuration of FENE chains was presented. The average velocity profile of the fully developed channel flow with 64 FENE chains was given and fitted well with the power-law curve. The change of chain cluster configuration was also simulated. These simulations show that DLM/FD method is capable of simulating the motion of not only FENE chain, but also FENE chain cluster.