One-dimensional Velocity Profiles Research Articles

Assessing The Potential Of PIV Data To Resolve Hidden Frequency Scales

The present work investigates the performance of an advection-based flow reconstruction model to increase the temporal resolution of PIV measurements recorded in a turbulent, planar jet. Using a semi-Lagrangian technique in combination with Rapid Distortion Theory, a modified trajectory tracking procedure is implemented. The method introduces a specification for a two-dimensional convective velocity based on the least squares minimization of the linearized advection equation, in contrast to the previously introduced one-dimensional mean velocity profile outlined in Vocke et al. (2023). The new implementation is based on local flow measurements, making it well-suited towards streamwise heterogeneity and spatially developing flows. With this, the flow at some unknown time can be estimated from the know flow measurements at the forward and backward time. Spectral analysis illustrates the model's proficiency in recovering spatiotemporal information far exceeding the Nyquist frequency, with spectral reconstruction errors of less than 5% for the most extreme case. It is demonstrated that the spectral content can be estimated at least two orders of magnitude beyond the sampling frequency of the original recording. Improved performance compared to alternative methods is demonstrated with only minor impact on computational time. This indicates the approach may be used as a tool for experimental researchers a) having no access to a high-speed PIV or b) with high-speed PIV to increase the spectral resolution even further.

Multiline molecular tagging velocimetry of nitric oxide at 100 kHz using an injection-seeded burst-mode OPO.

An injection-seeded, burst-mode optical parametric oscillator (OPO) operating at a repetition rate of 100 kHz is used to demonstrate the multiline molecular tagging velocimetry of an underexpanded jet using nitric oxide fluorescence. The very narrow linewidth of the OPO system, along with the relatively high pulse energies of the burst-mode system, enables efficient single-photon excitation of nitric oxide along multiple laser beam lines at a high repetition rate. Simultaneous one-dimensional velocity profile measurements were obtained of an underexpanded jet system at six different locations using a reference initial image and single-shot delayed images. A methodology for calculating the uncertainty of single-shot velocity is also described. Mean and root-mean-square velocity profiles are obtained at multiple locations simultaneously over a sampling time of 1 ms. The high-repetition-rate velocity measurements also appear to capture the onset of velocity oscillations and has the potential to reveal velocity frequency content occurring in the tens of kHz. The demonstrated velocimetry technique could be paired with other emerging burst-mode laser capabilities for a quantitative multiparameter gas property or multicomponent gas velocity measurements for supersonic and hypersonic flows, especially within ground test facilities that are limited to very short run durations.

NO measurements in high temperature hydrogen flames: The crucial role of the hydrogen oxidation chemistry for accurate NO predictions

The current work investigates the formation of Nitric Oxide (NO) in hydrogen–air flames, over a wide range of flame temperatures. The use of hydrogen allows improved focus on the thermal-NO pathway by removing the complexity introduced by the prompt-NO pathway, which has been shown to be an important contributor to inaccurate predictions of absolute post-flame NO concentrations in hydrocarbon flames. This experimental study is conducted at atmospheric pressure using stoichiometric, premixed, laminar stagnation flames. Adiabatic flame temperatures ranging from 1600K to 2300K are achieved by varying the argon concentration in air. One-dimensional velocity, temperature, and NO concentration profiles are measured using non-intrusive laser diagnostics: Particle Tracking Velocimetry (PTV), NO multiline thermometry, and NO Laser Induced Fluorescence (NO-LIF), respectively. Results show that the experimental velocity profiles are incorrectly captured by the studied mechanisms, especially at low and high temperatures. This suggests that major inaccuracies are present in the hydrogen oxidation chemistry of the thermochemical models, regardless of their optimisation methodology. Furthermore, NO-LIF profiles show major discrepancies between all the studied mechanisms and the experiments, especially at elevated temperatures. The disagreement stems from an inaccurate description of the base chemistry of the models. These inaccuracies arise specifically from the description of the radical pool driving the flame behaviour and NO formation. This study demonstrates the need for model optimisation on experimental measurements using pure hydrogen 1D flames to obtain an accurate description of the hydrogen oxidation chemistry at play. This would lead to an improved description of the NOx sub-chemistry of any hydrogen, or hydrocarbon, combustion system.

How well do we know thermal-NO?An investigation of NO formation in flames over a wide temperature range

