首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 27 毫秒
1.
2.
A spatially two-dimensional sixth order PDE describing the evolution of a growing crystalline surface h(xyt) that undergoes faceting is considered with periodic boundary conditions, as well as its reduced one-dimensional version. These equations are expressed in terms of the slopes \(u_1=h_{x}\) and \(u_2=h_y\) to establish the existence of global, connected attractors for both equations. Since unique solutions are guaranteed for initial conditions in \(\dot{H}^2_{per}\), we consider the solution operator \(S(t): \dot{H}^2_{per} \rightarrow \dot{H}^2_{per}\), to gain our results. We prove the necessary continuity, dissipation and compactness properties.  相似文献   

3.
Given bounded vector field \({b : {\mathbb{R}^{d}} \to {\mathbb{R}^{d}}}\), scalar field \({u : {\mathbb{R}^{d}} \to {\mathbb{R}}}\), and a smooth function \({\beta : {\mathbb{R}} \to {\mathbb{R}}}\), we study the characterization of the distribution \({{\rm div}(\beta(u)b)}\) in terms of div b and div(ub). In the case of BV vector fields b (and under some further assumptions), such characterization was obtained by L. Ambrosio, C. De Lellis and J. Malý, up to an error term which is a measure concentrated on the so-called tangential set of b. We answer some questions posed in their paper concerning the properties of this term. In particular, we construct a nearly incompressible BV vector field b and a bounded function u for which this term is nonzero. For steady nearly incompressible vector fields b (and under some further assumptions), in the case when d = 2, we provide complete characterization of div(\({\beta(u)b}\)) in terms of div b and div(ub). Our approach relies on the structure of level sets of Lipschitz functions on \({{\mathbb{R}^{2}}}\) obtained by G. Alberti, S. Bianchini and G. Crippa. Extending our technique, we obtain new sufficient conditions when any bounded weak solution u of \({\partial_t u + b \cdot \nabla u=0}\) is renormalized, that is when it also solves \({\partial_t \beta(u) + b \cdot \nabla \beta(u)=0}\) for any smooth function \({\beta \colon{\mathbb{R}} \to {\mathbb{R}}}\). As a consequence, we obtain new a uniqueness result for this equation.  相似文献   

4.
We construct a Sobolev homeomorphism in dimension \({n \geqq 4,\,f \in W^{1,1}((0, 1)^n,\mathbb{R}^n)}\) such that \({J_f = {\rm det} Df > 0}\) on a set of positive measure and J f  < 0 on a set of positive measure. It follows that there are no diffeomorphisms (or piecewise affine homeomorphisms) f k such that \({f_k\to f}\) in \({W^{1,1}_{\rm loc}}\).  相似文献   

5.
This study considers the quasilinear elliptic equation with a damping term,
$$\begin{aligned} \text {div}(D(u)\nabla u) + \frac{k(|{\mathbf {x}}|)}{|{\mathbf {x}}|}\,{\mathbf {x}}\cdot (D(u)\nabla u) + \omega ^2\big (|u|^{p-2}u + |u|^{q-2}u\big ) = 0, \end{aligned}$$
where \({\mathbf {x}}\) is an N-dimensional vector in \(\big \{{\mathbf {x}} \in \mathbb {R}^N: |{\mathbf {x}}| \ge \alpha \big \}\) for some \(\alpha > 0\) and \(N \in {\mathbb {N}}\setminus \{1\}\); \(D(u) = |\nabla u|^{p-2} + |\nabla u|^{q-2}\) with \(1 < q \le p\); k is a nonnegative and locally integrable function on \([\alpha ,\infty )\); and \(\omega \) is a positive constant. A necessary and sufficient condition is given for all radially symmetric solutions to converge to zero as \(|{\mathbf {x}}|\rightarrow \infty \). Our necessary and sufficient condition is expressed by an improper integral related to the damping coefficient k. The case that k is a power function is explained in detail.
  相似文献   

