In this paper, we study the stationary compressible nonisothermal nematic liquid crystal flows affected by the external force of general form in three-dimensional space. By using the contraction mapping principle, we prove the existence and uniqueness of strong solution around the constant state in some suitable function space.

1. Introduction and Main Result

We consider the problem of identifying the stationary motion of compressible nonisothermal nematic liquid crystal flows effected by the external force of general form:

Here, the unknown functions , , and are the density, velocity, absolute temperature, and macroscopic average of the nematic liquid crystal orientation field, respectively. The pressure is a smooth function of satisfying . The constants are the shear and bulk viscosity coefficients of the fluids, respectively, which are assumed to satisfy the physical restrictions and . The constant is the ratio of the heat conductivity coefficient over the heat capacity. Moreover, is the dissipation function:

In addition, , , , and are the given mass source, external force, director source, and energy source, respectively, which are given by the following form where and

As the space variable tends to infinity, we assume where and are two given constants, and is a unit constant vector.

The flow of nematic liquid crystals can be regarded as slow moving particles, in which the fluid velocity and the particles arrangement influence each other. The continuum theory of the nematic liquid crystals was first proposed by Ericksen [1] and Leslie [2] during the period between 1958 and 1968, see also the book by De Gennes [3]. The hydrodynamic flow of nematic liquid crystal system has attracted great interest and generated many important developments. Here, we only mention some of related to our study. When the temperature is absent, i.e., . In [4, 5], Huang et al. addressed some issues on the strong solutions in three-dimensional space. Specifically, in [4], the authors established a blow-up criterion, while the local existence of a unique strong solution if the initial data are sufficiently regular and satisfy a natural compatibility condition are studied in [5]. Later, Jiang et al. [6] considered the global existence of weak solutions with large initial energy and without any smallness condition on the initial density and velocity in a bounded domain of the multidimensional space. By using the domain expansion technique and the rigidity theorem, these authors [7] also proved the global existence of large weak solutions in two dimensions, provided that the second component of initial data of the direction field satisfies some geometric angle condition. One may also see [8] for some recent progress on the existence, regularity, uniqueness, and large time asymptotic of the nematic liquid crystal flows.

Recently, Feireisl et al. and Feireisl et al. [9, 10] established nonisothermal models of incompressible nematic liquid crystals and obtained the global existence of weak solutions. For further study about the well-posedness of solutions to the incompressible nematic liquid crystal flows, we refer to [1115] and references therein. For the compressible nonisothermal nematic liquid crystal flows, Guo et al. [16] shown the existence of global weak solutions by a three-level approximation and weak convergence when the adiabatic exponent . Moreover, they [17] also considered the global existence and decay rates of small smooth solutions in the whole space . Later, Fan et al. [18] proved the local well-posedness of strong solutions to the initial boundary value problem in . By some delicate energy estimates and the critical Sobolev inequalities of logarithmic type, Zhong [19] obtained the singularity formation of strong solutions to the initial boundary problem in . Then, Liu and Zhong [20] established the global existence and uniqueness of strong solutions with vacuum as far field density in three-dimensional space. However, a lot of physical and mathematical important problems are still open due to the lack of a smoothing mechanism and the strong nonlinearity. Up to now, there are no result available on the existence of stationary solutions to the three-dimensional compressible nonisothermal nematic liquid crystal flows with external force.

Inspired by the work of [21, 22] for the compressible Navier-Stokes equations, the aim of this paper is to study the existence and uniqueness of stationary solution, which is a small strong solution around the constant state to problem (1)–(4). It is worth mentioning here that the nonisothermal nematic liquid crystal flow system (1) adds the strong nonlinear terms, which seems more complicated than that of the Navier-Stokes equations, and we should carefully deal with the direction field of liquid crystals in the angular momentum equation. Moreover, we think that discussing the stationary solutions to (1) is of significance to some extent.

As in [21], we choose as the independent variables and regard as a smooth function of . To this end, set and denote Then, (1) can be reformulated as

As a preparation for stating our main result, we introduce some notations and conventions to be used throughout this paper. Here, denotes a generic positive constant which may vary depending on the estimate. with any integer stands for all derivatives up to -th order of the function with respect to the spatial variable . stands for the usual spaces with norm and for any integer , stands for the usual Sobolev spaces with norm . Let

Finally, we define some function spaces. For any , denote with

Moreover, we put where

