Acessibilidade / Reportar erro

A Note on the Well-Posedness of Control Complex Ginzburg-Landau Equations in Zhidkov Spaces

ABSTRACT

In this note, we consider the Complex Ginzburg-Landau equations with a bilinear control term in the real line. We prove well-posedness results concerned with the initial value problem for these equations in Zhidkov spaces using splitting methods.

Keywords:
well-posedness; Zhidkov spaces; Lie-Trotter method

1 INTRODUCTION

In this note, we deal with the 1-dimensional system

t u = ( α + i β ) x x u + γ u + ( c + i d ) | u | 2 u + ( a + i b ) v ( x , t ) u , u ( 0 ) = u 0 (1.1)

where u(x,t) is a complex valued function with x, t>0, α>0, β>0, γ0, a,b,c,d>0 and v is a bounded control function. The linear term represented by (α+iβ)xx characterizes the Complex Ginzburg-Landau equation (CGLe). For β = 0 (1.1) reduces to a nonlinear heat equation and for α = 0 to a nonlinear Schro¨dinger equation. The cubic CGLe is one of the most important nonlinear equations with applications in physics. It describes a large number of linear and nonlinear phenomena from superconductivity, superfluidity and Bose-Einstein condensation to liquid crystals 11 I.S. Aranson & L. Kramer. The world of the complex Ginzburg-Landau equation. Reviews of modern physics, 74(1) (2002), 99-143.. Well-posedness of (1.1) has been studied with different nonlinearities and in different spaces (see for instance, 44 A. Besteiro. A note on dark solitons in nonlinear complex Ginzburg-Landau equations. Mathematica, 62-85(1) (2020), 11-15.),(1111 J. Ginibre & G. Velo. The Cauchy problem in local spaces for the complex Ginzburg-Landau equation Compactness methods. Physica D: Nonlinear Phenomena , 95(3-4) (1996), 191-228.),(1212 J. Ginibre & G. Velo . The Cauchy Problem in Local Spaces for the Complex Ginzburg-Landau Equation II. Contraction Methods. Communications in mathematical physics, 187(1) (1997), 45-79.). Our aim is to study the well-posedness of the Complex Ginzburg-Landau equation with a bilinear control term, in Zhidkov spaces, using splitting methods. Controllability problems in parabolic equations were studied with different control alternatives and nonlinearities 22 D. Battogtokh & A. Mikhailov. Controlling turbulence in the complex Ginzburg-Landau equation. Physica D: Nonlinear Phenomena, 90(1-2) (1996), 84-95.), (33 D. Battogtokh , A. Preusser & A. Mikhailov . Controlling turbulence in the complex Ginzburg-Landau equation II. Two-dimensional systems. Physica D: Nonlinear Phenomena , 106(3-4) (1997), 327-362.), (1515 M. Ouzahra. Approximate and exact controllability of a reaction-diffusion equation governed by bilinear control. European Journal of Control, 32 (2016), 32-38.), (1616 J. Xiao, G. Hu, J. Yang & J. Gao. Controlling turbulence in the complex Ginzburg-Landau equation. Physical review letters, 81(25) (1998), 5552.. Zhidkov spaces were introduced by P. Zhidkov in 1717 P. Zhidkov. The Cauchy problem for the nonlinear Schrödinger equation. Technical report, Joint Inst. for Nuclear Research (1987). defined as bounded and uniformly continuous functions, with derivatives up to k order in L 2. Many applications were found for these spaces, for instance, in nonlinear optics, Zhidkov functions are used to model dark solitons. In 88 N. Efremidis, K. Hizanidis, H. Nistazakis, D. Frantzeskakis & B. Malomed. Stabilization of dark solitons in the cubic Ginzburg-Landau equation. Physical Review E, 62(5) (2000), 7410., dark soliton solutions are described for a special case of the complex Ginzburg-Landau equation. A typical example of a function in Zhidkov spaces is described in 1010 C. Gallo et al. Schrödinger group on Zhidkov spaces. Advances in Differential Equations, 9(5-6) (2004), 509-538.), (1313 Y.S. Kivshar & B. Luther-Davies. Dark optical solitons: physics and applications. Physics reports, 298(2-3) (1998), 81-197.. These are solutions of the form u(x,t)=uv(x-vt), in particular for the one dimensional case we have:

u v ( x ) = 1 - v 2 2 tanh 1 - v 2 2 x 2 + i v 2

The goal of this article is to prove well-posedness of (1.1) for Zhidkov spaces in the real line, using splitting the result in 44 A. Besteiro. A note on dark solitons in nonlinear complex Ginzburg-Landau equations. Mathematica, 62-85(1) (2020), 11-15.. These are numerical methods that split the flow of the equation, to approximate solutions. The time interval is divided into equal parts, and in each one, the equation evolves alternating the linear and nonlinear flows. This establishes an advantage for us, it exchanges a complicated problem (1.1), for two simpler equations. In this case, we extend the method to a “triple splitting”, dividing the equation into three parts. It is important to remark that this same method can be applied to prove well-posedness for other well known equations such as, reaction-diffusion and Schro¨dinger control equations in L p spaces. The splitting method is based on a Lie-Trotter method developed recently for numerical purposes 66 J.P. Borgna, M.D. Leo, D. Rial & C.S. de la Vega. General Splitting methods for abstract semilinear evolution equations. Commun. Math. Sci, Int. Press Boston, Inc., 13 (2015), 83-101.), (1414 M.D. Leo , D. Rial & C.F.S. de la Vega. High-order time-splitting methods for irreversible equations. IMA J. Numer. Anal., (2015), 1842-1866..

The paper is organized as follows: In Section 2 we set notations and state some preliminary results. In section 3 we analyze the nonlinear problem. Finally, in section 4 and using splitting methods, we combine results from sections 2 and 3 to show that the solution of (1.1) is in a Zhidkov space.

2 NOTATIONS AND PRELIMINARIES.

We introduce some definitions and preliminary results.

Definition 2.1.We define Cu(ℝ) as the set of uniformly continuous and bounded functions on.

Definition 2.2. For k > d / 2 , we define the Zhidkov space as,

X k ( d ) = { u L ( d ) C u ( d ) : j u L 2 ( d ) , 1 | j | k }

equipped with the norm:

u X k = u L + 1 | j | k j u L 2 (2.1)

Remark 2.1.Zhidkov spaces are closed for the norm defined in (2.1). (See1010 C. Gallo et al. Schrödinger group on Zhidkov spaces. Advances in Differential Equations, 9(5-6) (2004), 509-538.)

The following definitions and proofs can be extended to xd (See 99 K.J. Engel & R. Nagel. “One-parameter semigroups for linear evolution equations”, volume 194. Springer Science & Business Media (1999).).

Definition 2.3.We denote U(t) as the one parameter semigroup that solves the underlying linear equation

t u = ( α + i β ) x x u + γ u (2.2)

The operator can be represented by the convolution in x

U ( t ) u 0 = ( 4 π t ( α + i β ) ) - 1 / 2 e ( - x 2 / [ 4 t ( α + i β ) ] ) + γ t * u 0 = G t ( x ) * u 0

and the kernel G t satisfies:

| G t ( x ) | = ( 4 π t ( α 2 + β 2 ) 1 / 2 ) - 1 / 2 e ( - α x 2 / [ 4 t ( α 2 + β 2 ) ] ) + γ t

Proposition 2.1.The one-parameter family{U(t)}t0of operators defined asU(t)u0=Gt*u0is a strongly continuous semigroup on Cu (ℝ).

Proof. The proof is similar to Proposition 2.2 in 55 A. Besteiro & D. Rial. Global existence for vector valued fractional reaction-diffusion equations. Publicacions Matemàtiques, 96 (2021), 653-680.. □

Lemma 2.1. If u 0 X 1 ( ) then U ( t ) u 0 C ( [ 0 , T * ( u 0 ) ) , X 1 ( ) ) for t > 0

Proof. As u0L() and Gt(x)L1() then using Young’s inequality we have

G t * u 0 L G t L 1 u 0 L

On the other hand, we obtain

x ( G t * u 0 ) L 2 = G t * x u 0 L 2 G t L 1 * x u 0 L 2

As GtL1() and xu0L2() we have the result. □

Remark 2.2.Similarly, ifxdand we have k derivatives of U(t)u0, a similar procedure proves thatU(t)u0C([0,T*(u0)), Xk(d)).

Next, we consider integral solutions of the problem (1.1). We say that uC([0,T], Cu()) is a mild solution of (1.1) if and only if u verifies

u ( t ) = U ( t ) u 0 + 0 t U ( t - t ' ) B ( x , t ' , u ( t ' ) ) d t ' . (2.3)

where B(x,t,u)=(c+id)|u|2u+(a+ib)v(x,t)u. If B is a locally Lipschitz map, for any z0Cu() there exists a unique solution of the equation

t z = B ( t , z ) , z ( 0 ) = z 0 , (2.4)

defined in the interval [0,T*(z0)). Moreover, there exists a nonincreasing function T¯ : [0,)[0,), such that T*(z0)T¯(|z0|). The solution of (2.4) is solution of the integral equation

z ( t ) = z 0 + 0 t B ( t ' , z ( t ' ) ) d t ' . (2.5)

Also, one of the following alternatives holds:

  • - T*(z0)= ;

  • - T*(z0)< and |z(t)| when tT*(z0) .

We denote by N:+×+×Cu()Cu() the flow generated by the ordinary equation, i.e.: for any x, N(t,t0,u0)(x) is the solution of the problem (2.4) with initial datum z 0 = u 0(x). Therefore, if u(t) = N(t, t 0 , u 0)

u ( x , t ) = u 0 ( x ) + 0 t B ( x , t ' , u ( x , t ' ) ) d t '

We recall a well-known local existence result for evolution equations.

Theorem 2.1.There exists a functionT*:Cu()+such that foru0Cu(), exists a uniqueuC([0,T*(u0)), Cu())mild solution of (1.1) with u(0) = u 0 . Moreover, one of the following alternatives holds:

  • T*(u0)= ;

  • T*(u0)<andlimtT*(u0) supx|u(t)|= .

Proof. See Theorem 4.3.4 in 77 T. Cazenave & A. Haraux. “An Introduction to Semilinear Evolution Equations”. Oxford Lecture Ser. Math. Appl., Clarendon Press, Rev ed. (1999).. □

Proposition 2.2. Under conditions of the theorem above, the following statements hold true:

  1. T * : C u ( ) + is lower semi-continuous;

  2. Ifu0,nu0in Cu (ℝ) and0<T<T*(u0), then u nu in the Banach spaceC([0,T], Cu()) .

Proof. See Proposition 4.3.7 in 77 T. Cazenave & A. Haraux. “An Introduction to Semilinear Evolution Equations”. Oxford Lecture Ser. Math. Appl., Clarendon Press, Rev ed. (1999).. □

3 NONLINEAR EQUATION

In this section, we first analyze the following control equation:

t u = ( a + i b ) v ( x , t ) u , u ( 0 ) = u 0 (3.1)

where vC([0,T*(u0)), L(Ω)) and Ω is a bounded interval of ℝ with supp(v(·,t))Ω. The following Lemma allows us to have well-posedness of the control equation in X 1(ℝ) which is essential to apply the splitting method.

Lemma 3.2.LetvC([0,T*(u0)), L(Ω))withxvC([0,T*(u0)), L(Ω))and letB(x,t,u)=(a+ib)v(x,t)u, thenB : +×+×X1()X1()is a well-defined operator and B(x,t,u) is a locally Lipschitz map in u.

Proof. Let uX1() then

B ( x , t , u ) X 1 = ( a + i b ) v ( x , t ) u X 1 K ( v ( x , t ) u L + x ( v ( x , t ) u ) L 2 ) = K ( v ( x , t ) u L + x v ( x , t ) u + v ( x , t ) x u L 2 ) K ( v ( x , t ) u L + x v ( x , t ) u L 2 + v ( x , t ) x u L 2 )

As v(x,t) is a bounded function in Ω then B(x,t,u)X1<. On the other hand, using the notation B(x,t,u)=B(u), we have

B ( u ) - B ( w ) X 1 = ( a + i b ) v ( x , t ) u - ( a + i b ) v ( x , t ) w X 1 = ( a + i b ) v ( x , t ) ( u - w ) X 1 K v ( x , t ) ( u - w ) X 1 K v ( x , t ) ( u - w ) L + K x ( v ( x , t ) ( u - w ) ) L 2 K v ( x , t ) ( u - w ) L + K ( x v ( x , t ) ) ( u - w ) L 2 + K ( v ( x , t ) ) ( x ( u - w ) ) L 2 K v ( x , t ) L ( u - w ) L + K ( x v ( x , t ) ) L ( u - w ) ) L 2 + K ( v ( x , t ) ) L ( x ( u - w ) ) L 2 = K ' ( u - w ) L + K ' ' ( u - w ) ) L 2 + K ' ( x ( u - w ) ) L 2 k u - w X 1

