Strong Density Gradients Research Articles

An improved stochastic weighted particle method for boundary driven flows

The stochastic weighted particle method (SWPM) is a generalization of the Direct Simulation Monte Carlo (DSMC) method where particle weights are variable and dynamic. SWPM is backed by a strong theoretical foundation but has not been critically evaluated for problems of practical interest. A thorough assessment of SWPM for boundary-driven flows reveals significant numerical artifacts near the boundary, notably a diverging heat flux. To correct the boundary heat flux, two modifications to SWPM are proposed: separated grouping and a spatially-dependent weight transfer function. To gauge the relative efficiency of SWPM in comparison to DSMC, a high-Mach-number wheel flow which forms a strong density gradient is also simulated.

Dynamics of the Polar Front in the southwestern area of Svalbard, Norway

We have investigated the dynamics of the Polar Front (PF) in the southwestern area of Svalbard, where the dynamics are defined by the interaction between the Sørkapp Current (transporting cold and fresh Arctic-type water) and the West Spitsbergen Current (carrying the warmer and more saline Atlantic Water from the Norwegian Sea). Our approach was based on the global circulation model predictions. In accord with the model, the cold and fresher Arctic water masses originated from the eastern side of Svalbard and were flowing along the southern-most point of Svalbard into the southwestern area of Svalbard. This cold and fresher water was spreading and pushing warmer and more saline Atlantic water offshore. The model currents indicated that cold and fresher water masses were mostly transported northward along the south-western coast of Svalbard. In accord with the model and observed density profiles, the PF is a density compensated front with the density gradient mostly in the top 50-75 m. Presence of strong density gradients in the upper 50-75 m of the PF leads to the development of the submesoscale processes as surface frontogenesis and nonlinear Ekman transport. We found strong presence of surface frontogenesis at the PF in the model. The second submesoscale process, the nonlinear Ekman transport, is due to the forcing interaction. During the considered time frame, the winds were blowing mostly from the north-east to south-west across the PF, and therefore, the along-front wind component was very small, resulting in weak non-linear Ekman transport at the PF front.

Turbulence-reduced high-performance scenarios in Wendelstein 7-X

In the Wendelstein 7-X (W7-X) stellarator, turbulence is the dominant transport mechanism in most discharges. This leads to a ‘clamping’ of ion temperature over a wide range of heating power, predominantly flat density profiles where hollow profiles driven by neoclassical thermo-diffusion would be expected and by rapid impurity transport in injection experiments. Significantly reduced turbulent transport is observed in the presence of strong core density gradients found transiently after core pellet injection and irregularly after boronisation or boron pellet injection. Density peaking is also achieved in a controlled manner in purely neutral beam heated discharges where particle transport analysis reveals an abrupt reduction in the main-ion particle flux leading to significant density profile peaking not explained by the NBI particle source alone. The plasmas exhibit a heat diffusivity of around χ=0.25±0.1m2 s−1 at mid radius, a factor of around 4 lower than ECRH dominated discharges. Despite the improved confinement, the achieved ion temperature is limited by broader heat deposition and the lower power-per-particle given the higher density. This is overcome with limited reintroduction of ECRH power, where the low heat diffusivity diffusivity is maintained, the density rise supressed and ion temperatures above the clamping limit are achieved. The applicability of these plasmas for a high performance scenario on transport relevant time scales is assessed, including initial predictions for planned heating upgrades of W7-X, based on a range of assumptions about particle transport.

Characterization of heated volume generation by nanosecond pulsed plasma actuator with various pressure environments

Nanosecond dielectric barrier discharge (NS-DBD) has emerged as a promising technique for controlling high-speed flows, generating a heated volume that generates strong density and viscosity gradients, thereby perturbing flow dynamics. Since its potential application in low-pressure, high-speed flows, understanding how the size and growth of the heated volume correlate with surrounding pressure is crucial. In this study, we employed typical schlieren and background-oriented schlieren (BOS) techniques to investigate the heated volume’s sensitivity to surrounding pressure in quiescent air. The observed heated volume’s size variations with surrounding pressure likely stemmed from the increase in thermal diffusivity at lower pressures. BOS findings unveiled a nearly linear decrease in heated volume’s core density with energy input. Meanwhile, the heated volume’s size augmented with energy input but exhibited gradual saturation, attributable possibly to shear stresses impeding volume expansion as temperature and viscosity rose, or to consumption of energy in vibration excitation and other reactions. In the cases of 100 and 50 kPa, the sensitivity of the heated volume’s size to the reduced electric field appeared to be similar. However, at 10 kPa, where the reduced electric field is higher compared to that of the 100 and 50 kPa cases due to the lower air density, the size sensitivity drastically decreased. This suggested a transition in discharge mode from filamentary to diffusive behavior at lower pressures.

