Shortening design time through multiplatform simulations with a portable OpenCL golden-model: the LDPC decoder case

Hardware designers and engineers typically need to explore a multi-parametric design space in order to find the best configuration for their designs using simulations that can take weeks to months to complete. For example, designers of special purpose chips need to explore parameters such as the optimal bitwidth and data representation. This is the case for the development of complex algorithms such as Low-Density Parity-Check (LDPC) decoders used in modern communication systems. Currently, high-performance computing offers a wide set of acceleration options, that range from multicore CPUs to graphics processing units (GPUs) and FPGAs. Depending on the simulation requirements, the ideal architecture to use can vary. In this paper we propose a new design flow based on OpenCL, a unified multiplatform programming model, which accelerates LDPC decoding simulations, thereby significantly reducing architectural exploration and design time. OpenCL-based parallel kernels are used without modifications or code tuning on multicore CPUs, GPUs and FPGAs. We use SOpenCL (Silicon to OpenCL), a tool that automatically converts OpenCL kernels to RTL for mapping the simulations into FPGAs. To the best of our knowledge, this is the first time that a single, unmodified OpenCL code is used to target those three different platforms. We show that, depending on the design parameters to be explored in the simulation, on the dimension and phase of the design, the GPU or the FPGA may suit different purposes more conveniently, providing different acceleration factors. For example, although simulations can typically execute more than 3x faster on FPGAs than on GPUs, the overhead of circuit synthesis often outweighs the benefits of FPGA-accelerated execution.

Laboratory formation of a scaled protostellar jet by coaligned poloidal magnetic field

Although bipolar jets are seen emerging from a wide variety of astrophysical systems, the issue of their formation and morphology beyond their launching is still under study. Our scaled laboratory experiments, representative of young stellar object outflows, reveal that stable and narrow collimation of the entire flow can result from the presence of a poloidal magnetic field whose strength is consistent with observations. The laboratory plasma becomes focused with an interior cavity. This gives rise to a standing conical shock from which the jet emerges. Following simulations of the process at the full astrophysical scale, we conclude that it can also explain recently discovered x-ray emission features observed in low-density regions at the base of protostellar jets, such as the well-studied jet HH 154.

Multi-Dimensional Simulations of the Expanding Supernova Remnant of SN 1987A

The expanding remnant from SN 1987A is an excellent laboratory for investigating the physics of supernovae explosions. There is still a large number of outstanding questions, such as the reason for the asymmetric radio morphology, the structure of the pre-supernova environment, and the efficiency of particle acceleration at the supernova shock. We explore these questions using three-dimensional simulations of the expanding remnant between days 820 and 10,000 after the supernova. We combine a hydrodynamical simulation with semi-analytic treatments of diffusive shock acceleration and magnetic field amplification to derive radio emission as part of an inverse problem. Simulations show that an asymmetric explosion, combined with magnetic field amplification at the expanding shock, is able to replicate the persistent one-sided radio morphology of the remnant. We use an asymmetric Truelove & McKee progenitor with an envelope mass of 10 M-circle dot and an energy of 1.5 x 10(44) J. A termination shock in the progenitor's stellar wind at a distance of 0 ''.43-0 ''.51 provides a good fit to the turn on of radio emission around day 1200. For the H II region, a minimum distance of 0 ''.63 +/- 0 ''.01 and maximum particle number density of (7.11 +/- 1.78) x 10(7) m(-3) produces a good fit to the evolving average radius and velocity of the expanding shocks from day 2000 to day 7000 after explosion. The model predicts a noticeable reduction, and possibly a temporary reversal, in the asymmetric radio morphology of the remnant after day 7000, when the forward shock left the eastern lobe of the equatorial ring.

Automatic extraction of Reduced order models from CFD simulations for building energy modeling