6.
7.
This paper is concerned with the following fractional Schrödinger equation
$$\begin{aligned} \left\{ \begin{array}{ll} (-\Delta )^{s} u+u= k(x)f(u)+h(x) \text{ in } \mathbb {R}^{N}\\ u\in H^{s}(\mathbb {R}^{N}), \, u>0 \text{ in } \mathbb {R}^{N}, \end{array} \right. \end{aligned}$$
where \(s\in (0,1),N> 2s, (-\Delta )^{s}\) is the fractional Laplacian, k is a bounded positive function, \(h\in L^{2}(\mathbb {R}^{N}), h\not \equiv 0\) is nonnegative and f is either asymptotically linear or superlinear at infinity. By using the s-harmonic extension technique and suitable variational methods, we prove the existence of at least two positive solutions for the problem under consideration, provided that \(|h|_{2}\) is sufficiently small.
  相似文献   

8.
This paper is concerned with time periodic traveling curved fronts for periodic Lotka–Volterra competition system with diffusion in two dimensional spatial space
$$\begin{aligned} {\left\{ \begin{array}{ll} \dfrac{\partial u_{1}}{\partial t}=\Delta u_{1} +u_{1}(x,y,t)\left( r_{1}(t)-a_{1}(t)u_{1}(x,y,t)-b_{1}(t)u_{2}(x,y,t)\right) ,\\ \dfrac{\partial u_{2}}{\partial t}=d\Delta u_{2} +u_{2}(x,y,t)\left( r_{2}(t)-a_{2}(t)u_{1}(x,y,t)-b_{2}(t)u_{2}(x,y,t)\right) , \end{array}\right. } \end{aligned}$$
where \(\Delta \) denotes \(\frac{\partial ^{2}}{\partial x^{2} }+ \frac{\partial ^{2}}{\partial y^{2} }\), \(x,y\in {\mathbb {R}}\) and \(d>0\) is a constant, the functions \(r_i(t),a_i(t)\) and \(b_i(t)\) are T-periodic and Hölder continuous. Under suitable assumptions that the corresponding kinetic system admits two stable periodic solutions (p(t), 0) and (0, q(t)), the existence, uniqueness and stability of one-dimensional traveling wave solution \(\left( \Phi _{1}(x+ct,t),\Phi _{2}(x+ct,t)\right) \) connecting two periodic solutions (p(t), 0) and (0, q(t)) have been established by Bao and Wang ( J Differ Equ 255:2402–2435, 2013) recently. In this paper we continue to investigate two-dimensional traveling wave solutions of the above system under the same assumptions. First, we establish the asymptotic behaviors of one-dimensional traveling wave solutions of the system at infinity. Using these asymptotic behaviors, we then construct appropriate super- and subsolutions and prove the existence and non-existence of two-dimensional time periodic traveling curved fronts. Finally, we show that the time periodic traveling curved front is asymptotically stable.
  相似文献   

9.
Flow development and degradation during Type B turbulent drag reduction by 0.10 to 10 wppm solutions of a partially-hydrolysed polyacrylamide B1120 of MW \(=\) 18x106 was studied in a smooth pipe of ID \(=\) 4.60 mm and L/D \(=\) 210 at Reynolds numbers from 10000 to 80000 and wall shear stresses Tw from 8 to 600 Pa. B1120 solutions exhibited facets of a Type B ladder, including segments roughly parallel to, but displaced upward from, the P-K line; those that attained asymptotic maximum drag reduction at low Re f but departed downwards into the polymeric regime at a higher retro-onset Re f; and segments at MDR for all Re f. Axial flow enhancement profiles of S\(^{\prime }\) vs L/D reflected a superposition of flow development and polymer degradation effects, the former increasing and the latter diminishing S\(^{\prime }\) with increasing distance downstream. Solutions that induced normalized flow enhancements S\(^{\prime }\)/S\(^{\prime }_{\mathrm {m}} <\) 0.4 developed akin to solvent, with Le,p/D \(=\) Le,n/D \(<\) 42.3, while those at maximum drag reduction showed entrance lengths Le,m/D \(\sim \) 117, roughly 3 times the solvent Le,n/D. Degradation kinetics were inferred by first detecting a falloff point (Ref, S\(^{{\prime }\wedge }\)), of maximum observed flow enhancement, for each polymer solution. A plot of S\(^{{\prime }\wedge }\)vs C revealed S\(^{{\prime }\wedge }\)linear in C at low C, with lower bound [S\(^{\prime }\)] \(=\) 5.0 wppm??1, and S\(^{{\prime }\wedge }\) independent of C at high C, with upper bound S\(^{\prime }_{\mathrm {m}} =\) 15.9. The ratio S\(^{\prime }\)/S\(^{{\prime }\wedge }\) in any pipe section was interpreted to be the undegraded fraction of original polymer therein. Semi-log plots of (S\(^{\prime }\)/S\(^{{\prime }\wedge }\)) at a section vs transit time from pipe entrance thereto revealed first order kinetics, from which apparent degradation rate constants kdeg s??1 and entrance severities ?ln(S\(^{\prime }\)/S\(^{{\prime }\wedge }\))0 were extracted. At constant C, kdeg increased linearly with increasing wall shear stress Tw, and at constant Tw, kdeg was independent of C, providing a B1120 degradation modulus (kdeg/Tw) \(=\) (0.012 \(\pm \) 0.001) (Pa s)??1 for 8 \(<\) Tw Pa \(<\) 600, 0.30 \(<\) C wppm \(<\) 10. Entrance severities were negligible below a threshold Twe \(\sim \) 30 Pa and increased linearly with increasing Tw for Tw \(>\) Twe. The foregoing methods were applied to Type A drag reduction by 0.10 to 10 wppm solutions of a polyethyleneoxide PEO P309, MW \(=\) 11x106, in a smooth pipe of ID \(=\) 7.77 mm and L/D \(=\) 220 at Re from 4000 to 115000. P309 solutions that induced S\(^{\prime }\)/S\(^{\prime }_{\mathrm {m}} <\) 0.4 developed akin to solvent, with Le,p/D \(=\) Le,n/D \(<\) 23, while those at MDR had entrance lengths Le,m/D \(\sim \) 93, roughly 4 times the solvent Le,n/D. P309 solutions described a Type A fan distorted by polymer degradation. A typical trajectory departed the P-K line at an onset point Re f* followed by ascending and descending polymeric regime segments separated by a falloff point Ref, of maximum flow enhancement; for all P309 solutions, onset Re f* = 550 \(\pm \) 100 and falloff Ref = 2550 \(\pm \) 250, the interval between them delineating Type A drag reduction unaffected by degradation. A plot of falloff S\(^{{\prime }\wedge }\) vs C for PEO P309 solutions bore a striking resemblance to the analogous S\(^{{\prime }\wedge }\) vs C plot for solutions of PAMH B1120, indicating that the initial Type A drag reduction by P309 after onset at Re f* had evolved to Type B drag reduction by falloff at Ref. Presuming that Type B behaviour persisted past falloff permitted inference of P309 degradation kinetics; kdeg was found to increase linearly with increasing Tw at constant C and was independent of C at constant Tw, providing a P309 degradation modulus (kdeg/Tw) \(=\) (0.011 \(\pm \) 0.002) (Pa s)??1 for 4 \(<\) Tw Pa \(<\) 400, 0.10 \(<\) C wppm < 5.0. Comparisons between the present degradation kinetics and previous literature showed (kdeg/Tw) data from laboratory pipes of D \(\sim \) 0.01 m to lie on a simple extension of (kdeg/Tw) data from pipelines of D \(\sim \) 0.1 m and 1.0 m, along a power-law relation (kdeg/Tw) \(=\) 10??5.4.D??1.6. Intrinsic slips derived from PAMH B1120 and PEO P309-at-falloff experiments were compared with previous examples from Type B drag reduction by polymers with vinylic and glycosidic backbones, showing: (i) For a given polymer, [S\(^{\prime }\)] was independent of Re f and pipe ID, implying insensitivity to both micro- and macro-scales of turbulence; and (ii) [S\(^{\prime }\)] increased linearly with increasing polymer chain contour length Lc, the proportionality constant \(\beta =\) 0.053 \(\pm \) 0.036 enabling estimation of flow enhancement S\(^{\prime } =\) C.Lc.β for all Type B drag reduction by polymers.  相似文献   

10.
We formulated a paradox in the theory of turbulent premixed flame in the flamelet regime: discrepancy between the Damköhler (1940) and Shelkin (1943) estimate of the turbulence flame speed \(U_{t} \sim {u}^{\prime }\) in the case of strong turbulence (\({u}^{\prime }>>S_{L} \)) and numerous experiments that show a strong dependence of Ut on the speed of the instantaneous flame SL. We name this discrepancy the Damköhler-Shelkin paradox. The first aim of the research is to validate and clarify this estimate, which is based on intuitive considerations, as the paradox must be a statement that seems contradictory to observations but is actually true. We analysed the turbulent flame in the context of the original hyperbolic combustion equation that directly describes the leading edge of the flame, which is a locus of the Zel’dovich “leading points” controlling the speed of the turbulent flame. Analysis of the corresponding characteristic equations results in the expression for speed on the steady-state turbulent flame \(U_{t} ={u}^{\prime }\sqrt {1+(S_{L} /{u}^{\prime })^{2}} \), which is the case when \({u}^{\prime }>>S_{L} \) becomes \(U_{t} \cong {u}^{\prime }\). This result confirms and improves the Damköhler-Shelkin estimate \(U_{t} \sim {u}^{\prime }\). The second aim is to resolve the Damköhler-Shelkin paradox. We explain the discrepancy with observations by the fact that turbulent flames are transient due to insufficient residence time in the real burners to reach statistical equilibrium of wrinkle structures of the random flame surface. We consider the transient flame in the intermediate asymptotic stage when the small-scales wrinkles are in statistical equilibrium, while at the same time the large-scale wrinkles are far from equilibrium. The expressions for the flame speed and width, which we deduce using the dimensional analysis and general properties of the ransom surface, \(U_{t} \sim ({u}^{\prime }S_{L})^{1/2}\) and \(\delta _{t} \sim ({u}^{\prime }Lt)^{1/2}\), show that this transient flame is in fact a turbulent mixing layer travelling with constant speed Ut depending on SL, the intermediate steady propagation (ISP) flame. Qualitative estimations of the times required for the small-scale and large-scale wrinkles to reach statistical equilibrium show that the turbulent Bunsen- and V-flames correspond to the intermediated asymptotic stage, and the turbulent flames with a complete equilibrium structure of the wrinkled flamelet surface are not attainable under laboratory conditions. We present the results of numerical simulations of the impingent flames, which count in favour of the belief that these flames are also transient.  相似文献   

11.
We prove the uniqueness of positive ground state solutions of the problem \({ {\frac {d^{2}u}{dr^{2}}} + {\frac {n-1}{r}}{\frac {du}{dr}} + u \ln(|u|) = 0}\), \({u(r) > 0~\forall r \ge 0}\), and \({(u(r),u'(r)) \to (0, 0)}\) as \({r \to \infty}\). This equation is derived from the logarithmic Schrödinger equation \({{\rm i}\psi_{t} = {\Delta} \psi + u \ln \left(|u|^{2}\right)}\), and also from the classical equation \({{\frac {\partial u}{\partial t}} = {\Delta} u +u \left(|u|^{p-1}\right) -u}\). For each \({n \ge 1}\), a positive ground state solution is \({ u_{0}(r) = \exp \left(-{\frac{r^2}{4}} + {\frac{n}{2}}\right),~0 \le r < \infty}\). We combine \({u_{0}(r)}\) with energy estimates and associated Ricatti equation estimates to prove that, for each \({n \in \left[1, 9 \right]}\), \({u_{0}(r)}\) is the only positive ground state. We also investigate the stability of \({u_{0}(r)}\). Several open problems are stated.  相似文献   

12.
Let \({S\subset\mathbb{R}^2}\) be a bounded Lipschitz domain and denote by \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\) the set of mappings \({u\in W^{2,2}(S;\mathbb{R}^3)}\) which satisfy \({(\nabla u)^T(\nabla u) = Id}\) almost everywhere. Under an additional regularity condition on the boundary \({\partial S}\) (which is satisfied if \({\partial S}\) is piecewise continuously differentiable), we prove that the strong W 2,2 closure of \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)\cap C^{\infty}(\overline{S};\mathbb{R}^3)}\) agrees with \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\).  相似文献   

