Nonlinear preconditioning for efficient and accurate interface capturing in simulation of multicomponent compressible flows

Single fluid schemes that rely on an interface function for phase identification in multicomponent compressible flows are widely used to study hydrodynamic flow phenomena in several diverse applications. Simulations based on standard numerical implementation of these schemes suffer from an artificial increase in the width of the interface function owing to the numerical dissipation introduced by an upwind discretization of the governing equations. In addition, monotonicity requirements which ensure that the sharp interface function remains bounded at all times necessitate use of low-order accurate discretization strategies. This results in a significant reduction in accuracy along with a loss of intricate flow features. In this paper we develop a nonlinear transformation based interface capturing method which achieves superior accuracy without compromising the simplicity, computational efficiency and robustness of the original flow solver. A nonlinear map from the signed distance function to the sigmoid type interface function is used to effectively couple a standard single fluid shock and interface capturing scheme with a high-order accurate constrained level set reinitialization method in a way that allows for oscillation-free transport of the sharp material interface. Imposition of a maximum principle, which ensures that the multidimensional preconditioned interface capturing method does not produce new maxima or minima even in the extreme events of interface merger or breakup, allows for an explicit determination of the interface thickness in terms of the grid spacing. A narrow band method is formulated in order to localize computations pertinent to the preconditioned interface capturing method. Numerical tests in one dimension reveal a significant improvement in accuracy and convergence; in stark contrast to the conventional scheme, the proposed method retains its accuracy and convergence characteristics in a shifted reference frame. Results from the test cases in two dimensions show that the nonlinear transformation based interface capturing method outperforms both the conventional method and an interface capturing method without nonlinear transformation in resolving intricate flow features such as sheet jetting in the shock-induced cavity collapse. The ability of the proposed method in accounting for the gravitational and surface tension forces besides compressibility is demonstrated through a model fully three-dimensional problem concerning droplet splash and formation of a crownlike feature. (C) 2014 Elsevier Inc. All rights reserved.

Spiral-wave dynamics in ionically realistic mathematical models for human ventricular tissue: the effects of periodic deformation

We carry out an extensive numerical study of the dynamics of spiral waves of electrical activation, in the presence of periodic deformation (PD) in two-dimensional simulation domains, in the biophysically realistic mathematical models of human ventricular tissue due to (a) ten-Tusscher and Panfilov (the TP06 model) and (b) ten-Tusscher, Noble, Noble, and Panfilov (the TNNPO4 model). We first consider simulations in cable-type domains, in which we calculate the conduction velocity theta and the wavelength lambda of a plane wave; we show that PD leads to a periodic, spatial modulation of theta and a temporally periodic modulation of lambda; both these modulations depend on the amplitude and frequency of the PD. We then examine three types of initial conditions for both TP06 and TNNPO4 models and show that the imposition of PD leads to a rich variety of spatiotemporal patterns in the transmembrane potential including states with a single rotating spiral (RS) wave, a spiral-turbulence (ST) state with a single meandering spiral, an ST state with multiple broken spirals, and a state SA in which all spirals are absorbed at the boundaries of our simulation domain. We find, for both TP06 and TNNPO4 models, that spiral-wave dynamics depends sensitively on the amplitude and frequency of PD and the initial condition. We examine how these different types of spiral-wave states can be eliminated in the presence of PD by the application of low-amplitude pulses by square- and rectangular-mesh suppression techniques. We suggest specific experiments that can test the results of our simulations.

Soft-spring wall based non-periodic boundary conditions for non-equilibrium molecular dynamics of dense fluids

Non-equilibrium molecular dynamics (MD) simulations require imposition of non-periodic boundary conditions (NPBCs) that seamlessly account for the effect of the truncated bulk region on the simulated MD region. Standard implementation of specular boundary conditions in such simulations results in spurious density and force fluctuations near the domain boundary and is therefore inappropriate for coupled atomistic-continuum calculations. In this work, we present a novel NPBC model that relies on boundary atoms attached to a simple cubic lattice with soft springs to account for interactions from particles which would have been present in an untruncated full domain treatment. We show that the proposed model suppresses the unphysical fluctuations in the density to less than 1% of the mean while simultaneously eliminating spurious oscillations in both mean and boundary forces. The model allows for an effective coupling of atomistic and continuum solvers as demonstrated through multiscale simulation of boundary driven singular flow in a cavity. The geometric flexibility of the model enables straightforward extension to nonplanar complex domains without any adverse effects on dynamic properties such as the diffusion coefficient. (c) 2015 AIP Publishing LLC.