Thermal comfort is defined as “that condition of mind which expresses satisfaction with the thermal environment’ [1] [2]. Field studies have been completed in order to establish the governing conditions for thermal comfort [3]. These studies showed that the internal climate of a room was the strongest factor in establishing thermal comfort. Direct manipulation of the internal climate is necessary to retain an acceptable level of thermal comfort. In order for Building Energy Management Systems (BEMS) strategies to be efficiently utilised it is necessary to have the ability to predict the effect that activating a heating/cooling source (radiators, windows and doors) will have on the room. The numerical modelling of the domain can be challenging due to necessity to capture temperature stratification and/or different heat sources (radiators, computers and human beings). Computational Fluid Dynamic (CFD) models are usually utilised for this function because they provide the level of details required. Although they provide the necessary level of accuracy these models tend to be highly computationally expensive especially when transient behaviour needs to be analysed. Consequently they cannot be integrated in BEMS. This paper presents and describes validation of a CFD-ROM method for real-time simulations of building thermal performance. The CFD-ROM method involves the automatic extraction and solution of reduced order models (ROMs) from validated CFD simulations. The test case used in this work is a room of the Environmental Research Institute (ERI) Building at the University College Cork (UCC). ROMs have shown that they are sufficiently accurate with a total error of less than 1% and successfully retain a satisfactory representation of the phenomena modelled. The number of zones in a ROM defines the size and complexity of that ROM. It has been observed that ROMs with a higher number of zones produce more accurate results. As each ROM has a time to solution of less than 20 seconds they can be integrated into the BEMS of a building which opens the potential to real time physics based building energy modelling.

Effects of resistive magnetic field on fast electron divergence measured in experiments

Transport of fast electrons driven by an ultraintense laser through a tracer layer buried in solid targets is studied by particle-in-cell simulations. It is found that intense resistive magnetic fields, having a magnitude of several thousand Tesla, are generated at the interfaces of the materials due to the steep resistivity gradient between the target and tracer layer. Such magnetic fields can significantly inhibit the fast electron propagation. The electrons that can penetrate the first interface are mostly confined in the buried layer by the magnetic fields and cause heating of the tracer layer. The lateral extent of the heated region can be significantly larger than that of the relativistic electron beam. This finding suggests that the relativistic electron divergence inferred from Ká x-ray emission in experiments might be overestimated.

A Review of Wave Makers for 3D Numerical Simulations

Numerical methods have enabled the simulation of complex problems in off-shore and marine engineering. A significant challenge in these simulations is the creation of a realistic wave field. A good numerical tank requires wave creation and absorption of waves at various locations. Several numerical wavemakers with these capabilities have been presented in the past. This paper reviews four different wave-maker methods and discusses limitations, computational efficiency, requirements on the mesh and preprocessing and complexity of implementation.

Alfvén Waves in Simulations of Solar Photospheric Vortices

Using advanced numerical magneto-hydrodynamic simulations of the magnetized solar photosphere, including non-gray radiative transport and a non-ideal equation of state, we analyze plasma motions in photospheric magnetic vortices. We demonstrate that apparent vortex-like motions in photospheric magnetic field concentrations do not exhibit "tornado"-like behavior or a "bath-tub" effect. While at each time instance the velocity field lines in the upper layers of the solar photosphere show swirls, the test particles moving with the time-dependent velocity field do not demonstrate such structures. Instead, they move in a wave-like fashion with rapidly changing and oscillating velocity field, determined mainly by magnetic tension in the magnetized intergranular downflows. Using time-distance diagrams, we identify horizontal motions in the magnetic flux tubes as torsional Alfvén perturbations propagating along the nearly vertical magnetic field lines with local Alfvén speed.

Field testing a full-scale tidal turbine Part 2: In-line Wake Effects

Recent research has shown that higher ambient turbulence leads to better wake recovery, so turbines could be installed in closer proximity in real tidal flows than might be assumed from typical towing tank tests that do not take into account turbulent inflow conditions. The standard tools to assess flow velocities in field conditions are Doppler based sonar devices, such as Acoustic Doppler Profilers (ADPs) or Acoustic Doppler Velocimeters (ADVs). The use of these devices poses some challenges when assessing the wake of a tidal turbine. While ADPs allow the three-dimensional measurement of a velocity profile over a distance, the data is calculated as a mean of three diverging beams and with low temporal resolution. ADVs can measure with higher sampling frequency but only at a single point in the flow. During the MaRINET testing of the SCHOTTELSIT turbine at the QUB tidal test site in Portaferry, Northern Ireland, ADP and ADV measurements were successfully tested.Two methods were employed for measuring the wake: firstly, with a rigidly mounted ADP and secondly, with a submerged ADV which was streamed behind the turbine. This paper presents the experimental set-up and results and discusses limitations and challenges of the two methods used.

