A numerical comparison between two unconstrained variational formulations

In an effort to relieve the often cumbersome burden of meeting the requirements of the end conditions and to unify the solution formulation for boundary- and initial-value problems, unconstrained variational statements have been introduced in conjunction with some approximate methods. In the case of a boundary value problem, it is shown in this paper that two different variational statements can be established: one is arrived at by the use of the Lagrange multipliers, the other by energy considerations. The numerical convergence of the solutions associated with finite element schemes involving use of one of these two different variational statements is compared with that of the other. In the case of an initial value problem, both formulations can again be established when the adjoint field variable and the adjoint variational statement are introduced. The numerical data presented here indicate that while both methods generate excellent convergent results for the boundary problem, the method of stiff springs yields results which show much better convergence for the initial value problem than those achieved by Lagrange multipliers.     

自由界面等离子体平衡的变分原理

在自由界面情况下用求位能汛函极值的变分原理推出了等离子体平衡所需满足的平衡方程和边界条件;其泛函的Euler方程即为具有自由界面条件的磁面函数方程。也给出了便于计算的变分泛函,它相当于等值面边界问题。对于导体壁为简单几何形状的情况,用Rit法进行了数值计算。 关键词：     

Variational Principle in the Theory of Partial Distribution Functions

We have constructed thermodynamically consistent theory of correlations in a classical equilibrium system. The theory is based on a variational principle for thermodynamic potential, as functional of correlation functions. Both the thermodynamic potential and correlation functions are determined simultaneously from a uniform variation problem for this functional. We have offered also methods of consecutive approximations system construction for the solution of the obtained variational problem.     

Stability of a pretwisted tapered cantilever beam subjected to dissipative and follower forces

Stability of a pretwisted tapered cantilever beam of rectangular cross-section subjected to a follower force at its free end is investigated. The effects of internal and external damping are included in the study. The non-self-adjoint boundary value problem is formulated with the Euler-Bernoulli theory and an associated adjoint boundary value problem is introduced. Approximate values of the critical load are calculated on the basis of a suitable adjoint variational principle for several values of the geometric and material parameters of the beam. The results are shown in graphs.     

里兹方法在光波导理论中的应用

本文讨论了在三层平板介质波导中TE波的传播问题.三层平板的芯区由Kerr介质构成,其折射率平方(n²)的剖面分布随x作抛物分布变化.本文将应用变分法中的里兹方法,求出亥姆霍兹方程的第一级近似的解析解.     

Variational approximation methods for long-range force transmission in biopolymer gels

The variational principle of minimum free energy (MFEVP) has been widely used in research of soft matter statics. The MFEVP can be used not only to derive equilibrium equations (including both bulk equations and boundary conditions), but also to develop direct variational methods (such as Ritz method) to find approximate solutions to these equilibrium equations. We apply these variational methods to study long-range force transmission in nonlinear elastic biopolymer gels. It is shown that the slow decay of cell-induced displacements measured experimentally for fibroblast spheroids in three-dimensional fibrin gels can be well explained by variational approximations based on the three-chain model of biopolymer gels.     

Effect of angle of attack on the stability of a rotating non-uniform cantilever with a tip mass subjected to dissipative and non-conservative forces

The influence of angle of attack on the stability of a rotating viscoelastic tapered cantilever beam of rectangular cross-section carrying a tip mass and subjected to a circulatory force at its free end is investigated. The effect of external damping is included in the study. The non-self-adjoint boundary value problem is formulated with the Euler-Bernoulli hypothesis and an appropriate adjoint boundary value problem is introduced. Approximate values of the critical flutter load are calculated on the basis of an adjoint variational principle for several values of geometric and material parameters of the system. The results are presented through a series of graphs.     

Unconventional Hamilton-type variational principles for analytical mechanics

