skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Arctic Curves of the Twenty-Vertex Model with Domain Wall Boundaries
We use the tangent method to compute the arctic curve of the Twenty-Vertex (20V) model with particular domain wall boundary conditions for a wide set of integrable weights. To this end, we extend to the finite geometry of domain wall boundary conditions the standard connection between the bulk 20V and 6V models via the Kagome lattice ice model. This allows to express refined partition functions of the 20V model in terms of their 6V counterparts, leading to explicit parametric expressions for the various portions of its arctic curve. The latter displays a large variety of shapes depending on the weights and separates a central liquid phase from up to six different frozen phases. A number of numerical simulations are also presented, which highlight the arctic curve phenomenon and corroborate perfectly the analytic predictions of the tangent method. We finally compute the arctic curve of the Quarter Turn symmetric Holey Aztec Domino Tiling (QTHADT) model, a problem closely related to the 20V model and whose asymptotics may be analyzed via a similar tangent method approach. Again results for the QTHADT model are found to be in perfect agreement with our numerical simulations.  more » « less
Award ID(s):
1802044 1937241
PAR ID:
10142527
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Statistical Physics
ISSN:
0022-4715
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Recent experimental and computational studies indicate that near-wall turbulent flows can be characterized by universal small-scale autonomous dynamics that is modulated by large-scale structures. We formulate numerical simulations of near-wall turbulence in a small domain localized to the boundary, whose size scales in viscous units. To mimic the environment in which the near-wall turbulence evolves, the formulation accounts for the flux of mean momentum through the upper boundary of the domain. Comparisons of the model's two-dimensional energy spectra and low-order single-point statistics with the corresponding quantities computed from direct numerical simulations indicate that it successfully captures the dynamics of the small-scale near-wall turbulence. 
    more » « less
  2. In this study, we conduct a parametric analysis to evaluate the sensitivities of wall-modeled large-eddy simulation (LES) with respect to subgrid-scale (SGS) models, mesh resolution, wall boundary conditions and mesh anisotropy. While such investigations have been conducted for attached/flat-plate flow configurations, systematic studies specifically targeting turbulent flows with separation are notably sparse. To bridge this gap, our study focuses on the flow over a two-dimensional Gaussian-shaped bump at a moderately high Reynolds number, which involves smooth-body separation of a turbulent boundary layer under pressure-gradient and surface- curvature effects. In the simulations, the no-slip condition at the wall is replaced by three different forms of boundary condition based on the thin boundary layer equations and the mean wall-shear stress from high-fidelity numerical simulation to avoid the additional complexity of modeling the wall-shear stress. Various statistics, including the mean separation bubble size, mean velocity profile, and dissipation from SGS model, are compared and analyzed. The results reveal that capturing the separation bubble strongly depends on the choice of SGS model. While simulations approach grid convergence with resolutions nearing those of wall-resolved LES meshes, above this limit, the LES predictions exhibit intricate sensitivities to mesh resolution. Furthermore, both wall boundary conditions and the anisotropy of mesh cells exert discernible impacts on the turbulent flow predictions, yet the magnitudes of these impacts vary based on the specific SGS model chosen for the simulation. 
    more » « less
  3. Bone healing has been traditionally described as a four-phase process: inflammatory response, soft callus formation, hard callus development, and remodeling. The remodeling phase has been largely neglected in most numerical mechanoregulation models of fracture repair in favor of capturing early healing using a pre-defined callus domain. However, in vivo evidence suggests that remodeling occurs concurrently with repair and causes changes in cortical bone adjacent to callus that are typically neglected in numerical models of bone healing. The objective of this study was to use image processing techniques to quantify this early-stage remodeling in ovine osteotomies. To accomplish this, we developed a numerical method for radiodensity profilometry with optimization-based curve fitting to mathematically model the bone density gradients in the radial direction across the cortical wall and callus. After assessing data from 26 sheep, we defined a dimensionless density fitting function that revealed significant remodeling occurring in the cortical wall adjacent to callus during early healing, a 23% average reduction in density compared to intact. This fitting function is robust for modeling radial density gradients in both intact bone and fracture repair scenarios and can capture a wide variety of the healing responses. The fitting function can also be scaled easily for comparison to numerical model predictions and may be useful for validating future mechanoregulatory models of coupled fracture repair and remodeling. 
    more » « less
  4. Abstract We consider the finite element approximation of a coupled fluid‐structure interaction (FSI) system, which comprises a three‐dimensional (3D) Stokes flow and a two‐dimensional (2D) fourth‐order Euler–Bernoulli or Kirchhoff plate. The interaction of these parabolic and hyperbolic partial differential equations (PDE) occurs at the boundary interface which is assumed to be fixed. The vertical displacement of the plate dynamics evolves on the flat portion of the boundary where the coupling conditions are implemented via the matching velocities of the plate and fluid flow, as well as the Dirichlet boundary trace of the pressure. This pressure term also acts as a coupling agent, since it appears as a forcing term on the flat, elastic plate domain. Our main focus in this work is to generate some numerical results concerning the approximate solutions to the FSI model. For this, we propose a numerical algorithm that sequentially solves the fluid and plate subsystems through an effective decoupling approach. Numerical results of test problems are presented to illustrate the performance of the proposed method. 
    more » « less
  5. Abstract A simple, flow-physics-based model of flat-plate, transitional boundary layer skin friction and heat transfer is presented. The model is based on the assumption of negligible time-, spanwise-, and streamwise-average wall-normal velocity at the top of the boundary layer. This results in a threefold increase in boundary layer thickness over the transition region. This simple velocity assumption and its boundary-layer growth implications seem to be reasonably consistent with more sophisticated (direct numerical simulation (DNS)) modeling simulations. Only two modeling parameters need to be assumed, the Reynolds numbers at the onset and at the completion of transition, for which there is guidance based on freestream turbulence intensity for smooth plates. Several experimental datasets for air are modeled. New criteria are proposed to help define the onset and completion of transition: zero net vertical (wall-normal) velocity or mass flux (integrated in time and space, spanwise and streamwise) at the top of the boundary layer, and tripling of boundary layer thickness. Also presented is a minor improvement to a previously published unheated starting length factor for flat-plate laminar boundary layers with uniform wall heat flux. 
    more » « less