This study investigates the large variability and uncertainty in the thermal-initiation reaction rates found in the literature. An experimental study is conducted at atmospheric pressure in lean, premixed, laminar methane-air stagnation flames. Flame temperatures ranging from 1900K to 2500K are achieved by varying the argon concentration, in 21% and 40% oxygen mixtures balanced with nitrogen. The conditions are selected to promote the thermal route in the overall post-flame NO formation. One-dimensional velocity, temperature, and NO concentration profiles are measured with Particle Tracking Velocimetry (PTV), NO multi-line thermometry, and NO Laser-Induced Fluorescence (NO-LIF), respectively. While velocity and temperature measurements are accurately predicted by different thermochemical models, the simulated NO-LIF signal profiles show significant disagreement and large variability, relative to the measurements. Results show that, across all conditions, none of the studied mechanisms are able to capture accurately both the NO concentration, and formation rate in the post-flame region, suggesting that the fundamental chemistry remains inaccurate. The discrepancy in the predictions appears to be linked to the chosen parameters of the Arrhenius rate, specifically the pre-exponential factor, and the activation energy. Sensitivity and Reaction Pathway Analyses suggest that the differences in the Arrhenius parameters could originate from different consideration of the base radical chemistry, as well as different relative contributions of the four NO-formation routes. As a result, some models can predict similar NO concentrations but using significantly different underlying base and NOx chemistry. This implies that the models could diverge significantly in conditions where other non-thermal routes are more important. This study demonstrates the need for spatially-resolved experimental data across a broad range of experimental conditions, promoting and considering a variety of routes, in order to optimize NOx chemistry models with reduced uncertainty.

Analytical steady-state model based on Fourier integral transforms for cylindrical heat pipes under axisymetric conditions

A thermohydraulic analytical model of a capillary cylindrical heat pipe in steady-state is proposed in this article. It is based on an original representation by thermal quadrupoles to describe heat transfer in the wall and in the porous wick, via the use of Fourier integral transforms. Thanks to a validation from literature results, this model provides two-dimensional axisymetric thermal fields and one-dimensional pressure and velocity profiles of both liquid and vapour flows. Different developments and solutions are introduced according to the kind of boundary conditions at evaporator and at condenser, and with a more or less strong thermohydraulic coupling at the liquid/vapour interface. For the simple case with imposed uniform heat fluxes, intrinsic properties of the heat pipe are originally defined. The introduced model offers a generalisation of analytical models of standard heat pipe as a design or optimisation tool. Wider developments of analytical models for more complex three-dimensional geometries of heat pipe and in transient regime can be expected.

Hypersonic N2 boundary layer flow velocity profile measurements using FLEET.

Femtosecond laser electronic excitation tagging (FLEET) velocimetry was used in the boundary layer of an ogive-cylinder model in a Mach-6 Ludwieg tube. One-dimensional velocity profiles were extracted from the FLEET signal in laminar boundary layers from pure N2 flows at unit Reynolds numbers ranging from 3.4×106/m to3.9×106/m. The effects of model tip bluntness and the unit Reynolds number on the velocity profiles were investigated. The challenges and strategies of applying FLEET for direct boundary layer velocity measurement are discussed. The potential of utilizing FLEET velocimetry for understanding the dynamics of laminar and turbulent boundary layers in hypersonic flows is demonstrated.

Site characterization in Ismailia, Egypt using seismic ambient vibration array

Seismic ambient vibration array measurements were carried out at seven sites in Ismailia, Egypt in order to obtain one-dimensional (1D) shear wave velocity (Vs) profiles and horizontal-to-vertical spectral ratios (HVSR). We applied the spatial autocorrelation (SPAC) method to the ambient vibration data for the estimating Rayleigh-wave dispersion curves and obtained 1D Vs profiles. The deeper parts (90–350 m) of the estimated models have Vs values ranging between 500 and 700 m/s, while the subsurface parts have Vs varying between 140 and 200 m/s. The averaged Vs in the upper 30 m based on the results indicate that site classes of all the sites are classified as C or D according to NEHRP soil classification. The computed site amplification factors derived from the estimated Vs profiles show that the resonant frequencies at sites on the east side of the Suez Canal are below 1 Hz, while those at sites on the west side vary between 1 and 2 Hz, indicating the difference of subsurface geology across the Suez Canal. The HVSR curves at each site are generally similar for all sensors, indicating nearly flat subsurface layers beneath each site. We also compared the dominant frequencies of HVSRs and those derived from 1D amplification factors. The comparisons show that both the frequencies are similar to each other, but there are slight differences possibly due to the lack of deeper sedimentary layer structures and lower resolution of HVSRs in the lower frequency range, indicating the importance of further ambient vibration survey with broadband sensors. The results show that the combination of SPAC and HVSR methods is advantageous to estimate proper Vs profiles to evaluate site effects for earthquake disaster mitigation.

