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

Dependence of convective boundary mixing on boundary properties and turbulence strength

Convective boundary mixing is one of the major uncertainties in stellar evolution. In order to study its dependence on boundary properties and turbulence strength in a controlled way, we computed a series of 3D hydrodynamical simulations of stellar convection during carbon burning with a varying boosting factor of the driving luminosity. Our 3D implicit large eddy simulations were computed with the prompi code. We performed a mean field analysis of the simulations within the Reynolds-averaged Navier–Stokes framework. Both the vertical rms velocity within the convective region and the bulk Richardson number of the boundaries are found to scale with the driving luminosity as expected from theory: $v$ ∝ L1/3 and RiB ∝ L−2/3, respectively. The positions of the convective boundaries were estimated through the composition profiles across them, and the strength of convective boundary mixing was determined by analysing the boundaries within the framework of the entrainment law. We find that the entrainment is approximately inversely proportional to the bulk Richardson number, RiB (\propto \textrmRi_\textrmB^-α , α \sim 0.75). Although the entrainment law does not encompass all the processes occurring at boundaries, our results support the use of the entrainment law to describe convective boundary mixing in 1D models, at least for the advanced phases. The next steps and challenges ahead are also discussed

Developing graduate attributes through participation in undergraduate research conferences

© 2016 Taylor & Francis. Abstract: Graduate attributes are a framework of skills, attitudes, values and knowledge that graduates should develop by the end of their degree programmes. Adopting a largely qualitative approach and using semi-structured interviews, this paper outlines students’ experiences at a national undergraduate research conference over three years and evidences the graduate attributes developed. The students demonstrated intellectual autonomy, repurposing their work for presentation to a multidisciplinary audience through conversation with and benchmarking against peers. They gained confidence in expressing their identity as researchers and moved towards self-authorship, consciously balancing the contextual nature of their disciplinary knowledge with intra-personally grounded goals and values

The SOLAS air-sea gas exchange experiment (SAGE) 2004

Author Posting. © The Author(s), 2010. This is the author's version of the work. It is posted here by permission of Elsevier B.V. for personal use, not for redistribution. The definitive version was published in Deep Sea Research Part II: Topical Studies in Oceanography 58 (2011): 753-763, doi:10.1016/j.dsr2.2010.10.015.The SOLAS air-sea gas exchange experiment (SAGE) was a multiple-objective study investigating gas-transfer processes and the influence of iron fertilisation on biologically driven gas exchange in high-nitrate low-silicic acid low-chlorophyll (HNLSiLC) Sub-Antarctic waters characteristic of the expansive Subpolar Zone of the southern oceans. This paper provides a general introduction and summary of the main experimental findings. The release site was selected from a pre-voyage desktop study of environmental parameters to be in the south-west Bounty Trough (46.5°S 172.5°E) to the south-east of New Zealand and the experiment conducted between mid-March and mid-April 2004. In common with other mesoscale iron addition experiments (FeAX’s), SAGE was designed as a Lagrangian study quantifying key biological and physical drivers influencing the air-sea gas exchange processes of CO2, DMS and other biogenic gases associated with an iron-induced phytoplankton bloom. A dual tracer SF6/3He release enabled quantification of both the lateral evolution of a labelled volume (patch) of ocean and the air-sea tracer exchange at the 10’s of km’s scale, in conjunction with the iron fertilisation. Estimates from the dual-tracer experiment found a quadratic dependency of the gas exchange coefficient on windspeed that is widely applicable and describes air-sea gas exchange in strong wind regimes. Within the patch, local and micrometeorological gas exchange process studies (100 m scale) and physical variables such as near-surface turbulence, temperature microstructure at the interface, wave properties, and wind speed were quantified to further assist the development of gas exchange models for high-wind environments. There was a significant increase in the photosynthetic competence (Fv/Fm) of resident phytoplankton within the first day following iron addition, but in contrast to other FeAX’s, rates of net primary production and column-integrated chlorophyll a concentrations had only doubled relative to the unfertilised surrounding waters by the end of the experiment. After 15 days and four iron additions totalling 1.1 tonne Fe2+, this was a very modest response compared to the other mesoscale iron enrichment experiments. An investigation of the factors limiting bloom development considered co- limitation by light and other nutrients, the phytoplankton seed-stock and grazing regulation. Whilst incident light levels and the initial Si:N ratio were the lowest recorded in all FeAX’s to date, there was only a small seed-stock of diatoms (less than 1% of biomass) and the main response to iron addition was by the picophytoplankton. A high rate of dilution of the fertilised patch relative to phytoplankton growth rate, the greater than expected depth of the surface mixed layer and microzooplankton grazing were all considered as factors that prevented significant biomass accumulation. In line with the limited response, the enhanced biological draw-down of pCO2 was small and masked by a general increase in pCO2 due to mixing with higher pCO2 waters. The DMS precursor DMSP was kept in check through grazing activity and in contrast to most FeAX’s dissolved dimethylsulfide (DMS) concentration declined through the experiment. SAGE is an important low-end member in the range of responses to iron addition in FeAX’s. In the context of iron fertilisation as a geoengineering tool for atmospheric CO2 removal, SAGE has clearly demonstrated that a significant proportion of the low iron ocean may not produce a phytoplankton bloom in response to iron addition.SAGE was jointly funded through the New Zealand Foundation for Research, Science and Technology (FRST) programs (C01X0204) "Drivers and Mitigation of Global Change" and (C01X0223) "Ocean Ecosystems: Their Contribution to NZ Marine Productivity." Funding was also provided for specific collaborations by the US National Science Foundation from grants OCE-0326814 (Ward), OCE-0327779 (Ho), and OCE 0327188 OCE-0326814 (Minnett) and the UK Natural Environment Research Council NER/B/S/2003/00282 (Archer). The New Zealand International Science and Technology (ISAT) linkages fund provided additional funding (Archer and Ziolkowski), and the many collaborator institutions also provided valuable support

