Computational fluid dynamic analysis of bioprinted self-supporting perfused tissue models

Natural tissues are incorporated with vasculature, which is further integrated with a cardiovascular system responsible for driving perfusion of nutrient‐rich oxygenated blood through the vasculature to support cell metabolism within most cell‐dense tissues. Since scaffold‐free biofabricated tissues being developed into clinical implants, research models, and pharmaceutical testing platforms should similarly exhibit perfused tissue‐like structures, we generated a generalizable biofabrication method resulting in self‐supporting perfused (SSuPer) tissue constructs incorporated with perfusible microchannels and integrated with the modular FABRICA perfusion bioreactor. As proof of concept, we perfused an MLO‐A5 osteoblast‐based SSuPer tissue in the FABRICA. Although our resulting SSuPer tissue replicated vascularization and perfusion observed in situ, supported its own weight, and stained positively for mineral using Von Kossa staining, our in vitro results indicated that computational fluid dynamics (CFD) should be used to drive future construct design and flow application before further tissue biofabrication and perfusion. We built a CFD model of the SSuPer tissue integrated in the FABRICA and analyzed flow characteristics (net force, pressure distribution, shear stress, and oxygen distribution) through five SSuPer tissue microchannel patterns in two flow directions and at increasing flow rates. Important flow parameters include flow direction, fully developed flow, and tissue microchannel diameters matched and aligned with bioreactor flow channels. We observed that the SSuPer tissue platform is capable of providing direct perfusion to tissue constructs and proper culture conditions (oxygenation, with controllable shear and flow rates), indicating that our approach can be used to biofabricate tissue representing primary tissues and that we can model the system in silico

Evolution of Magnetism in Single-Crystal Honeycomb Iridates

We report the successful synthesis of single-crystals of the layered iridate, (Na$_{1-x}$Li$_{x}$)$_2$IrO$_3$, $0\leq x \leq 0.9$, and a thorough study of its structural, magnetic, thermal and transport properties. The new compound allows a controlled interpolation between Na$_2$IrO$_3$ and Li$_2$IrO$_3$, while maintaing the novel quantum magnetism of the honeycomb Ir$^{4+}$ planes. The measured phase diagram demonstrates a dramatic suppression of the N\'eel temperature, $T_N$, at intermediate $x$ suggesting that the magnetic order in Na$_2$IrO$_3$ and Li$_2$IrO$_3$ are distinct, and that at $x\approx 0.7$, the compound is close to a magnetically disordered phase that has been sought after in Na$_2$IrO$_3$ and Li$_2$IrO$_3$. By analyzing our magnetic data with a simple theoretical model we also show that the trigonal splitting, on the Ir$^{4+}$ ions changes sign from Na$_2$IrO$_3$ and Li$_2$IrO$_3$, and the honeycomb iridates are in the strong spin-orbit coupling regime, controlled by \jeff=1/2 moments.Comment: updated version with more dat

KKT conditions satisfied using adaptive neighboring in hybrid cellular automata for topology optimization

The hybrid cellular automaton (HCA) method is a biologically inspired algorithm capable of topology synthesis that was developed to simulate the behavior of the bone functional adaptation process. In this algorithm, the design domain is divided into cells with some communication property among neighbors. Local evolutionary rules, obtained from classical control theory, iteratively establish the value of the design variables in order to minimize the local error between a field variable and a corresponding target value. Karush-Kuhn-Tucker (KKT) optimality conditions have been derived to determine the expression for the field variable and its target. While averaging techniques mimicking intercellular communication have been used to mitigate numerical instabilities such as checkerboard patterns and mesh dependency, some questions have been raised whether KKT conditions are fully satisfied in the final topologies. Furthermore, the averaging procedure might result in cancellation or attenuation of the error between the field variable and its target. Several examples are presented showing that HCA converges to different final designs for different neighborhood configurations or averaging schemes. Although it has been claimed that these final designs are optimal, this might not be true in a precise mathematical sense—the use of the averaging procedure induces a mathematical incorrectness that has to be addressed. In this work, a new adaptive neighboring scheme will be employed that utilizes a weighting function for the influence of a cell’s neighbors that decreases to zero over time. When the weighting function reaches zero, the algorithm satisfies the aforementioned optimality criterion. Thus, the HCA algorithm will retain the benefits that result from utilizing neighborhood information, as well as obtain an optimal solution

Convergence analysis of hybrid cellular automata for topology optimization

The hybrid cellular automaton (HCA) algorithm was inspired by the structural adaptation of bones to their ever changing mechanical environment. This methodology has been shown to be an effective topology synthesis tool. In previous work, it has been observed that the convergence of the HCA methodology is affected by parameters of the algorithm. As a result, questions have been raised regarding the conditions by which HCA converges to an optimal design. The objective of this investigation is to examine the conditions that guarantee convergence to a Karush-Kuhn-Tucker (KKT) point. In this paper, it is shown that the HCA algorithm is a fixed point iterative scheme and the previously reported KKT optimality conditions are corrected. To demonstrate the convergence properties of the HCA algorithm, a simple cantilevered beam example is utilized. Plots of the spectral radius for projections of the design space are used to show regions of guaranteed convergence

Spin waves in ultrathin ferromagnetic overlayers

The influence of a non-magnetic metallic substrate on the spin wave excitations in ultrathin ferromagnetic overlayers is investigated for different crystalline orientations. We show that spin wave dumping in these systems occur due to the tunneling of holes from the substrate into the overlayer, and that the spin wave energies may be considerably affected by the exchange coupling mediated by the substrate.Comment: RevTeX 4, 7 pages, 5 figures; submitted to Phys. Rev.