Submitted:
15 November 2023
Posted:
15 November 2023
Read the latest preprint version here
Abstract
Keywords:
MSC: 49N15; 35A15; 49J40
1. Introduction
2. A General Duality Principle Non-Convex Optimization
3. Another Duality Principle for a Simpler Related Model in Phase Transition with a Respective Numerical Example
- Set and and
- Choose such that and
- Set
- Calculate solution of the system of equations:andthat isandso thatand
- Calculate by solving the system of equations:andthat isand
- If , then stop, else set and go to item 4.
-
clear allm8=300;d=1/m8;K=0.1;K1=120;for i=1:m8vo(i,1)=i*d/10;yo(i,1)=sin(i*d*pi)/2;end;k=1;b12=1.0;while andk=k+1;for i=1:m8-1duo(i,1)=(uo(i+1,1)-uo(i,1))/d;dvo(i,1)=(vo(i+1,1)-vo(i,1))/d;end;m9=zeros(2,2);m9(1,1)=1;i=1;m80(1,1,i)=-f1-K;m80(1,2,i)=-f1;m80(2,1,i)=-f1;m80(2,2,i)=-f1-K1;m50(:,:,i)=m80(:,:,i)*inv(m12);z(:,i)=inv(m12)*y11(:,i)*;for i=2:m8-1;m80(1,1,i)=-f1-K;m80(1,2,i)=-f1;m80(2,1,i)=-f1;m80(2,2,i)=-f1-K1;m50(:,:,i)=inv(m12)*m80(:,:,i);end;U(1,m8)=1/2;U(2,m8)=0.0;for i=1:m8-1U(:,m8-i)=m50(:,:,m8-i)*U(:,m8-i+1)+z(:,m8-i);end;for i=1:m8u(i,1)=U(1,i);v(i,1)=U(2,i);end;b12=max(abs(u-uo))uo=u;vo=v;u(m8/2,1)end;for i=1:m8y(i)=i*d;end;plot(y,uo)**************************************
3.1. A General Proposal for Relaxation
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)
5. An Improvement of the Convexity Conditions for a Non-Convex Related Model through an Approximate Primal Formulation
5.1. A Duality Principle for the Concerning Quasi-Convex Envelope
6. A Duality Principle for a Related Relaxed Formulation Concerning the Vectorial Approach in the Calculus of Variations
6.1. An Example in Finite Elasticity
7. An Exact Convex Dual Variational Formulation for a Non-Convex Primal One
8. Another Primal Dual Formulation for a Related Model
9. A Third Primal Dual Formulation for a Related Model
10. An Algorithm for a Related Model in Shape Optimization
10.1. Introduction
10.2. Mathematical Formulation of the Topology Optimization Problem
10.3. About a Concerning Algorithm and Related Numerical Method
- Set and .
- Calculate such that
- Calculate such that
- If or then stop, else set and go to item 2.
-
clear allglobal P m8 d w u v Ea Eb Lo d1 z1 m9 du1 du2 dv1 dv2 c3m8=27;m9=24;c3=0.95;d=1.0/m8;d1=0.5/m9;Ea=; (stronger material)Eb=1000; (softer material simulating voids)w=0.30;P=-42000000;z1=(m8-1)*(m9-1);A3=zeros(z1,z1);for i=1:z1A3(1,i)=1.0;end;b=zeros(z1,1);uo=0.000001*ones(z1,1);u1=ones(z1,1);b(1,1)=c3*z1;for i=1:m9-1for j=1:m8-1Lo(i,j)=c3;end; end;for i=1:z1x1(i)=c3*z1;end;for i=1:2*m8*m9xo(i)=0.000;end;xw=xo;xv=Lo;for k2=1:24c3=0.98*c3;b(1,1)=c3*z1;k2b14=1.0;k3=0;while andk3=k3+1;b12=1.0;k=0;while andk=k+1;k2k3kX=fminunc(’funbeam’,xo);xo=X;b12=max(abs(xw-xo));xw=X;end;for i=1:m9-1for j=1:m8-1ex=du1(i,j);ey=dv2(i,j);exy=1/2*(dv1(i,j)+du2(i,j));Sxy=E1/(2*(1+w))*exy;dc3(i,j)=-(Sx*ex+Sy*ey+2*Sxy*exy);end;end;for i=1:m9-1for j=1:m8-1f(j+(i-1)*(m8-1))=dc3(i,j);end;end;for k1=1:1k1X1=linprog(f,,,A3,b,uo,u1,x1);x1=X1;end;for i=1:m9-1for j=1:m8-1Lo(i,j)=X1(j+(m8-1)*(i-1));end;end;b14=max(max(abs(Lo-xv)))xv=Lo;colormap(gray); imagesc(-Lo); axis equal; axis tight; axis off;pause(1e-6)end;end;
- function S=funbeam(x)
- For a two dimensional beam of dimensions and we have obtained the following results:
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
15. A Related Restricted Problem in Phase Transition
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
- Similar constraints are valid corresponding to the charge of a single proton.
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
- 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)
- (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
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
- At this point, we also recall that the (iron) atom has 26 protons, 26 electrons and 30 neutrons.
- On the other hand a atom has 6 protons and this same number of electrons and neutrons.
- Here we define the density function , representing the Austenite phase, where:
- Similarly, we define the density function for the Martensite phase, which is denoted by , where:
-
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
27. A Dual Variational Formulation for a Related Model
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. Conclusion
- (1)
- Conflict of interest declaration: The author declares no conflict of interest concerning this article.
- (2)
-
Data Avaliability: Details on the software for numerical results avaialable upon request.e-mail:fabio.botelho@ufsc.br
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. Bulletin of the Polish Academy of Sciences, Technical Sciences 1988, 38. [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, Warszawa. [Google Scholar]
- 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-1317 (1989).
- 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.
- 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, 2003. [Google Scholar]
- Botelho, F. Functional Analysis and Applied Optimization in Banach Spaces; Springer: Switzerland, 2014. [Google Scholar]
- Botelho, F.S. Functional Analysis, Calculus of Variations and Numerical Methods in Physics and Engineering; CRC Taylor and Francis: Florida, 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, 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. Applied Mathematics and Computation 2012, 218(11976) 11989. [Google Scholar] [CrossRef]
- Rockafellar, R.T. Convex Analysis; Princeton Univ. Press, 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, Reprint, 2010).
- L.D. Landau and E.M. Lifschits, Course of Theoretical Physics, Vol. 5- Statistical Physics, part 1. (Butterworth-Heinemann, Elsevier, reprint 2008).
- Botelho, F.S. Advanced Calculus and its Applications in Variational Quantum Mechanics and Relativity Theory; CRC Taylor and Francis: Florida, 2021. [Google Scholar]
- Ciarlet, P. Mathematical Elasticity; Vol. II – Theory of Plates; North Holland Elsevier, 1997. [Google Scholar]
- Attouch, H.; Buttazzo, G.; Michaille, G. Variational Analysis in Sobolev and BV Spaces; MPS-SIAM Series in Optimization: Philadelphia, 2006. [Google Scholar]
- Strikwerda, J.C. Finite Difference Schemes and Partial Differential Equations, 2nd ed.; SIAM: Philadelphia, 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, USA, 1959; Volume III. [Google Scholar]
- Folland, G.B. Quantum Field Theory, A Tourist Guide for Mathematicians, Mathematical Surveys and Monographs; Americam Mathematical Society, AMS, Providence: Rhode Island, 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/).