Problems on control for a steady-state magnetic-hydrodynamic model of a viscous heat-conducting fluid under mixed boundary conditions

The control problems for steady-state equations of magnetic hydrodynamics (MHD) for a viscous heat-conducting fluid considered under mixed boundary conditions for the magnetic field and temperature are investigated. Their solvability is proved, the optimality systems describing the necessary conditions of an extremum are derived, and the theorems of local uniqueness and stability of the optimum solutions for explicit quality functionals are formulated.     

Initial boundary value problems for the equations of motion of compressible viscous and heat-conductive fluids

The equations of motion of compressible viscous and heat-conductive fluids are investigated for initial boundary value problems on the half space and on the exterior domain of any bounded region. The global solution in time is proved to exist uniquely and approach the stationary state ast→∞, provided the prescribed initial data and the external force are sufficiently small.     

Extremal points for fractional boundary value problems

This article is concerned with characterizing the first extremal point, b₀, for a Riemann–Liouville fractional boundary value problem, D^α₀₊y + p(t)y = 0, 0 < t < b, y(0) = y′(0) = y″(b) = 0, 2 < α ≤ 3, by applying the theory of u₀-positive operators with respect to a suitable cone in a Banach space. The key argument is that a mapping, which maps a linear, compact operator, depending on b to its spectral radius, is continuous and strictly increasing as a function of b. Furthermore, an application to a nonlinear case is given.     

Solvability of the inhomogeneous mixed boundary value problem for stationary magnetohydrodynamic equations

Global solvability of the boundary value problem for stationary magnetohydrodynamic equations considered under inhomogeneous mixed boundary conditions for a magnetic field and the inhomogeneous Dirichlet condition for the velocity is proved.     

A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

A minimum dissipation principle for the Rayleigh-Taylor problem in viscous magnetohydrodynamics

A variational principle for magnetohydrodynamic gravitational modes, including isotropic viscosity and the shear of magnetic field lines, is given. It is shown to imply a minimum dissipation requirement for the modes. A sufficient condition for stability given in the literature for the case of ordinary fluids is corrected here.     

A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in 2D

We present a new method for the evolution of inextensible vesicles immersed in a Stokesian fluid. We use a boundary integral formulation for the fluid that results in a set of nonlinear integro-differential equations for the vesicle dynamics. The motion of the vesicles is determined by balancing the non-local hydrodynamic forces with the elastic forces due to bending and tension. Numerical simulations of such vesicle motions are quite challenging. On one hand, explicit time-stepping schemes suffer from a severe stability constraint due to the stiffness related to high-order spatial derivatives and a milder constraint due to a transport-like stability condition. On the other hand, an implicit scheme can be expensive because it requires the solution of a set of nonlinear equations at each time step. We present two semi-implicit schemes that circumvent the severe stability constraints on the time step and whose computational cost per time step is comparable to that of an explicit scheme. We discretize the equations by using a spectral method in space, and a multistep third-order accurate scheme in time. We use the fast multipole method (FMM) to efficiently compute vesicle–vesicle interaction forces in a suspension with a large number of vesicles. We report results from numerical experiments that demonstrate the convergence and algorithmic complexity properties of our scheme.     

Dynamic model and stationary shapes of fluid vesicles

A phase-field model that takes into account the bending energy of fluid vesicles is presented. The Canham-Helfrich model is derived in the sharp-interface limit. A dynamic equation for the phase-field has been solved numerically to find stationary shapes of vesicles with different topologies and the dynamic evolution towards them. The results are in agreement with those found by minimization of the Canham-Helfrich free energy. This fact shows that our phase-field model could be applied to more complex problems of instabilities.     

On some boundary value problems on a strip in the complex plane

The Sommerfeld integral inversion method for the Helmholtz equation in an angular region with different boundary values leads to boundary value problems in an infinite strip of the complex plane. We investigate a generic system for such boundary value problems and give the existence and uniqueness results with optimal growth estimates on the solution. We also give the solutions of the Dirichlet problem in a strip when the boundary functions grow exponentially.     

A variational principle for boundary value problems in kinetic theory

A variational principle which applies directly to the integrodifferential form of the linearized Boltzmann equation is introduced. Extremely general boundary conditions and collision terms are allowed. For a class of interesting problems, the value of the functional to be varied is shown to be closely related to quantities of great physical interest. The formalism is applied to the treatment of plane Couette flow for different forms of the collision term (BGK model, rigid spheres, Maxwell's molecules).Research sponsored by the Air Force Office of Scientific Research under contract F 61(052)-68-C-0020, through the European Office of Aerospace Research, OAR, United States Air Force.     

Method for solving moving boundary value problems for linear evolution equations

We introduce a method of solving initial boundary value problems for linear evolution equations in a time-dependent domain, and we apply it to an equation with dispersion relation omega(k), in the domain l(t)相似文献