Laser-driven pointed acceleration of electrons with preformed plasma lens

The simultaneous laser-driven acceleration and angular manipulation of the fast electron beam are experimentally demonstrated. The bunch of multi-MeV energy charged particles is generated during the propagation of the femtosecond laser pulse through the near-critical plasma slab accompanied by plasma channeling. Plasma is formed by the controlled breakdown of a thin-tape target by a powerful nanosecond prepulse. The electron beam pointing approach is based on the refraction of a laser pulse in the presence of a strong radial density gradient in the breakdown of the tape with a small displacement of the femtosecond laser beam relative to the breakdown symmetry axis. A shift of several micrometers makes it possible to achieve beam deflection by an angle up to 10° with acceptable beam charge and spectrum conservation. This opens up opportunities for applications for scanning objects with an electron beam and the multistage electron beam energy gain in consecutive laser accelerators without bulk magnetic optics for particles. Experimental findings are supported by numerical particle-in-cell calculations of laser-plasma acceleration and hydrodynamic simulations. Published by the American Physical Society 2024

Scintillation modeling with random phase gradient screens

Multiple studies of scintillation phenomena have shown that, in certain situations, the intense phase fluctuations of trans-ionospheric radio signals are associated with the scattering on strong electron density gradients. The present study provides a theoretical framework for modeling such types of phase fluctuation events. Using the geometrical-optics approximation and retaining the second-order smallness correction in the expansion of the eikonal function, we relate the phase of the transmitted wave not only to the total electron content (TEC) of the ionosphere but also to the spatial gradient of the TEC. The considered correction term is related to the random refraction of signal rays on large-scale ionospheric structures, an effect, that becomes significant in the presence of strong electron density gradients. To conveniently simulate the wave propagation under such conditions, we propose the random phase gradient screen algorithm. For this purpose, we use the novel spatial electron density gradient product (NeGIX) based on in-situ observations of the Swarm Langmuir probe and ground-based TEC and TEC gradient observations. To illustrate the performance of the algorithm, we apply it to simulate a scintillation event over Europe and in the low-latitude region and compare the simulation results with scintillation indices, measured from GNSS ground observations. We show that in regions of the ionosphere where spatial ionospheric gradients are large, the phase gradient method shows better agreement with the observed scintillation levels than the conventional phase screen approach.

BONE REMODELLING PREDICTION FOR IMPROVED POROUS IMPLANT DESIGN

AbstractObjectivesYoung patients receiving metallic bone implants after surgical resection of bone cancer require implants that last into adulthood, and ideally life-long. Porous implants with similar stiffness to bone can promote bone ingrowth and thus beneficial clinical outcomes. A mechanical remodelling stimulus, strain energy density (SED), is thought to be the primary control variable of the process of bone growth into porous implants. The sequential process of bone growth needs to be taken into account to develop an accurate and validated bone remodelling algorithm, which can be employed to improve porous implant design and achieve better clinical outcomes.MethodsA bone remodelling algorithm was developed, incorporating the concept of bone connectivity (sequential growth of bone from existing bone) to make the algorithm more physiologically relevant. The algorithm includes adaptive elastic modulus based on apparent bone density, using a node-based model to simulate local remodelling variations while alleviating numerical checkerboard problems. Strain energy density (SED) incorporating stress and strain effects in all directions was used as the primary stimulus for bone remodelling. The simulations were developed to run in MATLAB interfacing with the commercial FEA software ABAQUS and Python. The algorithm was applied to predict bone ingrowth into a porous implant for comparison against data from a sheep model.ResultsThe accuracy of the predicted bone remodelling was verified for standard loading cases (bending, torsion) against analytical calculations. Good convergence was achieved. The algorithm predicted good bone remodelling and growth into the investigated porous implant. Using the standard algorithm without connectivity, bone started to remodel at locations unconnected to any bone, which is physiologically implausible. The implementation of bone connectivity ensures the gradual process of bone growth into the implant pores from the sides. The bone connectivity algorithm predicted that the full remodelling required more time (approximately 50% longer), which should be considered when developing post-surgical rehabilitation strategies for patients. Both algorithms with and without bone connectivity implementation converged to same final stiffness (less than 0.01% difference). Almost all nodes reached the same density value, with only a limited number of nodes (less than 1%) in transition areas with a strong density gradient having noticeable differences.ConclusionsAn improved bone remodelling algorithm based on strain energy density that modelled the sequential process of bone growth has been developed and tested. For a porous metallic bone implant the same final bone density distribution as for the original adaptive elasticity theory was predicted, with a slower and more fidelic process of growth from existing surrounding bone into the porous implant.Declaration of Interest(b) declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported:I declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research project.

