Morsingortiz7933

We propose a procedure to implement Dirichlet velocity boundary conditions for complex shapes that use data from a single node only, in the context of the lattice Boltzmann method. Two ideas are at the base of this approach. The first is to generalize the geometrical description of boundary conditions combining bounce-back rule with interpolations. The second is to enhance them by limiting the interpolation extension to the proximity of the boundary. Despite its local nature, the resulting method exhibits second-order convergence for the velocity field and shows similar or better accuracy than the well-established Bouzidi's scheme for curved walls [M. Bouzidi, M. Firdaouss, and P. Lallemand, Phys. Fluids 13, 3452 (2001)]PHFLE61070-663110.1063/1.1399290. Among the infinite number of possibilities, we identify several meaningful variants of the method, discerned by their approximation of the second-order nonequilibrium terms and their interpolation coefficients. For each one, we provide two parametrized versions that produce viscosity independent accuracy at steady state. The method proves to be suitable to simulate moving rigid objects or surfaces moving following either the rigid body dynamics or a prescribed kinematic. Also, it applies uniformly and without modifications in the whole domain for any shape, including corners, narrow gaps, or any other singular geometry.The lattice-Boltzmann method is convenient for simulating flow fields in porous media. Siponimod in vivo However, due to its lattice characteristics, the velocity near a solid surface is not accurate, which results in significant errors when simulating colloid transport in porous media. Based on the general properties of a flow field close to a solid surface, we propose an alternative velocity interpolation method in which the velocity at a solid surface is strictly zero. Numerical simulation results show that the proposed method can give more accurate results than the usual bilinear interpolation. In addition, we use this method to simulate the contact efficiency of colloids in porous media and obtain a new power-law form of the contact efficiency.We study the wetting critical behavior of the three-state (s=±1,0) Blume-Emery-Griffiths model using numerical simulations. This model provides a suitable scenario for the study of the role of vacancies on the wetting behavior of a thin magnetic film. To this aim we study a system confined between parallel walls with competitive short-range surface magnetic fields (h_L=-|h_1|). We locate relevant critical curves for different values of the biquadratic interaction and use a thermodynamic integration method to calculate the surface tension as well as the interfacial excess energy and determine the wetting transition. Furthermore, we also calculate the local position of the interface along the film and its fluctuations (capillary waves), which are a measure of the interface width. To characterize the role played by vacancies on the interfacial behavior we evaluate the excess density of vacancies, i.e., the density difference between a system with and without interface. We also show that the temperature dependence of both the local position of the interface and its width can be rationalized in term of a finite-size scaling description, and we propose and successfully test the same scaling behavior for the average position of the center of mass of the vacancies and its fluctuations. This shows that the excess of vacancies can be associated to the presence of the interface that causes the observed segregation. This segregation phenomena is also evidenced by explicitly evaluating the interfacial free energy.Chimera states refer to the dynamical states in which the inherent symmetry of the system is broken. The system composed of two interacting identical subpopulations of phase oscillators provides a platform to study chimera states. In this system, different types of chimera states have been identified and the transitions between them have been investigated. However, the parameter space is not fully explored in this system. In this work, we study a system comprised of two interacting subpopulations of nonidentical phase oscillators. Through numerical simulations and theoretical analyses, we find three symmetry-reserving states, including incoherent state, in-phase synchronous state, and antiphase synchronous state, and three types of symmetry-breaking states, including in-phase chimera states, antiphase chimera states, and weak chimera states. The stability diagrams of these dynamical states are explored on different parameter planes and transition scenarios amongst these states are investigated. We find that the weak chimera states act as the bridge between in-phase and antiphase chimera states. We also observe the existence of a period-two chimera state, chaotic chimera state, and drifting chimera states.We study the relation between stochastic thermodynamics and nonequilibrium thermodynamics by evaluating the entropy production and the relation between fluxes and forces in a harmonic system with N particles in contact with N different reservoirs. We suppose that the system is in a nonequilibrium stationary state in a first phase and we study the relaxation to equilibrium in a second phase. During this relaxation, we can identify the linear relation between fluxes and forces satisfying the Onsager reciprocity and we obtain a nonlinear expression for the entropy production. Only when forces and fluxes are small does the entropic production turn into a quadratic form in the forces, as predicted by the Onsager theory.Three-dimensional (3D) instabilities on a (potentially turbulent) two-dimensional (2D) flow are still incompletely understood, despite recent progress. Here, based on known physical properties of such 3D instabilities, we propose a simple, energy-conserving model describing this situation. It consists of a regularized 2D point-vortex flow coupled to localized 3D perturbations ("ergophages"), such that ergophages can gain energy by altering vortex-vortex distances through an induced divergent velocity field, thus decreasing point-vortex energy. We investigate the model in three distinct stages of evolution (i) The linear regime, where the amplitude of the ergophages grows or decays exponentially on average, with an instantaneous growth rate that fluctuates randomly in time. The instantaneous growth rate has a small auto-correlation time, and a probability distribution featuring a power-law tail with exponent between -2 and -5/3 (up to a cutoff) depending on the point-vortex base flow. Consequently, the logarithm of the ergophage amplitude performs a Lévy flight.