13.
We consider the temporal homogenization of linear ODEs of the form \({\dot{x}=Ax+\epsilon P(t)x+f(t)}\), where P(t) is periodic and \({\epsilon}\) is small. Using a 2-scale expansion approach, we obtain the long-time approximation \({x(t)\approx {\rm exp}(At) \left( \Omega(t)+\int_0^t {\rm exp}(-A \tau) f(\tau) {\rm d}\tau \right)}\), where \({\Omega}\) solves the cell problem \({\dot{\Omega}=\epsilon B \Omega + \epsilon F(t)}\) with an effective matrix B and an explicitly-known F(t). We provide necessary and sufficient conditions for the accuracy of the approximation (over a \({{\mathcal{O}}(\epsilon^{-1})}\) time-scale), and show how B can be computed (at a cost independent of \({\epsilon}\)). As a direct application, we investigate the possibility of using RLC circuits to harvest the energy contained in small scale oscillations of ambient electromagnetic fields (such as Schumann resonances). Although a RLC circuit parametrically coupled to the field may achieve such energy extraction via parametric resonance, its resistance R needs to be smaller than a threshold \({\kappa}\) proportional to the fluctuations of the field, thereby limiting practical applications. We show that if n RLC circuits are appropriately coupled via mutual capacitances or inductances, then energy extraction can be achieved when the resistance of each circuit is smaller than \({n\kappa}\). Hence, if the resistance of each circuit has a non-zero fixed value, energy extraction can be made possible through the coupling of a sufficiently large number n of circuits (\({n\approx 1000}\) for the first mode of Schumann resonances and contemporary values of capacitances, inductances and resistances). The theory is also applied to the control of the oscillation amplitude of a (damped) oscillator.  相似文献   