Stirring across the Antarctic Circumpolar Current's southern boundary at the prime meridian, Weddell Sea

Abstract. At the southern boundary of the Antarctic Circumpolar Current (ACC), relatively warm ACC waters encounter the colder waters surrounding Antarctica. Strong density gradients across the southern boundary indicate the presence of a frontal jet and are thought to modulate the southward heat transport across the front. In this study, the southern boundary in the Weddell Sea sector at the prime meridian is surveyed for the first time in high resolution over 2 months during an austral summer with underwater gliders occupying a transect across the front on five occasions. The five transects show that the frontal structure (i.e. hydrography, velocities and lateral density gradients) varies temporally. The results demonstrate significant, transient (a few weeks) variability of the southern boundary and its frontal jet in location, strength and width. A mesoscale cold-core eddy is identified to disrupt the southern boundary’s frontal structure and strengthen lateral density gradients across the front. The front's barrier properties are assessed using mixing length scales and potential vorticity to establish the cross-frontal exchange of properties between the ACC and the Weddell Gyre. The results show that stronger lateral density gradients caused by the mesoscale eddy strengthen the barrier-like properties of the front through reduced mixing length scales and pronounced gradients of potential vorticity. In contrast, the barrier-like properties of the southern boundary are reduced when no mesoscale eddy is influencing the density gradients across the front. Using satellite altimetry, we further demonstrate that the barrier properties over the past decade have strengthened as a result of increased meridional gradients of absolute dynamic topography and increased frontal jet speeds in comparison to previous decades. Our results emphasise that locally and rapidly changing barrier properties of the southern boundary are important to quantify the cross-frontal exchange, which is particularly relevant in regions where the southern boundary is located near the Antarctic shelf break (e.g. in the West Antarctic sector).

Numerical Study of Hypersonic Near-Continuum Flow by Improved Slip Boundary Condition

Appropriate numerical models and boundary conditions are required to accurately model hypersonic flows involving local nonequilibrium effects, which can be induced by the wall effect and large gradients of flow properties in the flowfield. Considering the strong density gradient near the wall for hypersonic flows, an improved Gökçen slip model is proposed in this paper by introducing a scaled mean free path to the conventional Gökçen’s model. Hypersonic argon flows over a circular cylinder and a blunt plate in the near-continuum regime are simulated by the proposed slip model with the conventional Navier–Stokes–Fourier equations and an extended continuum model, which was recently developed to model the shear nonequilibrium effect in high-speed flows. The accuracy of the new boundary model is evaluated by comparing the results with the direct simulation Monte Carlo method. The effects of Knudsen number, Mach number, and wall temperature on the nonequilibrium degree and aerodynamic predictions are also investigated. The numerical results show that the improved model can well predict the surface friction and heat transfer coefficients, exhibiting much better performance than the existing continuum models.

Electromagnetic drift wave instability in tokamak plasmas with strong pedestal gradient

