Submitted:
16 December 2023
Posted:
18 December 2023
Read the latest preprint version here
Abstract
Keywords:
1. Introduction
2. A general duality principle non-convex optimization
4. A convex dual variational formulation for a third similar model
4.1. The algorithm through which we have obtained the numerical results
7. An exact convex dual variational formulation for a non-convex primal one
11. A duality principle for a general vectorial case in the calculus of variations
12. A note on the Galerkin Functional
13. A note on the Legendre-Galerkin functional
13.1. Numerical examples
14. A general concave dual variational formulation for global optimization
16. One more dual variational formulation
17. A model in superconductivity through an eigenvalue approach
18. A simplified qualitative many body model for the hydrogen nuclear fusion
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- 8.
- 1.
- 2.
- 3.
- 4.
- 5.
19. A more detailed mathematical description of the hydrogen nuclear fusion
- 1.
- For the Deuterium field
- 2.
- For the Tritium field
- 3.
- For the Helium field
- 4.
- For the Neutron field
- 5.
- For the electronic field resulting from the ionization
- 1.
- 2.
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- so that
- 8.
- 9.
- 10.
- 11.
- 12.
- 2.
- Mass conservation equation:
- 3.
-
Energy equationwhere we assume the Fourier lawwhere is the scalar field of temperature.Also,and
- 4.
- for an appropriate scalar function .
- 5.
- Mass relations
- (a)
- (b)
- (c)
- (d)
- (e)
- (f)
- 6.
- Other mass constraints
- (a)
- (b)
- (c)
- (d)
- (e)
- 7.
- For the induced electric field, we must havewhere and are appropriate real constants related to the respective charges.
- 8.
- A Maxwell equation:where
- 9.
- Another Maxwell equation:where the total electric field stands forand where generically denotingwe have also
- 1.
- For the Deuterium field
- 2.
- For the Tritium field
- 3.
- For the Helium field
- 4.
- For the Neutron field
- 5.
- For the electronic field resulting from the ionization
- 1.
- 2.
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- 8.
- 9.
- so that
- 10.
- 11.
- 12.
20. A final mathematical description of the hydrogen nuclear fusion
- 1.
- For a single Deuterium atom indexed by s:
- 2.
- For a single Tritium atom indexed by s:
- 3.
- For a single Helium atom indexed by s:
- 4.
- For the Neutron field:
- 5.
- For the electronic field resulting from the ionization
- 1.
- 2.
- (a)
- (b)
- (c)
- (d)
- (e)
- (f)
- (g)
- so that
- (h)
- (i)
- 6.
- Other mass constraints
- (a)
- (b)
- (c)
- (d)
- (e)
- (f)
- (g)
- (h)
- 7.
- 8.
- For the induced electric field, we must havewhere and are appropriate real constants related to the respective charges.
- 9.
- A Maxwell equation:where
- 10.
- Another Maxwell equation:where the total electric field stands forand where generically denotingwe have also
- 1.
- For the Deuterium field
- 2.
- For the Tritium field
- 3.
- For the Helium field
- 4.
- For the Neutron field
- 5.
- For the electronic field resulting from the ionization
- 1.
- 2.
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- 8.
- 9.
- so that
- 10.
- 11.
- 12.
21. A qualitative modeling for a general phase transition process
22. A mathematical description of a hydrogen molecule in a quantum mechanics context
- 1.
- 2.
- 3.
- (1)
- : mass of electron in the atom , where
- (2)
- : mass of proton in the atom , where
23. A mathematical model for the water hydrolysis
- 1.
- molecule generically corresponds to wave function .
- 2.
- molecule corresponds to wave function
- 3.
- hydrogen atom corresponds to wave function
- (1)
- For the water density (for charges), denoted by , we havewhere is the mass of a single water molecule and generically refers to the hydrogen proton at the hydrogen atom concerning the molecular density and so on.
- 2.
- For the density, denoted by , we havewhere is the mass of a single molecule of .
- 3.
- For the ionized hydrogen atom have
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- 8.
- 9.
- 10.
- 11.
- 12.
- 13.
- 14.
- 15.
- 16.
- 17.
24. A mathematical model for the Austenite and Martensite phase transition
- 1.
- For the Austenite phase:
- (a)
- (b)
- (c)
- (d)
- (e)
- (f)
- (g)
- (h)
- 2.
- For the Martensite phase:
- (a)
- (b)
- (c)
- (d)
- (e)
- (f)
- (g)
- (h)
- 3.
- For the total (iron) mass,
- 4.
- For the total Carbon mass
25. A note on classical free fields through a variational perspective
25.1. The Angular-Momentum tensor
25.2. A note on the solution of the Klein-Gordon equation
25.3. A note on the Dirac equation
26. A note on quantum field operators
11.1. An application concerning the harmonic oscillator operator in quantum mechanics
28. The generalized method of lines applied to fourth order differential equations
28.1. A numerical example
29. A note on hyper-finite differences for the generalized method of lines
30. Applications to the optimal shape design for a beam model
- 1.
- Set and
- 2.
- Calculate solution of equationwhere
- 3.
- Calculate such thatwhere
- 4.
- Set and go to step 2 until an appropriate convergence criterion is satisfied.
- (1)
-
clear allglobal m8 d d2wo H e1 ho h1 xo b5m8=100;d=1.0/m8;b5=0.1;e1=210*;ho=0.18;A=zeros(m8-1,m8-1);for i=1:m8-1A(1,i)=1.0;xo(i,1)=0.55;x3(i,1)=0.55;end;lb=0.4*ones(m8-1,1);ub=ones(m8-1,1);b=zeros(m8-1,1);b(1,1)=0.65*(m8-1);for i=1:m8f(i,1)=1.0;L(i,1)=1/2;P(i,1)=36.0*;end;i=1;m12=2;m50(i)=1/m12;z(i)=1/m50(i)*(-P(i,1)*);for i=2:m8-1m12=2-m50(i-1);m50(i)=1/m12;z(i)=m50(i)*(-P(i,1)*+z(i-1));end;v(m8,1)=0;for i=1:m8-1v(m8-i,1)=m50(m8-i)*v(m8-i+1,1)+z(m8-i);end;k=1;b12=1.0;while andkk=k+1;for i=1:m8-1H(i,1)=b5*/12*e1;f1(i,1)=v(i,1)/H(i,1);end;i=1;m12=2;m70(i)=1/m12;z1(i)=m70(i)*(-f1(i,1)*);for i=2:m8-1m12=2-m70(i-1);m70(i)=1/m12;z1(i)=m70(i)*(-f1(i,1)*+z1(i-1));end;w(m8,1)=0;for i=1:m8-1w(m8-i,1)=m70(m8-i)*w(m8-i+1,1)+z1(m8-i);end;d2wo(1,1)=(-2*w(1,1)+w(2,1))/;for i=2:m8-1d2wo(i,1)=(w(i+1,1)-2*w(i,1)+w(i-1,1))/;end;k9=1;b14=1.0;whilek9k9=k9+1;X=fmincon(’beamNov2023’,xo,A,b,,lb,ub);b14=max(abs(xo-X))xo=X;end;b12=max(abs(xo-x3))x3=xo;for i=1:m8-1L(i,1)=xo(i,1);end;end;***************
- (1)
-
function S=beamNov2023(x)global m8 d d2wo H e1 ho h1 xo b5S=0;for i=1:m8-1S=S+1//b5/e1*(H(i,1)**12;end;*****************************
- (1)
-
clear allglobal m8 d d2wo H e1 ho h1 xo b5m8=100;d=1.0/m8;b5=0.1;e1=210*;ho=0.18;A=zeros(m8-1,m8-1);for i=1:m8-1A(1,i)=1.0;xo(i,1)=0.55;x3(i,1)=0.55;end;lb=0.4*ones(m8-1,1);ub=ones(m8-1,1);b=zeros(m8-1,1);b(1,1)=0.65*(m8-1);for i=1:m8f(i,1)=1.0;L(i,1)=1/2;P(i,1)=36.0*;end;i=1;m12=2;m50(i)=1/m12;z(i)=1/m50(i)*(-P(i,1)*);for i=2:m8-1m12=2-m50(i-1);m50(i)=1/m12;z(i)=m50(i)*(-P(i,1)*+z(i-1));end;v(m8,1)=0;for i=1:m8-1v(m8-i,1)=m50(m8-i)*v(m8-i+1,1)+z(m8-i);end;k=1;b12=1.0;whilekk=k+1;for i=1:m8-1H(i,1)=b5*/12*e1;f1(i,1)=v(i,1)/H(i,1);f2(i,1)=i*d/H(i,1);f3(i,1)=1/H(i,1);end;i=1;m12=2;m70(i)=1/m12;z1(i)=m70(i)*(-f1(i,1)*);z2(i)=m70(i)*(-f2(i,1)*);z3(i)=m70(i)*(-f3(i,1)*);for i=2:m8-1m12=2-m70(i-1);m70(i)=1/m12;z1(i)=m70(i)*(-f1(i,1)*+z1(i-1));z2(i)=m70(i)*(-f2(i,1)*+z2(i-1));z3(i)=m70(i)*(-f3(i,1)*+z3(i-1));end;w1(m8,1)=0;w2(m8,1)=0;w3(m8,1)=0;for i=1:m8-1w1(m8-i,1)=m70(m8-i)*w1(m8-i+1,1)+z1(m8-i);w2(m8-i,1)=m70(m8-i)*w2(m8-i+1,1)+z2(m8-i);w3(m8-i,1)=m70(m8-i)*w3(m8-i+1,1)+z3(m8-i);end;m3(1,1)=w2(1,1);m3(1,2)=w3(1,1);m3(2,1)=w2(m8-1,1);m3(2,2)=w3(m8-1,1);h3(1,1)=-w1(1,1);h3(2,1)=-w1(m8-1,1);h5(:,1)=inv(m3)*h3;for i=1:m8wo(i,1)=w1(i,1)+h5(1,1)*w2(i,1)+h5(2,1)*w3(i,1);end;d2wo(1,1)=(-2*wo(1,1)+wo(2,1))/;for i=2:m8-1d2wo(i,1)=(wo(i+1,1)-2*wo(i,1)+wo(i-1,1))/;end;k9=1;b14=1.0;whilek9k9=k9+1;X=fmincon(’beamNov2023’,xo,A,b,,lb,ub);b14=max(abs(xo-X))xo=X;end;b12=max(abs(xo-x3))x3=xo;for i=1:m8-1L(i,1)=xo(i,1);end;end;*****************************
31. Applications to the optimal shape design for a plate model
- (1)
- Set and
- (2)
- Calculate solution of equationwhere
- (3)
- Calculate such thatwhere
- (4)
- Set and go to step 2 until an appropriate convergence criterion is satisfied.

