Numerical derivation of constitutive models for unbonded flexible risers

This is the post-print version of the final paper published in International Journal of Mechanical Sciences. The published article is available from the link below. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication.In this paper a new constitutive model for flexible risers is proposed and a procedure for the identification of the related input parameters is developed using a multi-scale approach. The constitutive model is formulated in the framework of an Euler–Bernoulli beam model, with the addition of suitable pressure terms to the generalized stresses to account for the internal and external pressures, and therefore can be efficiently used for large-scale analyses. The developed non-linear relationship between generalized stresses and strains in the beam is based on the analogy between frictional slipping between different layers of a flexible riser and frictional slipping between micro-planes of a continuum medium in non-associative elasto-plasticity. Hence, a linear elastic relationship is used for the initial response in which no-slip occurs; an onset-slip function is introduced to define the ‘no-slip’ domain, i.e. the set of generalized stresses for which no slip occurs; a non-associative rule with linear kinematic hardening is used to model the full-slip phase. The results of several numerical simulations for a riser of small-length, obtained with a very detailed (small-scale) non-linear finite-element model, are used to identify the parameters of the constitutive law, bridging in this way the small scale of the detailed finite-element simulations with the large scale of the beam model. The effectiveness of the proposed method is validated by the satisfactory agreement between the results of various detailed finite-element simulations for a short riser, subject to internal and external uniform pressure and uniform cyclic bending loading, with those given by the proposed constitutive law.Lloyds Register EME

Nonlinear multi-scale homogenization with different structural models at different scales

We present an extension of the computational homogenization theory to cases where different structural models are used at different scales and no energy potential can be defined at the small scale. We observe that volumetric averaging, which is not applicable in such cases unless similarities exist in the macro-scale and micro-scale models, is not a necessary prerequisite to carry out computational homogenization. At each material point of the macro-model, we replace the conventional representative volume element with a representative domain element (RDE). To link the large-scale and small-scale problems, we then introduce a linear operator, mapping the smooth part of the small-scale displacement field of each RDE to the large-scale strain field and a trace operator to impose boundary conditions in the RDE. The latter is defined on the basis of engineering judgement, analogously to the conventional theory. A generalized Hill’s condition, rather than being invoked, is derived from duality principles and is used to recover the stress measures at the large scale. For the implementation in a nonlinear finite-element analysis, ‘control nodes’ and constraint equationsare used. The effectiveness of the procedure is demonstrated for three beam-to-truss example problems, for which multi-scale convergence is numerically analysed.Lloyd’s Register EME

Influence of natural convection on gold nanorods-assisted photothermal treatment of bladder cancer in mice

Background: The thermally-induced urine flow can generate cooling that may alter the treatment outcome during hyperthermic treatments of bladder cancer. This paper investigates the effects of natural convection inside the bladder and at skin surface during gold nanorods (GNR) - assisted photothermal therapy (PTT) of bladder cancer in mice. Methods: 3D models of mouse bladder at orientations corresponding to the mouse positioned on its back, its side and its abdomen were examined. Numerical simulations were carried out for GNR volume fractions of 0.001, 0.005 and 0.01% and laser power of 0.2 and 0.3 W. Results: The obtained results showed that cooling due to natural convection inside the bladder and above the skin depends on the mouse orientation. For a mouse positioned on its back, on its side or on its abdomen, the maximum temperature achieved inside the tumour at 0.001% GNR volume fraction and 0.2 W laser power was 55.2°C, 50.0°C and 52.2°C, respectively compared to 56.8°C when natural convection was not considered. The average thermal gradients when natural convection was considered were also lower, suggesting a more homogenous temperature distribution. Conclusions: Natural convection inside the bladder can be beneficial but also detrimental to GNR-assisted PTT depending on the level of heating. At low levels of heating due to low GNR volume fraction and/or laser power, flow inside the bladder may dissipate heat from the targeted tissue; making the treatment ineffective. At high levels of heating due to high GNR volume fraction and/or laser power, cooling may prevent excessive thermal damage to surrounding tissues

Time-gated transillumination and reflection by biological tissues and tissuelike phantoms: simulation versus experiment

A numerical method is presented to solve exactly the time-dependent diffusion equation that describes light transport in turbid media. The simulation takes into account spatial variations of the scattering and absorption factors of the medium and the objects as well as random fluctuations of these quantities. The technique is employed to explore the possibility of locating millimeter-sized objects immersed in turbid media from time-gated measurements of the transmitted or reflected (near-infrared) light. The simulation results for tissue-like phantoms are compared with experimental transillumination data, and excellent agreement is found. Simulations of time-gated reflection experiments indicate that it may be possible to detect objects of 1-mm radius.

Extended Powell–Sabin finite element scheme for linear elastic fracture mechanics