Precision of VS30 values derived from noninvasive surface wave methods at 31 sites in California

We study the inter- and intra-method variability of VS30 results by inverting/forward-modeling individual dispersion data for 31 seismographic stations located in California where combinations of surface-wave methods were applied and the minimum recorded wavelength from each method satisfies the 30-m depth criteria. These methods consist of noninvasive geophysical (active and passive surface-wave techniques) multi-station approaches, including the Multi-channel Analysis of Surface Waves (MASW; Rayleigh and Love waves), Spectral Analysis of Surface Waves (SASW), Microtremor Array [using Extended Spatial Autocorrelation (ESAC) processing methods], and Refraction Microtremor (ReMi) methods. Depending on the apparent geologic or seismic complexity of the site, field crews applied one or a combination of these methods whenever economically feasible to estimate the one-dimensional shear-wave velocity (VS) profile and calculate VS30, the time-averaged VS to a depth of 30 m. For each of the 31 sites, we find both types of variability in VS30 estimates generally remain insignificant (arithmetic mean of 5% difference). We also find similar results (3%) when we evaluate individual-method based VS30 estimates against composite-method based estimates. We note that VS30 values vary insignificantly when using a combination of complementary methods, e.g., active MASW data combined with passive MAM data, and that the most reliable results are also based on close fitting of the theoretical dispersion data to the representative (experimental) dispersion data.

Ionized and Molecular Gas Kinematics in a z = 1.4 Star-forming Galaxy*

Abstract We present deep observations of a z = 1.4 massive, star-forming galaxy (SFG) in molecular and ionized gas at comparable spatial resolution (CO 3–2, NOrthern Extended Millimeter Array (NOEMA); Hα, Large Binocular Telescope (LBT)). The kinematic tracers agree well, indicating that both gas phases are subject to the same gravitational potential and physical processes affecting the gas dynamics. We combine the one-dimensional velocity and velocity dispersion profiles in CO and Hα to forward-model the galaxy in a Bayesian framework, combining a thick exponential disk, a bulge, and a dark matter halo. We determine the dynamical support due to baryons and dark matter, and find a dark matter fraction within one effective radius of f DM ( ≤ R e ) = 0.18 − 0.04 + 0.06 . Our result strengthens the evidence for strong baryon-dominance on galactic scales of massive z ∼ 1–3 SFGs recently found based on ionized gas kinematics alone.

Sporadic low-velocity volumes spatially correlate with shallow very low frequency earthquake clusters

A low-velocity zone (LVZ) has been detected by seismic exploration surveys within the Nankai accretionary prism toe off the Kii Peninsula, southwestern Japan, and is considered to be a mechanically weak volume at depth. Such mechanical heterogeneities potentially influence seismic and tsunamigenic slips on megathrust earthquakes in the subduction zone. However, the spatial distribution of the LVZ along the trough-parallel direction is still elusive. Here we show sporadic LVZs in the prism toe from one-dimensional shear wave velocity (Vs) profiles obtained at 49 cabled ocean bottom stations, which were estimated by a nonlinear inversion technique, simulated annealing, using the displacement–pressure ratios of the Rayleigh wave. The estimated distribution of LVZs along the trough widely correlates with the epicentres of shallow very low frequency earthquakes (sVLFEs), which suggests that sVLFEs are activated in the sporadically distributed low-velocity and mechanically weak volumes where fluids significantly reduce the shear strength of faults.

Microtremor exploration for shallow S-wave velocity profiles at stations in local strong motion network in Bursa, Yalova, and Kocaeli in north-western Turkey

