Effects of instabilities and coherent structures on the performance of a thermocline based thermal energy storage

We investigate the transient dynamics of disturbances inside a thermocline based molten salt thermal energy storage (TES). Numerical simulations were conducted with four inlet flow configurations. The disturbances introduced at the inlet grow via Rayleigh Taylor instability. The formed vortical motions inside the tank propagate downstream and destroy the thermocline. The vortex-thermocline interaction upsets the stratification inside the TES. The disturbance growth rate, penetration length and vortex Reynolds number are measured. The growth of penetration length prior to the vortex-thermocline interaction is quadratic. The vortex Reynolds number of the eddy which causes thermocline breakdown increases with increase in Atwood number. The impingement of vortex on thermocline is studied. (C) 2015 Elsevier Ltd. All rights reserved.

Die Filling Behaviour of Semi Solid A356 Al Alloy Slurry During Rheo Pressure Die Casting

The present work discusses the findings obtained from simulations of semi solid die filling of a steering knuckle, prior to actual component development using in-house developed rheo pressure die casting system. Die filling capability of A356 Al alloy at semi-solid state has been investigated using commercial software Flow-3Dcast to optimise the pouring temperature of semi-solid slurry into the die cavity, while all other variables such as gating design, die preheat temperature and injection velocity are kept constant based on the prior knowledge obtained from trial numerical simulations and experimentation. Efforts have been made to nullify the essence of costly, time consuming experiments towards obtaining high-quality castings out of the findings obtained from numerical simulations. The optimum pouring temperature identified in the present study is 610 A degrees C, which facilitates smoother slurry flow, minimum surface defect concentration, uniform temperature field and solid fraction distribution within the component cavity.

Rayleigh-Taylor instability in a thermocline based thermal storage tank

Numerical simulations are performed to study the stability characteristics of a molten salt thermocline storage unit. Perturbations are introduced into a stable flow field in such a way as to make the top-fluid heavier than the fluid at the bottom, thereby causing a possible instability in the system. The evolution pattern of the various disturbances are examined in detail. Disturbances applied for short duration get decayed before they could reach the thermocline, whereas medium and long duration disturbances evolve into a ``falling spike'' or ``stalactite-like'' structure and destabilize the thermocline. Rayleigh Taylor instability is observed inside the storage tank. The effect of the duration, velocity and temperature of the disturbance on thermocline thickness and penetration length are studied. A quadratic time dependence of penetration length was observed. New perspectives on thermocline breakdown phenomena are obtained from the numerical flow field. (C) 2015 Elsevier Masson SAS. All rights reserved.

DESIGN REQUIREMENTS FOR DIRECT SUPERCRITICAL CARBON DIOXIDE RECEIVER DEVELOPMENT AND TESTING

This paper establishes the design requirements for the development and testing of direct supercritical carbon dioxide (sCO2) solar receivers. Current design considerations are based on the ASME Boiler and Pressure Vessel Code (BPVC). Section I (BPVC) considers typical boilers/superheaters (i.e. fired pressure vessels) which work under a constant low heat flux. Section VIII (BPVC) considers pressure vessels with operating pressures above 15 psig 2 bar] (i.e. unfired pressure vessels). Section III, Division I - Subsection NH (BPVC) considers a more detailed stress calculation, compared to Section I and Section VIII, and requires a creep-fatigue analysis. The main drawback from using the BPVC exclusively is the large safety requirements developed for nuclear power applications. As a result, a new set of requirements is needed to perform detailed thermal-structural analyses of solar thermal receivers subjected to a spatially-varying, high-intensity heat flux. The last design requirements document of this kind was an interim Sandia report developed in 1979 (SAND79-8183), but it only addresses some of the technical challenges in early-stage steam and molten-salt solar receivers but not the use of sCO2 receivers. This paper presents a combination of the ASME BPVC and ASME B31.1 Code modified appropriately to achieve the reliability requirements in sCO(2) solar power systems. There are five main categories in this requirements document: Operation and Safety, Materials and Manufacturing, Instrumentation, Maintenance and Environmental, and General requirements. This paper also includes the modeling guidelines and input parameters required in computational fluid dynamics and structural analyses utilizing ANSYS Fluent, ANSYS Mechanical, and nCode Design Life. The main purpose of this document is to serve as a reference and guideline for design and testing requirements, as well as to address the technical challenges and provide initial parameters for the computational models that will be employed for the development of sCO(2) receivers.