The linear eigenmode characterizations and the nonlinear turbulence energy spreading of the drift waves in a tokamak plasma with strong pedestal gradient are numerically investigated based on an electromagnetic Landau fluid model. By the linear eigenmode analysis, it is found that the dominant instability in the low β regime is the ion-temperature-gradient (ITGc) mode and the electron drift wave instability (eDWI p ) in the core and edge region with strong density gradient, respectively. Multiple eigenstates of the eDWI p with different peak locations in the poloidal direction can be obtained by the eigenvalue problem solver. The dominant one is the high order eDWI p corresponding to the unconventional ballooning mode structure with multiple peaks in the poloidal position, in contrast to the conventional modes that peak at the outboard mid-plane, and has been verified through initial value simulation. In the high β regime, the dominant eigenmodes in the core and edge region are the conventional and unconventional kinetic ballooning modes respectively. In the nonlinear simulation, an inward turbulence spreading phenomenon during the quasi-saturation phase of the edge turbulence is clearly observed. The inward speed of the turbulence energy front in the high β regime is much faster than that in the low β regime. It is interestingly found that the speed of the turbulence energy front increases with the increase of the plasma β in the low β regime, while it is almost unchanged in the high β regime. It is identified that the turbulence spreading in the low and high β regimes are determined by the nonlinear dynamics and the linear toroidal coupling respectively.

Heat and carbon uptake in the Southern Ocean: the state of the art and future priorities.

Heat and carbon uptake in the Southern Ocean: the state of the art and future priorities.

Warming waters of the Gulf of Maine: The role of Shelf, Slope and Gulf Stream Water masses

We present an analysis of a 19-year time series of mooring data (T, S,) collected since 2003 in Jordan Basin in the eastern Gulf of Maine and beginning in 2016, in situ nitrate; we also include hydrographic station data collected by NOAA’s ECOMON Program. Analyses support, in part, conclusions drawn by others that the Gulf of Maine has warmed and continues to warm as a result of long-term and shorter term air-sea heat fluxes. We show that warming is to a large extent a result of changes in advective heat fluxes via lateral exchanges with water masses in neighboring shelf and slope regions (Warm Slope Water and Scotian Shelf Water). Water properties of the source waters and their inflows to the Gulf have been shown to have been changing over time, and now reveal a mix that includes warm, salty, low nutrient Gulf Stream Water. The result is a new baseline in the Gulf of Maine, of warmer temperatures and higher salinities, which appears to have become established after 2010. Annual cycles of spring-to summer water column warming and fall-to-winter cooling and convective mixing are presented and discussed, showing that the maximum depth of convective mixing in the eastern interior Gulf of Maine does not usually exceed 150 m, a depth constrained by the presence a strong vertical density gradient created by layer of high salinity bottom waters. Approximately-six months after the annual late-summer influx to the Gulf of warm and salty water masses at depth and along the bottom, a mass of cold, low salinity Scotian Shelf Water arrives in the upper water layers in winter – these two major water mass fluxes are a half year out of phase with one another, every year. We suggest that these two water mass intrusion events are linked – that the winter-spring inflow of Scotian Shelf Water creates a barotropic pressure gradient in the Gulf that persists into the summer, thereby impeding the inflow of warm and salty Slope, and/or Gulf Stream Water, until that Scotian Shelf Water mass has become sufficiently dispersed in late summer, thereby effectively opening a valve that allows warm and salty deep waters from outside the Gulf to penetrate.

Surface viscosity in simple liquids.

The response of Newtonian liquids to small perturbations is usually considered to be fully described by homogeneous transport coefficients like shear and dilatational viscosity. However, the presence of strong density gradients at the liquid/vapor boundary of fluids hints at the possible existence of an inhomogeneous viscosity. Here, we show that a surface viscosity emerges from the collective dynamics of interfacial layers in molecular simulations of simple liquids. We estimate the surface viscosity to be 8-16 times smaller than that of the bulk fluid at the thermodynamic point considered. This result can have important implications for reactions at liquid surfaces in atmospheric chemistry and catalysis.

Statistical Properties of Lower Hybrid Waves in the Magnetopause Reconnection Exhaust Region