14.
A result by Pogorelov asserts that C 1 isometric immersions u of a bounded domain \({S \subset \mathbb R^2}\) into \({\mathbb {R}^3}\) whose normal takes values in a set of zero area enjoy the following regularity property: the gradient \({f := \nabla u}\) is ‘developable’ in the sense that the nondegenerate level sets of f consist of straight line segments intersecting the boundary of S at both endpoints. Motivated by applications in nonlinear elasticity, we study the level set structure of such f when S is an arbitrary bounded Lipschitz domain. We show that f can be approximated by uniformly bounded maps with a simplified level set structure. We also show that the domain S can be decomposed (up to a controlled remainder) into finitely many subdomains, each of which admits a global line of curvature parametrization.  相似文献   

15.
We consider the system Δu ? W u (u) = 0, where \({u : \mathbb{R}^n \to \mathbb{R}^n}\) , for a class of potentials \({W : \mathbb{R}^n \to \mathbb{R}}\) that possess several global minima and are invariant under a general finite reflection group G. We establish existence of nontrivial G-equivariant entire solutions connecting the global minima of W along certain directions at infinity.  相似文献   

16.
In this paper we focused our study on derived from Anosov diffeomorphisms (DA diffeomorphisms ) of the torus \(\mathbb {T}^3,\) it is, an absolute partially hyperbolic diffeomorphism on \(\mathbb {T}^3\) homotopic to a linear Anosov automorphism of the \(\mathbb {T}^3.\) We can prove that if \(f: \mathbb {T}^3 \rightarrow \mathbb {T}^3 \) is a volume preserving DA diffeomorphism homotopic to a linear Anosov A,  such that the center Lyapunov exponent satisfies \(\lambda ^c_f(x) > \lambda ^c_A > 0,\) with x belongs to a positive volume set, then the center foliation of f is non absolutely continuous. We construct a new open class U of non Anosov and volume preserving DA diffeomorphisms, satisfying the property \(\lambda ^c_f(x) > \lambda ^c_A > 0\) for \(m-\)almost everywhere \(x \in \mathbb {T}^3.\) Particularly for every \(f \in U,\) the center foliation of f is non absolutely continuous.  相似文献   