连续移动三维瞬态激光熔池温度场数值模拟

详细介绍了在ANSYS软件平台上，建立连续移动三维瞬态激光熔池温度场计算模型的方法，计算模型中考虑了材料表面温度对激光吸收率的影响及材料相变过程对激光熔池温度场的影响。系统分析了连续移动三维激光熔池温度场随时间的变化规律。通过该计算模型，可以掌握激光加工过程中连续移动激光熔池的加热和冷却规律。计算结果表明，当激光沿45#钢基板表面由一端向另一端沿直线扫描时，由于热传导的作用，激光熔池温度随时间增加而升高，同时连续移动熔池表面温度最高点不在激光束中心，而是稍稍偏后于激光束中心。在相同激光工艺参数下，计算熔池横截面尺寸与实验所测熔池横截面尺寸相吻合，表明所建立的连续移动熔池温度场计算模型是正确和可靠的。

Feasibility of Laminar Plasma-Jet Hardening of Cast Iron Surface

Nontransferred DC laminar plasma jets of stable flow and low impinging pressure acting on the substrate were used to heat W–Mo–Cu cast iron for phase transfer hardening of the surface layer. Substrates were heated in multipass with or without overlapping or heated with only single-pass. Surface morphologies of the molten trace and microstructure of the cross-section were observed, and the hardness distribution of the treated surface layer was examined. The surface layer of single-pass-heated specimen has an average hardness of about 900 HV0.1, while the specimen treated with multipass shows an average hardness of about 700 HV0.1, because of the heat effect from the neighboring pass treating, compared with the substrate hardness of about 300 HV0.1. The results demonstrate the stable and favorably controlled heating of the laminar plasma jet on the substrate surface and feasibility of using it as a tool for surface hardening of cast iron.

Three-Dimensional Modeling of Heat Transfer and Fluid Flow in Laminar-Plasma Material Re-Melting Processing

Modeling study is performed concerning the heat transfer and fluid flow for a laminar argon plasma jet impinging normally upon a flat workpiece exposed to the ambient air. The diffusion of the air into the plasma jet is handled by using the combined-diffusion-coefficient approach. The heat flux density and jet shear stress distributions at the workpiece surface obtained from the plasma jet modeling are then used to study the re-melting process of a carbon steel workpiece. Besides the heat conduction within the workpiece, the effects of the plasma-jet inlet parameters (temperature and velocity), workpiece moving speed, Marangoni convection, natural convection etc. on the re-melting process are considered. The modeling results demonstrate that the shapes and sizes of the molten pool in the workpiece are influenced appreciably by the plasma-jet inlet parameters, workpiece moving speed and Marangoni convection. The jet shear stress manifests its effect at higher plasma-jet inlet velocities, while the natural convection effect can be ignored. The modeling results of the molten pool sizes agree reasonably with available experimental data.

Microstructure and mechanical properties at TiCp/Ni-alloy interfaces in laser-synthesized coatings

