Domain modeling and grid generation for multi-block structured grids with application to aerodynamic and hydrodynamic configurations

About five years ago, a joint development was started of a flow simulation system for engine-airframe integration studies on propeller as well as jet aircraft. The initial system was based on the Euler equations and made operational for industrial aerodynamic design work. The system consists of three major components: a domain modeller, for the graphical interactive subdivision of flow domains into an unstructured collection of blocks; a grid generator, for the graphical interactive computation of structured grids in blocks; and a flow solver, for the computation of flows on multi-block grids. The industrial partners of the collaboration and NLR have demonstrated that the domain modeller, grid generator and flow solver can be applied to simulate Euler flows around complete aircraft, including propulsion system simulation. Extension to Navier-Stokes flows is in progress. Delft Hydraulics has shown that both the domain modeller and grid generator can also be applied successfully for hydrodynamic configurations. An overview is given about the main aspects of both domain modelling and grid generation

Elliptic surface grid generation on minimal and parmetrized surfaces

An elliptic grid generation method is presented which generates excellent boundary conforming grids in domains in 2D physical space. The method is based on the composition of an algebraic and elliptic transformation. The composite mapping obeys the familiar Poisson grid generation system with control functions specified by the algebraic transformation. New expressions are given for the control functions. Grid orthogonality at the boundary is achieved by modification of the algebraic transformation. It is shown that grid generation on a minimal surface in 3D physical space is in fact equivalent to grid generation in a domain in 2D physical space. A second elliptic grid generation method is presented which generates excellent boundary conforming grids on smooth surfaces. It is assumed that the surfaces are parametrized and that the grid only depends on the shape of the surface and is independent of the parametrization. Concerning surface modeling, it is shown that bicubic Hermite interpolation is an excellent method to generate a smooth surface which is passing through a given discrete set of control points. In contrast to bicubic spline interpolation, there is extra freedom to model the tangent and twist vectors such that spurious oscillations are prevented

Choice of implicit and explicit operators for the upwind differencing method

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76274/1/AIAA-1988-624-513.pd

A non-linear multigrid method for the steady Euler equations

Higher-order accurate Euler-flow solutions are presented for some airfoil test cases. Second-order accurate solutions are computed by an Iterative Defect Correction process. For two test cases even higher accuracy is obtained by the additional use of a ~xtrapolation technique. Finite volume Osher-type discretizations are applied throughout. Two interpolation schemes (one with and one w~hout a flux limiter) are used for the computation of the second-order defect. In each Defect Correction cycle, the solution is computed by non-linear mu~igrid iteration, in which Collective Symmetric Gauss-Seidel relaxation is used as the smoothing procedure. The computational method does not require tuning of parameters. The solutions show a good resolution of discontinuities, and they are obtained at low computational costs. The rate of convergence seems to be grid-independent

Modelling the Innate Immune Response against Avian Influenza Virus in Chicken

At present there is limited understanding of the host immune response to (low pathogenic) avian influenza virus infections in poultry. Here we develop a mathematical model for the innate immune response to avian influenza virus in chicken lung, describing the dynamics of viral load, interferon-α, -β and -γ, lung (i.e. pulmonary) cells and Natural Killer cells. We use recent results from experimentally infected chickens to validate some of the model predictions. The model includes an initial exponential increase of the viral load, which we show to be consistent with experimental data. Using this exponential growth model we show that the duration until a given viral load is reached in experiments with different inoculation doses is consistent with a model assuming a linear relationship between initial viral load and inoculation dose. Subsequent to the exponential-growth phase, the model results show a decline in viral load caused by both target-cell limitation as well as the innate immune response. The model results suggest that the temporal viral load pattern in the lungs displayed in experimental data cannot be explained by target-cell limitation alone. For biologically plausible parameter values the model is able to qualitatively match to data on viral load in chicken lungs up until approximately 4 days post infection. Comparison of model predictions with data on CD107-mediated degranulation of Natural Killer cells yields some discrepancy also for earlier days post infection

Intrinsic gain modulation and adaptive neural coding

In many cases, the computation of a neural system can be reduced to a receptive field, or a set of linear filters, and a thresholding function, or gain curve, which determines the firing probability; this is known as a linear/nonlinear model. In some forms of sensory adaptation, these linear filters and gain curve adjust very rapidly to changes in the variance of a randomly varying driving input. An apparently similar but previously unrelated issue is the observation of gain control by background noise in cortical neurons: the slope of the firing rate vs current (f-I) curve changes with the variance of background random input. Here, we show a direct correspondence between these two observations by relating variance-dependent changes in the gain of f-I curves to characteristics of the changing empirical linear/nonlinear model obtained by sampling. In the case that the underlying system is fixed, we derive relationships relating the change of the gain with respect to both mean and variance with the receptive fields derived from reverse correlation on a white noise stimulus. Using two conductance-based model neurons that display distinct gain modulation properties through a simple change in parameters, we show that coding properties of both these models quantitatively satisfy the predicted relationships. Our results describe how both variance-dependent gain modulation and adaptive neural computation result from intrinsic nonlinearity.Comment: 24 pages, 4 figures, 1 supporting informatio

Compressive properties of min-mod-type limiters in modelling shockwave-containing flows

The long-ignored compressive properties of Min-mod-type limiter is investigated in this manuscript by demonstrating its potential in numerically modelling shockwave-containing flows, especially in shock wave/boundary layer interaction (SWBLI) problems. Theoretical studies were firstly performed based on Sweby’s total variation diminishing (TVD) limiter region and Spekreijse’s monotonicity-preserving limiter region to indicate Min-mod-type limiters’ compressive properties. The influence of limiters on the solution accuracy was evaluated using a hybrid-order analysis method based on the grid-independent study in three typical shockwave-containing flows. The conclusions are that, Min-mod-type limiter can be utilized as a dissipative and/or compressive limiter, but depending on the reasonable value of the compression parameter. The compressive Min-mod limiter tends to be more attractive in modelling shockwave-containing flows as compared to other commonly preferred limiters because of its stable computational process and its high-resolution predictions. However, the compressive Min-mod limiter may suffer from its slightly poor convergence, as that observed in other commonly accepted smooth limiters in modelling SWBLI problems. © 2020, The Author(s)

Modeling convergent ON and OFF pathways in the early visual system

For understanding the computation and function of single neurons in sensory systems, one needs to investigate how sensory stimuli are related to a neuron’s response and which biological mechanisms underlie this relationship. Mathematical models of the stimulus–response relationship have proved very useful in approaching these issues in a systematic, quantitative way. A starting point for many such analyses has been provided by phenomenological “linear–nonlinear” (LN) models, which comprise a linear filter followed by a static nonlinear transformation. The linear filter is often associated with the neuron’s receptive field. However, the structure of the receptive field is generally a result of inputs from many presynaptic neurons, which may form parallel signal processing pathways. In the retina, for example, certain ganglion cells receive excitatory inputs from ON-type as well as OFF-type bipolar cells. Recent experiments have shown that the convergence of these pathways leads to intriguing response characteristics that cannot be captured by a single linear filter. One approach to adjust the LN model to the biological circuit structure is to use multiple parallel filters that capture ON and OFF bipolar inputs. Here, we review these new developments in modeling neuronal responses in the early visual system and provide details about one particular technique for obtaining the required sets of parallel filters from experimental data