Turbulent and Transitional Modeling of Drag on Oceanographic Measurement Devices

Computational fluid dynamic techniques have been applied to the determination of drag on oceanographic devices (expendable bathythermographs). Such devices, which are used to monitor changes in ocean heat content, provide information that is dependent on their drag coefficient. Inaccuracies in drag calculations can impact the estimation of ocean heating associated with global warming. Traditionally, ocean-heating information was based on experimental correlations which related the depth of the device to the fall time. The relation of time-depth is provided by a fall-rate equation (FRE). It is known that FRE depths are reasonably accurate for ocean environments that match the experiments from which the correlations were developed. For other situations, use of the FRE may lead to depth errors that preclude XBTs as accurate oceanographic devices. Here, a CFD approach has been taken which provides drag coefficients that are used to predict depths independent of an FRE

Temporal blocking of finite-difference stencil operators with sparse "off-the-grid" sources

Stencil kernels dominate a range of scientific applications, including seismic and medical imaging, image processing, and neural networks. Temporal blocking is a performance optimization that aims to reduce the required memory bandwidth of stencil computations by re-using data from the cache for multiple time steps. It has already been shown to be beneficial for this class of algorithms. However, applying temporal blocking to practical applications' stencils remains challenging. These computations often consist of sparsely located operators not aligned with the computational grid (“off-the-grid”). Our work is motivated by modelling problems in which source injections result in wavefields that must then be measured at receivers by interpolation from the grided wavefield. The resulting data dependencies make the adoption of temporal blocking much more challenging. We propose a methodology to inspect these data dependencies and reorder the computation, leading to performance gains in stencil codes where temporal blocking has not been applicable. We implement this novel scheme in the Devito domain-specific compiler toolchain. Devito implements a domain-specific language embedded in Python to generate optimized partial differential equation solvers using the finite-difference method from high-level symbolic problem definitions. We evaluate our scheme using isotropic acoustic, anisotropic acoustic, and isotropic elastic wave propagators of industrial significance. After auto-tuning, performance evaluation shows that this enables substantial performance improvement through temporal blocking over highly-optimized vectorized spatially-blocked code of up to 1.6x

Rituximab monitoring and redosing in pediatric neuromyelitis optica spectrum disorder.

Abstract OBJECTIVE: To study rituximab in pediatric neuromyelitis optica (NMO)/NMO spectrum disorders (NMOSD) and the relationship between rituximab, B cell repopulation, and relapses in order to improve rituximab monitoring and redosing. METHODS: Multicenter retrospective study of 16 children with NMO/NMOSD receiving 652 rituximab courses. According to CD19 counts, events during rituximab were categorized as "repopulation," "depletion," or "depletion failure" relapses (repopulation threshold CD19 6510 7 10(6) cells/L). RESULTS: The 16 patients (14 girls; mean age 9.6 years, range 1.8-15.3) had a mean of 6.1 events (range 1-11) during a mean follow-up of 6.1 years (range 1.6-13.6) and received a total of 76 rituximab courses (mean 4.7, range 2-9) in 42.6-year cohort treatment. Before rituximab, 62.5% had received azathioprine, mycophenolate mofetil, or cyclophosphamide. Mean time from rituximab to last documented B cell depletion and first repopulation was 4.5 and 6.8 months, respectively, with large interpatient variability. Earliest repopulations occurred with the lowest doses. Significant reduction between pre- and post-rituximab annualized relapse rate (ARR) was observed (p = 0.003). During rituximab, 6 patients were relapse-free, although 21 relapses occurred in 10 patients, including 13 "repopulation," 3 "depletion," and 4 "depletion failure" relapses. Of the 13 "repopulation" relapses, 4 had CD19 10-50 7 10(6) cells/L, 10 had inadequate monitoring ( 641 CD19 in the 4 months before relapses), and 5 had delayed redosing after repopulation detection. CONCLUSION: Rituximab is effective in relapse prevention, but B cell repopulation creates a risk of relapse. Redosing before B cell repopulation could reduce the relapse risk further. CLASSIFICATION OF EVIDENCE: This study provides Class IV evidence that rituximab significantly reduces ARR in pediatric NMO/NMOSD. This study also demonstrates a relationship between B cell repopulation and relapses

Combining checkpointing and data compression for large scale seismic inversion

Seismic inversion and imaging are adjoint-based optimization problems that processes up to terabytes of data, regularly exceeding the memory capacity of available computers. Data compression is an effective strategy to reduce this memory requirement by a certain factor, particularly if some loss in accuracy is acceptable. A popular alternative is checkpointing, where data is stored at selected points in time, and values at other times are recomputed as needed from the last stored state. This allows arbitrarily large adjoint computations with limited memory, at the cost of additional recomputations. In this paper we combine compression and checkpointing for the first time to compute a realistic seismic inversion. The combination of checkpointing and compression allows larger adjoint computations compared to using only compression, and reduces the recomputation overhead significantly compared to using only checkpointing

Full-waveform inversion, Part 3: Optimization

