Numerical Description of Dilute Particle-Laden FLows by a Quadrature-Based Moment Method

The numerical simulation of gas-particle flows is divided into two families of methods. In Euler-Lagrange methods individual particle trajectories are computed, whereas in Euler-Euler methods particles are characterized by statistical descriptors. Lagrangian methods are very precise but their computational cost increases with instationarity and particle volume fraction. In Eulerian methods (also called moment methods) the particle-phase computational cost is comparable to that of the fluid phase but requires strong simplificaions. Existing Eulerian models consider unimodal or close-to-equilibrium particle velocity distributions and then fail when the actual distribution is far from equilibrium. Quadrature-based Eulerian methods introduce a new reconstruction of the velocity distribution, written as a sum of delta functions in phase space constrained to give the right values for selected low-order moments. Two of the quadrature-based Eulerian methods, differing by the reconstruction algorithm, are the focus of this work. Computational results for two academic cases (crossing jets, Taylor-Green flow) are compared to those of a Lagrangian method (considered as the reference solution) and of an existing second-order moment method. With the quadrature-based Eulerian methods, significant qualitative improvement is noticed compared to the second-order moment method in the two test cases

Swept-Wing Ice Accretion Characterization and Aerodynamics

NASA, FAA, ONERA, the University of Illinois and Boeing have embarked on a significant, collaborative research effort to address the technical challenges associated with icing on large-scale, three-dimensional swept wings. The overall goal is to improve the fidelity of experimental and computational simulation methods for swept-wing ice accretion formation and resulting aerodynamic effect. A seven-phase research effort has been designed that incorporates ice-accretion and aerodynamic experiments and computational simulations. As the baseline, full-scale, swept-wing-reference geometry, this research will utilize the 65% scale Common Research Model configuration. Ice-accretion testing will be conducted in the NASA Icing Research Tunnel for three hybrid swept-wing models representing the 20%, 64% and 83% semispan stations of the baseline-reference wing. Three-dimensional measurement techniques are being developed and validated to document the experimental ice-accretion geometries. Artificial ice shapes of varying geometric fidelity will be developed for aerodynamic testing over a large Reynolds number range in the ONERA F1 pressurized wind tunnel and in a smaller-scale atmospheric wind tunnel. Concurrent research will be conducted to explore and further develop the use of computational simulation tools for ice accretion and aerodynamics on swept wings. The combined results of this research effort will result in an improved understanding of the ice formation and aerodynamic effects on swept wings. The purpose of this paper is to describe this research effort in more detail and report on the current results and status to date.

A 3D finite-volume integral boundary layer method for icing applications

A three-dimensional integral boundary layer code was developed to allow fast computations of boundary layer flows for the purpose of ice accretion modelling. The model is derived in this paper. It is based on a surface Finite-Volume approach. The unsteady equations of momentum deficit and kinetic energy deficit are solved until convergence is reached, preventing from specifying explicitly the stagnation point or separation line. A validation of the code is also presented in the present article. First, the 3D solver is cross-checked against a 2D solver on test cases of self-similar flows and on a NACA0012 configuration. The modelling of the effects of three-dimensionality is also assessed on a self-similar flow test-case. Moreover, the use of unstructured grids is also validated. Finally, an example of the use of the code for the computation of ice accretion is presented

Dynamics and statistics of heavy particles in turbulent flows

We present the results of Direct Numerical Simulations (DNS) of turbulent flows seeded with millions of passive inertial particles. The maximum Taylor's Reynolds number is around 200. We consider particles much heavier than the carrier flow in the limit when the Stokes drag force dominates their dynamical evolution. We discuss both the transient and the stationary regimes. In the transient regime, we study the growt of inhomogeneities in the particle spatial distribution driven by the preferential concentration out of intense vortex filaments. In the stationary regime, we study the acceleration fluctuations as a function of the Stokes number in the range [0.16:3.3]. We also compare our results with those of pure fluid tracers (St=0) and we find a critical behavior of inertia for small Stokes values. Starting from the pure monodisperse statistics we also characterize polydisperse suspensions with a given mean Stokes.Comment: 13 pages, 10 figures, 2 table

NOXA-Induced Alterations in the Bax/Smac Axis Enhance Sensitivity of Ovarian Cancer Cells to Cisplatin

Ovarian cancer is the most common cause of death from gynecologic malignancy. Deregulation of p53 and/or p73-associated apoptotic pathways contribute to the platinum-based resistance in ovarian cancer. NOXA, a pro-apoptotic BH3-only protein, is identified as a transcription target of p53 and/or p73. In this study, we found that genetic variants of Bcl-2 proteins exist among cisplatin-sensitive and -resistant ovarian cancer cells, and the responses of NOXA and Bax to cisplatin are regulated mainly by p53. We further evaluated the effect of NOXA on cisplatin. NOXA induced apoptosis and sensitized A2780s and SKOV3 cells to cisplatin in vitro and in vivo. The effects were mediated by elevated Bax expression, enhanced caspase activation, release of Cyt C and Smac into the cytosol. Furthermore, gene silencing of Bax or Smac significantly attenuated NOXA and/or cisplatin-induced apoptosis in chemosensitive A2780s cells, whereas overexpression of Bax or addition of Smac-N7 peptide significantly increased NOXA and/or cisplatin-induced apoptosis in chemoresistant SKOV3 cells. To our knowledge, these data suggest a new mechanism by which NOXA chemosensitized ovarian cancer cells to cisplatin by inducing alterations in the Bax/Smac axis. Taken together, our findings show that NOXA is potentially useful as a chemosensitizer in ovarian cancer therapy

Velocity-space sensitivity of the time-of-flight neutron spectrometer at JET

The velocity-space sensitivities of fast-ion diagnostics are often described by so-called weight functions. Recently, we formulated weight functions showing the velocity-space sensitivity of the often dominant beam-target part of neutron energy spectra. These weight functions for neutron emission spectrometry (NES) are independent of the particular NES diagnostic. Here we apply these NES weight functions to the time-of-flight spectrometer TOFOR at JET. By taking the instrumental response function of TOFOR into account, we calculate time-of-flight NES weight functions that enable us to directly determine the velocity-space sensitivity of a given part of a measured time-of-flight spectrum from TOFOR

Relationship of edge localized mode burst times with divertor flux loop signal phase in JET

A phase relationship is identified between sequential edge localized modes (ELMs) occurrence times in a set of H-mode tokamak plasmas to the voltage measured in full flux azimuthal loops in the divertor region. We focus on plasmas in the Joint European Torus where a steady H-mode is sustained over several seconds, during which ELMs are observed in the Be II emission at the divertor. The ELMs analysed arise from intrinsic ELMing, in that there is no deliberate intent to control the ELMing process by external means. We use ELM timings derived from the Be II signal to perform direct time domain analysis of the full flux loop VLD2 and VLD3 signals, which provide a high cadence global measurement proportional to the voltage induced by changes in poloidal magnetic flux. Specifically, we examine how the time interval between pairs of successive ELMs is linked to the time-evolving phase of the full flux loop signals. Each ELM produces a clear early pulse in the full flux loop signals, whose peak time is used to condition our analysis. The arrival time of the following ELM, relative to this pulse, is found to fall into one of two categories: (i) prompt ELMs, which are directly paced by the initial response seen in the flux loop signals; and (ii) all other ELMs, which occur after the initial response of the full flux loop signals has decayed in amplitude. The times at which ELMs in category (ii) occur, relative to the first ELM of the pair, are clustered at times when the instantaneous phase of the full flux loop signal is close to its value at the time of the first ELM