According to the basic idea of classical yin-yang complementarity and modern dual-complementarity, in a simple and unified new way proposed by Luo, the un-conventional Hamilton-type variational principles of holonomic conservative system in analytical mechanics can be established systematically. This unconventional Hamilton-type variational principle can fully characterize the initial-value problem of analytical mechanics, so that it is an important innovation for the Hamilton-type variational principle. In this paper, an important integral relation is given, which can be considered as the expression of the generalized principle of virtual work for analytical mechanics in mechanics. Based on this relation, it is possible not only to obtain the principle of virtual work of holonomic conservative system in analytical mechanics, but also to derive systematically the complementary functionals for three-field and two-field unconventional variational principles, and the functional for the one-field one by the generalized Legendre transformation given in this paper. Further, with this new approach, the intrinsic relationship among various principles can be explained clearly. Meanwhile, the unconventional Hamilton-type variational principles of nonholonomic conservative system in analytical mechanics can also be established systematically in this paper.     

Vibration of a rectangular plate with a central power-law profiled groove by the Rayleigh–Ritz method

The prediction of natural frequencies of rectangular plates where a profiled indentation is present is made using a Rayleigh–Ritz variational energy method. Panels with holes are often found for access cables and access gaps, and it is shown that the application of damping to the profile leads to a more efficient method of reducing vibration than covering a whole rectangular plate, which is advantageous where weight saving is necessary.     

基于二阶离散变分方法的非线性映射参数估计

提出一种估计非线性映射未知参数的二阶离散变分方法.首先针对非线性离散混沌系统, 利用变分方法导出了伴随方程和目标泛函梯度, 以此为基础利用二阶离散变分方法给出了二阶伴随方程和精确计算Hessian矩阵-向量乘积的显式表达式; 其次设计了估计非线性映射未知参数的新算法, 并以此对Hyperhenón映射和二维抛物映射中的未知参数进行了精确的估计. 数值仿真结果表明了该方法的有效性和优点. 关键词： 非线性映射 参数估计 二阶离散变分方法 伴随方程     

Unconventional Hamilton-type variational principles for electromagnetic elastodynamics

According to the basic idea of classical yin-yang complementarity and modern dual-complementarity, in a simple and unified new way proposed by Luo, the unconventional Hamilton-type variational principles for electromagnetic elastodynamics can be established systematically. This new variational principles can fully characterize the initial-boundary-value problem of this dynamics. In this paper, the expression of the generalized principle of virtual work for electromagnetic dynamics is given. Based on this equation, it is possible not only to obtain the principle of virtual work in electromagnetic dynamics, but also to derive systematically the complementary functionals for eleven-field, nine-field and six-field unconventional Hamilton-type variational principles for electromagnetic elastodynamics, and the potential energy functionals for four-field and three-field ones by the generalized Legendre transformation given in this paper. Furthermore, with this approach, the intrinsic relationship among various principles can be explained clearly.     

Two variational approaches for domain decomposition problems solved by SGBEM with nonconforming discretizations

Two original approaches for the solution of elastic boundary value problems with domain decomposition (DDBVP) using the symmetric Galerkin boundary element method (SGBEM) are presented. Each approach is based on a variational principle, a difference between them consisting in the treatment of the coupling conditions which connect the solutions through an interface. The computer codes developed are able to deal with curved interfaces in a domain decomposition problem discretized by nonmatching meshes of linear elements along the interfaces. The effectiveness of the methods is documented by a numerically solved example. The text was submitted by the authors in English.     

基于变分方法的混沌系统参数估计

提出一种基于变分原理的估计混沌系统未知参数的方法,对以x= F(x,θ) 为控制方程的所有混沌系统具有普适性.首先将混沌系统方程引入到目标泛函中;接着利用变分原理导出了混沌系统的伴随方程和待辨识参数泛函梯度的通用公式;然后设计了估计混沌系统未知参数的算法;最后对典型的Lorenz混沌系统和超混沌Chen系统的未知参数进行了估计.数值仿真结果表明该方法是一种非常有效的估计混沌系统未知参数的方法. 关键词： 混沌系统 参数估计 变分方法 伴随方程     

厄尔尼诺和南方涛动海气耦合模型中参数估计的变分方法

提出一种基于变分原理估计厄尔尼诺和南方涛动海气耦合模型中未知参数的方法. 首先将所研究的非线性海气耦合动力方程引入到目标泛函中; 接着利用变分方法导出伴随方程和待辨识参数泛函梯度的公式; 然后设计了估计未知参数的算法.数值试验结果表明变分方法是一 种能有效估计海气耦合非线性系统未知参数的方法.     

