Jacobian-free Newton-Krylov methods with GPU acceleration for computing nonlinear ship wave patterns

The nonlinear problem of steady free-surface flow past a submerged source is considered as a case study for three-dimensional ship wave problems. Of particular interest is the distinctive wedge-shaped wave pattern that forms on the surface of the fluid. By reformulating the governing equations with a standard boundary-integral method, we derive a system of nonlinear algebraic equations that enforce a singular integro-differential equation at each midpoint on a two-dimensional mesh. Our contribution is to solve the system of equations with a Jacobian-free Newton-Krylov method together with a banded preconditioner that is carefully constructed with entries taken from the Jacobian of the linearised problem. Further, we are able to utilise graphics processing unit acceleration to significantly increase the grid refinement and decrease the run-time of our solutions in comparison to schemes that are presently employed in the literature. Our approach provides opportunities to explore the nonlinear features of three-dimensional ship wave patterns, such as the shape of steep waves close to their limiting configuration, in a manner that has been possible in the two-dimensional analogue for some time.

A finite volume method for solving the two-sided time-space fractional advection-dispersion equation

We present a finite volume method to solve the time-space two-sided fractional advection-dispersion equation on a one-dimensional domain. The spatial discretisation employs fractionally-shifted Grünwald formulas to discretise the Riemann-Liouville fractional derivatives at control volume faces in terms of function values at the nodes. We demonstrate how the finite volume formulation provides a natural, convenient and accurate means of discretising this equation in conservative form, compared to using a conventional finite difference approach. Results of numerical experiments are presented to demonstrate the effectiveness of the approach.

The resistance of PRESENT-80 against related-key differential attacks

We examine the security of the 64-bit lightweight block cipher PRESENT-80 against related-key differential attacks. With a computer search we are able to prove that for any related-key differential characteristic on full-round PRESENT-80, the probability of the characteristic only in the 64-bit state is not higher than 2−64. To overcome the exponential (in the state and key sizes) computational complexity of the search we use truncated differences, however as the key schedule is not nibble oriented, we switch to actual differences and apply early abort techniques to prune the tree-based search. With a new method called extended split approach we are able to make the whole search feasible and we implement and run it in real time. Our approach targets the PRESENT-80 cipher however,with small modifications can be reused for other lightweight ciphers as well.

Truncated differential analysis of reduced-round LBlock

In this paper we present truncated differential analysis of reduced-round LBlock by computing the differential distribution of every nibble of the state. LLR statistical test is used as a tool to apply the distinguishing and key-recovery attacks. To build the distinguisher, all possible differences are traced through the cipher and the truncated differential probability distribution is determined for every output nibble. We concatenate additional rounds to the beginning and end of the truncated differential distribution to apply the key-recovery attack. By exploiting properties of the key schedule, we obtain a large overlap of key bits used in the beginning and final rounds. This allows us to significantly increase the differential probabilities and hence reduce the attack complexity. We validate the analysis by implementing the attack on LBlock reduced to 12 rounds. Finally, we apply single-key and related-key attacks on 18 and 21-round LBlock, respectively.

Poly(2-oxazoline) hydrogels as next generation three dimensional cell supports

Synthetic hydrogels selectively decorated with cell adhesion motifs are rapidly emerging as promising substrates for 3D cell culture. When cells are grown in 3D they experience potentially more physiologically relevant cell-cell interactions and physical cues compared with traditional 2D cell culture on stiff surfaces. A newly developed polymer based on poly(2-oxazoline)s has been used for the first time to control attachment of fibroblast cells and is discussed here for its potential use in 3D cell culture with particular focus on cancer cells towards the ultimate aim of high throughput screening of anti-cancer therapies. Advantages and limitations of using poly(2-oxazoline) hydrogels are discussed and compared with more established polymers, especially polyethylene glycol (PEG).

Gelatine methacrylamide-based hydrogels : an alternative three-dimensional cancer cell culture system

