Modeling of subsurface pulsed radar for nondestructive testing of structures

The use of pulsed radar for investigating the integrity of structural elements is gaining popularity and becoming firmly established as a nondestructive test method in civil engineering. Difficulties can often arise in the interpretation of results obtained, particularly where internal details are relatively complex. One approach that can be used to understand and evaluate radar results is through numerical modeling of signal propagation and reflection. By comparing the results of a numerical modeling with those from field measurements, engineers can gain valuable insight into the probable features embedded beneath the surface of a structural element. This paper discusses a series of numerical techniques for modeling subsurface radar and compares the precision of the results with those taken from real field data. It is found that more complex problems require more sophisticated analysis techniques to obtain realistic results, with a consequential increase in the computational resources to carry out the modeling.

Numerical assessment of potential impacts of hydraulically fractured Bowland Shale on overlying aquifers

Natural gas extracted from hydraulically fractured shale formations potentially has a big impact on the global energy landscape. However, there are concerns of potential environmental impacts of hydraulic fracturing of the shale formations, particularly those related to water quality. To evaluate the potential impact of hydraulically fractured shale on overlying aquifers, we conduct realizations of numerical modeling simulations to assess fluid flow and chloride transport from a synthetic Bowland Shale over a period of 11,000 years. The synthetic fractured shale was represented by a three-dimensional discrete fracture model that was developed by using the data from a Bowland Shale gas exploration in Lancashire, UK. Chloride mass exchange between fractures and the rock matrix was fully accounted for in the model. The assessment was carried out to investigate fluid and chloride mass fluxes before, during, and after hydraulic fracturing of the Bowland Shale. Impacts of the upward fracture height and aperture, as well as hydraulic conductivity of the multilayered bedrock system, are also included this assessment. This modeling revealed that the hydraulically fractured Bowland Shale is unlikely to pose a risk to its overlying groundwater quality when the induced fracture aperture is ≤200 µm. With the fracture aperture ≥1000 µm, the upward chloride flux becomes very sensitive to the upward fracture height growth and hydraulic conductivity of the multilayered bedrock system. In the extremely unlikely event of the upward fracture growth directly connecting the shale formation to the overlying Sherwood Sandstone aquifer with the fracture aperture ≥1000 µm, the upward chloride mass flux could potentially pose risks to the overlying aquifer in 100 years. The model study also revealed that the upward mass flux is significantly intercepted by the horizontal mass flux within a high permeable layer between the Bowland Shale and its overlying aquifers, reducing further upward flux toward the overlying aquifers.

Meanline Modeling of Inlet Recirculation in Automotive Turbocharger Centrifugal Compressors

This study provides a novel meanline modeling approach for centrifugal compressors. All compressors analyzed are of the automotive turbocharger variety and have typical upstream geometry with no casing treatments or preswirl vanes. Past experience dictates that inducer recirculation is prevalent toward surge in designs with high inlet shroud to outlet radius ratios; such designs are found in turbocharger compressors due to the demand for operating range. The aim of the paper is to provide further understanding of impeller inducer flow paths when operating with significant inducer recirculation. Using three-dimensional (3D) computational fluid dynamics (CFD) and a single-passage model, the flow coefficient at which the recirculating flow begins to develop and the rate at which it grows are used to assess and correlate work and angular momentum delivered to the incoming flow. All numerical modeling has been fully validated using measurements taken from hot gas stand tests for all compressor stages. The new modeling approach links the inlet recirculating flow and the pressure ratio characteristic of the compressor. Typically for a fixed rotational speed, between choke and the onset of impeller inlet recirculation the pressure ratio rises gradually at a rate dominated by the aerodynamic losses. However, in modern automotive turbocharger compressors where operating range is paramount, the pressure ratio no longer changes significantly between the onset of recirculation and surge. Instead the pressure ratio remains relatively constant for reducing mass flow rates until surge occurs. Existing meanline modeling techniques predict that the pressure ratio continues to gradually rise toward surge, which when compared to test data is not accurate. A new meanline method is presented here which tackles this issue by modeling the direct effects of the recirculation. The result is a meanline model that better represents the actual fluid flow seen in the CFD results and more accurately predicts the pressure ratio and efficiency characteristics in the region of the compressor map affected by inlet recirculation.

Quasi-phasematching of harmonic generation via multimode beating in waveguides

A new scheme for quasi-phasematching high harmonic generation (HHG) in gases is proposed. In this, the rapid variation of the axial intensity resulting from excitation of more than one mode of a waveguide is used to achieve quasi phasematching. Numerical modeling demonstrates enhancement of the harmonic signal over that achieved for a single coherence length by factors > 10(4). (C) 2007 Optical Society of America

Energy conserving schemes for the simulation of musical instrument contact dynamics