Titanium carbide particle (TiCp) reinforced Ni alloy composite coatings were synthesized by laser cladding using a cw 3 kW CO2 laser. Two kinds of coatings were present in terms of TiCp origins, i.e. undissolved and in situ reacted TiCp, respectively. The former came from the TiCp pre-coated on the sample, whereas the latter from in situ reaction between titanium and graphite in the molten pool during laser irradiation. Conventional and high-resolution transmission electron microscope observations showed the epitaxial growth of TiC, the precipitation of CrB, and the chemical reaction between Ti and B elements around phase interfaces of undissolved TiCp. The hardness, H, and elastic modulus, E, were measured by nanoindentation of the matrix near the TiCp interface. For undissolved TiCp, the loading curve revealed pop-in phenomena caused by the plastic deformation of the crack formation or debounding of TiCp from the matrix. As for in situ generated TiCp, no pop-in mark appears. On the other hand, in situ reacted TiCp led to much higher hardness and modulus than that in the case of undissolved TiCp. The coating reinforced by in situ generated TiCp displayed the highest impact wear resistance at both low and high impact conditions, as compared with coatings with undissolved TiCp and without TiCp. The impact wear resistance of the coating reinforced by undissolved TiCp increases at a low impact work but decreases at a high impact work, as compared with the single Ni alloy coating. The degree of wear for the composite coating depends primarily on the debonding removal of TiCp.

Plasma-Arc Technology For The Thermal Treatment Of Chemical Wastes

Plasma-arc technology was developed to dispose of chemical wastes from a chemical plant by the Institute of Mechanics, Chinese Academy of Sciences (CAS-IMECH). A pilot plant system with this technology was constructed to destroy two types of chemical wastes. The system included shredding, mixing, and feeding subsystems, a plasma-arc reactor of 150 kW, an off-gas burning subsystem, and a scrubbing subsystem. The additives (CaO, SiO2, and Fe) were added into the reactor to form vitrified slag and capture the hazardous elements. The molten slag was quickly quenched to form an amorphous glassy structure. A direct current (DC) experimental facility of 30kW with plasma-arc technology was also set up to study the pyrolysis process in the laboratory, and the experimental results showed the cooling speed is the most important factor for good vitrified structure of the slag. According to previous tests, the destruction and removal efficiency (DRE) for these chemical wastes was more than 99.999%, and the polychlorinated biphenyls (PCBs) concentration in the solid residues was in the range of 1.28 to 12.9mg/kg, which is far below the Chinese national emission limit for the hazardous wastes. A simplified electromagneto model for numerical simulation was developed to predict the temperature and velocity fields. This model can make satisfactory maximum temperature and velocity distributions in the arc region, as well as the results by the magneto hydrodynamic approach.

THE FORMATION MECHANISM OF THE VITRIFIED SLAG IN A PLASMA ARC REACTOR FOR HAZARDOUS WASTE TREATMENT

Various hazardous wastes with additives have been vitrified to investigate the formation mechanism of the glassy slag by a 30 kW DC plasma-arc reactor developed by the Institute of Mechanics, Chinese Academy of Sciences. The average temperature in the reaction area is controlled at 1500°C. The chemical compositions of three sorts of fly ashes are analyzed by XRF (X-Ray Fluorescence). Fly ashes with vitrifying additives can be vitrified to form glassy slag, which show that the ratio of the whole oxygen ions to the whole network former ions in glass (R) is appropriate in the range of 2~3 to form durable vitrified slag. In this experiment, the arc power is controlled below 5 kW to inhibit waste evaporation. To enhance the effects of heat transfer to wastes, ferrous powder has been added into the graphite crucible, which aggregates as ingot below the molten silicate after vitrification. The slag fails to form glass if the quenching rate is less than 1 K/min. Therefore, the slag will break into small chips due to the sharp quenching rate, which is more than 100 K/sec.

Liquid silicate equation of state : using shock waves to understand the properties of the deep Earth