Modern cancer research requires physiological, three-dimensional (3-D) cell culture platforms, wherein the physical and chemical characteristics of the extracellular matrix (ECM) can be modified. In this study, gelatine methacrylamide (GelMA)-based hydrogels were characterized and established as in vitro and in vivo spheroid-based models for ovarian cancer, reflecting the advanced disease stage of patients, with accumulation of multicellular spheroids in the tumour fluid (ascites). Polymer concentration (2.5-7% w/v) strongly influenced hydrogel stiffness (0.5±0.2kPa to 9.0±1.8kPa) but had little effect on solute diffusion. The diffusion coefficient of 70kDa fluorescein isothiocyanate (FITC)-labelled dextran in 7% GelMA-based hydrogels was only 2.3 times slower compared to water. Hydrogels of medium concentration (5% w/v GelMA) and stiffness (3.4kPa) allowed spheroid formation and high proliferation and metabolic rates. The inhibition of matrix metalloproteinases and consequently ECM degradability reduced spheroid formation and proliferation rates. The incorporation of the ECM components laminin-411 and hyaluronic acid further stimulated spheroid growth within GelMA-based hydrogels. The feasibility of pre-cultured GelMA-based hydrogels as spheroid carriers within an ovarian cancer animal model was proven and led to tumour development and metastasis. These tumours were sensitive to treatment with the anti-cancer drug paclitaxel, but not the integrin antagonist ATN-161. While paclitaxel and its combination with ATN-161 resulted in a treatment response of 33-37.8%, ATN-161 alone had no effect on tumour growth and peritoneal spread. The semi-synthetic biomaterial GelMA combines relevant natural cues with tunable properties, providing an alternative, bioengineered 3-D cancer cell culture in in vitro and in vivo model systems.

Known and chosen key differential distinguishers for block ciphers

In this paper we investigate the differential properties of block ciphers in hash function modes of operation. First we show the impact of differential trails for block ciphers on collision attacks for various hash function constructions based on block ciphers. Further, we prove the lower bound for finding a pair that follows some truncated differential in case of a random permutation. Then we present open-key differential distinguishers for some well known round-reduced block ciphers.

Chemiluminescence imaging in supersonic combustors operating in radical-farming mode

Emission spectroscopy was used to investigate ignition and combustion characteristics of supersonic combustion ramjet engines. Two-dimensional scramjet models with inlet injection, fuelled with hydrogen gas, were used in the study. The scramjet engines were configured to operate in radical farming mode, where combustion radicals are formed behind shock waves reflected at the walls. The chemiluminescence emission signals were recorded in a two-dimensional, time-integrated fashion to give information on the location and distribution of the radical farms in the combustors. High signal levels were detected in localised regions immediately downstream of shock reflections, an indication of localised hydroxyl formation supporting the concept of radical farming. Results are presented for a symmetric as well as an asymmetric scramjet geometry. These data represent the first successful visualisation of radical farms in the hot pockets of a supersonic combustor. Spectrally resolved measurements have been obtained in the ultraviolet wavelength range between 300 and 400 nm. This data shows that the OH! chemiluminescence signal around 306nm is not the most dominant source of radiation observed in the radical farms.

Two-scale computational modelling of unsaturated flow in soils exhibiting small scale heterogeneities

Unsaturated water flow in soil is commonly modelled using Richards’ equation, which requires the hydraulic properties of the soil (e.g., porosity, hydraulic conductivity, etc.) to be characterised. Naturally occurring soils, however, are heterogeneous in nature, that is, they are composed of a number of interwoven homogeneous soils each with their own set of hydraulic properties. When the length scale of these soil heterogeneities is small, numerical solution of Richards’ equation is computationally impractical due to the immense effort and refinement required to mesh the actual heterogeneous geometry. A classic way forward is to use a macroscopic model, where the heterogeneous medium is replaced with a fictitious homogeneous medium, which attempts to give the average flow behaviour at the macroscopic scale (i.e., at a scale much larger than the scale of the heterogeneities). Using the homogenisation theory, a macroscopic equation can be derived that takes the form of Richards’ equation with effective parameters. A disadvantage of the macroscopic approach, however, is that it fails in cases when the assumption of local equilibrium does not hold. This limitation has seen the introduction of two-scale models that include at each point in the macroscopic domain an additional flow equation at the scale of the heterogeneities (microscopic scale). This report outlines a well-known two-scale model and contributes to the literature a number of important advances in its numerical implementation. These include the use of an unstructured control volume finite element method and image-based meshing techniques, that allow for irregular micro-scale geometries to be treated, and the use of an exponential time integration scheme that permits both scales to be resolved simultaneously in a completely coupled manner. Numerical comparisons against a classical macroscopic model confirm that only the two-scale model correctly captures the important features of the flow for a range of parameter values.