- (1)
-
clear allglobal m8 d d2xwo d2ywo H e1 ho xo b5m8=40;d=1.0/m8;w5=0.3;e1=200*/;ho=0.12;A=zeros();for i=1:A(1,i)=1.0;xo(i,1)=0.55;x3(i,1)=0.55;end;lb=0.45*ones(,1);ub=ones(,1);b=zeros(,1);b(1,1)=0.75*for i=1:(m8-1)for j=1:m8-1f(i,j,1)=1.0;L(i,j,1)=1/2;P(i,j,1)=2*; end;end;for i=1:m8wo(:,i)=0.001*ones(m8-1,1);end;m2=zeros(m8-1,m8-1);for i=2:m8-2m2(i,i)=-2.0;m2(i,i-1)=1.0;m2(i,i+1)=1.0;end;m2(1,1)=-2.0;m2(1,2)=1.0;m2(m8-1,m8-1)=-2.0;m2(m8-1,m8-2)=1.0;Id=eye(m8-1);i=1;m12=2*Id-m2*; m50(:,:,i)=inv(m12);z(:,i)=m50(:,:,i)*(-P(:,i,1)*);for i=2:m8-1m12=2*Id-m2*-m50(:,:,i-1);m50(:,:,i)=inv(m12);z(:,i)=m50(:,:,i)*(-P(:,i,1)*+z(:,i-1));end; v(:,m8)=zeros(m8-1,1);for i=1:m8-1v(:,m8-i)=m50(:,:,m8-i)*v(:,m8-i+1)+z(:,m8-i);end;k=1;b12=1.0;while () and ()kk=k+1;for i=1:m8-1for j=1:m8-1H(j,i,1)=/12*e1;f1(j,i,1)=v(j,i)/H(j,i,1);end;end;i=1;m12=2*Id-m2*;m70(:,:,i)=inv(m12);z1(:,i)=m70(:,:,i)*(-f1(:,i,1)*);for i=2:m8-1m12=2*Id-m2*-m70(:,:,i-1);m70(:,:,i)=inv(m12);z1(:,i)=m70(:,:,i)*(-f1(:,i,1)*+z1(:,i-1));end;w(:,m8)=zeros(m8-1,1);for i=1:m8-1w(:,m8-i)=m70(:,:,m8-i)*w(:,m8-i+1)+z1(:,m8-i);end;d2xwo(:,1)=(-2*w(:,1)+w(:,2))/;for i=2:m8-1d2xwo(:,i)=(w(:,i+1)-2*w(:,i)+w(:,i-1))/;end;for i=1:m8-1d2ywo(:,i)=m2*w(:,i)/;end;k9=1; b14=1.0;while () and ()k9k9=k9+1;X=fmincon(’beamNov2023A3’,xo,A,b,,lb,ub);b14=max(abs(xo-X))xo=X;end;b12=max(max(abs(w-wo)))wo=w;x3=xo;for i=1:m8-1for j=1:m8-1L(j,i,1)=xo((i-1)*(m8-1)+j,1);end;end;end;for i=1:m8-1x8(i,1)=i*d;end;mesh(x8,x8,L);*********************
- (1)
-
function S=beamNov2023A3(x)global m8 d d2xwo d2ywo H e1 ho xo b5S=0;for i=1:m8-1for j=1:m8-1x1(j,i)=x((m8-1)*(i-1)+j,1);end;end;for i=1:m8-1for j=1:m8-1S=S+;end;end;********************************
32. A note on the first Maxwell equation of electromagnetism
33. A note on relaxation for a general model in the vectorial calculus of variations
18.1. Some related numerical results
- (1)
-
clear allglobal m8 d u e3m8=100;d=1/m8;e3=0.0005;for i=1:2*m8+1xo(i,1)=0.36;end;b12=1.0;k=1;whilekk=k+1;X=fminunc(’funDecember2023’,xo);b12=max(abs(xo-X))xo=X;u(m8/2)end;for i=1:m8x(i,1)=i*d;end;plot(x,u);***********************
- (1)
-
function S=funDecember2023(x)global m8 d u e3for i=1:m8u(i,1)=x(i,1);v(i,1)=x(i+m8,1);yo(i,1)=sin(pi*i*d)/2;end;L=(1+sin(x(2*m8+1,1)))/2;u(m8,1)=1/2;v(m8,1)=0.0;du(1,1)=u(1,1)/d;dv(1,1)=v(1,1)/d;for i=2:m8du(i,1)=(u(i,1)-u(i-1,1))/d;dv(i,1)=(v(i,1)-v(i-1,1))/d;end;d2u(1,1)=(-2*u(1,1)+u(2,1))/;for i=2:m8-1d2u(i,1)=(u(i-1,1)-2*u(i,1)+u(i+1,1))/;end;S=0;for i=1:m8S=S+;S=S+;S=S+;end;for i=1:m8-1S=S+e3*;end;*******************
18.2. A related duality principle and concerning convex dual formulation
34. A general convex primal dual formulation for an originally non-convex primal one
35. A note on the special relativistic physics
36. Conclusion
Data Availability Statement
Conflicts of Interest
References
- W.R. Bielski, A. W.R. Bielski, A. Galka, J.J. Telega, The Complementary Energy Principle and Duality for Geometrically Nonlinear Elastic Shells. I. Simple case of moderate rotations around a tangent to the middle surface. Bulletin of the Polish Academy of Sciences, Technical Sciences, Vol. 38, No. 7-9, 1988.
- W.R. Bielski and J.J. Telega,A Contribution to Contact Problems for a Class of Solids and Structures, Arch. Mech., 37, 4-5, pp. 303-320, Warszawa 1985.
- J.J. Telega, On the complementary energy principle in non-linear elasticity. Part I: Von Karman plates and three dimensional solids, C.R. Acad. Sci. Paris, Serie II, 308, 1193-1198; Part II: Linear elastic solid and non-convex boundary condition. Minimax approach, ibid, pp. 1313.
- A.Galka and J.J.Telega, Duality and the complementary energy principle for a class of geometrically non-linear structures. Part I. Five parameter shell model; Part II. Anomalous dual variational priciples for compressed elastic beams, Arch. Mech. 47 (1995) 677-698, 699-724.
- J.F. Toland, A duality principle for non-convex optimisation and the calculus of variations, Arch. Rat. Mech. Anal., 71, No. 1 (1979), 41-61.
- R.A. Adams and J.F. Fournier, Sobolev Spaces, 2nd edn. (Elsevier, New York, 2003).
- F. Botelho, Functional Analysis and Applied Optimization in Banach Spaces, Springer Switzerland, 2014.
- F.S. Botelho, Functional Analysis, Calculus of Variations and Numerical Methods in Physics and Engineering, CRC Taylor and Francis, Florida, 2020.
- F.S. Botelho, Variational Convex Analysis, Ph.D. thesis, Virginia Tech, Blacksburg, VA -USA, (2009).
- F. Botelho, Topics on Functional Analysis, Calculus of Variations and Duality, Academic Publications, Sofia, (2011).
- F. Botelho, Existence of solution for the Ginzburg-Landau system, a related optimal control problem and its computation by the generalized method of lines, Applied Mathematics and Computation, 218, 11976-11989, (2012).
- R.T. Rockafellar, Convex Analysis, Princeton Univ. Press, (1970).
- F.S. Dual Variational Formulations for a Large Class of Non-Convex Models in the Calculus of Variations, 2023. [CrossRef]
- J.F. Annet, Superconductivity, Superfluids and Condensates, 2nd edn. 2010.
- L.D. Landau and E.M. Lifschits, Course of Theoretical Physics, Vol. 5- Statistical Physics, part 1. (Butterworth-Heinemann, Elsevier, reprint 2008).
- F.S. Botelho, Advanced Calculus and its Applications in Variational Quantum Mechanics and Relativity Theory, CRC Taylor and Francis, Florida, 2021.
- P.Ciarlet, Mathematical Elasticity, Vol. II – Theory of Plates, North Holland Elsevier (1997).
- H. Attouch, G. H. Attouch, G. Buttazzo and G. Michaille, Variational Analysis in Sobolev and BV Spaces, MPS-SIAM Series in Optimization, Philadelphia, 2006.
- J.C. Strikwerda, Finite Difference Schemes and Partial Differential Equations, SIAM, second edition (Philadelphia, 2004).
- N.N. Bogoliubov and D.V. Shirkov, Introduction of the Theory of Quantized Fields, Monographs and Texts in Physics and Astronomy, Vol. III, Intercience Publishers INC., New York, USA, 1959.
- G.B. Folland, Quantum Field Theory, A Tourist Guide for Mathematicians, Mathematical Surveys and Monographs, Vol. 149, Americam Mathematical Society, AMS, Providence, Rhode Island, USA, 2008.



























Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).