Conservation laws of nonconservative nonholonomic system based on Herglotz variational problem

The generalized variational principle of Herglotz type provides a variational method for describing nonconservative or dissipative processes. The purpose of this letter is to extend this variational principle to a first order linear nonholonomic system and study the conservation laws of the nonconservative nonholonomic system based on Herglotz variational problem. A new differential variational principle of the nonconservative nonholonomic system is proposed, which is based on Herglotz variational problem. And the differential equations of motion of the system are also obtained. Then, according to the condition for the invariance of the differential variational principle, the conservation theorem based on Herglotz variational problem for the nonconservative nonholonomic system are obtained. The theorem contains the conservation theorem of the nonconservative holonomic system as its special case, which can be reduced to the first Noether's theorem based on Herglotz variational problem under proper conditions. The inverse theorem of the conservation theorem is also provided and proved. An example is given to illustrate the application at the end of this letter.     

Improvement of variational approach in an interacting two-fermion system

A more reasonable trial ground state wave function is constructed for the relative motion of an interacting two-fermion system in a one-dimensional(1D) harmonic potential. At the boundaries both the wave function and its first derivative are continuous and the quasi-momentum is determined by a more practical constraint condition which associates two variational parameters. The upper bound of the ground state energy is obtained by applying the variational principle to the expectation value of the Hamiltonian of relative motion on the trial wave function. The resulted energy and wave function show better agreement with the analytical solution than the original proposal.     

Surface Tension of Multi-phase Flow with Multiple Junctions Governed by the Variational Principle

We explore a computational model of an incompressible fluid with a multi-phase field in three-dimensional Euclidean space. By investigating an incompressible fluid with a two-phase field geometrically, we reformulate the expression of the surface tension for the two-phase field found by Lafaurie et al. (J Comput Phys 113:134–147, 1994) as a variational problem related to an infinite dimensional Lie group, the volume-preserving diffeomorphism. The variational principle to the action integral with the surface energy reproduces their Euler equation of the two-phase field with the surface tension. Since the surface energy of multiple interfaces even with singularities is not difficult to be evaluated in general and the variational formulation works for every action integral, the new formulation enables us to extend their expression to that of a multi-phase (N-phase, N\geqslant2N\geqslant2) flow and to obtain a novel Euler equation with the surface tension of the multi-phase field. The obtained Euler equation governs the equation for motion of the multi-phase field with different surface tension coefficients without any difficulties for the singularities at multiple junctions. In other words, we unify the theory of multi-phase fields which express low dimensional interface geometry and the theory of the incompressible fluid dynamics on the infinite dimensional geometry as a variational problem. We apply the equation to the contact angle problems at triple junctions. We computed the fluid dynamics for a two-phase field with a wall numerically and show the numerical computational results that for given surface tension coefficients, the contact angles are generated by the surface tension as results of balances of the kinematic energy and the surface energy.     

Application of variational calculus to propagation of coupled pulses in optical fibers

Two nonlinear Schrödinger equations, linked by cross-modulation terms, are used to study the nature of coupled pulse propagation in an optical fiber. The problem is formulated within the framework of variational calculus. Using Gaussian trial functions for the propagating pulses, an expression is constructed for the effective Lagrangian of the system. It is shown that this Lagrangian, via the Ritz optimization procedure, provides a basis to construct approximate solutions of the problem. Some judicious approximations are invoked to investigate how the cross modulation affects the behavior of pulse propagation. Conditions are derived under which both pulses can propagate without distortion.     

Data assimilation for the Navier–Stokes- equations

The well-posedness of the data assimilation problem for the Navier–Stokes-α equations on a bounded three-dimensional domain is investigated. The data assimilation procedures under consideration are the adjoint method of variational data assimilation (4D-Var) and the method of continuous data assimilation. Concerning the adjoint method the existence of optimal initial conditions with respect to an observation-dependent cost functional is proven, the optimizers are characterized by a first-order necessary condition involving the adjoint linearized Navier–Stokes-α equations and conditions for the uniqueness of the initial conditions are given. Well-posedness of the continuous data assimilation problem is proven and convergence rates in terms of observational resolution are provided.