Are in vitro estimates of cell diffusivity and cell proliferation rate sensitive to assay geometry?

Cells respond to various biochemical and physical cues during wound–healing and tumour progression. In vitro assays used to study these processes are typically conducted in one particular geometry and it is unclear how the assay geometry affects the capacity of cell populations to spread, or whether the relevant mechanisms, such as cell motility and cell proliferation, are somehow sensitive to the geometry of the assay. In this work we use a circular barrier assay to characterise the spreading of cell populations in two different geometries. Assay 1 describes a tumour–like geometry where a cell population spreads outwards into an open space. Assay 2 describes a wound–like geometry where a cell population spreads inwards to close a void. We use a combination of discrete and continuum mathematical models and automated image processing methods to obtain independent estimates of the effective cell diffusivity, D, and the effective cell proliferation rate, λ. Using our parameterised mathematical model we confirm that our estimates of D and λ accurately predict the time–evolution of the location of the leading edge and the cell density profiles for both assay 1 and assay 2. Our work suggests that the effective cell diffusivity is up to 50% lower for assay 2 compared to assay 1, whereas the effective cell proliferation rate is up to 30% lower for assay 2 compared to assay 1.

Numerical simulation of idealised three-dimensional downburst wind fields

Convective downburst wind storms generate the peak annual gust wind speed for many parts of the non-cyclonic world at return periods of importance for ultimate limit state design. Despite this there is little clear understanding of how to appropriately design for these wind events given their significant dissimilarities to boundary layer winds upon which most design is based. To enhance the understanding of wind fields associated with these storms a three-dimensional numerical model was developed to simulate a multitude of idealised downburst scenarios and to investigate their near-ground wind characteristics. Stationary and translating downdraft wind events in still and sheared environments were simulated with baseline results showing good agreement with previous numerical work and full-scale observational data. Significant differences are shown in the normalised peak wind speed velocity profiles depending on the environmental wind conditions in the vicinity of the simulated event. When integrated over the height of mid- to high rise structures, all simulated profiles are shown to produce wind loads smaller than an equivalent 10 m height matched open terrain boundary layer profile. This suggests that for these structures the current design approach is conservative from an ultimate loading standpoint. Investigating the influence of topography on the structure of the simulated near-ground downburst wind fields, it is shown that these features amplify wind speeds in a manner similar to that expected for boundary layer winds, but the extent of amplification is reduced. The level of reduction is shown to be dependent on the depth of the simulated downburst outflow.

Joint identification of infinite-frequency added mass and fluid-memory models of marine structures

This paper addresses the problem of joint identification of infinite-frequency added mass and fluid memory models of marine structures from finite frequency data. This problem is relevant for cases where the code used to compute the hydrodynamic coefficients of the marine structure does not give the infinite-frequency added mass. This case is typical of codes based on 2D-potential theory since most 3D-potential-theory codes solve the boundary value associated with the infinite frequency. The method proposed in this paper presents a simpler alternative approach to other methods previously presented in the literature. The advantage of the proposed method is that the same identification procedure can be used to identify the fluid-memory models with or without having access to the infinite-frequency added mass coefficient. Therefore, it provides an extension that puts the two identification problems into the same framework. The method also exploits the constraints related to relative degree and low-frequency asymptotic values of the hydrodynamic coefficients derived from the physics of the problem, which are used as prior information to refine the obtained models.

Open access to infinite content (or "in praise of law reviews")

This Article is about legal scholarly publication in a time of plenitude. It is an attempt to explain why the most pressing questions in legal scholarly publishing are about how we ensure access to an infinity of content. It explains why standard assumptions about resource scarcity in publication are wrong in general, and how the changes in the modality of publication affect legal scholarship. It talks about the economics of open access to legal material, and how this connects to a future where there is infinite content. And because student-edited law reviews fit this future better than their commercially-produced, peer-refereed cousins, this Article is, in part, a defense of the crazy-beautiful institution that is the American law review.