We conducted microtremor array surveys for shallow S-wave velocity profiles at 20 sites in Bursa, Yalova and Kocaeli provinces in the north-western part of Turkey to provide fundamental data to assess the seismic hazard in the area. All of the measurement sites were positioned very close to strong motion stations belonging to the Disaster and Emergency Management Presidency of Turkey (AFAD) in order to further understand site amplification factors in strong motion records. Of the 20 study sites, two were located in Yalova, four in Bursa and 14 in Kocaeli.We temporarily installed two small arrays to obtain simultaneous records of vertical microtremors. Then, the spatial autocorrelation method was applied to retrieve Rayleigh wave phase velocity curves in a frequency range from 1 to 30 Hz from the array records. The phase velocities in the western part of the Kocaeli area are low across a wide frequency range, while relatively high phase velocities are found in the eastern part of the Kocaeli province. The phase velocities in the Yalova and Bursa provinces are widely distributed suggesting large variations in soil conditions. The observed phase velocity curve at each site was inverted to a one-dimensional (1D) S-wave velocity profile to a depth of 100 m, using a hybrid heuristic inversion method. All the S-wave velocity profiles in the eastern Kocaeli area are similar; however, the sites in the western Kocaeli and Yalova-Bursa areas have profiles with different features from the others.Finally, we discuss amplification factors for S-waves using the inverted profiles. The dominant fundamental periods of the amplification factors were distributed in a frequency range from 0.7 to 5 Hz. The profiles obtained are also used to map average S-wave velocities in the study area, with an addition of existing data at strong motion stations of the AFAD.

Shear Wave Velocity Profiles of Roadway Substructures from Multichannel Analysis of Surface Waves and Waveform Tomography

Assessment of roadway subsidence caused by embedded low-velocity anomalies is critical to the health and safety of the traveling public. Surface-based seismic techniques are often used to assess roadways because of data acquisition convenience and large depths of characterization. To mitigate the negative impact of closing a traffic lane under traditional seismic testing, a new test system that uses a land streamer is presented. The main advantages of the system are the elimination of the need to couple the geophones to the roadway, the use of only one source at the end of the geophone array, and the movement of the whole test system along the roadway quickly. For demonstration, experimental data were collected on asphalt pavement overlying a backfilled sinkhole that was experiencing further subsidence. For the study, a 24-channel land streamer and a propelled energy generator to generate seismic energy were used. The test system was pulled by a pickup truck along the roadway and the data were collected with 81 shots at every 3 m for a road segment of 277.5 m, with a total data acquisition time of about 1 h. The measured seismic data set was analyzed by the standard multichannel analysis of surface waves (MASW) and advanced two-dimensional (2-D) waveform tomography methods. Eighty-one one-dimensional shear wave velocity (VS) profiles from the MASW were combined to obtain a single 2-D profile. The waveform tomography method was able to characterize subsurface structures at a high resolution (1.5- × 1.5-m cells) along the test length to a depth of 22.5 m. Very low S-wave velocity was obtained at the repaired sinkhole location. The 2-D VS profiles from the MASW and waveform tomography methods are consistent. Both methods were able to delineate high- and low-velocity soil layers and variable bedrock.

Observations of the influence of Taylor-Couette geometry on the onset of shear-banding in surfactant wormlike micelles

One-dimensional velocity profiles have been acquired using nuclear magnetic resonance (NMR) techniques to observe the onset of shear-banding in the cetylpyridinium chloride (CPCl)/sodium salicylate (NaSal) in brine wormlike micelle system during steady shear rate controlled experiments. In particular, velocity profiles acquired with NMR were used to calculate the local shear rate next to the fixed wall of concentric cylinder (Taylor-Couette or TC) shear geometry. The change in this local shear rate, as a function of the applied shear rate, was then used to infer when the fluid began to shear-band. Using this technique, measurements were made on a series of TC geometries to explore the influence of the inherent shear rate inhomogeneity between the curved surfaces of the shear geometry. These experiments showed that the onset of shear-banding in the CPCl/NaSal/brine system was dependent on the TC geometry dimensions. Specifically, greater differences in curvature between the two shear surfaces caused shear-...

Transdimensional Love-wave tomography of the British Isles and shear-velocity structure of the East Irish Sea Basin from ambient-noise interferometry