AbstractLower hybrid waves (LHWs) are frequently observed associated with magnetic reconnection. A general picture of the statistical properties of LHWs in magnetopause reconnection, however, has not yet been established. Here we report such a statistical study of LHWs using magnetospheric multiscale mission data in the magnetopause reconnection exhaust region. LHWs are identified in reconnection outflow regions over a broad range of magnetic local time. Most of the LHWs are near strong density gradients located on the low‐density magnetospheric side of the reconnection layer. The wave amplitude is typically several mV/m but can occasionally reach a few tens of mV/m. Compared with intervals that do not exhibit LHWs, intervals with LHWs show a ∼50% increase in parallel electron temperature, providing statistical evidence that LHWs contribute to parallel electron heating in reconnection exhausts and the dayside low‐latitude boundary layer.

Vertical Motions Driven by Geostrophic Stress in Surface Boundary Layer under a Finite Ekman Number Regime

Abstract Geostrophic stress caused by a strong horizontal density gradient embedded in the surface boundary layer plays an important role in generating vertical motion and associated tracer transport. However, dependence of this frictionally driven vertical velocity on the Ekman number (Ek), a key dimensionless parameter for frictional flows in a rotating reference frame, has not been systematically analyzed, especially for a finite Ek. In this study, we theoretically demonstrate that the geostrophic stress always induces an ageostrophic stress acting to offset itself, and such an offsetting effect becomes more evident with increasing Ek. When Ek approaches unity or larger, vertical motion driven by geostrophic stress is much weaker than that derived by Garrett and Loder (GL81), who neglect effects of ageostrophic stress and predict a vertical velocity magnitude scaled with curl of geostrophic stress. Although the cancellation tendency between geostrophic and ageostrophic stress is universal, its underlying dynamics depends on vertical structures of turbulent viscosity and geostrophic flows. A realistic simulation in the winter Kuroshio Extension is conducted to validate the theoretical results and examine which regime, a small versus finite Ek, is more relevant in this region. It is found that the characteristic vertical scale involved in the definition of Ek is primarily determined by the vertical structure of turbulent viscosity and evidently smaller than that of geostrophic flow. The value of Ek in the winter Kuroshio Extension is generally larger than unity. Correspondingly, the GL81 model results in severe overestimation of the geostrophic stress-driven vertical velocity and tracer transport.

Impacts of climate change on estuarine stratification and implications for hypoxia within a shallow subtropical system

Vertical density stratification often plays an important role in the formation and expansion of coastal hypoxic zones through its effect on near-bed circulation and vertical oxygen flux. However, the impact of future climate change on estuarine circulation is widely unknown. Here, we developed and calibrated a three-dimensional hydrodynamic model for Pensacola Bay, a shallow subtropical estuary in the northeastern Gulf of Mexico. Model simulations based on years 2013–2017 were applied to examine changes in salinity, temperature, and density under future climate scenarios, including increased radiative forcing (IR) and temperature (T), increased freshwater discharge (D), sea level rise (SLR), and wind intensification (W). Simulations showed that the impacts of climate change on modeled state variables varied over time with external forcing conditions. The model demonstrated the potential for sea level rise and increased freshwater discharge to episodically increase vertical density gradients in the Bay. However, increased wind forcing destabilized vertical gradients, at times reducing the spatial extent and duration of stable stratification. For time periods with low freshwater discharge, moderate increases in wind speed (10%) can destabilize density gradients strengthened by increased discharge (10%) and sea level rise (0.48 m). In contrast, destruction of strong density gradients that form near the mid-Bay channel following flood events requires stronger wind forcing. These results highlight the importance of considering natural variability in freshwater and wind forcing, as well as local phenomena that are generally unresolved by global climate models.

Spatially heterogeneous dynamics and locally arrested density fluctuations from first principles

We present a first-principles formalism for studying dynamical heterogeneities in glass-forming liquids. Based on the non-equilibrium self-consistent generalized Langevin equation theory, we were able to describe the time-dependent local density profile during the particle interchange among small regions of the fluid. The final form of the diffusion equation contains both the contribution of the chemical potential gradient written in terms of a coarse-grained density and a collective diffusion coefficient as well as the effect of a history-dependent mobility factor. With this diffusion equation, we captured interesting phenomena in glass-forming liquids such as the cases when a strong density gradient is accompanied by a very low mobility factor attributable to the denser part: in such circumstances, the density profile falls into an arrested state even in the presence of a density gradient. On the other hand, we also show that above a certain critical temperature, which depends on the volume fraction, any density heterogeneity relaxes to a uniform state in a finite time, known as equilibration time. We further show that such equilibration time varies little with the temperature in diluted systems but can change drastically with temperature in concentrated systems.

