Lax Pairs for Equations Describing Compatible Nonlocal Poisson Brackets of Hydrodynamic Type and Integrable Reductions of the Lamé Equations

We solve the problem of describing compatible nonlocal Poisson brackets of hydrodynamic type. We prove that for nonsingular pairs of compatible nonlocal Poisson brackets of hydrodynamic type, there exist special local coordinates such that the metrics and the Weingarten operators of both brackets are diagonal. The nonlinear evolution equations describing all nonsingular pairs of compatible nonlocal Poisson brackets of hydrodynamic type are derived in these special coordinates, and the integrability of these equations is proved using the inverse scattering transform. The Lax pairs with a spectral parameter for these equations are found. We construct various classes of integrable reductions of the derived equations. These classes of reductions are of an independent differential-geometric and applied interest. In particular, if one of the compatible Poisson brackets is local, we obtain integrable reductions of the classical Lamé equations describing all orthogonal curvilinear coordinate systems in a flat space; if one of the compatible brackets is generated by a constant-curvature metric, the corresponding equations describe integrable reductions of the equations for orthogonal curvilinear coordinate systems in a space of constant curvature.     

Thick grids with circular apertures: A comparison of theoretical and experimental performance

In this paper a comparison between a rigorous electromagnetic model for transmission through a hexagonal array of circular waveguides in a series of thick, metallic screens and experimental measurements in the far infrared is made. It is found that there is excellent agreement between theory and experiment when the frequency is below that where any diffracted orders propagate. The agreement is still very good above this frequency. Below a frequency approximately equal to the cut-off frequency of the circular waveguides little power is transmitted. As the thickness of the screen is increased, this decrease in transmission becomes more abrupt. Also, for thick screens, resonances appear in the transmission spectrum which are analogous to those which appear in the spectra of two grids separated by a distance comparable to the wavelength of the radiation being used.     

Moving grid method for numerical simulation of stratified flows

We develop a second‐order accurate Navier–Stokes solver based on r‐adaptivity of the underlying numerical discretization. The motion of the mesh is based on the fluid velocity field; however, certain adjustments to the Lagrangian velocities are introduced to maintain quality of the mesh. The adjustments are based on the variational approach of energy minimization to redistribute grid points closer to the areas of rapid solution variation. To quantify the numerical diffusion inherent to each method, we monitor changes in the background potential energy, computation of which is based on the density field. We demonstrate on a standing interfacial gravity wave simulation how using our method of grid evolution decreases the rate of increase of the background potential energy compared with using the same advection scheme on the stationary grid. To further highlight the benefit of the proposed moving grid method, we apply it to the nonhydrostatic lock‐exchange flow where the evolution of the interface is more complex than in the standing wave test case. Naive grid evolution based on the fluid velocities in the lock‐exchange flow leads to grid tangling as Kelvin–Helmholtz billows develop at the interface. This is remedied by grid refinement using the variational approach. Copyright © 2012 John Wiley & Sons, Ltd.     

A class of finite difference schemes for interface problems with an HOC approach

In this paper, we propose a new methodology for numerically solving elliptic and parabolic equations with discontinuous coefficients and singular source terms. This new scheme is obtained by clubbing a recently developed higher‐order compact methodology with special interface treatment for the points just next to the points of discontinuity. The overall order of accuracy of the scheme is at least second. We first formulate the scheme for one‐dimensional (1D) problems, and then extend it directly to two‐dimensional (2D) problems in polar coordinates. In the process, we also perform convergence and related analysis for both the cases. Finally, we show a new direction of implementing the methodology to 2D problems in cartesian coordinates. We then conduct numerous numerical studies on a number of problems, both for 1D and 2D cases, including the flow past circular cylinder governed by the incompressible Navier–Stokes equations. We compare our results with existing numerical and experimental results. In all the cases, our formulation is found to produce better results on coarser grids. For the circular cylinder problem, the scheme used is seen to capture all the flow characteristics including the famous von Kármán vortex street. Copyright © 2016 John Wiley & Sons, Ltd.     

Generalization of tangential trace spaces of H(curl,Ω) for curvilinear Lipschitz polyhedral domains Ω

For curvilinear Lipschitz polyhedral domains Ω, explicit characterizations of the tangential trace spaces of H ¹(Ω) are presented. These extend the original characterizations given by Buffa and Ciarlet that hold on Lipschitz polyhedral domains with plane faces. The tangential trace spaces of H ¹(Ω) are fundamental for the definition, analysis and intuitive understanding of the trace spaces of H ( curl ,Ω) and therefore, more general characterizations of the latter are obtained at the same time. Copyright © 2013 John Wiley & Sons, Ltd.     

A computational method for geometric optimization of enhanced heat transfer devices based upon entropy generation minimization

A computational fluid dynamics‐based optimization methodology is developed, appropriate for the geometric optimization of enhanced heat transfer devices based upon the principle of entropy generation minimization, in which the objective function is evaluated from a flow field obtained by computational simulation. A quasi‐Newton optimization procedure is employed, with computation of the objective function gradients based upon a finite difference approach. The optimization procedure is developed to be general with regard to the choice of objective function, the details of the problem under consideration, and the computational methodology employed in solving the fluid flow and heat transfer problems. A novel implementation of a Taylor series‐based procedure for the fast solution of nearby problems is presented, which is found to greatly benefit the efficiency of the present methodology. Finally, a numerical experiment is presented, illustrating the use of the present method in the geometric optimization of a practical enhanced heat transfer device on the basis of the criterion of entropy generation minimization. The optimization of the fin spacing of a simple plate fin heat sink is considered, and a comparison of the computational results with results obtained by analytical optimization based upon empirical friction factor and Nusselt number correlations is given. Copyright © 2012 John Wiley & Sons, Ltd.     

A note on compressive limiting for two‐material flows

In this short note we describe a simple extension to the multi‐material shock‐capturing algorithm presented in (J. Comput. Phys. 2007; 223 :262–297) that can be used to maintain sharp material interfaces. The method takes the form of an artificial compression which is designed so that the material indicator jumps across only a few cells but which does not excite physical instabilities in the flow. The advantages of the approach include its simplicity and flexibility in that it provides a parameter that effectively determines the captured interface thickness. Copyright © 2009 John Wiley & Sons, Ltd.     

Magnetic Coupling Constants of Self‐Assembled CuII [3×3] Grids: Alternative Spin Model from Theoretical Calculations

This paper reports a theoretical analysis of the electronic structure and magnetic properties of a ferromagnetic Cu^II [3×3] grid. A two‐step strategy, combining calculations on the whole grid and on binuclear fragments, has been employed to evaluate all the magnetic interactions in the grid. The calculations confirm an S=7/2 ground state, which is in accordance with the magnetisation versus field curve and the thermal dependence of the magnetic moment data. Only the first‐neighbour coupling terms present non‐negligible amplitudes, all of them in agreement with the structure and arrangement of the Cu 3d magnetic orbitals. The results indicate that the dominant interaction in the system is the antiferromagnetic coupling between the ring and the central Cu sites (J₃=J₄≈?31 cm^?1). In the ring two different interactions can be distinguished, J₁=4.6 cm^?1 and J₂=?0.1 cm^?1, in contrast to the single J model employed in the magnetic data fit. The calculated J values have been used to determine the energy level distribution of the Heisenberg magnetic states. The effective magnetic moment versus temperature plot resulting from this ab initio energy profile is in good agreement with the experimental curve and the fitting obtained with the simplified spin model, despite the differences between these two spin models. This study underlines the role that the theoretical evaluations of the coupling constants can play on the rationalisation of the magnetic properties of these complex polynuclear systems.