where in the last step we used Hölder’s inequality. □

We have the same result for the solution for the nonlinear problem associated with the term |u|2 u, that is the equation

t u = - ( c + i d ) | u | 2 u , u ( 0 ) = u 0 ,

Lemma 3.3.Ifu0X1()then the solution of the problem (3.2),u(t)C([0,T*(u0)), X1()).

Proof. See 44 A. Besteiro. A note on dark solitons in nonlinear complex Ginzburg-Landau equations. Mathematica, 62-85(1) (2020), 11-15. Lemma 3.1. □

4 SPLITTING METHOD

This section is based on the splitting method developed in 55 A. Besteiro & D. Rial. Global existence for vector valued fractional reaction-diffusion equations. Publicacions Matemàtiques, 96 (2021), 653-680.), (66 J.P. Borgna, M.D. Leo, D. Rial & C.S. de la Vega. General Splitting methods for abstract semilinear evolution equations. Commun. Math. Sci, Int. Press Boston, Inc., 13 (2015), 83-101.), (1414 M.D. Leo , D. Rial & C.F.S. de la Vega. High-order time-splitting methods for irreversible equations. IMA J. Numer. Anal., (2015), 1842-1866.. We apply the Lie-Trotter method to the linear and nonlinear problem. The temporal variable must be broken down into regular intervals and the evolution of the linear, nonlinear and control problems are considered alternately. This is described by three sequences: {V 0,k } for the linear equation and {V 1,k }, {W 2,k } for the nonlinearity and the control term, respectively. Using Theorem 3.9 from 1414 M.D. Leo , D. Rial & C.F.S. de la Vega. High-order time-splitting methods for irreversible equations. IMA J. Numer. Anal., (2015), 1842-1866., the approximate solution converges to the solution of problem (1.1), when the time intervals h=t/n0.

Let X be a Banach space and we define αj : a periodic function of period 1 as:

α j ( t ) = 3 , if k + j / 3 t < k + ( j + 1 ) / 3 , 0 , if k + ( j + 1 ) / 3 t < k + j + 1 ,

for k and j = 0, 1, 2.

Given h>0, we define the function αh : as αjh(t)=αj(t/h). Clearly 0αjh3, α jh is h-periodic and its mean value is 1.

We consider τh : 2 given by

τ h ( t , t ' ) = t ' t α 0 h ( t ' ' ) d t ' ' ,

We define ω={(t,t')2:0t't} and Uh:ωB(X) given by Uh(t,t')=U(τh(t,t')).

We consider the system,

t u h + α 0 h ( t ) σ ( - x x ) u h ( x , t ) = α 1 h ( t ) B 1 ( x , t , u h ( x , t ) ) + α 2 h ( t ) B 2 ( x , t , u h ( x , t ) ) , u h ( x , 0 ) = u h 0 ( x )

where u(x,t)X, t>0, σ and Bj:××XX is a continuous function with j = 1, 2. Similarly, we define the integral equation:

u h ( t ) = U h ( t , 0 ) u h 0 + 0 t U h ( t , t ' ) ( α 1 h ( t ' ) B 1 ( x , t , u h ( t ' ) ) + α 2 h ( t ' ) B 2 ( x , t , u h ( t ' ) ) ) d t ' (4.1)

The following two theorems are developed similar to the results of section 4 of 55 A. Besteiro & D. Rial. Global existence for vector valued fractional reaction-diffusion equations. Publicacions Matemàtiques, 96 (2021), 653-680., where all results are proved for one nonlinearity. We extend the proofs to consider two nonlinearities. The following Theorem is similar to Propostion 4.3 of 55 A. Besteiro & D. Rial. Global existence for vector valued fractional reaction-diffusion equations. Publicacions Matemàtiques, 96 (2021), 653-680..

Theorem 4.2.Let uhbe the solution of (4.1), ifW2,k=uh(kh+h), V0,k=uh(kh+h/3)andV1,k=uh(kh+2h/3)then

V 0 , k + 1 = U ( h ) W 2 , k , (4.2a)

V 1 , k + 1 = N 1 ( k h + h / 3 , k h + 2 h / 3 , V 0 , k + 1 ) , (4.2b)

W 2 , k + 1 = N 2 ( k h + 2 h / 3 , k h + h , V 1 , k + 1 ) , (4.2c)

where Nj, (j = 1, 2) is the flux associated to:

t w = α j h ( t ' ) B j ( t , w ( t ) ) , w ( 0 ) = w 0 ,

Proof. For t1(0,t) it is verified

u h ( t ) = U h ( t , t 1 ) u h 0 ( t 1 ) + t 1 t U h ( t , t ' ) ( α 1 h ( t ' ) B 1 ( x , t , u h ( t ' ) ) + α 2 h ( t ' ) B 2 ( x , t , u h ( t ' ) ) ) d t '

taking that t 1 = kh and t = kh+h/3, we have

V 0 , k + 1 = U h ( k h + h / 3 , k h ) W 2 , k + k h k h + h / 3 U h ( t , t ' ) ( α 1 h ( t ' ) B 1 ( x , t , u h ( t ' ) ) + α 2 h ( t ' ) B 2 ( x , t , u h ( t ' ) ) ) d t ' ,

since α0h(t')=3 for t'[kh,kh+h/3), we have τh(kh+h/3,kh)=h and we get (4.2a). Similarly, α0h(t')=α2h(t')=0 for t[kh+h/3,kh+2h/3), then τh(t,kh+h/3)=0 and therefore

u h ( t ) = V 0 , k + 1 + 3 k h + h / 3 t B 1 ( x , t , u h ( t ' ) ) d t ' ,

evaluating in t=kh+2h/3, we obtain (4.2b). The same process can be done to obtain (4.2c). □

Theorem 4.3. Let u C ( [ 0 , T * ) , X ) the be solution of the integral problem

u ( t ) = U ( t ) u 0 + 0 t U ( t - t ' ) B ( x , t , u ( t ' ) ) d t ' ,

withB=B1+B2both defined as in (4.1). Let alsoT(0,T*)andε>0. Then there existsh*>0such that if0<h<h*, then u h the solution of (4.1) withuh(x,0)=u0(x), is defined in the interval [0, T] and verifiesu(t)-uh(t)Xεfort[0,T].

Proof. The proof is similar to Theorem 4.4 from 55 A. Besteiro & D. Rial. Global existence for vector valued fractional reaction-diffusion equations. Publicacions Matemàtiques, 96 (2021), 653-680., considering two distinct nonlinear Lipschitz terms. □

Now, we apply Lemma 2.1 from Section 2 related to linear equation and Lemmas 3.2 and 3.3 from Section 3 related to the control equation. In order to obtain the well-posedness result for the solution u(t) of equation (1.1), we use Theorem 4.3. We denote by N1:+×+×Cu()Cu() the flow generated by the equation (3.2) as u(t) = N 1(t,t 0 ,u 0), and similarly u(t) = N 2(t,t 0 ,u 0) the flow generated by the equation (3.1) defined for t0t<T*(t0,u0).

Theorem 4.4.Letu0X1(), then the solution of satisfies (1.1) u(t)C([0,T*(u0)), X1())fort(0,T*(u0)).

Proof. For t[0,T*(u0)), let n, h=t/n and {W2,k}0kn, {V0,k}1kn, {V1,k}1kn be the sequences given by W 2,0 = u 0,

V 0 , k + 1 = U ( h ) W 2 , k , V 1 , k + 1 = N 1 ( k h + 2 h / 3 , k h + h / 3 , V 0 , k + 1 ) , W 2 , k + 1 = N 2 ( k h + h , k h + 2 h / 3 , V 1 , k + 1 ) , k = 0 , , n - 1 .

We claim that W2,kC([0,T*(u0)), X1()) for k = 0, . . . , n. Clearly, the assertion is true for k = 0. If W2,kC([0,T*(u0)), X1()), from Lemma 2.1, we have U(h)W2,kC([0,T*(u0)), X1()) and from Lemma 3.3 we can see that

V 1 , k + 1 = N 1 ( h , V 0 , k + 1 ) C ( [ 0 , T * ( u 0 ) ) , X 1 ( ) ) .

Similarly, using Lemma 3.2, we have that

W 2 , k + 1 = N 2 ( h , V 1 , k + 1 ) C ( [ 0 , T * ( u 0 ) ) , X 1 ( ) ) .

By Theorem 4.3 we have that W2,nu(t) when n. As X 1(ℝ) is closed, we obtain the result.

Aknowledgments

This work was supported by Universidad Abierta Interamericana (UAI) and CONICET- Argentina.

REFERENCES

  • 1
    I.S. Aranson & L. Kramer. The world of the complex Ginzburg-Landau equation. Reviews of modern physics, 74(1) (2002), 99-143.
  • 2
    D. Battogtokh & A. Mikhailov. Controlling turbulence in the complex Ginzburg-Landau equation. Physica D: Nonlinear Phenomena, 90(1-2) (1996), 84-95.
  • 3
    D. Battogtokh , A. Preusser & A. Mikhailov . Controlling turbulence in the complex Ginzburg-Landau equation II. Two-dimensional systems. Physica D: Nonlinear Phenomena , 106(3-4) (1997), 327-362.
  • 4
    A. Besteiro. A note on dark solitons in nonlinear complex Ginzburg-Landau equations. Mathematica, 62-85(1) (2020), 11-15.
  • 5
    A. Besteiro & D. Rial. Global existence for vector valued fractional reaction-diffusion equations. Publicacions Matemàtiques, 96 (2021), 653-680.
  • 6
    J.P. Borgna, M.D. Leo, D. Rial & C.S. de la Vega. General Splitting methods for abstract semilinear evolution equations. Commun. Math. Sci, Int. Press Boston, Inc., 13 (2015), 83-101.
  • 7
    T. Cazenave & A. Haraux. “An Introduction to Semilinear Evolution Equations”. Oxford Lecture Ser. Math. Appl., Clarendon Press, Rev ed. (1999).
  • 8
    N. Efremidis, K. Hizanidis, H. Nistazakis, D. Frantzeskakis & B. Malomed. Stabilization of dark solitons in the cubic Ginzburg-Landau equation. Physical Review E, 62(5) (2000), 7410.
  • 9
    K.J. Engel & R. Nagel. “One-parameter semigroups for linear evolution equations”, volume 194. Springer Science & Business Media (1999).
  • 10
    C. Gallo et al. Schrödinger group on Zhidkov spaces. Advances in Differential Equations, 9(5-6) (2004), 509-538.
  • 11
    J. Ginibre & G. Velo. The Cauchy problem in local spaces for the complex Ginzburg-Landau equation Compactness methods. Physica D: Nonlinear Phenomena , 95(3-4) (1996), 191-228.
  • 12
    J. Ginibre & G. Velo . The Cauchy Problem in Local Spaces for the Complex Ginzburg-Landau Equation II. Contraction Methods. Communications in mathematical physics, 187(1) (1997), 45-79.
  • 13
    Y.S. Kivshar & B. Luther-Davies. Dark optical solitons: physics and applications. Physics reports, 298(2-3) (1998), 81-197.
  • 14
    M.D. Leo , D. Rial & C.F.S. de la Vega. High-order time-splitting methods for irreversible equations. IMA J. Numer. Anal., (2015), 1842-1866.
  • 15
    M. Ouzahra. Approximate and exact controllability of a reaction-diffusion equation governed by bilinear control. European Journal of Control, 32 (2016), 32-38.
  • 16
    J. Xiao, G. Hu, J. Yang & J. Gao. Controlling turbulence in the complex Ginzburg-Landau equation. Physical review letters, 81(25) (1998), 5552.
  • 17
    P. Zhidkov. The Cauchy problem for the nonlinear Schrödinger equation. Technical report, Joint Inst. for Nuclear Research (1987).

Publication Dates

  • Publication in this collection
    05 Sept 2022
  • Date of issue
    Jul-Sep 2022

History

  • Received
    26 Aug 2021
  • Accepted
    24 Mar 2022
Sociedade Brasileira de Matemática Aplicada e Computacional - SBMAC Rua Maestro João Seppe, nº. 900, 16º. andar - Sala 163, Cep: 13561-120 - SP / São Carlos - Brasil, +55 (16) 3412-9752 - São Carlos - SP - Brazil
E-mail: sbmac@sbmac.org.br