We present the first Love-wave group velocity and shear velocity maps of the British Isles obtained from ambient noise interferometry and fully non-linear inversion. We computed interferometric inter-station Green's functions by cross-correlating the transverse component of ambient noise records retrieved by 61 seismic stations across the UK and Ireland. Group velocity measurements along each possible inter-station path were obtained using frequency-time analysis and converted into a series of inter-station traveltime datasets between 4 and 15 seconds period. Traveltime uncertainties estimated from the standard deviation of dispersion curves constructed by stacking randomly-selected subsets of daily cross-correlations, were observed to be too low to allow reasonable data fits to be obtained during tomography. Data uncertainties were therefore estimated again during the inversion as distance-dependent functionals. We produced Love-wave group velocity maps within 8 different period bands using a fully non-linear tomography method which combines the transdimensional reversible-jump Markov chain Monte Carlo (rj-McMC) algorithm with an eikonal raytracer. By modelling exact raypaths at each step of the Markov chain we ensured that the non-linear character of the inverse problem was fully and correctly accounted for. Between 4 and 10 seconds period, the group velocity maps show remarkable agreement with the known geology of the British Isles and correctly identify a number of low-velocity sedimentary basins and high-velocity features. Longer period maps, in which most sedimentary basins are not visible, are instead mainly representative of basement rocks. In a second stage of our study we used the results of tomography to produce a series of Love-wave group velocity dispersion curves across a grid of geographical points focussed around the East Irish Sea sedimentary basin. We then independently inverted each curve using a similar rj-McMC algorithm to obtain a series of one-dimensional shear velocity profiles. By merging all 1D profiles, we created a fully three-dimensional model of the crust beneath the East Irish Sea. The depth to basement in this model compares well with that averaged from seismic reflection profiles. This result is the first 3-dimensional model in the UK with fully quantified uncertainties: it shows basin depths and basement structures, and their concomitant uncertainties.

Mathematical modeling and simulation of flow velocity profile for rectangular open channels

The accuracy of discharge measurement in open channels using propeller current meters can be improved considerably by optimizing the number and locations of the current meters to be deployed for this purpose. The optimization is possible if the nature of the flow profile in the channel can be assessed. This paper reports development of a mathematical model of the flow velocity profile of rectangular open channels. One-dimensional velocity profile equation is taken as a sum of logarithmic and parabolic terms and used to determine the unit-width discharge profile across the width of the channel. The unit-width discharge equation comprises two logarithmic terms along with a beta function. The two equations together involve 10 mathematical parameters, which are linked to the physical conditions, geometry, and configuration of the channel. MATLAB code has been developed to account for these parameters and obtain the velocity profile of flow in any rectangular open channel. Flow under similar conditions is simulated using computational fluid dynamics (CFD) and the results for six different cases are used for validating the mathematical model. The CFD based simulation, in turn, is validated by carrying out actual measurements in a rectangular open channel using an acoustic Doppler current profiler.

One-dimensional velocity profiles in open-channel flow with intense transport of coarse sediment

The paper deals with laboratory experiments in open-channel flows with intense transport of model sediment (coarse plastic particles) in our new tilting flume. The major objectives of the paper are: 1. to discuss applied measuring methods, 2. to analyze measured velocity profiles. Ad 1. A profile of the longitudinal component of local velocity was measured across the vertical axis of symmetry of a flume cross section using three independent measuring methods (Prandtl tube, Ultrasonic Velocity Profiler, Acoustic Doppler Velocity Profiler). Due to strong stratification of the flow in the flume, parts of the profile are measured in regions of very different local concentrations of sediment (from virtually zero concentration to the maximum concentration of bed packing). This makes measurements complicated, particularly for ultrasonic measuring techniques. Profiles measured using the different techniques are evaluated and mutually compared. Ad 2. The layered character of the flow causes that shapes of velocity profiles tend to be different in the transport layer (rich on transported particles) above the bed and in the solids-free region between the top of the transport layer and the water surface. Shapes of the profiles are analyzed. Particular attention is paid to the logarithmic profile in the solids-free region of the flow cross section. The profile can be handled using the law of the hydraulically-rough wall. In the law, the eroded top of the bed with the transport layer is supposed to be the rough boundary and appropriate values are sought for its variables.

ICONE23-1907 EFFECT OF COLD CAP BOUNDARY CONDITIONS ON JOULE-HEATING FLOW IN THE SLOPING BOTTOM CAVITY