Data availability: No data was used for the research described in the article.Powell–Sabin B-splines, which are based on triangles, are employed in the framework of the extended finite element method (XFEM) for fracture analysis. This avoids the necessity of remeshing in discrete fracture models and increases the solution accuracy around the crack tip. Powell–Sabin B-splines are -continuous throughout the whole domain. The stresses around crack tips are captured more accurately than when using elements with a standard Lagrangian interpolation. Although Powell–Sabin B-splines do not hold the Kronecker-delta property, the Heaviside function and the tip enrichment function are confined to the cracked elements only, similar to the traditional XFEM but different from the extended isogeometric method. In addition, Powell–Sabin B-splines still hold -continuous throughout cracked elements. There is no need to lower the continuity at element boundaries, to confine basis function support in cracked elements. Shifting is used to ensure compatibility with the surrounding discretization. The sub-triangle technique is employed for the numerical integration over crack elements. The versatility and accuracy of the approach to simulate crack problems are assessed in case studies, featuring mode-I and mixed-mode crack problems.Engineering and Physical Sciences Research Council (EPSRC: EP/T033940/1)

Mechanical Response of a Lined Pipe Under Dynamic Impact

An experimental and numerical investigation on the mechanical response of a lined pipe (compound pipe) under a dynamic impact is presented. The influence of the impact energy has been studied in terms of the depth of the dent formed, and of the strains and residual stresses. To this end, a three-dimensional explicit dynamic non-linear finite element model has been developed and successfully validated against the results of impact-test experiments conducted on pipes made of AISI 10305 steel, with and without the AISI304 stainless steel liner. The validation was made by comparing numerically computed strains with those measured by strain gauges, as well as in terms of permanent deformation. The model is then utilized to evaluate the residual stresses, the amount of energy dissipation and the velocity of impact process as a function of different pipes (i.e. with or without liner) and of the free drop heights.Design for Resource and Energy efficiency in cerAMic kiln

Midinfrared intersubband absorption in ZnxCd1−xSe/Znx′Cdy′Mg1−x′−y′Se multiple quantum well structures

The authors report the observation of intersubband absorption in ZnxCd(1−x)Se/Znx′Cdy′Mg(1−x′−y′)Semultiple quantum wells. Lattice-matched samples were grown by molecular beam epitaxy on InP (001) substrates. Photoluminescence measurements indicate that the samples have excellent material quality. The peak absorption wavelengths measured by Fourier transform infrared spectroscopy are 3.99 and 5.35μm for two samples with ZnxCd(1−x)Se well widths of 28 and 42Å, respectively. These values fall within the 3–5μm wavelength range, which is of interest for midinfrared intersubband devices, such as quantum cascade lasers and quantum well infrared photodetectors. Their experimental results fit well with theoretical predictions based on the envelope function approximation. The results indicate that these wide band gap II-VI materials are very promising for midinfrared intersubband device applications

Meat freshness revealed by visible to near-infrared spectroscopy and principal component analysis

Increasing concerns about adulterated meat encouraged industry looking for new non-invasive methods for rapid accurate meat quality assessment. Main meat chromophores (myoglobin, oxy-myoglobin, fat, water, collagen) are characterized by close comparable absorption in visible to near-infrared (NIR) spectral region. Therefore, structural and compositional variations in meat may lead to relative differences in the absorption of light. Utilizing typical fiber-optic probes and integrating sphere, a degradation of pork samples freshness was observed at room temperature referring to the relative changes in absorbance of main meat chromophores. The application of principal component analysis (PCA) used for examination of measured absorbance spectra revealed more detailed sub-stages of freshness, which are not observed by the conventional analysis of the reflectance spectra. The results show a great potential of the combined application of optical-NIR spectroscopy with complementary use of PCA approach for assessing meat quality and monitoring relative absorbance alternation of oxymyoglobin and myoglobin in visible, and fat, water, collagen in NIR spectral ranges

The Static Failure of Adhesively Bonded Metal Laminate Structures: A Cohesive Zone Approach

Data on distribution, ecology, biomass, recruitment, growth, mortality and productivity of the West African bloody cockle Anadara senilis were collected at the Banc d'Aguuin, Mauritania, in early 1985 and 1986. Ash-free dry weight appeared to be correlated best with shell height. A. senilis was abundant on the tidal flats of landlocked coastal bays, but nearly absent on the tidal flats bordering the open sea. The average biomass for the entire area of tidal flats was estimated at 5.5 g·m−2 ash-free dry weight. The A. senilis population appeared to consist mainly of 10 to 20-year-old individuals, showing a very slow growth and a production: biomass ratio of about 0.02 y−1. Recruitment appeared negligible and mortality was estimated to be about 10% per year. Oystercatchers (Haematopus ostralegus), the gastropod Cymbium cymbium and unknown fish species were responsible for a large share of this. The distinction of annual growth marks permitted the assessment of year-class strength, which appeared to be correlated with the average discharge of the river Senegal. This may be explained by assuming that year-class strength and river discharge both are correlated with rainfall at the Banc d'Arguin.