Submitted:
23 November 2023
Posted:
23 November 2023
Read the latest preprint version here
Abstract
Keywords:
MSC: 49N15; 35A15; 49J40
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
- clear all
- (number of nodes)
-
(we have fixed the number of iterations)
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
19. A more detailed mathematical description of the hydrogen nuclear fusion
- For the Deuterium field
- For the Tritium field
- For the Helium field
- For the Neutron field
- For the electronic field resulting from the ionization
-
so that
-
The momentum conservation equation for the fluid motionHere is the total density and P is the fluid pressure field.Furthermore,and
- Mass conservation equation:
-
Energy equationwhere we assume the Fourier lawwhere is the scalar field of temperature.Also,and
-
for an appropriate scalar function .
-
Mass relations
- (a)
- (b)
- (c)
- (d)
- (e)
- (f)
- (g)
- (h)
-
so that
where,- (a)
- (b)
- (c)
- (d)
- (e)
- (f)
-
Other mass constraints
- (a)
- (b)
- (c)
- (d)
- (e)
- For the induced electric field, we must havewhere and are appropriate real constants related to the respective charges.
- A Maxwell equation:where
- Another Maxwell equation:where the total electric field stands forand where generically denotingwe have also
- For the Deuterium field
- For the Tritium field
- For the Helium field
- For the Neutron field
- For the electronic field resulting from the ionization
-
so that
20. A final mathematical description of the hydrogen nuclear fusion
- For a single Deuterium atom indexed by s:
- For a single Tritium atom indexed by s:
- For a single Helium atom indexed by s:
- For the Neutron field:
- For the electronic field resulting from the ionization
-
so that
- Here denotes the outward normal vectorial fields to the concerning surfaces.
-
The momentum conservation equation for the fluid motionHere is the total density and P is the fluid pressure field.Furthermore,and
- Mass conservation equation:
-
Energy equationwhere we assume the Fourier lawwhere is the scalar field of temperature.Also,and
-
for an appropriate scalar function .
-
Mass relations
- (a)
- (b)
- (c)
- (d)
- (e)
where,- (a)
- (b)
- (c)
- (d)
- (e)
- (f)
- (g)
-
so that
- (h)
- (i)
-
Other mass constraints
- (a)
- (b)
- (c)
- (d)
- (e)
- (f)
- (g)
- (h)
- For the induced electric field, we must havewhere and are appropriate real constants related to the respective charges.
- A Maxwell equation:where
- Another Maxwell equation:where the total electric field stands forand where generically denotingwe have also
- For the Deuterium field
- For the Tritium field
- For the Helium field
- For the Neutron field
- For the electronic field resulting from the ionization
-
so that
21. A qualitative modeling for a general phase transition process
-
We start with corresponding to and in .
-
We end the process with corresponding to and in .
22. A mathematical description of a hydrogen molecule in a quantum mechanics context
- : mass of electron in the atom , where
- : mass of proton in the atom , where
- From the proton in the atom :
- For the proton in the atom :
- For the atom :
- For the atom :
- For the electrons and , concerning the physical electronic link between the atoms:
- For the total molecular density:
23. A mathematical model for the water hydrolysis
- molecule generically corresponds to wave function .
- molecule corresponds to wave function
- hydrogen atom corresponds to wave function
- 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.
- For the density, denoted by , we havewhere is the mass of a single molecule of .
- For the ionized hydrogen atom have
24. A mathematical model for the Austenite and Martensite phase transition
-
For the Austenite phase:
- (a)
- (b)
- (c)
- (d)
- (e)
- (f)
- (g)
- (h)
-
For the Martensite phase:
- (a)
- (b)
- (c)
- (d)
- (e)
- (f)
- (g)
- (h)
- For the total (iron) mass,
- For the total Carbon mass
25. A note on classical free fields through a variational perspective
25.1. The Angular-Momentum tensor
26. A note on quantum field operators
26.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
-
clear allm8=100;d=1/m8;e1=1.0;for i=1:m8f(i,1)=1.0;end;a(1)=2/3;b(1)=-1/6;c(1)=f(1,1)*/(6e1);m12=(6-4*a(1));a(2)=(4*b(1)+4)/m12;b(2)=-1/m12;c(2)=1/m12*(4*c(1)+f(2,1)*/e1);for i=3:m8-2m12=(a(i-2)*a(i-1)+b(i-2)-4*a(i-1)+6);a(i)=-1/m12*(a(i-2)*b(i-1)-4*b(i-1)-4);b(i)=-1/m12;c(i)=1/m12*(f(i,1)*/e1-c(i-2)-a(i-2)*c(i-1)+4*c(i-1));end;u(m8,1)=0;u(m8-1,1)=0;for i=2:m8-1;u(m8-i,1)=a(m8-i)*u(m8-i+1,1)+b(m8-i)*u(m8-i+2,1)+c(m8-i);end;for i=1:m8x(i)=i*d;end;plot(x,u)******************
29. Applications to the optimal shape design for a beam model
- Set and
- Calculate solution of equationwhere
- Calculate such thatwhere
- Set and go to step 2 until an appropriate convergence criterion is satisfied.
-
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;***************
-
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;*****************************
-
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;*****************************
30. Applications to the optimal shape design for a plate model
- Set and
- Calculate solution of equationwhere
- Calculate such thatwhere
- Set and go to step 2 until an appropriate convergence criterion is satisfied.
-
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);*********************
-
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;********************************
31. A note on the first Maxwell equation of electromagnetism
32. Conclusion
Data Avaliability
Conflict of interest declaration
References
- Bielski, W.R.; Galka, A.; Telega, J.J. The Complementary Energy Principle and Duality for Geometrically Nonlinear Elastic Shells. I. Simple case of moderate rotations around a tangent to the middle surface. Bull. Pol. Acad. Sci. Tech. Sci. 1988, 38, 7–9. [Google Scholar]
- Bielski, W.R.; Telega, J.J. A Contribution to Contact Problems for a Class of Solids and Structures. Arch. Mech. 1985, 37, 303–320. [Google Scholar]
- Telega, J.J. On the complementary energy principle in non-linear elasticity. Part I: Von Karman plates and three dimensional solids. Part II: Linear elastic solid and non-convex boundary condition. Minimax approach, ibid. C.R. Acad. Sci. Paris Serie II 1989, 308, 1193–1198. [Google Scholar]
- Galka, A.; Telega, J.J. 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. 1995, 47, 677–698. [Google Scholar]
- Toland, J.F. A duality principle for non-convex optimisation and the calculus of variations. Arch. Rat. Mech. Anal. 1979, 71, 41–61. [Google Scholar] [CrossRef]
- Adams, R.A.; Fournier, J.F. Sobolev Spaces, 2nd ed.; Elsevier: New York, NY, USA, 2003. [Google Scholar]
- Botelho, F. Functional Analysis and Applied Optimization in Banach Spaces; Springer: Cham, Switzerland, 2014. [Google Scholar]
- Botelho, F.S. Functional Analysis, Calculus of Variations and Numerical Methods in Physics and Engineering; CRC Taylor and Francis: Boca Raton, FL, USA, 2020. [Google Scholar]
- Botelho, F.S. Variational Convex Analysis. Ph.D. Thesis, Virginia Tech, Blacksburg, VA, USA, 2009. [Google Scholar]
- Botelho, F. Topics on Functional Analysis, Calculus of Variations and Duality; Academic Publications: Sofia, Bulgaria, 2011. [Google Scholar]
- Botelho, F. Existence of solution for the Ginzburg-Landau system, a related optimal control problem and its computation by the generalized method of lines. Appl. Math. Comput. 2012, 218, 11976–11989. [Google Scholar] [CrossRef]
- Rockafellar, R.T. Convex Analysis; Princeton University Press: Princeton, NJ, USA, 1970. [Google Scholar]
- Botelho, F.S. Dual Variational Formulations for a Large Class of Non-Convex Models in the Calculus of Variations. Mathematics 2023, 11, 63. [Google Scholar] [CrossRef]
- Annet, J.F. Superconductivity, Superfluids and Condensates, 2nd ed.; Oxford Master Series in Condensed Matter Physics; Oxford University Press: Oxford, UK, 2010. [Google Scholar]
- Landau, L.D.; Lifschits, E.M. Course of Theoretical Physics, Vol. 5-Statistical Physics, Part 1; Elsevier Butterworth-Heinemann: Oxford, UK, 2008. [Google Scholar]
- Botelho, F.S. Advanced Calculus and its Applications in Variational Quantum Mechanics and Relativity Theory; CRC Taylor and Francis: Boca Raton, FL, USA, 2021. [Google Scholar]
- Ciarlet, P. Mathematical Elasticity, Vol. II—Theory of Plates; Elsevier: North Holland, The Netherlands, 1997. [Google Scholar]
- Attouch, H.; Buttazzo, G.; Michaille, G. Variational Analysis in Sobolev and BV Spaces; MPS-SIAM Series in Optimization; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2006. [Google Scholar]
- Strikwerda, J.C. Finite Difference Schemes and Partial Differential Equations, 2nd ed.; SIAM: Philadelphia, PA, USA, 2004. [Google Scholar]
- Bogoliubov, N.N.; Shirkov, D.V. Introduction of the Theory of Quantized Fields, Monographs and Texts in Physics and Astronomy; Intercience Publishers INC.: New York, NY, USA, 1959; Volume 3. [Google Scholar]
- Folland, G.B. Quantum Field Theory, A Tourist Guide for Mathematicians, Mathematical Surveys and Monographs; Americam Mathematical Society, AMS: Providence, RI, USA, 2008; Volume 149. [Google Scholar]

























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/).