Now, we state our main result about the existence and uniqueness of stationary solutions to the stationary problem (1).

Theorem 1. There exist small constants and depending only on , such that if for some positive constant , then the problem (6) admits a solution .

Moreover, the solution is unique in the following sense: if there is another solution satisfying (6) with the same , and , then

The rest of this paper is organized as follows. In Section 2, we establish the weighted theory for the corresponding linearized problem of (6). The proof for the existence and uniqueness of the stationary solutions to (6) will be considered in Section 3.

2. Weighted Theory for Linearized Problem

In this section, we study the weighted theory for the following linearized system of (6):

Here, and are given as

Moreover, we fix or and always assume that

We would like to point out that equations (13)1–(13)3, which are independent of , have been well studied in [21], while (13)4 is easy to be handled. Hence, by using the same arguments as in [21], that is by using the Banach closed range theorem and some weighted– estimates on the linearized problem (13), we get the following result.

Theorem 2. There exists a positive constants depending only on and such that if in (15) satisfies , then the problem (13) has a solution satisfying where is a constant depending only on .

3. The Proof of Theorem 1

In this section, we prove the main Theorem 1 for the existence of stationary solution of (6) by using the contraction mapping principle in . For this purpose, we study the following iterated equations where is given, , etc., and

3.1. Construction of Solution Map for (18)

Now, we devote ourselves to establishing a solution map for (18). To begin with, we derive the weighted estimates on the solution of (18) by applying Theorem 2. In fact, let

and in Theorem 2 be denoted as in (19). By the Sobolev inequality, we can choose sufficiently small such that and in (15) satisfy . It is straightforward to check that (15) holds for , and with defined by

Thus, one has the desired result by applying Theorem 2 with for (18).

Lemma 3. Let satisfy (22). Then, there exists a positive constants such that if , the problem (18) with has a solution satisfying where is a constant depending only on .

To continue, we cite a lemma in [22] which will be used to establish the norm of the solution to (18).

Lemma 4 (see [22]). Let be a scalar function satisfying (i)If is a smooth scalar function in the form: and satisfiesthen for any multi-index with , one has (ii)If is a smooth scalar function in the form: and satisfiesthen for any multi-index with one has Here, the constant depends only on .
According to Lemma 3, we can introduce the solution map for (18): by . In fact, in the following proposition, we can see that maps into itself. That is, .

Proposition 5. There exists such that for sufficiently small , if satisfies then (18) with has a solution .

Proof. We prove this result by two steps.

Step 1. We begin with the -norm of the solution to (18). From the Helmholtz decomposition and Fourier transform, the solution of (18) can be written of the form and compare [22]: where Now, we apply Lemma 4 to get the estimate for . As a start, we study more details of in order to estimate and . Since , there exists and such that This in turn yields In view of (32), the Sobolev inequality and the fact that , a direct computation gives that with defined by By Lemma 4(i), we thus obtain Next, we calculate the estimate for . Notice that Since , it follows from Lemma 3 that Using Lemma 4(ii), we immediately find Hence, according to (30), (36), and (39), we can conclude the estimate for and as To continue, we treat the estimate for . To this end, we first write in the form This together with and (32) that with defined by By Lemma 4(i), it holds that As to the estimate for , we note that can be rewritten as By the Sobolev inequality and the fact that , we obtain with defined by With this help and Lemma 4(i), one has Next, for the case of , in terms of Lemma 3 and the Sobolev inequality, we immediately deduce that Consequently, combining (40), (44), (48), and (49), one has

Step 2. By Lemma 3 and (50), it follows that (18) admits a solution satisfying where the constants depend only on . To complete our proof, we need to study more details on . To this end, let us define and as Thus, from (18)1, the solution can be written as In addition, by and (50), we can get The proof of this proposition is completed.

3.2. Contraction of the Solution Map

Now, we are in a position to prove the solution map for (18) is contractive. Assume that and for . By (18), a straightforward calculation yields that where

By the Sobolev inequality, it holds that with defined in (22). By applying Theorem 2 with to (55), we obtain

With the same computations as used in (50), one can arrive at where are functions satisfying

In addition, we denote and as

We can derived from that which combined with (58) and (59) yields

With the above preparation in hand, we then have the following proposition.

Proposition 6. Assume that satisfies the estimate (for defined in Proposition 5): for some and sufficiently small . Then, for and , one has