Flow behavior in a sloping bottom cavity is observed to study the effect of cavity shape on flow behavior for Joule-heating flow. In the former study, a simple cubic cavity is applied to study the chaotic flow behavior of Joule-heating convection due to simplification as the real melter case is complicated. In this study, a sloping bottom cavity of the dimension one-fifth of the actual melter is applied to study the detail flow behavior. Carbon electrodes and top cooling surface are placed to make Joule-heating and the chaotic flow behavior. The working fluid is 80%wt Glycerol-water solution with LiCl as electrolyte. To observe the chaotic flow behavior spatio-temporally, Ultrasonic Velocity Profiler (UVP) is applied in this experiment to obtain the one-dimensional continuous velocity profiles in the center line of cavity. Particle Image Velocity (PIV) method is also applied to observe the two-dimensional flow behavior and to examine the cross-check between UVP and PIV for the chaotic flow behavior with temperature distribution. The flow profiles of the former cubic cavity and the sloping bottom cavity are compared changing voltage magnitude and cooling temperature of the electrodes side to analyze the effect of cavity shape under Joule-heating condition. The flow behavior in the upper part of the sloping bottom cavity is similar to that in the cubic cavity in the experiment in whole cavity, the range down-flow achieved is larger than the cubic cavity.

Aircraft-Based Aerosol Sampling in Clouds: Performance Characterization of Flow-Restriction Aerosol Inlets

AbstractInteraction of liquid cloud droplets and ice particles with aircraft aerosol inlets can result in the generation of a large number of secondary particles and contaminate aerosol measurements. Recent studies have shown that a sampler designed with a perpendicular subsampling tube located within a flow-through conduit (i.e., a flow-restriction inlet) was best suited for in-cloud sampling. Analysis of field data obtained from different flow-restriction inlets shows that their critical cloud droplet breakup diameters are strongly dependent on design details and operating conditions. Using computational fluid dynamics (CFD) simulations, in-cloud sampling performance of a selected inlet can be predicted reasonably accurately for known operating conditions. To understand the relation between inlet design parameters and its sampling performance, however, CFD calculations are impractical. Here, using a simple, representative one-dimensional velocity profile and a validated empirical droplet breakup criteria, a parametric study is conducted to understand the relationship between different inlet design features and operating conditions on its critical breakup diameters. The results of this study suggest that an optimal inlet for in-cloud aerosol sampling should have a combination of a restriction nozzle at the aft end of the flow-through conduit to minimize wall-impaction shatter artifacts and a blunt leading edge to minimize shatter artifact generation from the aerodynamic breakup of cloud droplets. Inlets for in-cloud aerosol sampling from aircraft will, therefore, differ significantly in design from those used for clear-air aerosol sampling.

General pure convection residence time distribution theory of fully developed laminar flows in straight planar and axisymmetric channels

In literature, the diffusion-free residence time distribution (RTD) of laminar flows – the so-called convection model – has been determined for various velocity profiles mostly on a case-by-case basis. In this analytical paper, we derive general mathematical relations which allow computing the diffusion-free differential and cumulative RTD in straight planar, circular and concentric annular channels for arbitrary monotonic and piece-wise monotonic one-dimensional velocity profiles. The theory is used to determine the RTD of plane Couette–Poiseuille flow with non-monotonic velocity profile, and the optimal value of the volumetric flow rate where the RTD becomes most narrow. It is shown that any velocity profile that depends in a sub-layer linearly on the distance from a stationary or moving no-slip wall has a differential RTD which follows a −3 power law as the residence time approaches its maximum. The variance of the RTD is directly associated with the asymptotic behavior of the RTD and can be finite or infinite.

Effect of surface wave inversion non-uniqueness on 1D seismic ground response analysis

Surface wave methods are increasingly being used for geotechnical site characterization. The methodology is based on the dispersive characteristic of Rayleigh waves in vertically heterogeneous medium. Experimental dispersion curve is inverted to obtain one-dimensional shear-wave velocity profile by inverse problem solution. Uncertainty associated with this surface wave inversion has drawn much attention. Inverse problem solution can provide different equivalent shear-wave velocity profiles, which may lead to different seismic site response analysis. In this study, a neighborhood algorithm has been used for inversion of dispersion curve to get a set of equivalent shear-wave velocity profiles. These equivalent velocity profiles are then used for 1D ground response analysis for different input motion record of the same earthquake at different epicentral distances. Results show significant variation in amplification spectrum in terms of maximum amplification as well as peak frequency. The extent of this uncertainty largely depends on the characteristics of the ground motion records at different epicentral distances. A linear variation is observed between mean coefficients of variation of amplification spectrum and epicentral distance of ground motion records. A gradual increase in mean value of peak frequency and peak amplification with the epicentral distance is also observed.