17.
Let \({f: U\rightarrow {\mathbb R}^2}\) be a continuous map, where U is an open subset of \({{\mathbb R}^2}\). We consider a fixed point p of f which is neither a sink nor a source and such that {p} is an isolated invariant set. Under these assumption we prove, using Conley index methods and Nielsen theory, that the sequence of fixed point indices of iterations \({\{{\rm ind}(f^n,p)\}_{n=1}^\infty}\) is periodic, bounded from above by 1, and has infinitely many non-positive terms, which is a generalization of Le Calvez and Yoccoz theorem (Annals of Math., 146, 241–293 (1997)) onto the class of non-injective maps. We apply our result to study the dynamics of continuous maps on 2-dimensional sphere.  相似文献   

18.
Consider a weakly nonlinear CGL equation on the torus \(\mathbb {T}^d\):
$$\begin{aligned} u_t+i\Delta u=\epsilon [\mu (-1)^{m-1}\Delta ^{m} u+b|u|^{2p}u+ ic|u|^{2q}u]. \end{aligned}$$
(*)
Here \(u=u(t,x)\), \(x\in \mathbb {T}^d\), \(0<\epsilon <<1\), \(\mu \geqslant 0\), \(b,c\in \mathbb {R}\) and \(m,p,q\in \mathbb {N}\). Define \(I(u)=(I_{\mathbf {k}},\mathbf {k}\in \mathbb {Z}^d)\), where \(I_{\mathbf {k}}=v_{\mathbf {k}}\bar{v}_{\mathbf {k}}/2\) and \(v_{\mathbf {k}}\), \(\mathbf {k}\in \mathbb {Z}^d\), are the Fourier coefficients of the function \(u\) we give. Assume that the equation \((*)\) is well posed on time intervals of order \(\epsilon ^{-1}\) and its solutions have there a-priori bounds, independent of the small parameter. Let \(u(t,x)\) solve the equation \((*)\). If \(\epsilon \) is small enough, then for \(t\lesssim {\epsilon ^{-1}}\), the quantity \(I(u(t,x))\) can be well described by solutions of an effective equation:
$$\begin{aligned} u_t=\epsilon [\mu (-1)^{m-1}\Delta ^m u+ F(u)], \end{aligned}$$
where the term \(F(u)\) can be constructed through a kind of resonant averaging of the nonlinearity \(b|u|^{2p}+ ic|u|^{2q}u\).
  相似文献   