Collisions are an innate part of the function of many musical instruments. Due to the nonlinear nature of contact forces, special care has to be taken in the construction of numerical schemes for simulation and sound synthesis. Finite difference schemes and other time-stepping algorithms used for musical instrument modelling purposes are normally arrived at by discretising a Newtonian description of the system. However because impact forces are non-analytic functions of the phase space variables, algorithm stability can rarely be established this way. This paper presents a systematic approach to deriving energy conserving schemes for frictionless impact modelling. The proposed numerical formulations follow from discretising Hamilton׳s equations of motion, generally leading to an implicit system of nonlinear equations that can be solved with Newton׳s method. The approach is first outlined for point mass collisions and then extended to distributed settings, such as vibrating strings and beams colliding with rigid obstacles. Stability and other relevant properties of the proposed approach are discussed and further demonstrated with simulation examples. The methodology is exemplified through a case study on tanpura string vibration, with the results confirming the main findings of previous studies on the role of the bridge in sound generation with this type of string instrument.

The radiated energy budget of Chromospheric Plasma in a major solar flare deduced from multi-wavelength observations

This paper presents measurements of the energy radiated by the lower solar atmosphere, at optical, UV, and EUV wavelengths, during an X-class solar flare (SOL2011-02-15T01:56) in response to an injection of energy assumed to be in the form of nonthermal electrons. Hard X-ray observations from RHESSI were used to track the evolution of the parameters of the nonthermal electron distribution to reveal the total power contained in flare accelerated electrons. By integrating over the duration of the impulsive phase, the total energy contained in the nonthermal electrons was found to be >2 × 1031 erg. The response of the lower solar atmosphere was measured in the free–bound EUV continua of H i (Lyman), He i, and He ii, plus the emission lines of He ii at 304 Å and H i (Lyα) at 1216 Å by SDO/EVE, the UV continua at 1600 Å and 1700 Å by SDO/AIA, and the white light continuum at 4504 Å, 5550 Å, and 6684 Å, along with the Ca ii H line at 3968 Å using Hinode/SOT. The summed energy detected by these instruments amounted to ~3 × 1030 erg; about 15% of the total nonthermal energy. The Lyα line was found to dominate the measured radiative losses. Parameters of both the driving electron distribution and the resulting chromospheric response are presented in detail to encourage the numerical modeling of flare heating for this event, to determine the depth of the solar atmosphere at which these line and continuum processes originate, and the mechanism(s) responsible for their generation.

Physical and numerical constraints in source modeling for finite difference simulation of room acoustics

In finite difference time domain simulation of room acoustics, source functions are subject to various constraints. These depend on the way sources are injected into the grid and on the chosen parameters of the numerical scheme being used. This paper addresses the issue of selecting and designing sources for finite difference simulation, by first reviewing associated aims and constraints, and evaluating existing source models against these criteria. The process of exciting a model is generalized by introducing a system of three cascaded filters, respectively, characterizing the driving pulse, the source mechanics, and the injection of the resulting source function into the grid. It is shown that hard, soft, and transparent sources can be seen as special cases within this unified approach. Starting from the mechanics of a small pulsating sphere, a parametric source model is formulated by specifying suitable filters. This physically constrained source model is numerically consistent, does not scatter incoming waves, and is free from zero- and low-frequency artifacts. Simulation results are employed for comparison with existing source formulations in terms of meeting the spectral and temporal requirements on the outward propagating wave.

Diagnostic based modeling for determining absolute atomic oxygen densities in atmospheric pressure helium-oxygen plasmas

Absolute atomic oxygen ground state densities in a radio-frequency driven atmospheric pressure plasma jet, operated in a helium-oxygen mixture, are determined using diagnostic based modeling. One-dimensional numerical simulations of the electron dynamics are combined with time integrated optical emission spectroscopy. The population dynamics of the upper O 3p 3P (l=844 nm) atomic oxygen state is governed by direct electron impact excitation, dissociative excitation, radiation losses, and collisional induced quenching. Absolute values for atomic oxygen densities are obtained through comparison with the upper Ar 2p1 (l=750.4 nm) state. Results for spatial profiles and power variations are presented and show excellent quantitative agreement with independent two-photon laser-induced fluorescence measurements.

Diagnostic-based modeling on a micro-scale atmospheric-pressure plasma jet

Diagnostic-based modeling (DBM) actively combines complementary advantages of numerical plasma simulations and relatively simple optical emission spectroscopy (OES). DBM is applied to determine spatial absolute atomic oxygen ground-state density profiles in a micro atmospheric-pressure plasma jet operated in He–O2. A 1D fluid model with semi-kinetic treatment of the electrons yields detailed information on the electron dynamics and the corresponding spatio-temporal electron energy distribution function. Benchmarking this time- and space-resolved simulation with phase-resolved OES (PROES) allows subsequent derivation of effective excitation rates as the basis for DBM. The population dynamics of the upper O(3p3P) oxygen state (? = 844 nm) is governed by direct electron impact excitation, dissociative excitation, radiation losses, and collisional induced quenching. Absolute values for atomic oxygen densities are obtained through tracer comparison with the upper Ar(2p1) state (? = 750.4 nm). The resulting spatial profile for the absolute atomic oxygen density shows an excellent quantitative agreement to a density profile obtained by two-photon absorption laser-induced fluorescence spectroscopy.