The equations of state (EOS) of several geologically important silicate liquids have been constrained via preheated shock wave techniques. Results on molten Fe₂SiO₄ (fayalite), Mg₂SiO₄ (forsterite), CaFeSi₂O₆ (hedenbergite), an equimolar mixture of CaAl₂Si₂O₈-CaFeSi₂O₆ (anorthite-hedenbergite), and an equimolar mixture of CaAl₂Si₂O₈-CaFeSi₂O₆-CaMgSi₂O₆(anorthite-hedenbergite-diopside) are presented. This work represents the first ever direct EOS measurements of an iron-bearing liquid or of a forsterite liquid at pressures relevant to the deep Earth (> 135 GPa). Additionally, revised EOS for molten CaMgSi₂O₆ (diopside), CaAl₂Si₂O₈ (anorthite), and MgSiO₃ (enstatite), which were previously determined by shock wave methods, are also presented.

The liquid EOS are incorporated into a model, which employs linear mixing of volumes to determine the density of compositionally intermediate liquids in the CaO-MgO-Al₂O₃-SiO₂-FeO major element space. Liquid volumes are calculated for temperature and pressure conditions that are currently present at the core-mantle boundary or that may have occurred during differentiation of a fully molten mantle magma ocean.

The most significant implications of our results include: (1) a magma ocean of either chondrite or peridotite composition is less dense than its first crystallizing solid, which is not conducive to the formation of a basal mantle magma ocean, (2) the ambient mantle cannot produce a partial melt and an equilibrium residue sufficiently dense to form an ultralow velocity zone mush, and (3) due to the compositional dependence of Fe2+ coordination, there is a threshold of Fe concentration (molar X_Fe ≤ 0.06) permitted in a liquid for which its density can still be approximated by linear mixing of end-member volumes.

Oxygen, carbon and hydrogen isotope studies of contact metamorphism

The O¹⁸/O¹⁶, C¹³/C¹², and D/H ratios have been determined for rocks and coexisting minerals from several granitic plutons and their contact metamorphic aureoles in northern Nevada, eastern California, central Colorado, and Texas, with emphasis on oxygen isotopes. A consistent order of O¹⁸/O¹⁶, C¹³/C¹², and D/H enrichment in coexisting minerals, and a correlation between isotopic fractionations among coexisting mineral pairs are in general observed, suggesting that mineral assemblages tend to approach isotopic equilibrium during contact metamorphism. In certain cases, a correlation is observed between oxygen isotopic fractionations of a mineral pair and sample distance from intrusive contacts. Isotopic temperatures generally show good agreement with heat flow considerations. Based on the experimentally determined quartz-muscovite O¹⁸/O¹⁶ fractionation calibration curve, temperatures are estimated to be 525 to 625°C at the contacts of the granitic stocks studied.

Small-scale oxygen isotope exchange effects between intrusive and country rock are observed over distances of 0.5 to 3 feet on both sides of the contacts; the isotopic gradients are typically 2 to 3 per mil per foot. The degree of oxygen isotopic exchange is essentially identical for different coexisting minerals. This presumably occurred through a diffusion-controlled recrystallization process. The size of the oxygen isotope equilibrium systems in the small-scale exchanged zones vary from about 1.5 cm to 30 cm. A xenolith and a re-entrant of country rock projecting into on intrusive hove both undergone much more extensive isotopic exchange (to hundreds of feet); they also show abnormally high isotopic temperatures. The marginal portions of most plutons have unusually high O¹⁸/O¹⁶ ratios compared to "normal" igneous rocks, presumably due to large-scale isotopic exchange with meta-sedimentary country rocks when the igneous rocks were essentially in a molten state. The isotopic data suggest that outward horizontal movement of H₂O into the contact metamorphic aureoles is almost negligible, but upward movement of H₂O may be important. Also, direct influx and absorption of water from the country rock may be significant in certain intrusive stocks.

Except in the exchanged zones, the O¹⁸/O¹⁶ ratios of pelitic rocks do not change appreciably during contact metamorphism, even in the cordierite and sillimanite grades; this is in contrast to regional metamorphic rocks which commonly decrease in O¹⁸ with increasing grade. Low O¹⁸/O¹⁶ and C¹³/C¹² ratios of the contact metamorphic marbles generally correlate well with the presence of calc-silicate minerals, indicating that the CO₂ liberated during metamorphic decarbonation reactions is enriched in both O¹⁸ and C¹³ relative to the carbonates.