Conceptual design of three-dimensional scaffolds of powder-based materials for bone tissue engineering applications

Purpose - The purpose of this paper is to investigate the possibility to construct tissue-engineered bone repair scaffolds with pore size distributions using rapid prototyping techniques. Design/methodology/approach - The fabrication of porous scaffolds with complex porous architectures represents a major challenge in tissue engineering and the design aspects to mimic complex pore shape as well as spatial distribution of pore sizes of natural hard tissue remain unexplored. In this context, this work aims to evaluate the three-dimensional printing process to study its potential for scaffold fabrication as well as some innovative design of homogeneously porous or gradient porous scaffolds is described and such design has wider implication in the field of bone tissue engineering. Findings - The present work discusses biomedically relevant various design strategies with spatial/radial gradient in pore sizes as well as with different pore sizes and with different pore geometries. Originality/value - One of the important implications of the proposed novel design scheme would be the development of porous bioactive/biodegradable composites with gradient pore size, porosity, composition and with spatially distributed biochemical stimuli so that stem cells loaded into scaffolds would develop into complex tissues such as those at the bone-cartilage interface.

Strength reliability and in vitro degradation of three-dimensional powder printed strontium-substituted magnesium phosphate scaffolds

Strontium ions (Sr2+) are known to prevent osteoporosis and also encourage bone formation. Such twin requirements have motivated researchers to develop Sr-substituted biomaterials for orthopaedic applications. The present study demonstrates a new concept of developing Sr-substituted Mg-3(PO4)(2) - based biodegradable scaffolds. In particular, this work reports the fabrication, mechanical properties with an emphasis on strength reliability as well as in vitro degradation of highly biodegradable strontium-incorporated magnesium phosphate cements. These implantable scaffolds were fabricated using three-dimensional powder printing, followed by high temperature sintering and/or chemical conversion, a technique adaptable to develop patient-specific implants. A moderate combination of strength properties of 36.7 MPa (compression), 242 MPa (bending) and 10.7 MPa (tension) were measured. A reasonably modest Weibull modulus of up to 8.8 was recorded after uniaxial compression or diametral tensile tests on 3D printed scaffolds. A comparison among scaffolds with varying compositions or among sintered or chemically hardened scaffolds reveals that the strength reliability is not compromised in Sr-substituted scaffolds compared to baseline Mg-3(PO4)(2). The micro-computed tomography analysis reveals the presence of highly interconnected porous architecture in three-dimension with lognormal pore size distribution having median in the range of 17.74-26.29 mu m for the investigated scaffolds. The results of extensive in vitro ion release study revealed passive degradation with a reduced Mg2+ release and slow but sustained release of Sr2+ from strontium-substituted magnesium phosphate scaffolds. Taken together, the present study unequivocally illustrates that the newly designed Sr-substituted magnesium phosphate scaffolds with good strength reliability could be used for biomedical applications requiring consistent Sr2+-release, while the scaffold degrades in physiological medium. Statement of significance The study investigates the additive manufacturing of scaffolds based on different strontium-substituted magnesium phosphate bone cements by means of three-dimensional powder printing technique (3DPP). Magnesium phosphates were chosen due to their higher biodegradability compared to calcium phosphates, which is due to both a higher solubility as well as the absence of phase changes (to low soluble hydroxyapatite) in vivo. Since strontium ions are known to promote bone formation by stimulating osteoblast growth, we aimed to establish such a highly degradable magnesium phosphate ceramic with an enhanced bioactivity for new bone ingrowth. After post-processing, mechanical strengths of up to 36.7 MPa (compression), 24.2 MPa (bending) and 10.7 MPa (tension) could be achieved. Simultaneously, the failure reliability of those bioceramic implant materials, measured by Weibull modulus calculations, were in the range of 4.3-8.8. Passive dissolution studies in vitro proved an ion release of Mg2+ and PO43- as well as Sr2+, which is fundamental for in vivo degradation and a bone growth promoting effect. In our opinion, this work broadens the range of bioceramic bone replacement materials suitable for additive manufacturing processing. The high biodegradability of MPC ceramics together with the anticipated promoting effect on osseointegration opens up the way for a patient-specific treatment with the prospect of a fast and complete healing of bone fractures. (C) 2015 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.