Modeling Frequency-Dependent Boundaries as Digital Impedance Filters in FDTD and K-DWM Room Acoustics Simulations

A new method for modeling-frequency-dependent boundaries in finite-difference time-domain (FDTD) and Kirchhoff variable digital waveguide mesh (K-DWM) room acoustics simulations is presented. The proposed approach allows the direct incorporation of a digital impedance filter (DIF) in the Multidimensional (2D or 3D) FDTD boundary model of a locally reacting surface. An explicit boundary update equation is obtained by carefully constructing a Suitable recursive formulation. The method is analyzed in terms of pressure wave reflectance for different wall impedance filters and angles of incidence. Results obtained from numerical experiments confirm the high accuracy of the proposed digital impedance filter boundary model, the reflectance of which matches locally reacting surface (LRS) theory closely. Furthermore a numerical boundary analysis (NBA) formula is provided as a technique for an analytic evaluation of the numerical reflectance of the proposed digital impedance filter boundary formulation.

A Phase Grating Approach to Modeling Surface Diffusion in FDTD Room Acoustics Simulations

In this paper, a method for modeling diffusive boundaries in finite difference time domain (FDTD) room acoustics simulations with the use of impedance filters is presented. The proposed technique is based on the concept of phase grating diffusers, and realized by designing boundary impedance filters from normal-incidence reflection filters with added delay. These added delays, that correspond to the diffuser well depths, are varied across the boundary surface, and implemented using Thiran allpass filters. The proposed method for simulating sound scattering is suitable for modeling high frequency diffusion caused by small variations in surface roughness and, more generally, diffusers characterized by narrow wells with infinitely thin separators. This concept is also applicable to other wave-based modeling techniques. The approach is validated by comparing numerical results for Schroeder diffusers to measured data. In addition, it is proposed that irregular surfaces are modeled by shaping them with Brownian noise, giving good control over the sound scattering properties of the simulated boundary through two parameters, namely the spectral density exponent and the maximum well depth.

Finite difference time domain modeling of phase grating diffusion

In this paper, a method for modeling diffusion caused by non-smooth boundary surfaces in simulations of room acoustics using finite difference time domain (FDTD) technique is investigated. The proposed approach adopts the well-known theory of phase grating diffusers to efficiently model sound scattering from rough surfaces. The variation of diffuser well-depths is attained by nesting allpass filters within the reflection filters from which the digital impedance filters used in the boundary implementation are obtained. The presented technique is appropriate for modeling diffusion at high frequencies caused by small surface roughness and generally diffusers that have narrow wells and infinitely thin separators. The diffusion coefficient was measured with numerical experiments for a range of fractional Brownian diffusers.

Hydraulic Subsurface Measurements and Hydrodynamic Modeling as Indicators for Groundwater Flow Systems in the Rotondo Granite, Central Alps (Switzerland)

Regional groundwater flow in high mountainous terrain is governed by a multitude of factors such as geology, topography, recharge conditions, structural elements such as fracturation and regional fault zones as well as man-made underground structures. By means of a numerical groundwater flow model, we consider the impact of deep underground tunnels and of an idealized major fault zone on the groundwater flow systems within the fractured Rotondo granite. The position of the free groundwater table as response to the above subsurface structures and, in particular, with regard to the influence of spatial distributed groundwater recharge rates is addressed. The model results show significant unsaturated zones below the mountain ridges in the study area with a thickness of up to several hundred metres. The subsurface galleries are shown to have a strong effect on the head distribution in the model domain, causing locally a reversal of natural head gradients. With respect to the position of the catchment areas to the tunnel and the corresponding type of recharge source for the tunnel inflows (i.e. glaciers or recent precipitation), as well as water table elevation, the influence of spatial distributed recharge rates is compared to uniform recharge rates. Water table elevations below the well exposed high-relief mountain ridges are observed to be more sensitive to changes in groundwater recharge rates and permeability than below ridges with less topographic relief. In the conceptual framework of the numerical simulations, the model fault zone has less influence on the groundwater table position, but more importantly acts as fast flow path for recharge from glaciated areas towards the subsurface galleries. This is in agreement with a previous study, where the imprint of glacial recharge was observed in the environmental isotope composition of groundwater sampled in the subsurface galleries. Copyright © 2012 John Wiley & Sons, Ltd.