Emma E. Davies, Hannah T. Rüdisser, Ute V. Amerstorfer, Christian Möstl, Maike Bauer, Eva Weiler, Tanja Amerstorfer, Satabdwa Majumdar, Phillip Hess, Andreas J. Weiss, Martin A. Reiss, Lucie M. Green, David M. Long, Teresa Nieves-Chinchilla, Domenico Trotta, Timothy S. Horbury, Helen O'Brien, Edward Fauchon-Jones, Jean Morris, Christopher J. Owen, et al (2) As both Parker Solar Probe (PSP) and Solar Orbiter (SolO) reach heliocentric distances closer to the Sun, they present an exciting opportunity to study the structure of CMEs in the inner heliosphere. We present an analysis of the global flux rope structure of the 2022 September 5 CME event that impacted PSP at a heliocentric distance of only 0.07 au and SolO at 0.69 au. We compare in situ measurements at PSP and SolO to determine global and local expansion measures, finding a good agreement between magnetic field relationships with heliocentric distance, but significant differences with respect to flux rope size. We use PSP/WISPR images as input to the ELEvoHI model, providing a direct link between remote and in situ observations; we find a large discrepancy between the resulting modeled arrival times, suggesting that the underlying model assumptions may not be suitable when using data obtained close to the Sun, where the drag regime is markedly different in comparison to larger heliocentric distances. Finally, we fit the SolO/MAG and PSP/FIELDS data independently with the 3DCORE model and find that many parameters are consistent between spacecraft, however, challenges are apparent when reconstructing a global 3D structure that aligns with arrival times at PSP and Solar Orbiter, likely due to the large radial and longitudinal separations between spacecraft. From our model results, it is clear the solar wind background speed and drag regime strongly affect the modeled expansion and propagation of CMEs and need to be taken into consideration.
Optimization of functionally graded metamaterial arrays with a high dimensional and continuous geometric design space is cumbersome and could be accelerated via machine learning tools. Mechanical metamaterials can manipulate acoustic or ultrasonic waves by introducing large dispersive and attenuative effects near their natural frequency. In this work functionally graded structures are designed and optimized to combine the energy attenuation performance of a number of unit cells with varying frequency responses and to reduce the interlayer mismatch effects. Optimization through genetic algorithm avoids the many local minima related to high dimensionality of the design space but requires many iterations. A reduced order model (ROM) is applied that can reproduce the transmission response that is traditionally calculated with FEM, in a fraction of the time. Pairing GA and the ROM together, an array of 6 unit cells (with a total of 18 independent geometric design variables) is optimized to have stop bands with extended width and sharper boundaries. Symmetric functionally graded structures are determined to be optimal geometric configurations. Measured 3D printed features are projected onto the ROM solutions to quantify the effect of printing uncertainty on array performance. Repeatability error of $\pm~20$ $\mu$m is determined to reduce the mean depth of the transmission stop band by a factor of $10^2$ and introduce small shifts in center frequency and band width. Proposed methods to improve the resolution of accessible points in the ROM space, reduce sensitivity to geometric uncertainty, and add design freedom include introducing out-of-plane perforations and varying constituent materials using tunable filled resin systems.
Although inertial particle-laden flows occur in a wide range of industrial and natural processes, there is both a lack of fundamental understanding of these flows and continuum-level governing equations needed to predict transport and particle distribution. Towards this effort, the Taylor-Couette flow (TCF) system has been used recently to study the flow behavior of particle-laden fluids under inertia. This article provides an overview of experimental, theoretical, and computational work related to the TCF of neutrally buoyant non-Brownian suspensions, with an emphasis on the effect of finite-sized particles on the series of flow transitions and flow structures. Particles, depending on their size and concentration, cause several significant deviations from Newtonian fluid behavior, including shifting the Reynolds number corresponding to transitions in flow structure and changing the possible structures present in the flow. Furthermore, particles may also migrate depending on the flow structure, leading to hysteretic effects that further complicate the flow behavior. The current state of theoretical and computational modeling efforts to describe the experimental observations is discussed, and suggestions for potential future directions to improve the fundamental understanding of inertial particle-laden flows are provided.
Scattering experiments can be leveraged to extract the effective properties of a heterogeneous metamaterial slab based on multi-point measurements in surrounding media. In this technique, two measurements are made in the ambient media on each side of a finite thickness micro-structured slab to decompose incoming and outgoing waves. The method is applied to an example with locally resonant micro-structured inclusions while paying close attention to parameters that influence the extracted material parameters. It is observed that the extracted overall parameters converge to limiting values when the number of unit cells across the slab thickness or ambient media modulus are increased. Dependence of extracted material parameters on the number of cells through thickness or ambient media properties are attributed to the different response of boundary (skin) and interior cells. A method is presented which represents a finite array with different effective properties for the skin regions vs. the interior regions with great success in reproducing the scattering for different slab thicknesses, and through which the interior regions properties become independent of ambient media properties. In stop bands with lower material loss, challenges arise that are associated with extremely small transmission. The assumption of continuity in transmission phase is enforceable to remove phase ambiguity, although in certain cases (associated with nearly lossless specimens) this assumption appears to fail. In such cases, interference from coupled shear modes may lead to apparent higher longitudinal transmission within the stop band. This work provides a proof of concept for the development of an experimental methodology capable of extracting the effective properties and dispersion behavior of heterogeneous mechanical metamaterials without any knowledge of the internal fields.
Fast forecasting of reservoir pressure distribution in geologic carbon storage (GCS) by assimilating monitoring data is a challenging problem. Due to high drilling cost, GCS projects usually have spatially sparse measurements from wells, leading to high uncertainties in reservoir pressure prediction. To address this challenge, we propose to use low-cost Interferometric Synthetic-Aperture Radar (InSAR) data as monitoring data to infer reservoir pressure build up. We develop a deep learning-accelerated workflow to assimilate surface displacement maps interpreted from InSAR and to forecast dynamic reservoir pressure. Employing an Ensemble Smoother Multiple Data Assimilation (ES-MDA) framework, the workflow updates three-dimensional (3D) geologic properties and predicts reservoir pressure with quantified uncertainties. We use a synthetic commercial-scale GCS model with bimodally distributed permeability and porosity to demonstrate the efficacy of the workflow. A two-step CNN-PCA approach is employed to parameterize the bimodal fields. The computational efficiency of the workflow is boosted by two residual U-Net based surrogate models for surface displacement and reservoir pressure predictions, respectively. The workflow can complete data assimilation and reservoir pressure forecasting in half an hour on a personal computer.
Gyujoon Hwang, Ahmed Abdulkadir, Guray Erus, Mohamad Habes, Raymond Pomponio, Haochang Shou, Jimit Doshi, Elizabeth Mamourian, Tanweer Rashid, Murat Bilgel, Yong Fan, Aristeidis Sotiras, Dhivya Srinivasan, John C. Morris, Daniel Marcus, Marilyn S. Albert, Nick R. Bryan, Susan M. Resnick, Ilya M. Nasrallah, Christos Davatzikos, et al (1) Neuroimaging biomarkers that distinguish between typical brain aging and Alzheimer's disease (AD) are valuable for determining how much each contributes to cognitive decline. Machine learning models can derive multi-variate brain change patterns related to the two processes, including the SPARE-AD (Spatial Patterns of Atrophy for Recognition of Alzheimer's Disease) and SPARE-BA (of Brain Aging) investigated herein. However, substantial overlap between brain regions affected in the two processes confounds measuring them independently. We present a methodology toward disentangling the two. T1-weighted MRI images of 4,054 participants (48-95 years) with AD, mild cognitive impairment (MCI), or cognitively normal (CN) diagnoses from the iSTAGING (Imaging-based coordinate SysTem for AGIng and NeurodeGenerative diseases) consortium were analyzed. First, a subset of AD patients and CN adults were selected based purely on clinical diagnoses to train SPARE-BA1 (regression of age using CN individuals) and SPARE-AD1 (classification of CN versus AD). Second, analogous groups were selected based on clinical and molecular markers to train SPARE-BA2 and SPARE-AD2: amyloid-positive (A+) AD continuum group (consisting of A+AD, A+MCI, and A+ and tau-positive CN individuals) and amyloid-negative (A-) CN group. Finally, the combined group of the AD continuum and A-/CN individuals was used to train SPARE-BA3, with the intention to estimate brain age regardless of AD-related brain changes. Disentangled SPARE models derived brain patterns that were more specific to the two types of the brain changes. Correlation between the SPARE-BA and SPARE-AD was significantly reduced. Correlation of disentangled SPARE-AD was non-inferior to the molecular measurements and to the number of APOE4 alleles, but was less to AD-related psychometric test scores, suggesting contribution of advanced brain aging to these scores.
Fast assimilation of monitoring data to update forecasts of pressure buildup and carbon dioxide (CO2) plume migration under geologic uncertainties is a challenging problem in geologic carbon storage. The high computational cost of data assimilation with a high-dimensional parameter space impedes fast decision-making for commercial-scale reservoir management. We propose to leverage physical understandings of porous medium flow behavior with deep learning techniques to develop a fast history matching-reservoir response forecasting workflow. Applying an Ensemble Smoother Multiple Data Assimilation framework, the workflow updates geologic properties and predicts reservoir performance with quantified uncertainty from pressure history and CO2 plumes interpreted through seismic inversion. As the most computationally expensive component in such a workflow is reservoir simulation, we developed surrogate models to predict dynamic pressure and CO2 plume extents under multi-well injection. The surrogate models employ deep convolutional neural networks, specifically, a wide residual network and a residual U-Net. The workflow is validated against a flat three-dimensional reservoir model representative of a clastic shelf depositional environment. Intelligent treatments are applied to bridge between quantities in a true-3D reservoir model and those in a single-layer reservoir model. The workflow can complete history matching and reservoir forecasting with uncertainty quantification in less than one hour on a mainstream personal workstation.
This paper describes an experimental study of filtration of a colloidal suspension using microfluidic devices. A suspension of micrometer-scale colloids flows through parallel slit-shaped pores at fixed pressure drop. Clogs and cakes are systematically observed at pore entrance, for variable applied pressure drop and ionic strength. Based on image analysis of the layer of colloids close to the device wall, global and local studies are performed to analyse in detail the near-wall layer microstructure. Whereas global porosity of this layer does not seem to be affected by ionic strength and applied pressure drop, a local study shows some heterogeneity: clogs are more porous at the vicinity of the pore than far away. An analysis of medium-range order using radial distribution function shows a slightly more organized state at high ionic strength. This is confirmed by a local analysis using two-dimension continuous wavelet decomposition: the typical size of crystals of colloids is larger for low ionic strength, and it increases with distance from the pores. We bring these results together in a phase diagram involving colloid-colloid repulsive interactions and fluid velocity.
Mark T. Gruneisen, Mark L. Eickhoff, Scott C. Newey, Kurt E. Stoltenberg, Jeffery F. Morris, Michael Bareian, Mark A. Harris, Denis W. Oesch, Michael D. Oliker, Michael B. Flanagan, Brian T. Kay, Jonathan D. Schiller, R. Nicholas Lanning Previous demonstrations of free-space quantum communication in daylight have been touted as significant for the development of global-scale quantum networks. Until now, no one has carefully tuned their atmospheric channel to reproduce the daytime sky radiance and slant-path turbulence conditions as they exist between space and Earth. In this article we report a quantum communication field experiment under conditions representative of daytime downlinks from space. Higher-order adaptive optics increased quantum channel efficiencies far beyond those possible with tip/tilt correction alone while spatial filtering at the diffraction limit rejected optical noise without the need for an ultra-narrow spectral filter. High signal-to-noise probabilities and low quantum-bit-error rates were demonstrated over a wide range of channel radiances and turbulence conditions associated with slant-path propagation in daytime. The benefits to satellite-based quantum key distribution are quantified and discussed.
Quasi-2D experiments of a submerged sediment layer creeping downward were performed, varying the channel tilt and a porous flow under the respective thresholds for yielding. Logarithmic decay rates of the deformation are observed, with the rate increasing with both control parameters. A new dimensionless parameter, $P^*$, accounting for both mean porous flow and gravity force effects on particle motion, allows a collapse of all the deformation results on a single curve. Two distinct creep regimes are identified, and correspond to a systematic change of the void size distribution as $P^*$ increases.
The phenomenon of shear-induced jamming is a factor in the complex rheological behavior of dense suspensions. Such shear-jammed states are fragile, i.e., they are not stable against applied stresses that are incompatible with the stress imposed to create them. This peculiar flow-history dependence of the stress response is due to flow-induced microstructures. To examine jammed states realized under constant shear stress, we perform dynamic simulations of non-Brownian particles with frictional contact forces and hydrodynamic lubrication forces. We find clear signatures that distinguish these fragile states from the more conventional isotropic jammed states.
Predicting rapid and slower soil evolution remains a scientific challenge. This process involves poorly understood aspects of disordered granular matter and dense suspension dynamics. This study presents a novel two-dimensional experiment on a small-scale chip structure; this allows the observation of the deformation at the particle scale of a large-grained sediment bed, under conditions where friction dominate over cohesive and thermal forces, and with an imposed fluid flow. Experiments are performed at conditions which span the particle resuspension criterion, and particle motion is detected and analyzed. The void size population and statistics of particle trajectories bring insight to the sediment dynamics near fluidization conditions. Specifically, particle rearrangement and net bed compaction are observed at flow rates significantly below the criterion for instability growth. Above a threshold, a large vertical channel through the bed forms. In the range of flow rates where channelization can occur, the coexistence of compacting and dilating bed scenarios is observed. The results of the study enhance our capacity for modeling of both slow dynamics and eventual rapid destabilization of sediment beds. Microfluidics channel soil-on-a-chip studies open avenues to new investigations including dissolution-precipitation, fine particles transport, or micro-organisms swimming and population growth, which may depend on mechanics of the porous media itself.
Blockage of pores by particles is found in many processes, including filtration and oil extraction. We present filtration experiments through a linear array of ten channels with one dimension which is sub-micron, through which a dilute dispersion of Brownian polystyrene spheres flows under the action of a fixed pressure drop. The growth rate of a clog formed by particles at a pore entrance systematically increases with the number of already saturated (entirely clogged) pores, indicating that there is an interaction or "cross-talk" between the pores. This observation is interpreted based on a phenomenological model, stating that a diffusive redistribution of particles occurs along the membrane, from clogged to free pores. This one-dimensional model could be extended to two-dimensional membranes.
PROCESS is a reactor systems code - it assesses the engineering and economic viability of a hypothetical fusion power station using simple models of all parts of a reactor system. PROCESS allows the user to choose which constraints to impose and which to ignore, so when evaluating the results it is vital to study the list of constraints used. New algorithms submitted by collaborators can be incorporated - for example safety, first wall erosion, and fatigue life will be crucial and are not yet taken into account. This paper describes algorithms relating to the engineering aspects of the plant. The toroidal field (TF) coils and the central solenoid are assumed by default to be wound from niobium-tin superconductor with the same properties as the ITER conductors. The winding temperature and induced voltage during a quench provide a limit on the current density in the TF coils. Upper limits are placed on the stresses in the structural materials of the TF coil, using a simple two-layer model of the inboard leg of the coil. The thermal efficiency of the plant can be estimated using the maximum coolant temperature, and the capacity factor is derived from estimates of the planned and unplanned downtime, and the duty cycle if the reactor is pulsed. An example of a pulsed power plant is given.
The dynamics of two liquid droplets coalescing in their saturated vapor phase are investigated by Lattice Boltzmann numerical simulations. Attention is paid to the effect of the vapor phase on the formation and growth dynamics of the liquid bridge in the viscous regime. We observe that the onset of the coalescence occurs earlier and the expansion of the bridge initially proceeds faster when the coalescence takes place in a saturated vapor compared to the coalescence in a non-condensable gas. We argue that the initially faster evolution of the coalescence in the saturated vapor is caused by the vapor transport through condensation during the early stages of the coalescence.
Accelerated Stokesian Dynamics (ASD) simulation and a microstructural theory are applied to study structure and the viscosity of hard-sphere Brownian suspensions in active microrheology (MR).We consider moderate to dense suspensions, from near to far from equilibrium conditions.The theory explicitly considers many-body hydrodynamic interactions (HIs) in active MR, and is compared with ASD.Two conditions of moving the probe with constant force (CF) and constant velocity (CV) are considered.The structure is quantified using the probability distribution of colloidal particles around the probe, g(r), which is computed as a solution to the pair Smoluchowski equation (SE) for 0.2<\phi<0.50, and a range of Peclet numbers (Pe), describing the ratio of external force on the probe to thermal forces.Results of ASD and theory demonstrate that a wake zone depleted of bath particles behind the moving probe forms at large Pe, while a boundary-layer accumulation develops upstream.The wake length saturates at Pe>>1 for CF while it continuously grows in CV.This contrast in behavior is related to the dispersion in the motion of the probe under CF conditions, while CV motion has no dispersion.This effect is incorporated in the theory as a force-induced hydrodynamic diffusion flux in the pair SE. We also demonstrate that, despite this difference of structure in CF and CV, g(r) near the probe is set by Pe, for both CF and CV resulting in similar values for their viscosity.Using the theory, the structural anisotropy and Brownian viscosity near equilibrium are shown to be quantitatively similar in both CF and CV motions, which is in contrast with the dilute theory which predict distortions and Brownian viscosities twice as large in CV, relative to CF.This difference arises due to the many-body interactions associated with the equilibrium structure in the moderate to dense regime.
Dynamic particle-scale numerical simulations are used to show that the shear thickening observed in dense colloidal, or Brownian, suspensions is of a similar nature to that observed in non-colloidal suspensions, i.e., a stress-induced transition from a flow of lubricated near-contacting particles to a flow of a frictionally contacting network of particles. Abrupt (or discontinuous) shear thickening is found to be a geometric rather than hydrodynamic phenomenon; it stems from the strong sensitivity of the jamming volume fraction to the nature of contact forces between suspended particles. The thickening obtained in a colloidal suspension of purely hard frictional spheres is qualitatively similar to experimental observations. However, the agreement cannot be made quantitative with only hydrodynamics, frictional contacts and Brownian forces. Therefore the role of a short-range repulsive potential mimicking the stabilization of actual suspensions on the thickening is studied. The effects of Brownian and repulsive forces on the onset stress can be combined in an additive manner. The simulations including Brownian and stabilizing forces show excellent agreement with experimental data for the viscosity $\eta$ and the second normal stress difference $N_2$.
A global signature of the build-up to an intrinsic ELM is found in the phase of signals measured in full flux azimuthal loops in the divertor region of JET. Full flux loop signals provide a global measurement proportional to the voltage induced by changes in poloidal magnetic flux; they are electromagnetically induced by the dynamics of spatially integrated current density. We perform direct time-domain analysis of the high time-resolution full flux loop signals VLD2 and VLD3. We analyze plasmas where a steady H-mode is sustained over several seconds, during which all the observed ELMs are intrinsic; there is no deliberate intent to pace the ELMing process by external means. ELM occurrence times are determined from the Be II emission at the divertor. We previously found that the occurrence times of intrinsic ELMs correlate with specific phases of the VLD2 and VLD3 signals. Here, we investigate how the VLD2 and VLD3 phases vary with time in advance of the ELM occurrence time. We identify a build-up to the ELM in which the VLD2 and VLD3 signals progressively align to the phase at which ELMs preferentially occur, on a ~ 2 -5ms timescale. At the same time, the VLD2 and VLD3 signals become phase synchronized with each other, consistent with the emergence of coherent global dynamics in the integrated current density. In a plasma that remains close to a global magnetic equilibrium, this can reflect bulk displacement or motion of the plasma. This build-up signature to an intrinsic ELM can be extracted from a time interval of data that does not extend beyond the ELM occurrence time, so that these full flux loop signals could assist in ELM prediction or mitigation.
The discontinuous shear thickening (DST) of dense suspensions is a remarkable phenomenon in which the viscosity can increase by several orders of magnitude at a critical shear rate. It has the appearance of a first order phase transition between two hypothetical "states" that we have recently identified as Stokes flows with lubricated or frictional contacts, respectively. Here we extend the analogy further by means of novel stress-controlled simulations and show the existence of a non-monotonic steady-state flow curve analogous to a non-monotonic equation of state. While we associate DST with an S-shaped flow curve, at volume fractions above the shear jamming transition the frictional state loses flowability and the flow curve reduces to an arch, permitting the system to flow only at small stresses. Whereas a thermodynamic transition leads to phase separation in the coexistence region, we observe a uniform shear flow all along the thickening transition. A stability analysis suggests that uniform shear may be mechanically stable for the small Reynolds numbers and system sizes in a rheometer.
A unique experiment in the Joint European Torus (JET) consecutively produced 120 almost identical plasma pulses, providing two orders of magnitude more data than is usually available. This allows the statistical detection of previously unobservable phenomena such as a sequence of resonant-like waiting times between edge-localised instabilities (ELMs). Here we investigate the causes of this phenomenon. By synchronising data to the 1000s of ELM times and averaging the results, random errors are reduced by a factor of 50, allowing unprecedentedly detailed behaviour to be described. A clear link can then be observed between plasma confinement, ELM occurrence, vertical plasma oscillations, and an otherwise unobservable oscillation in a control coil current that is not usually associated with ELM occurrence. The results suggest a strong and unanticipated edge-plasma dependence on control system behaviour.
The microstructure and rheological properties of suspensions of neutrally buoyant hard spherical particles in Newtonian fluid under conditions of finite inertia are studied using the lattice-Boltzmann method (LBM), which is based on a discrete Boltzmann model for the fluid and Newtonian dynamics for the particles. The suspensions are subjected to simple-shear flow and the properties are studied as a function of Reynolds number and volume fraction, $\phi$. The inertia is characterized by the particle-scale shear flow Reynolds number $Re = \frac{\rho \dot{\gamma}a^{2}}{\mu}$, where $a$ is the particle radius, $\dot{\gamma}$ is the shear rate and $\rho$ and $\mu$ are the density and viscosity of the fluid, respectively. The influences of inertia and of the volume fraction are studied for $0.005\leqslant Re \leqslant 5$ and $0.1\leqslant \phi \leqslant 0.35$. The flow-induced microstructure is studied using the pair distribution function $g(\boldsymbol{r})$. Different stress mechanisms, including those due to surface tractions (stresslet), acceleration, and the Reynolds stress due to velocity fluctuations are computed and their influence on the first and second normal stress differences, the particle pressure and the viscosity of the suspensions are detailed. The probability density functions of particle force and torque are also presented.
Particles suspended in a Newtonian fluid raise the viscosity and also generally give rise to a shear-rate dependent rheology. In particular, pronounced shear thickening may be observed at large solid volume fractions. In a recent article (R. Seto, R. Mari, J. F. Morris, and M. M. Denn., Phys. Rev. Lett., 111:218301, 2013) we have considered the minimum set of components to reproduce the experimentally observed shear thickening behavior, including Discontinuous Shear Thickening (DST). We have found frictional contact forces to be essential, and were able to reproduce the experimental behavior by a simulation including this physical ingredient along with viscous lubrication. In the present article, we thoroughly investigate the effect of friction and express it in the framework of the jamming transition. The viscosity divergence at the jamming transition has been a well known phenomenon in suspension rheology, as reflected in many empirical laws for the viscosity. Friction can affect this divergence, and in particular the jamming packing fraction is reduced if particles are frictional. Within the physical description proposed here, shear thickening is a direct consequence of this effect: as the shear rate increases, friction is increasingly incorporated as more contacts form, leading to a transition from a mostly frictionless to a mostly frictional rheology. This result is significant because it shifts the emphasis from lubrication hydrodynamics and detailed microscopic interactions to geometry and steric constraints close to the jamming transition.
A.J. Webster, R.O. Dendy, F.A. Calderon, S.C. Chapman, E. Delabie, D. Dodt, R. Felton, T.N. Todd, F. Maviglia, J. Morris, V. Riccardo, B. Alper, S Brezinsek, P. Coad, J. Likonen, M. Rubel, JET EFDA Contributors For a two week period during the Joint European Torus (JET) 2012 experimental campaign, the same high confinement plasma was repeated 151 times. The dataset was analysed to produce a probability density function (pdf) for the waiting times between edge-localised plasma instabilities ("ELMS"). The result was entirely unexpected. Instead of a smooth single peaked pdf, a succession of 4-5 sharp maxima and minima uniformly separated by 7-8 millisecond intervals was found. Here we explore the causes of this newly observed phenomenon, and conclude that it is either due to a self-organised plasma phenomenon or an interaction between the plasma and a real-time control system. If the maxima are a result of "resonant" frequencies at which ELMs can be triggered more easily, then future ELM control techniques can, and probably will, use them. Either way, these results demand a deeper understanding of the ELMing process.
D.Allan, C.Andreopoulos, C.Angelsen, G.J.Barker, G.Barr, S.Bentham, I.Bertram, S.Boyd, K.Briggs, R.G.Calland, J.Carroll, S.L.Cartwright, A.Carver, C.Chavez, G.Christodoulou, J.Coleman, P.Cooke, G.Davies, C.Densham, F.Di Lodovico, et al (89) The T2K experiment studies oscillations of an off-axis muon neutrino beam between the J-PARC accelerator complex and the Super-Kamiokande detector. Special emphasis is placed on measuring the mixing angle theta_13 by observing electron neutrino appearance via the sub-dominant muon neutrino to electron neutrino oscillation, and searching for CP violation in the lepton sector. The experiment includes a sophisticated, off-axis, near detector, the ND280, situated 280 m downstream of the neutrino production target in order to measure the properties of the neutrino beam and to understand better neutrino interactions at the energy scale below a few GeV. The data collected with the ND280 are used to study charged- and neutral-current neutrino interaction rates and kinematics prior to oscillation, in order to reduce uncertainties in the oscillation measurements by the far detector. A key element of the near detector is the ND280 electromagnetic calorimeter (ECal), consisting of active scintillator bars sandwiched between lead sheets and read outwith multi-pixel photon counters (MPPCs). The ECal is vital to the reconstruction of neutral particles, and the identification of charged particle species. The ECal surrounds the Pi-0 detector (P0D) and the tracking region of the ND280, and is enclosed in the former UA1/NOMAD dipole magnet. This paper describes the design, construction and assembly of the ECal, as well as the materials from which it is composed. The electronic and data acquisition (DAQ) systems are discussed, and performance of the ECal modules, as deduced from measurements with particle beams, cosmic rays, the calibration system, and T2K data, is described.
Discontinuous shear thickening (DST) observed in many dense athermal suspensions has proven difficult to understand and to reproduce by numerical simulation. By introducing a numerical scheme including both relevant hydrodynamic interactions and granularlike contacts, we show that contact friction is essential for having DST. Above a critical volume fraction, we observe the existence of two states: a low viscosity, contactless (hence, frictionless) state, and a high viscosity frictional shear jammed state. These two states are separated by a critical shear stress, associated with a critical shear rate where DST occurs. The shear jammed state is reminiscent of the jamming phase of granular matter. Continuous shear thickening is seen as a lower volume fraction vestige of the jamming transition.
SuperB Collaboration, M. Baszczyk, P. Dorosz, J. Kolodziej, W. Kucewicz, M. Sapor, A. Jeremie, E. Grauges Pous, G. E. Bruno, G. De Robertis, D. Diacono, G. Donvito, P. Fusco, F. Gargano, F. Giordano, F. Loddo, F. Loparco, G. P. Maggi, V. Manzari, M. N. Mazziotta, et al (389) In this Technical Design Report (TDR) we describe the SuperB detector that was to be installed on the SuperB e+e- high luminosity collider. The SuperB asymmetric collider, which was to be constructed on the Tor Vergata campus near the INFN Frascati National Laboratory, was designed to operate both at the Upsilon(4S) center-of-mass energy with a luminosity of 10^36 cm^-2s^-1 and at the tau/charm production threshold with a luminosity of 10^35 cm^-2s^-1. This high luminosity, producing a data sample about a factor 100 larger than present B Factories, would allow investigation of new physics effects in rare decays, CP Violation and Lepton Flavour Violation. This document details the detector design presented in the Conceptual Design Report (CDR) in 2007. The R&D and engineering studies performed to arrive at the full detector design are described, and an updated cost estimate is presented. A combination of a more realistic cost estimates and the unavailability of funds due of the global economic climate led to a formal cancelation of the project on Nov 27, 2012.
The adsorption dynamics of a colloidal particle at a fluid interface is studied theoretically and numerically, documenting distinctly different relaxation regimes. The adsorption of a perfectly smooth particle is characterized by a fast exponential relaxation to thermodynamic equilibrium where the interfacial free energy has a minimum. The short relaxation time is given by the ratio of viscous damping to capillary forces. Physical and/or chemical heterogeneities in a colloidal system, however, can result in multiple minima of the free energy giving rise to metastability. In the presence of metastable states we observe a crossover to a slow logarithmic relaxation reminiscent of physical aging in glassy systems. The long relaxation time is determined by the thermally-activated escape rate from metastable states. Analytical expressions derived in this work yield quantitative agreement with molecular dynamics simulations and recent experimental observations. This work provides new insights on the adsorption dynamics of colloidal particles at fluid interfaces.
M. Adinolfi, G. Aglieri Rinella, E. Albrecht, T. Bellunato, S. Benson, T. Blake, C. Blanks, S. Brisbane, N.H. Brook, M. Calvi, B. Cameron, R. Cardinale, L. Carson, A. Contu, M. Coombes, C. D'Ambrosio, S. Easo, U. Egede, S. Eisenhardt, E. Fanchini, et al (77) The LHCb experiment has been taking data at the Large Hadron Collider (LHC) at CERN since the end of 2009. One of its key detector components is the Ring-Imaging Cherenkov (RICH) system. This provides charged particle identification over a wide momentum range, from 2-100 GeV/c. The operation and control software, and online monitoring of the RICH system are described. The particle identification performance is presented, as measured using data from the LHC. Excellent separation of hadronic particle types (pion, kaon and proton) is achieved.
We study numerically the hydrodynamics of dip coating from a suspension and report a mechanism for colloidal assembly and pattern formation on smooth and uniform substrates. Below a critical withdrawal speed of the substrate, capillary forces required to deform the meniscus prevent colloidal particles from entering the coating film. Capillary forces are overcome by hydrodynamic drag only after a minimum number of particles organize in a close-packed formation within the meniscus. Once within the film, the formed assembly moves at nearly the withdrawal speed and rapidly separates from the next assembly. The interplay between hydrodynamic and capillary forces can thus produce periodic and regular structures within the curved meniscus that extends below the withdrawn film. The hydrodynamically-driven assembly documented here is consistent with stripe pattern formations observed experimentally in the so-called thin-film entrainment regime.
In order to maintain efficient tracking in the forward region of H1 after the luminosity upgrade of the HERA machine, the H1 Forward Track Detector was also upgraded. While much of the original software and techniques used for the HERA I phase could be reused, the software for pattern recognition was completely rewritten. This, along with several other improvements in hit finding and high-level track reconstruction, are described in detail together with a summary of the performance of the detector.
This sequence of fluid dynamics videos illustrates the behavior of a suspension of noncolloidal particles flowing past various bluff body obstacles within a microfluidic device. The polystyrene particles, of 7 \mum diameter and volume fraction of 8.4%, are carefully made neutrally buoyant with the suspending liquid composed of a mixture of water and a small fraction of glycerol. The channel depth is 60 \mum and the typical length of the obstacles normal to the flow direction is 200 \mum. The flow rate is varied to generate Reynolds numbers based on the scale of the obstacle in the approximate range 60 < Re < 500; the narrow dimension in the depth direction suppresses onset of unsteadiness and vortex shedding, so that the flows studied are found to be steady (aside from particle-scale fluctuations). Particles are observed to be depleted in the wake region of the obstacle. In certain cases, the entire wake is clear of particles; in other cases there is a portion of the wake in which particles recirculate while a portion of the wake is completely devoid of particles. Experimental observations reveal that if particles are forced into an initially particle-depleted region, they will eventually leave and will bring the wake to its original state, implying these are steady-state distributions.
We investigate experimentally the influence of suspended particles on the transition to turbulence. The particles are monodisperse and neutrally-buoyant with the liquid. The role of the particles on the transition depends both upon the pipe to particle diameter ratios and the concentration. For large pipe-to-particle diameter ratios the transition is delayed while it is lowered for small ratios. A scaling is proposed to collapse the departure from the critical Reynolds number for pure fluid as a function of concentration into a single master curve.
The Relativistic Heavy Ion Collider (RHIC) at Brookhaven has completed nearly two years of successful commissioning and operation. The control system is briefly reviewed and its contribution to the RHIC effort is analyzed, with emphasis on special challenges and innovative design: continuing efforts are also discussed.
Reliable and reproducible performance of the more than 800 Relativistic Heavy Ion Collider (RHIC) magnet power supplies is essential to successful RHIC operation. In order to support power supply commissioning, a system was developed to capture detailed power supply measurements from all the RHIC ring power supplies during acceleration ramps. Diagnostic tools were developed to allow experts to assess ramp reproducibility and rapidly identify problems. The system has now become a routine part of RHIC operations, with data captured for every acceleration ramp. This paper describes the RHIC power supply ramp diagnostic system and considers its impact on RHIC operations.
The movement of the Relativistic Heavy Ion Collider (RHIC) through its various states (eg. injection, acceleration, storage, collisions) is controlled by an application called the Sequencer. This program orchestrates most magnet and instrumentation systems and is responsible for the coordinated acquisition and saving of data from various systems. The Sequencer system, its software infrastructure, support programs, and the language used to drive it are discussed in this paper. Initial operational experience is also described.