Mechanisms for mhd poynting flux generation in simulations of solar photospheric magnetoconvection

We investigate the generation mechanisms of MHD Poynting flux in the magnetized solar photosphere. Using radiative MHD modeling of the solar photosphere with initial magnetic configurations that differ in their field strength and geometry, we show the presence of two different mechanisms for MHD Poynting flux generation in simulations of solar photospheric magnetoconvection. The weaker mechanism is connected to vertical transport of weak horizontal magnetic fields in the convectively stable layers of the upper photosphere, while the stronger is the production of Poynting flux in strongly magnetized intergranular lanes experiencing horizontal vortex motions. These mechanisms may be responsible for the energy transport from the solar convection zone to the higher layers of the solar atmosphere.

A Far-field Pattern Separation Approach for the Synthesis of Directional Modulation Transmitter Arrays (Invited)

In this paper a far-field power pattern separation approach is proposed for the synthesis of directional modulation (DM) transmitter arrays. Separation into information pattern and interference patterns is enabled by far-field pattern null steering. Compared with other DM synthesis methods, e.g., BER-driven DM optimization and orthogonal vector injection, this approach facilitates manipulation of artificial interference spatial distributions. With such capability more interference power can be projected into those most vulnerable to eavesdropping spatial directions in free space, i.e., the information sidelobes. In such a fashion information leaked through radiation sidelobes can be effectively mitigated in a transmitter power efficient manner. The proposed synthesis approach is further validated via bit error rate (BER) simulations.

Shocks in unmagnetized plasma with a shear flow: Stability and magnetic field generation

A pair of curved shocks in a collisionless plasma is examined with a two-dimensional particle-in-cell simulation. The shocks are created by the collision of two electron-ion clouds at a speed that exceeds everywhere the threshold speed for shock formation. A variation of the collision speed along the initially planar collision boundary, which is comparable to the ion acoustic speed, yields a curvature of the shock that increases with time. The spatially varying Mach number of the shocks results in a variation of the downstream density in the direction along the shock boundary. This variation is eventually equilibrated by the thermal diffusion of ions. The pair of shocks is stable for tens of inverse ion plasma frequencies. The angle between the mean flow velocity vector of the inflowing upstream plasma and the shock's electrostatic field increases steadily during this time. The disalignment of both vectors gives rise to a rotational electron flow, which yields the growth of magnetic field patches that are coherent over tens of electron skin depths.

Resistivity and Induced Polarization Monitoring of Biogas combined with Microbial Ecology on a Brown field Site

The accumulation of biogenic greenhouse gases (methane, carbon dioxide) in organic sediments is an important factor in the redevelopment and risk management of many brownfield sites. Good practice with brownfield site characterization requires the identification of free-gas phases and pathways that allow its migration and release at the ground surface. Gas pockets trapped in the subsurface have contrasting properties with the surrounding porous media that favor their detection using geophysical methods. We have developed a case study in which pockets of gas were intercepted with multilevel monitoring wells, and their lateral continuity was monitored over time using resistivity. We have developed a novel interpretation procedure based on Archie’s law to evaluate changes in water and gas content with respect to a mean background medium. We have used induced polarization data to account for errors in applying Archie’s law due to the contribution of surface conductivity effects. Mosaics defined by changes in water saturation allowed the recognition of gas migration and groundwater infiltration routes and the association of gas and groundwater fluxes. The inference on flux patterns was analyzed by taking into account pressure measurements in trapped gas reservoirs and by metagenomic analysis of the microbiological content, which was retrieved from suspended sediments in groundwater sampled in multilevel monitoring wells. A conceptual model combining physical and microbiological subsurface processes suggested that biogas trapped at depth may have the ability to quickly travel to the surface.