The D/H ratios of biotites in the contact metamorphic rocks and their associated intrusions show a geographic correlation that is similar to that shown by the D/H ratios of meteoric surface waters, perhaps indicating that meteoric waters were present in the rocks during crystallization of the biotites.

High pressure states in condensed matter: I. High pressure behavior of the iron-sulphur system with applications to the earth's core. II. Empirical equation of state for organic compounds at high pressures

Part I:

The earth's core is generally accepted to be composed primarily of iron, with an admixture of other elements. Because the outer core is observed not to transmit shear waves at seismic frequencies, it is known to be liquid or primarily liquid. A new equation of state is presented for liquid iron, in the form of parameters for the 4th order Birch-Murnaghan and Mie-Grüneisen equations of state. The parameters were constrained by a set of values for numerous properties compiled from the literature. A detailed theoretical model is used to constrain the P-T behavior of the heat capacity, based on recent advances in the understanding of the interatomic potentials for transition metals. At the reference pressure of 10⁵ Pa and temperature of 1811 K (the normal melting point of Fe), the parameters are: ρ = 7037 kg/m³, K_S0 = 110 GPa, K_S' = 4.53, K_S" = -.0337 GPa-1, and γ = 2.8, with γ α ρ^-1.17. Comparison of the properties predicted by this model with the earth model PREM indicates that the outer core is 8 to 10 % less dense than pure liquid Fe at the same conditions. The inner core is also found to be 3 to 5% less dense than pure liquid Fe, supporting the idea of a partially molten inner core. The density deficit of the outer core implies that the elements dissolved in the liquid Fe are predominantly of lower atomic weight than Fe. Of the candidate light elements favored by researchers, only sulfur readily dissolves into Fe at low pressure, which means that this element was almost certainly concentrated in the core at early times. New melting data are presented for FeS and FeS₂ which indicate that the FeS₂ is the S-hearing liquidus solid phase at inner core pressures. Consideration of the requirement that the inner core boundary be observable by seismological means and the freezing behavior of solutions leads to the possibility that the outer core may contain a significant fraction of solid material. It is found that convection in the outer core is not hindered if the solid particles are entrained in the fluid flow. This model for a core of Fe and S admits temperatures in the range 3450K to 4200K at the top of the core. An all liquid Fe-S outer core would require a temperature of about 4900 K at the top of the core.

Part II.

The abundance of uses for organic compounds in the modern world results in many applications in which these materials are subjected to high pressures. This leads to the desire to be able to describe the behavior of these materials under such conditions. Unfortunately, the number of compounds is much greater than the number of experimental data available for many of the important properties. In the past, one approach that has worked well is the calculation of appropriate properties by summing the contributions from the organic functional groups making up molecules of the compounds in question. A new set of group contributions for the molar volume, volume thermal expansivity, heat capacity, and the Rao function is presented for functional groups containing C, H, and O. This set is, in most cases, limited in application to low molecular liquids. A new technique for the calculation of the pressure derivative of the bulk modulus is also presented. Comparison with data indicates that the presented technique works very well for most low molecular hydrocarbon liquids and somewhat less well for oxygen-bearing compounds. A similar comparison of previous results for polymers indicates that the existing tabulations of group contributions for this class of materials is in need of revision. There is also evidence that the Rao function contributions for polymers and low molecular compounds are somewhat different.

The Equation of State and Petrogenesis of Komatiite

(1) Equation of State of Komatiite