Apple skin patterning is associated with differential expression of MYb10

Background Some apple (Malus × domestica Borkh.) varieties have attractive striping patterns, a quality attribute that is important for determining apple fruit market acceptance. Most apple cultivars (e.g. 'Royal Gala') produce fruit with a defined fruit pigment pattern, but in the case of 'Honeycrisp' apple, trees can produce fruits of two different kinds: striped and blushed. The causes of this phenomenon are unknown. Results Here we show that striped areas of 'Honeycrisp' and 'Royal Gala' are due to sectorial increases in anthocyanin concentration. Transcript levels of the major biosynthetic genes and MYB10, a transcription factor that upregulates apple anthocyanin production, correlated with increased anthocyanin concentration in stripes. However, nucleotide changes in the promoter and coding sequence of MYB10 do not correlate with skin pattern in 'Honeycrisp' and other cultivars differing in peel pigmentation patterns. A survey of methylation levels throughout the coding region of MYB10 and a 2.5 Kb region 5' of the ATG translation start site indicated that an area 900 bp long, starting 1400 bp upstream of the translation start site, is highly methylated. Cytosine methylation was present in all three contexts, with higher methylation levels observed for CHH and CHG (where H is A, C or T) than for CG. Comparisons of methylation levels of the MYB10 promoter in 'Honeycrisp' red and green stripes indicated that they correlate with peel phenotypes, with an enrichment of methylation observed in green stripes. Conclusions Differences in anthocyanin levels between red and green stripes can be explained by differential transcript accumulation of MYB10. Different levels of MYB10 transcript in red versus green stripes are inversely associated with methylation levels in the promoter region. Although observed methylation differences are modest, trends are consistent across years and differences are statistically significant. Methylation may be associated with the presence of a TRIM retrotransposon within the promoter region, but the presence of the TRIM element alone cannot explain the phenotypic variability observed in 'Honeycrisp'. We suggest that methylation in the MYB10 promoter is more variable in 'Honeycrisp' than in 'Royal Gala', leading to more variable color patterns in the peel of this cultivar.

Three-dimensional hydrogen-bonded structures in the hydrated proton-transfer salts of isonipecotamide with the dicarboxylic oxalic and adipic acid analogues

The structures of the 1:1 hydrated proton-transfer compounds of isonipecotamide (piperidine-4-carboxamide) with oxalic acid, 4-carbamoylpiperidinium hydrogen oxalate dihydrate, C6H13N2O+·C2HO4-·2H2O, (I), and with adipic acid, bis­(4-car­bam­oylpiperidinium) adipate dihydrate, 2C6H13N2O+·C6H8O42-·2H2O, (II), are three-dimensional hydrogen-bonded constructs involving several different types of enlarged water-bridged cyclic associations. In the structure of (I), the oxalate monoanions give head-to-tail carb­oxy­lic acid O-HOcarboxyl hydrogen-bonding inter­actions, forming C(5) chain substructures which extend along a. The isonipecotamide cations also give parallel chain substructures through amide N-HO hydrogen bonds, the chains being linked across b and down c by alternating water bridges involving both carboxyl and amide O-atom acceptors and amide and piperidinium N-HOcarboxyl hydrogen bonds, generating cyclic R43(10) and R32(11) motifs. In the structure of (II), the asymmetric unit comprises a piperidinium cation, half an adipate dianion, which lies across a crystallographic inversion centre, and a solvent water mol­ecule. In the crystal structure, the two inversion-related cations are inter­linked through the two water mol­ecules, which act as acceptors in dual amide N-HOwater hydrogen bonds, to give a cyclic R42(8) association which is conjoined with an R44(12) motif. Further N-HOwater, water O-HOamide and piperidinium N-HOcarbox­yl hydrogen bonds give the overall three-dimensional structure. The structures reported here further demonstrate the utility of the isonipecotamide cation as a synthon for the generation of stable hydrogen-bonded structures. The presence of solvent water mol­ecules in these structures is largely responsible for the non-occurrence of the common hydrogen-bonded amide-amide dimer, promoting instead various expanded cyclic hydrogen-bonding motifs.