Environmental impacts of tidal power schemes

This paper describes the potential environmental changes caused by tidal power installations with illustration for schemes in the eastern Irish Sea, focusing mainly on major estuarine barrages. The generic impacts in the near-field and far-field are discussed. Results from a zero-dimensional and a two-dimensional model are presented: the former allows rapid calculations to be made for a large range of options while the latter allows the full effect on two-dimensional hydrodynamics to be investigated. It is shown that there may be a significant change in tidal amplitude at the coast of Northern Ireland. The bed stress in the Bristol Channel will be significantly reduced if a Severn barrage is constructed. Some effects on the tidal mixing are expected although the location of tidal fronts in the Irish and Celtic Seas will not be changed significantly. The largest environmental impact is expected to be on the amount of inter-tidal area retained after construction of an estuarine barrage. It is shown that the loss of mudflat scan be substantially reduced by using a dual-mode (ebb and flood generation) scheme with an increased number of turbines over the lowest-cost option

A robust and computationally efficient model of a two-dimensional coastal polynya

A two-dimensional polynya model incorporating a new computationally efficient parameterisation for the collection thickness of frazil ice at the polynya edge, H, is presented. The new parameterisation is designed to overcome the lack of robustness associated with the parameterisation of [Biggs, N.R.T., Morales Maqueda, M.A., Willmott, A.J. 2000. Polynya flux model solutions incorporating a parameterisation for the collection thickness of consolidated new ice. J. Fluid. Mech. 408, 179-204], in which H depends on the relative orientation of the polynya edge with respect to the flow of ice. As a result of this dependence, steady state polynya edge solutions exist that support seemingly unrealistic corners and that cannot be attained by evolving the polynya equations from an initial state in which the polynya is closed. The new parameterisation, removes the dependence of H on the orientation of the polynya edge by averaging the collection thickness of Biggs et al. (2000) along the polynya edge over length scales much smaller than the length of the polynya edge but much larger than the diameter of individual ice floes. The new parameterisation is robust in that steady state two-dimensional polynya solutions can now always be attained starting from a closed polynya stat

The opening of wind-driven polynyas

The opening of wind-driven coastal polynyas has often been investigated using idealised flux models. Polynya flux models postulate that the boundary separating the region of thin ice adjacent to the coast within the polynya from the thicker ice piling up downstream is a mathematical shock. To conserve mass, any divergence of the ice flux across the shock translates into a change in the shock’s position or, in other words, a change in the width of the thin-ice region of the polynya. Polynya flux models are physically incomplete in that, while they conserve ice mass, they do not conserve linear momentum. In this paper, we investigate the improvements that can be achieved in the simulation of polynyas by imposing conservation of momentum as well as mass. We start by adopting a mathematically solid formulation of the ice mass and momentum balances throughout the polynya region, from the coast to the pack ice. Hydrostatic and plastic versions of the ice internal forces are used in the model. Two different approaches are then explored. We first postulate the existence of a shock at the seaward edge of the thin-ice region of the polynya and derive jump conditions for the conservation of ice mass and momentum at the shock which are consistent with the continuous model physics. Polynyas simulated by this mass- and momentum-conserving shock model always reach a steady state if the polynya forcing is uniform in space and constant in time. This is also true for all polynya flux models presented previously in the literature, but the location of the steady-state polynya edge and the time required to reach it can greatly differ between shock formulations and more simplistic flux ones. We next relax the assumption that a shock exists and let the boundary between thin ice and piling up ice emerge naturally as part of the full solution of the continuous model equations. Polynyas simulated in this way are very different from those simulated by either shock or flux models. Most notably, we find that steady-state polynya solutions are not always attainable in the continuous model. We determine under which conditions this is so and explain how such unsteady solutions come about. We also show that, in those cases when a steady-state solution exists in the continuous model, the steady-state polynya width is considerably larger than in the shock model, and the time required to attain it is accordingly longer. The occurrence of such significant differences between the polynya solutions calculated with flux and shock models, on the one hand, and with more sophisticated continuous formulations, on the other hand, suggests that the former are, at best, incomplete, and should be used with caution

Appraising the extractable tidal energy resource of the UK's western coastal waters

A two-dimensional west coast tidal model, built on the ADCIRC platform (an unstructured grid two-dimensional depth-integrated shallow water model), has been developed to examine the scope for reliable and fully predictable electricity generation from UK coastal waters using an ambitious combination of estuary barrages, tidal lagoons and tidal stream generator arrays. The main emphasis has been towards conjunctive operation of major estuary barrages, initially including the presence of pilot-scale tidal stream developments, though ambitious exploitation of extensive tidal streams has also been explored