The equation of state (EOS) of a molten komatiite (27 wt% MgO) was detennined in the 5 to 36 GPa pressure range via shock wave compression from 1550°C and 0 bar. Shock wave velocity, U_S, and particle velocity, U_P, in km/s follow the linear relationship U_S = 3.13(±0.03) + 1.47(±0.03) U_P. Based on a calculated density at 1550°C, 0 bar of 2.745±0.005 glee, this U_S-U_P relationship gives the isentropic bulk modulus K_S = 27.0 ± 0.6 GPa, and its first and second isentropic pressure derivatives, K'_S = 4.9 ± 0.1 and K"_S = -0.109 ± 0.003 GPa^-1.

The calculated liquidus compression curve agrees within error with the static compression results of Agee and Walker [1988a] to 6 GPa. We detennine that olivine (FO₉₄) will be neutrally buoyant in komatiitic melt of the composition we studied near 8.2 GPa. Clinopyroxene would also be neutrally buoyant near this pressure. Liquidus garnet-majorite may be less dense than this komatiitic liquid in the 20-24 GPa interval, however pyropic-garnet and perovskite phases are denser than this komatiitic liquid in their respective liquidus pressure intervals to 36 GPa. Liquidus perovskite may be neutrally buoyant near 70 GPa.

At 40 GPa, the density of shock-compressed molten komatiite would be approximately equal to the calculated density of an equivalent mixture of dense solid oxide components. This observation supports the model of Rigden et al. [1989] for compressibilities of liquid oxide components. Using their theoretical EOS for liquid forsterite and fayalite, we calculate the densities of a spectrum of melts from basaltic through peridotitic that are related to the experimentally studied komatiitic liquid by addition or subtraction of olivine. At low pressure, olivine fractionation lowers the density of basic magmas, but above 14 GPa this trend is reversed. All of these basic to ultrabasic liquids are predicted to have similar densities at 14 GPa, and this density is approximately equal to the bulk (PREM) mantle. This suggests that melts derived from a peridotitic mantle may be inhibited from ascending from depths greater than 400 km.

The EOS of ultrabasic magmas was used to model adiabatic melting in a peridotitic mantle. If komatiites are formed by >15% partial melting of a peridotitic mantle, then komatiites generated by adiabatic melting come from source regions in the lower transition zone (≈500-670 km) or the lower mantle (>670 km). The great depth of incipient melting implied by this model, and the melt density constraint mentioned above, suggest that komatiitic volcanism may be gravitationally hindered. Although komatiitic magmas are thought to separate from their coexisting crystals at a temperature =200°C greater than that for modern MORBs, their ultimate sources are predicted to be diapirs that, if adiabatically decompressed from initially solid mantle, were more than 700°C hotter than the sources of MORBs and derived from great depth.

We considered the evolution of an initially molten mantle, i.e., a magma ocean. Our model considers the thermal structure of the magma ocean, density constraints on crystal segregation, and approximate phase relationships for a nominally chondritic mantle. Crystallization will begin at the core-mantle boundary. Perovskite buoyancy at > 70 GPa may lead to a compositionally stratified lower mantle with iron-enriched mangesiowiistite content increasing with depth. The upper mantle may be depleted in perovskite components. Olivine neutral buoyancy may lead to the formation of a dunite septum in the upper mantle, partitioning the ocean into upper and lower reservoirs, but this septum must be permeable.

(2) Viscosity Measurement with Shock Waves

We have examined in detail the analytical method for measuring shear viscosity from the decay of perturbations on a corrugated shock front The relevance of initial conditions, finite shock amplitude, bulk viscosity, and the sensitivity of the measurements to the shock boundary conditions are discussed. The validity of the viscous perturbation approach is examined by numerically solving the second-order Navier-Stokes equations. These numerical experiments indicate that shock instabilities may occur even when the Kontorovich-D'yakov stability criteria are satisfied. The experimental results for water at 15 GPa are discussed, and it is suggested that the large effective viscosity determined by this method may reflect the existence of ice VII on the Rayleigh path of the Hugoniot This interpretation reconciles the experimental results with estimates and measurements obtained by other means, and is consistent with the relationship of the Hugoniot with the phase diagram for water. Sound waves are generated at 4.8 MHz at in the water experiments at 15 GPa. The existence of anelastic absorption modes near this frequency would also lead to large effective viscosity estimates.

