Bioprinting can be defined as the art of combining materials and cells to fabricate designed, hierarchical 3D hybrid constructs. Suitable materials, so called bioinks, have to comply with challenging rheological processing demands and rapidly form a stable hydrogel postprinting in a cytocompatible manner. Gelatin is often adopted for this purpose, usually modified with (meth-)acryloyl functionalities for postfabrication curing by free radical photopolymerization, resulting in a hydrogel that is cross-linked via nondegradable polymer chains of uncontrolled length. The application of allylated gelatin (GelAGE) as a thiol-ene clickable bioink for distinct biofabrication applications is reported. Curing of this system occurs via dimerization and yields a network with flexible properties that offer a wider biofabrication window than (meth-)acryloyl chemistry, and without additional nondegradable components. An in-depth analysis of GelAGE synthesis is conducted, and standard UV-initiation is further compared with a recently described visible-light-initiator system for GelAGE hydrogel formation. It is demonstrated that GelAGE may serve as a platform bioink for several biofabrication technologies by fabricating constructs with high shape fidelity via lithography-based (digital light processing) 3D printing and extrusion-based 3D bioprinting, the latter supporting long-term viability postprinting of encapsulated chondrocytes.

In the framework of linear stability theory, we analyze how a liquid-gas mixing layer is affected by several parameters: viscosity ratio, density ratio, and several length scales. These scales reflect the presence of a velocity defect induced by the wake behind the splitter plate and the presence of boundary layers which develop ahead of the plate trailing edge. Incorporating such effects, we compute the various temporal and spatial instability modes and identify their driving instability mechanism based on their Reynolds number dependence, spatial structure, and energy budget. It is examined how the velocity defect modifies the temporal and the spatial stability properties. In addition, the transition from convective to absolute instability occurs at lower velocity contrast between gas and liquid free streams when a defect is present. This transition is also promoted by surface tension. Compared to inviscid stability computations, our spatial stability analysis displays a better agreement with measured growth rates obtained in two recent air-water experiments.

We present different simulations of primary atomization using an adaptive Volumeof-Fluid method based on octree meshes. The use of accurate numerical schemes for mesh adaptation, Volume-of-Fluid advection and balanced force surface tension calculation implemented in Gerris, the code used to perform the simulations included in this work, has made possible to carry out accurate simulations with characteristic scales spreading over several orders of magnitude. The code is validated by comparisons with the temporal linear theory for moderate density and viscosity ratios, which basically corresponds to atomization processes in high pressure chambers. In order to show the potential of the code in different scenarios related to atomization, preliminary results are shown in relation with the study of the two-dimensional and 3D temporal and spatial problem, the influence of the injector and the vortex generated inside the chamber, and the effect of swirling at high Reynolds numbers.

Direct numerical simulation is applied to investigate instability and transition to turbulence in the flow of an electrically conducting incompressible fluid between two parallel unbounded insulating walls affected by a wall-normal magnetic field (the Hartmann flow). The linear stability analysis of this flow provided unrealistically high critical Reynolds numbers, about two orders of magnitude higher than those observed in experiments. We propose an explanation based on the streak growth and breakdown mechanism described earlier for other shear flows. The mechanism is investigated using a two-step procedure that includes transient growth of two-dimensional optimal perturbations and the subsequent three-dimensional instability of the modulated streaky flow. In agreement with recent experimental investigations the calculations produce a critical range between 350 and 400 for the Hartmann thickness based Reynolds number, where the transition occurs at realistic amplitudes of two-and three-dimensional perturbations.

The effect of a uniform spanwise magnetic field on a turbulent channel flow is investigated for the case of a low magnetic Reynolds number. Direct numerical simulation ͑DNS͒ and large eddy simulation ͑LES͒ computations are performed for two values of the hydrodynamic Reynolds number ͑10 4 and 2 ϫ 10 4 ͒ and with the Hartmann number varying in a wide range. It is shown that the main effect of the magnetic field is the suppression of turbulent velocity fluctuations and momentum transfer in the wall-normal direction. This leads to drag reduction and transformation of the mean flow profile. The centerline velocity grows, the mean velocity gradients near the wall decrease, and the typical horizontal dimensions of the coherent structures enlarge upon increasing the Hartmann number. Comparison between LES and DNS results shows that the dynamic Smagorinsky model accurately reproduces the flow transformation.

The convective instability in a plane liquid layer with time-dependent temperature profile is investigated by means of a general method suitable for linear stability analysis of an unsteady basic flow. The method is based on a non-normal approach, and predicts the onset of instability, critical wavenumber and time. The method is applied to transient Rayleigh-Bénard-Marangoni convection due to cooling by evaporation. Numerical results as well as theoretical scalings for the critical parameters as function of the Biot number are presented for the limiting cases of purely buoyancy-driven and purely surface-tensiondriven convection. Critical parameters from calculations are in good agreement with those from experiments on drying polymer solutions, where the surface cooling is induced by solvent evaporation.

No abstract

Local dissipation-scale distributions and high-order statistics of the energy dissipation rate are examined in turbulent channel flow using very high-resolution direct numerical simulations at Reynolds numbers Re τ = 180, 381 and 590. For sufficiently large Re τ , the dissipation-scale distributions and energy dissipation moments in the channel bulk flow agree with those in homogeneous isotropic turbulence, including only a weak Reynolds-number dependence of both the finest and largest scales. Systematic, but Re τ -independent, variations in the distributions and moments arise as the wall is approached for y + 100. In the range 100 < y + < 200, there are substantial differences in the moments between the lowest and the two larger values of Re τ . This is most likely caused by coherent vortices from the near-wall region, which fill the whole channel for low Re τ .