Binaural Reproduction of Finite Difference Simulations Using Spherical Array Processing

Due to its efficiency and simplicity, the finite-difference time-domain method is becoming a popular choice for solving wideband, transient problems in various fields of acoustics. So far, the issue of extracting a binaural response from finite difference simulations has only been discussed in the context of embedding a listener geometry in the grid. In this paper, we propose and study a method for binaural response rendering based on a spatial decomposition of the sound field. The finite difference grid is locally sampled using a volumetric array of receivers, from which a plane wave density function is computed and integrated with free-field head related transfer functions, in the spherical harmonics domain. The volumetric array is studied in terms of numerical robustness and spatial aliasing. Analytic formulas that predict the performance of the array are developed, facilitating spatial resolution analysis and numerical binaural response analysis for a number of finite difference schemes. Particular emphasis is placed on the effects of numerical dispersion on array processing and on the resulting binaural responses. Our method is compared to a binaural simulation based on the image method. Results indicate good spatial and temporal agreement between the two methods.

Impact of age, various forms of cataract, and visual acuity on whole-field scotopic sensitivity screening for glaucoma in rural Taiwan.

OBJECTIVE: To evaluate the impact of age, various forms of cataract, and visual acuity on whole-field scotopic sensitivity screening for glaucoma in a rural population. DESIGN: Clinic-based study with population-based recruitment. SETTING: Jin Shan Township near Taipei, Taiwan. SUBJECTS: Three hundred forty-six residents (ages, > or = 40 years) of Jin Shan Township. INTERVENTIONS: Whole-field scotopic testing, ophthalmoscopy with dilation of the pupils, cataract grading against photographic standards, and screening visual field testing in a random one-third subsample. MAIN OUTCOME MEASURES: Whole-field scotopic sensitivity (in decibels) and diagnostic status as a case of glaucoma, glaucoma suspect, or normal. RESULTS: Participants in Jin Shan Township did not differ significantly in the rate of blindness, low visual acuity, or family history of glaucoma from a random sample of nonrespondents. Scotopic sensitivity testing detected 100% (6/6) of subjects with open-angle glaucoma at a specificity of 80.2%. The mean +/- SE scotopic sensitivity for six subjects with open-angle glaucoma (32.78 +/- 1.51 dB) differed significantly from that of 315 normal individuals (38.51 +/- 0.22 dB), when adjusted for age and visual acuity (P = .05, t test). With linear regression modeling, factors that correlated significantly with scotopic sensitivity were intraocular pressure, screening visual field, best corrected visual acuity, presence of cortical cataract, and increasing age. CONCLUSIONS: Although cataract affects the whole-field scotopic threshold, it appears that scotopic testing may be of value in field-based screening for glaucoma.

Hα Line Profile Asymmetries and the Chromospheric Flare Velocity Field

The asymmetries observed in the line profiles of solar flares can provide important diagnostics of the properties and dynamics of the flaring atmosphere. In this paper the evolution of the Hα and Ca ii λ8542 lines are studied using high spatial, temporal, and spectral resolution ground-based observations of an M1.1 flare obtained with the Swedish 1 m Solar Telescope. The temporal evolution of the Hα line profiles from the flare kernel shows excess emission in the red wing (red asymmetry) before flare maximum and excess in the blue wing (blue asymmetry) after maximum. However, the Ca ii λ8542 line does not follow the same pattern, showing only a weak red asymmetry during the flare. RADYN simulations are used to synthesize spectral line profiles for the flaring atmosphere, and good agreement is found with the observations. We show that the red asymmetry observed in Hα is not necessarily associated with plasma downflows, and the blue asymmetry may not be related to plasma upflows. Indeed, we conclude that the steep velocity gradients in the flaring chromosphere modify the wavelength of the central reversal in the Hα line profile. The shift in the wavelength of maximum opacity to shorter and longer wavelengths generates the red and blue asymmetries, respectively.