(3) Equation of State of Molybdenum at 1400°C

Shock compression data to 96 GPa for pure molybdenum, initially heated to 1400°C, are presented. Finite strain analysis of the data gives a bulk modulus at 1400°C, K'_S. of 244±2 GPa and its pressure derivative, K'_OS of 4. A fit of shock velocity to particle velocity gives the coefficients of U_S = C_O+S U_P to be C_O = 4.77±0.06 km/s and S = 1.43±0.05. From the zero pressure sound speed, C_O, a bulk modulus of 232±6 GPa is calculated that is consistent with extrapolation of ultrasonic elasticity measurements. The temperature derivative of the bulk modulus at zero pressure, θK_OSθT|_P, is approximately -0.012 GPa/K. A thermodynamic model is used to show that the thermodynamic Grüneisen parameter is proportional to the density and independent of temperature. The Mie-Grüneisen equation of state adequately describes the high temperature behavior of molybdenum under the present range of shock loading conditions.

Imaging the earth with ambient noise and earthquakes

In this thesis, I develop the velocity and structure models for the Los Angeles Basin and Southern Peru. The ultimate goal is to better understand the geological processes involved in the basin and subduction zone dynamics. The results are obtained from seismic interferometry using ambient noise and receiver functions using earthquake- generated waves. Some unusual signals specific to the local structures are also studied. The main findings are summarized as follows:

(1) Los Angeles Basin

The shear wave velocities range from 0.5 to 3.0 km/s in the sediments, with lateral gradients at the Newport-Inglewood, Compton-Los Alamitos, and Whittier Faults. The basin is a maximum of 8 km deep along the profile, and the Moho rises to a depth of 17 km under the basin. The basin has a stretch factor of 2.6 in the center decreasing to 1.3 at the edges, and is in approximate isostatic equilibrium. This "high-density" (~1 km spacing) "short-duration" (~1.5 month) experiment may serve as a prototype experiment that will allow basins to be covered by this type of low-cost survey.

(2) Peruvian subduction zone

Two prominent mid-crust structures are revealed in the 70 km thick crust under the Central Andes: a low-velocity zone interpreted as partially molten rocks beneath the Western Cordillera – Altiplano Plateau, and the underthrusting Brazilian Shield beneath the Eastern Cordillera. The low-velocity zone is oblique to the present trench, and possibly indicates the location of the volcanic arcs formed during the steepening of the Oligocene flat slab beneath the Altiplano Plateau.

The Nazca slab changes from normal dipping (~25 degrees) subduction in the southeast to flat subduction in the northwest of the study area. In the flat subduction regime, the slab subducts to ~100 km depth and then remains flat for ~300 km distance before it resumes a normal dipping geometry. The flat part closely follows the topography of the continental Moho above, indicating a strong suction force between the slab and the overriding plate. A high-velocity mantle wedge exists above the western half of the flat slab, which indicates the lack of melting and thus explains the cessation of the volcanism above. The velocity turns to normal values before the slab steepens again, indicating possible resumption of dehydration and ecologitization.

(3) Some unusual signals

Strong higher-mode Rayleigh waves due to the basin structure are observed in the periods less than 5 s. The particle motions provide a good test for distinguishing between the fundamental and higher mode. The precursor and coda waves relative to the interstation Rayleigh waves are observed, and modeled with a strong scatterer located in the active volcanic area in Southern Peru. In contrast with the usual receiver function analysis, multiples are extensively involved in this thesis. In the LA Basin, a good image is only from PpPs multiples, while in Peru, PpPp multiples contribute significantly to the final results.