This tutorial is the third part of a full-waveform inversion (FWI) tutorial series with a step-by-step walkthrough of setting up forward and adjoint wave equations and building a basic FWI inversion framework. For discretizing and solving wave equations, we use Devito (http://www.opesci.org/devito-public), a Python-based domain-specific language for automated generation of finite-difference code (Lange et al., 2016). The first two parts of this tutorial (Louboutin et al., 2017, 2018) demonstrated how to solve the acoustic wave equation for modeling seismic shot records and how to compute the gradient of the FWI objective function using the adjoint-state method. With these two key ingredients, we will now build an inversion framework that can be used to minimize the FWI least-squares objective function

A computational method for determining XBT depths

Abstract. A new technique for determining the depth of expendable bathythermographs (XBTs) is developed. This new method uses a forward-stepping calculation which incorporates all of the forces on the XBT devices during their descent. Of particular note are drag forces which are calculated using a new drag coefficient expression. That expression, obtained entirely from computational fluid dynamic modeling, accounts for local variations in the ocean environment. Consequently, the method allows for accurate determination of depths for any local temperature environment. The results, which are entirely based on numerical simulation, are compared with the experiments of LM Sippican T-5 XBT probes. It is found that the calculated depths differ by less than 3% from depth estimates using the standard fall-rate equation (FRE). Furthermore, the differences decrease with depth. The computational model allows an investigation of the fluid flow patterns along the outer surface of the probe as well as in the interior channel. The simulations take account of complex flow phenomena such as laminar-turbulent transition and flow separation

Scenarios of domain pattern formation in a reaction-diffusion system

We performed an extensive numerical study of a two-dimensional reaction-diffusion system of the activator-inhibitor type in which domain patterns can form. We showed that both multidomain and labyrinthine patterns may form spontaneously as a result of Turing instability. In the stable homogeneous system with the fast inhibitor one can excite both localized and extended patterns by applying a localized stimulus. Depending on the parameters and the excitation level of the system stripes, spots, wriggled stripes, or labyrinthine patterns form. The labyrinthine patterns may be both connected and disconnected. In the the stable homogeneous system with the slow inhibitor one can excite self-replicating spots, breathing patterns, autowaves and turbulence. The parameter regions in which different types of patterns are realized are explained on the basis of the asymptotic theory of instabilities for patterns with sharp interfaces developed by us in Phys. Rev. E. 53, 3101 (1996). The dynamics of the patterns observed in our simulations is very similar to that of the patterns forming in the ferrocyanide-iodate-sulfite reaction.Comment: 15 pages (REVTeX), 15 figures (postscript and gif), submitted to Phys. Rev.

Intra-Arterial Combination Therapy for Experimental Acute Ischemic Stroke

Acute ischemic stroke continues to devastate millions of individuals worldwide. Current treatments work to restore blood flow but not rescue affected tissue. Our goal was to develop a combination of neuroprotective agents administered intra-arterially following recanalization to target ischemic tissue. Using C57Bl/6J male mice, we performed tandem transient ipsilateral middle cerebral/common carotid artery occlusion, followed by immediate intra-arterial pharmacotherapy administration through a standardized protocol. Two pharmacotherapy agents, verapamil and lubeluzole, were selected based on their potential to modulate different aspects of the ischemic cascade; verapamil, a calcium channel blocker, works in an acute fashion blocking L-type calcium channels, whereas lubeluzole, an N-methyl-D-aspartate modulator, works in a delayed fashion blocking intracellular glutamate trafficking. We hypothesized that combination therapy would provide complimentary and potentially synergistic benefit treating brain tissue undergoing various stages of injury. Physiological measurements for heart rate and pulse distention (blood pressure) demonstrated no detrimental effects between groups, suggesting that the combination drug administration is safe. Tissue analysis demonstrated a significant difference between combination and control (saline) groups in infarct volume, neuronal health, and astrogliosis. Although a significant difference in functional outcome was not observed, we did note that the combination treatment group had a greater percent change from baseline in forced motor movement as compared with controls. This study demonstrates the safety and feasibility of intra-arterial combination therapy following successful recanalization and warrants further study

Towards self-verification in finite difference code generation

Code generation from domain-specific languages is becoming increasingly popular as a method to obtain optimised low-level code that performs well on a given platform and for a given problem instance. Ensuring the correctness of generated codes is crucial. At the same time, testing or manual inspection of the code is problematic, as the generated code can be complex and hard to read. Moreover, the generated code may change depending on the problem type, domain size, or target platform, making conventional code review or testing methods impractical. As a solution, we propose the integration of formal verification tools into the code generation process. We present a case study in which the CIVL verification tool is combined with the Devito finite difference framework that generates optimised stencil code for PDE solvers from symbolic equations. We show a selection of properties of the generated code that can be automatically specified and verified during the code generation process. Our approach allowed us to detect a previously unknown bug in the Devito code generation tool

Convective Fingering of an Autocatalytic Reaction Front

We report experimental observations of the convection-driven fingering instability of an iodate-arsenous acid chemical reaction front. The front propagated upward in a vertical slab; the thickness of the slab was varied to control the degree of instability. We observed the onset and subsequent nonlinear evolution of the fingers, which were made visible by a {\it p}H indicator. We measured the spacing of the fingers during their initial stages and compared this to the wavelength of the fastest growing linear mode predicted by the stability analysis of Huang {\it et. al.} [{\it Phys. Rev. E}, {\bf 48}, 4378 (1993), and unpublished]. We find agreement with the thickness dependence predicted by the theory.Comment: 11 pages, RevTex with 3 eps figures. To be published in Phys Rev E, [email protected], [email protected], [email protected]