Frontal Dynamics in the Alboran Sea: 1. Coherent 3D Pathways at the Almeria‐Oran Front Using Underwater Glider Observations

AbstractOcean fronts are areas that can support phytoplankton production through fertilization in the sunlit layer and the subduction of biogeochemical properties from the surface to the interior of the ocean. The Almeria‐Oran (AO) front is formed from the juxtaposition of fresh inflowing Atlantic waters and more saline re‐circulating Mediterranean waters. A fleet of three gliders flying in parallel lines was deployed across the AO to obtain observations in the CALYPSO project. These observations were combined with remote sensing and modeling simulations, thus providing a novel approach to identifying the three‐dimensional transport and the submesoscale across‐front circulation. The resulting 33 cross‐front sections reveal spatial and temporal changes in the frontal boundary, with isopycnals steepening and/or relaxing. The observations revealed strong horizontal density gradients (up to ∼1.4 kg m−3) and the spatial variability was observed over different length scales (∼10–45 km). The potential vorticity decreased across the front due to the vorticity component in the horizontal density gradient direction. The predominant cyclonic relative vorticity on the dense side of the AO is associated with downwelling processes. The biogeochemical observations also suggest vertical transport along coherent pathways through baroclinic instability. Phytoplankton biomass enhancement occurs as a result, and is subducted below the euphotic layer. The observed oxygen filaments show upwelling and downwelling, providing a mechanism for oxygenating deeper layers and reducing the ventilation of deep low‐oxygenated waters. Understanding the mechanisms of vertical transport can help us evaluate the dynamics of ocean fronts and their impacts on biological carbon storage.

Turbulence mitigation in maximum-J stellarators with electron-density gradient

In fusion devices, the geometry of the confining magnetic field has a significant impact on the instabilities that drive turbulent heat loss. This is especially true of stellarators, where the density-gradient-driven branch of the ‘trapped electron mode’ (TEM) is predicted to be linearly stable if the magnetic field has the maximum-J property, as is very approximately the case in certain magnetic configurations of the Wendelstein 7-X experiment (W7-X). Here we show, using both analytical theory and simulations, that the benefits of the optimisation of W7-X also serve to mitigate ion-temperature-gradient (ITG) modes as long as an electron density gradient is present. We find that the effect indeed carries over to nonlinear numerical simulations, where W7-X has low TEM-driven transport, and reduced ITG turbulence in the presence of a density gradient, giving theoretical support for the existence of enhanced confinement regimes, in the presence of strong density gradients (e.g. hydrogen pellet or neutral beam injection).

Ambient high-frequency seismic surface waves in the firn column of central west Antarctica

AbstractFirn is the pervasive surface material across Antarctica, and its structures reflect its formation and history in response to environmental perturbations. In addition to the role of firn in thermally isolating underlying glacial ice, it defines near-surface elastic and density structure and strongly influences high-frequency (> 5 Hz) seismic phenomena observed near the surface. We investigate high-frequency seismic data collected with an array of seismographs deployed on the West Antarctic Ice Sheet (WAIS) near WAIS Divide camp in January 2019. Cross-correlations of anthropogenic noise originating from the approximately 5 km-distant camp were constructed using a 1 km-diameter circular array of 22 seismographs. We distinguish three Rayleigh (elastic surface) wave modes at frequencies up to 50 Hz that exhibit systematic spatially varying particle motion characteristics. The horizontal-to-vertical ratio for the second mode shows a spatial pattern of peak frequencies that matches particle motion transitions for both the fundamental and second Rayleigh modes. This pattern is further evident in the appearance of narrow band spectral peaks. We find that shallow lateral structural variations are consistent with these observations, and model spectral peaks as Rayleigh wave amplifications within similarly scaled shallow basin-like structures delineated by the strong velocity and density gradients typical of Antarctic firn.