An axisymmetric time-domain spectral-element method for full-wave simulations: Application to ocean acoustics

The numerical simulation of acoustic waves in complex 3D media is a key topic in many branches of science, from exploration geophysics to non-destructive testing and medical imaging. With the drastic increase in computing capabilities this field has dramatically grown in the last twenty years. However many 3D computations, especially at high frequency and/or long range, are still far beyond current reach and force researchers to resort to approximations, for example by working in 2D (plane strain) or by using a paraxial approximation. This article presents and validates a numerical technique based on an axisymmetric formulation of a spectral finite-element method in the time domain for heterogeneous fluid-solid media. Taking advantage of axisymmetry enables the study of relevant 3D configurations at a very moderate computational cost. The axisymmetric spectral-element formulation is first introduced, and validation tests are then performed. A typical application of interest in ocean acoustics showing upslope propagation above a dipping viscoelastic ocean bottom is then presented. The method correctly models backscattered waves and explains the transmission losses discrepancies pointed out in Jensen et al. (2007). Finally, a realistic application to a double seamount problem is considered.Comment: Added a reference, and fixed a typo (cylindrical versus spherical

On a Cahn--Hilliard--Darcy system for tumour growth with solution dependent source terms

We study the existence of weak solutions to a mixture model for tumour growth that consists of a Cahn--Hilliard--Darcy system coupled with an elliptic reaction-diffusion equation. The Darcy law gives rise to an elliptic equation for the pressure that is coupled to the convective Cahn--Hilliard equation through convective and source terms. Both Dirichlet and Robin boundary conditions are considered for the pressure variable, which allows for the source terms to be dependent on the solution variables.Comment: 18 pages, changed proof from fixed point argument to Galerkin approximatio

Relative Importance of Convective Uncertainties in Massive Stars

In this work, we investigate the impact of uncertainties due to convective boundary mixing (CBM), commonly called ‘overshoot’, namely the boundary location and the amount of mixing at the convective boundary, on stellar structure and evolution. For this we calculated two grids of stellar evolution models with the MESA code, each with the Ledoux and the Schwarzschild boundary criterion, and vary the amount of CBM. We calculate each grid with the initial masses 15, 20 and $25\, \rm M_ødot$. We present the stellar structure of the models during the hydrogen and helium burning phases. In the latter, we examine the impact on the nucleosynthesis. We find a broadening of the main-sequence with more CBM, which is more in agreement with observations. Furthermore during the core hydrogen burning phase there is a convergence of the convective boundary location due to CBM. The uncertainties of the intermediate convective zone remove this convergence. The behaviour of this convective zone strongly affects the surface evolution of the model, i.e. how fast it evolves red-wards. The amount of CBM impacts the size of the convective cores and the nucleosynthesis, e.g. the 12C to 16O ratio and the weak s-process. Lastly, we determine the uncertainty that the range of parameter values investigated introduce and we find differences of up to $70$ for the core masses and the total mass of the star

3D hydrodynamic simulations of carbon burning in massive stars

We present the first detailed 3D hydrodynamic implicit large eddy simulations of turbulent convection of carbon burning in massive stars. Simulations begin with radial profiles mapped from a carbon-burning shell within a 15 M⊙ 1D stellar evolution model. We consider models with 1283, 2563, 5123, and 10243 zones. The turbulent flow properties of these carbon-burning simulations are very similar to the oxygen-burning case. We performed a mean field analysis of the kinetic energy budgets within the Reynolds-averaged Navier–Stokes framework. For the upper convective boundary region, we find that the numerical dissipation is insensitive to resolution for linear mesh resolutions above 512 grid points. For the stiffer, more stratified lower boundary, our highest resolution model still shows signs of decreasing sub-grid dissipation suggesting it is not yet numerically converged. We find that the widths of the upper and lower boundaries are roughly 30 per cent and 10 per cent of the local pressure scaleheights, respectively. The shape of the boundaries is significantly different from those used in stellar evolution models. As in past oxygen-shell-burning simulations, we observe entrainment at both boundaries in our carbon-shell-burning simulations. In the large Péclet number regime found in the advanced phases, the entrainment rate is roughly inversely proportional to the bulk Richardson number, RiB (∝RiB−α, 0.5 ≲ α ≲ 1.0). We thus suggest the use of RiB as a means to take into account the results of 3D hydrodynamics simulations in new 1D prescriptions of convective boundary mixing

On a diffuse interface model for tumour growth with non-local interactions and degenerate mobilities

We study a non-local variant of a diffuse interface model proposed by Hawkins--Darrud et al. (2012) for tumour growth in the presence of a chemical species acting as nutrient. The system consists of a Cahn--Hilliard equation coupled to a reaction-diffusion equation. For non-degenerate mobilities and smooth potentials, we derive well-posedness results, which are the non-local analogue of those obtained in Frigeri et al. (European J. Appl. Math. 2015). Furthermore, we establish existence of weak solutions for the case of degenerate mobilities and singular potentials, which serves to confine the order parameter to its physically relevant interval. Due to the non-local nature of the equations, under additional assumptions continuous dependence on initial data can also be shown.Comment: 28 page

The effect of interstitial pressure on therapeutic agent transport : coupling with the tumor blood and lymphatic vascular systems

Vascularized tumor growth is characterized by both abnormal interstitial fluid flow and the associated interstitial fluid pressure (IFP). Here, we study the effect that these conditions have on the transport of therapeutic agents during chemotherapy. We apply our recently developed vascular tumor growth model which couples a continuous growth component with a discrete angiogenesis model to show that hypertensive IFP is a physical barrier that may hinder vascular extravasation of agents through transvascular fluid flux convection, which drives the agents away from the tumor. This result is consistent with previous work using simpler models without blood flow or lymphatic drainage. We consider the vascular/interstitial/lymphatic fluid dynamics to show that tumors with larger lymphatic resistance increase the agent concentration more rapidly while also experiencing faster washout. In contrast, tumors with smaller lymphatic resistance accumulate less agents but are able to retain them for a longer time. The agent availability (area-under-the curve, or AUC) increases for less permeable agents as lymphatic resistance increases, and correspondingly decreases for more permeable agents. We also investigate the effect of vascular pathologies on agent transport. We show that elevated vascular hydraulic conductivity contributes to the highest AUC when the agent is less permeable, but to lower AUC when the agent is more permeable. We find that elevated interstitial hydraulic conductivity contributes to low AUC in general regardless of the transvascular agent transport capability. We also couple the agent transport with the tumor dynamics to simulate chemotherapy with the same vascularized tumor under different vascular pathologies. We show that tumors with an elevated interstitial hydraulic conductivity alone require the strongest dosage to shrink. We further show that tumors with elevated vascular hydraulic conductivity are more hypoxic during therapy and that the response slows down as the tumor shrinks due to the heterogeneity and low concentration of agents in the tumor interior compared with the cases where other pathological effects may combine to flatten the IFP and thus reduce the heterogeneity. We conclude that dual normalizations of the micronevironment ? both the vasculature and the interstitium ? are needed to maximize the effects of chemotherapy, while normalization of only one of these may be insufficient to overcome the physical resistance and may thus lead to sub-optimal outcomes.PostprintPeer reviewe

Miscarriage following dengue virus 3 infection in the first six weeks of pregnancy of a dengue virus-naive traveller returning from Bali to Italy, April 2016

We report miscarriage following dengue virus (DENV)-3 infection in a pregnant woman returning from Bali to Italy in April 2016. On her arrival, the woman had fever, rash, asthenia and headache. DENV RNA was detected in plasma and urine samples collected the following day. Six days after symptom onset, she had a miscarriage. DENV RNA was detected in fetal material, but in utero fetal infection cannot be demonstrated due to possible contamination by maternal blood