19.
Let (XG) be a G-action topological dynamical system (t.d.s. for short), where G is a countably infinite discrete amenable group. In this paper, we study the topological pressure of the sets of generic points. We show that when the system satisfies the almost specification property, for any G-invariant measure \(\mu \) and any continuous map \(\varphi \),
$$\begin{aligned} P\left( X_{\mu },\varphi ,\{F_n\}\right) = h_{\mu }(X)+\int \varphi d\mu , \end{aligned}$$
where \(\{F_n\}\) is a Følner sequence, \(X_{\mu }\) is the set of generic points of \(\mu \) with respect to (w.r.t. for short) \(\{F_n\}\), \(P(X_{\mu },\varphi ,\{F_n\})\) is the topological pressure of \(X_{\mu }\) for \(\varphi \) w.r.t. \(\{F_n\}\) and \(h_{\mu }(X)\) is the measure-theoretic entropy.
  相似文献   

20.
V. I. Fabrikant 《Meccanica》2018,53(11-12):2709-2723
The term of generalized material is introduced here as the material, whose state is described by the second order homogeneous differential equations with constant coefficients. The generalized point sources are described by homogeneous expressions, containing the first order derivatives with constant coefficients. The Green’s function for a half-space, made of generalized material, subjected to the action of generalized point sources, is derived in the form of a single integral over a circle. Some of the components of the surface Green’s function are presented in finite form, no computation of any integral is needed. The bonded contact problem is described as mixed–mixed boundary value problem for a half-space, with normal and tangential displacements prescribed inside domain S and vanishing tractions outside this domain on the plane \(x_{3}=0\). A set of governing integral equations for bonded contact problem was derived, with some of the kernels defined in finite form. The general crack problem is defined as that of a flat crack of arbitrary shape S in the plane \(x_{3}=0\), with normal tractions \(\sigma _{33}\) applied symmetrically to the crack faces and tangential tractions \(\sigma _{31}\) and \(\sigma _{23}\) applied anti-symmetrically to the crack faces. A set of 3 governing integral equations is derived, with some of the kernels presented in the finite form. The relationship between the integrands in Fourier transform of the kernels of the sets of the integral equations of both problems is established: they are related in such a way, as if they were the coefficients in the schematic sets of algebraic equations. This relationship can also be presented as a general relationship between matrices of arbitrary rank n, with their components being determinants of certain matrices of rank q, with \(q\ge n\).  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号