Science.gov

Sample records for plateaus numerical instabilities

  1. Peaks, plateaus, numerical instabilities, and achievable accuracy in Galerkin and norm minimizing procedures for solving Ax=b

    SciTech Connect

    Cullum, J.

    1994-12-31

    Plots of the residual norms generated by Galerkin procedures for solving Ax = b often exhibit strings of irregular peaks. At seemingly erratic stages in the iterations, peaks appear in the residual norm plot, intervals of iterations over which the norms initially increase and then decrease. Plots of the residual norms generated by related norm minimizing procedures often exhibit long plateaus, sequences of iterations over which reductions in the size of the residual norm are unacceptably small. In an earlier paper the author discussed and derived relationships between such peaks and plateaus within corresponding Galerkin/Norm Minimizing pairs of such methods. In this paper, through a set of numerical experiments, the author examines connections between peaks, plateaus, numerical instabilities, and the achievable accuracy for such pairs of iterative methods. Three pairs of methods, GMRES/Arnoldi, QMR/BCG, and two bidiagonalization methods are studied.

  2. NUMERICAL STUDY ON GROUNDWATER POLLUTION MECHANISM IN THE TWO SHIRASU PLATEAUS IN THE SAME CATCHMENT AREA

    NASA Astrophysics Data System (ADS)

    Nakagawa, Kei; Saito, Masahiko; Momii, Kazuro

    Groundwater pollution by nitrate due to agricultural activities has become serious environmental issues in the two shirasu plateaus of Kasanohara and Kanoyabaru of Osumi peninsula, Kagoshima prefecture. We developed a groundwater recharge model and examined the water budget characteristics. Furthermore, two-dimensional numerical simulations were carried out to evaluate the effects of countermeasures for the contaminated source. The results of water balance analysis showed that runoff is 10-15%, groundwater recharge is 65%, and evapotranspiration is 20-25% of precipitation in the catchment area. The numerical results indicate that it takes 50 years for Kasanohara and 80 years for Kanoyabaru to display the effects of countermeasure of 80% reduction in concentration of the source zone. The long-term numerical simulation suggested that detailed recharge rate fluctuation rather than constant recharge rate should be applied.

  3. Numerical analysis of engine instability

    NASA Astrophysics Data System (ADS)

    Habiballah, M.; Dubois, I.

    Following a literature review on numerical analyses of combustion instability, to give the state of the art in the area, the paper describes the ONERA methodology used to analyze the combustion instability in liquid propellant engines. Attention is also given to a model (named Phedre) which describes the unsteady turbulent two-phase reacting flow in a liquid rocket engine combustion chamber. The model formulation includes axial or radial propellant injection, baffles, and acoustic resonators modeling, and makes it possible to treat different engine types. A numerical analysis of a cryogenic engine stability is presented, and the results of the analysis are compared with results of tests of the Viking engine and the gas generator of the Vulcain engine, showing good qualitative agreement and some general trends between experiments and numerical analysis.

  4. Mountain growth, orographic precipitation, and the formation of high-plateaus. Insights from numerical modeling experiments

    NASA Astrophysics Data System (ADS)

    Garcia-Castellanos, D.

    2011-12-01

    Orographic precipitation is thought to exert a significant control on the growth of orogens, namely on the distribution of tectonic uplift and deformation in the multi-million-year time-scale. Computer tectonic models based on simple physical laws predict that the concentration of erosional surface mass removal in either side of an orogen can have a significant impact on the how tectonic deformation propagates within the orogenic building. Proxies or simple formulations of climatic processes such as precipitation and evaporation are key to such coupled tectonic/climatic models. Intramountain high plateaus are extreme cases of climatic control on sediment flow, because they often consist of internally-drained (endorheic) basins that trap most erosional products within the orogen. Such intramountain sediment accumulation leads to high-plateaus characterized by high topography and, in spite of the compressional tectonics, low relief. The development of high-plateaus has been extensively studied in relation to the tectonic setting and to inherited structures in the crust, but in this presentation I will concentrate on the potential climatic controls. Computer modeling results will be shown suggesting that a preexisting dry climate, through its effects on surface transport and orogenic evolution, may significantly modify tectonic deformation and contribute to the later formation of a high plateau. This is due to the following feedback: (1) dry climatic conditions at the early stages of orogenesis favor the tectonic defeat of rivers draining the orogen, promoting internal drainage (endorheism); (2) endorheism extends the life of intramountain basins maximizing the mass trapped within the orogen, rising lithostatic stress underneath and expelling fault deformation towards the flanks of the orogen; and (3) this propagation of deformation further isolates the central parts of the orogen from incoming precipitation, reinforcing sediment trapping and topographic leveling of

  5. Oceanic Plateaus

    NASA Astrophysics Data System (ADS)

    Kerr, A. C.

    2003-12-01

    Although the existence of large continental flood basalt provinces has been known for some considerable time, e.g., Holmes (1918), the recognition that similar flood basalt provinces also exist below the oceans is relatively recent. In the early 1970s increasing amounts of evidence from seismic reflection and refraction studies revealed that the crust in several large portions of the ocean floor is significantly thicker than "normal" oceanic crust, which is 6-7 km thick. One of the first areas of such over-thickened crust to be identified was the Caribbean plate ( Edgar et al., 1971) which Donnelly (1973) proposed to be an "oceanic flood basalt province". The term oceanic plateau was coined by Kroenke (1974), and was prompted by the discovery of a large area of thickened crust (>30 km) in the western Pacific known as the Ontong Java plateau (OJP). As our knowledge of the ocean basins has improved over the last 25 years, many more oceanic plateaus have been identified ( Figure 1). Coffin and Eldholm (1992) introduced the term "large igneous provinces" (LIPs) as a generic term encompassing oceanic plateaus, continental flood basalt provinces, and those provinces which form at the continent-ocean boundary (volcanic rifted margins). (22K)Figure 1. Map showing all major oceanic plateaus, and other large igneous provinces discussed in the text (after Saunders et al., 1992). LIPs are generally believed to be formed by decompression melting of upwelling hotter mantle, known as mantle plumes. Although ideas about hotpots and mantle plumes have been around for almost 40 years (Wilson, 1963), it is only in the past 15 years that LIPs have become the focus of major research. One of the main reasons for the increased research activity into LIPs is the realization that significant proportions of these LIPs erupted over a relatively short time, often less than 2-3 Myr (see review in Coffin, 1994). This has important implications for mantle processes and source regions ( Hart et

  6. A Numerical Study of Feathering Instability

    NASA Astrophysics Data System (ADS)

    Lee, Wing-Kit; Wang, Hsiang-Hsu

    2016-06-01

    The stability of a spiral shock of self-gravitating, magnetized interstellar medium is studied by performing two-dimensional numerical simulations of a local patch of tight-winding spiral arm. As previously suggested by the linear studies, two types of instabilities are identified, namely, wiggle instability and feathering instability. The former instability occurs in the hydrodynamics limit and results in short wavelength perturbations. On the other hand, the feathering instability requires both self-gravitating and magnetic fields and results in wider structures.

  7. The development of topographic plateaus in an India-Asia-like collision zone using 3D numerical simulations

    NASA Astrophysics Data System (ADS)

    Pusok, Adina E.; Kaus, Boris; Popov, Anton

    2014-05-01

    . We also address the question of how large topographic plateaus, such as the Tibetan Plateau, can form in an integrated lithospheric and upper-mantle scale model. Acknowledgements: Funding was provided by the European Research Council under the European Community's Seventh Framework Program (FP7/2007-2013) / ERC Grant agreement #258830. Numerical computations have been performed on MOGON (ZDV Mainz computing center) and JUQUEEN (Jülich high-performance computing center).

  8. Understanding Etna flank instability through numerical models

    NASA Astrophysics Data System (ADS)

    Apuani, Tiziana; Corazzato, Claudia; Merri, Andrea; Tibaldi, Alessandro

    2013-02-01

    As many active volcanoes, Mount Etna shows clear evidence of flank instability, and different mechanisms were suggested to explain this flank dynamics, based on the recorded deformation pattern and character. Shallow and deep deformations, mainly associated with both eruptive and seismic events, are concentrated along recognised fracture and fault systems, mobilising the eastern and south-eastern flank of the volcano. Several interacting causes were postulated to control the phenomenon, including gravity force, magma ascent along the feeding system, and a very complex local and/or regional tectonic activity. Nevertheless, the complexity of such dynamics is still an open subject of research and being the volcano flanks heavily urbanised, the comprehension of the gravitative dynamics is a major issue for public safety and civil protection. The present research explores the effects of the main geological features (in particular the role of the subetnean clays, interposed between the Apennine-Maghrebian flysch and the volcanic products) and the role of weakness zones, identified by fracture and fault systems, on the slope instability process. The effects of magma intrusions are also investigated. The problem is addressed by integrating field data, laboratory tests and numerical modelling. A bi- and tri-dimensional stress-strain analysis was performed by a finite difference numerical code (FLAC and FLAC3D), mainly aimed at evaluating the relationship among geological features, volcano-tectonic structures and magmatic activity in controlling the deformation processes. The analyses are well supported by dedicated structural-mechanical field surveys, which allowed to estimate the rock mass strength and deformability parameters. To take into account the uncertainties which inevitably occur in a so complicated model, many efforts were done in performing a sensitivity analysis along a WNW-ESE section crossing the volcano summit and the Valle del Bove depression. This was

  9. Numerical Modelling Of Pumpkin Balloon Instability

    NASA Astrophysics Data System (ADS)

    Wakefield, D.

    Tensys have been involved in the numerical formfinding and load analysis of architectural stressed membrane structures for 15 years. They have recently broadened this range of activities into the `lighter than air' field with significant involvement in aerostat and heavy-lift hybrid airship design. Since early 2004 they have been investigating pumpkin balloon instability on behalf of the NASA ULDB programme. These studies are undertaken using inTENS, an in-house finite element program suite based upon the Dynamic Relaxation solution method and developed especially for the non-linear analysis and patterning of membrane structures. The paper describes the current state of an investigation that started with a numerical simulation of the lobed cylinder problem first studied by Calladine. The influence of material properties and local geometric deformation on stability is demonstrated. A number of models of complete pumpkin balloons have then been established, including a 64-gore balloon with geometry based upon Julian Nott's Endeavour. This latter clefted dramatically upon initial inflation, a phenomenon that has been reproduced in the numerical model. Ongoing investigations include the introduction of membrane contact modelling into inTENS and correlation studies with the series of large-scale ULDB models currently in preparation.

  10. Numerical simulation of surface waves instability on a homogeneous grid

    NASA Astrophysics Data System (ADS)

    Korotkevich, Alexander O.; Dyachenko, Alexander I.; Zakharov, Vladimir E.

    2016-05-01

    We performed full-scale numerical simulation of instability of weakly nonlinear waves on the surface of deep fluid. We show that the instability development leads to chaotization and formation of wave turbulence. Instability of both propagating and standing waves was studied. We separately studied pure capillary wave, that was unstable due to three-wave interactions and pure gravity waves, that were unstable due to four-wave interactions. The theoretical description of instabilities in all cases is included in the article. The numerical algorithm used in these and many other previous simulations performed by the authors is described in detail.

  11. Tube Feeding Transition Plateaus

    ERIC Educational Resources Information Center

    Klein, Marsha Dunn

    2007-01-01

    The journey children make from tube feeding to oral feeding is personal for each child and family. There is a sequence of predictable plateaus that children climb as they move toward orally eating. By better understanding this sequence, parents and children can maximize the development, learning, enjoyment and confidence at each plateau. The…

  12. Numerical Simulations of Instabilities in Single-Hole Office Elements

    NASA Technical Reports Server (NTRS)

    Ahuja, Vineet; Hosangadi, Ashvin; Hitt, Matthew A.; Lineberry, David M.

    2013-01-01

    An orifice element is commonly used in liquid rocket engine test facilities either as a flow metering device, a damper for acoustic resonance or to provide a large reduction in pressure over a very small distance in the piping system. While the orifice as a device is largely effective in stepping down pressure, it is also susceptible to a wake-vortex type instability that generates pressure fluctuations that propagate downstream and interact with other elements of the test facility resulting in structural vibrations. Furthermore in piping systems an unstable feedback loop can exist between the vortex shedding and acoustic perturbations from upstream components resulting in an amplification of the modes convecting downstream. Such was the case in several tests conducted at NASA as well as in the Ariane 5 strap-on P230 engine in a static firing test where pressure oscillations of 0.5% resulted in 5% thrust oscillations. Exacerbating the situation in cryogenic test facilities, is the possibility of the formation of vapor clouds when the pressure in the wake falls below the vapor pressure leading to a cavitation instability that has a lower frequency than the primary wake-vortex instability. The cavitation instability has the potential for high amplitude fluctuations that can cause catastrophic damage in the facility. In this paper high-fidelity multi-phase numerical simulations of an orifice element are used to characterize the different instabilities, understand the dominant instability mechanisms and identify the tonal content of the instabilities.

  13. Numerical modeling of the Parker instability in a rotating plasma

    NASA Astrophysics Data System (ADS)

    Khalzov, Ivan; Brown, Ben; Katz, Noam; Forest, Cary

    2011-10-01

    We study numerically the analogue of the Parker (magnetic buoyancy) instability in a rotating plasma screw pinch confined in a bounded cylinder. The goal of the study is to show the possibility of reaching the Parker instability for the plasma parameters achievable in the Madison Plasma Couette Experiment (MPCX). Simulations are performed using the extended magnetohydrodynamic (MHD) code NIMROD for an isothermal compressible plasma model. Both linear and nonlinear regimes of the instability are studied, and the obtained results are compared with analytic results for a slab geometry. It is shown that the effect of plasma rotation in a cylindrical geometry is two-fold: first, centrifugal acceleration acts as analogue of gravity and provides the equilibrium density stratification; second, the presence of Coriolis force results in increase of critical gradient of magnetic field required for the onset of instability.

  14. Numerical simulation of the hydrodynamic instability experiments and flow mixing

    NASA Astrophysics Data System (ADS)

    Bai, Jingsong; Wang, Tao; Li, Ping; Zou, Liyong; Liu, Cangli

    2009-12-01

    Based on the numerical methods of volume of fluid (VOF) and piecewise parabolic method (PPM) and parallel circumstance of Message Passing Interface (MPI), a parallel multi-viscosity-fluid hydrodynamic code MVPPM (Multi-Viscosity-Fluid Piecewise Parabolic Method) is developed and performed to study the hydrodynamic instability and flow mixing. Firstly, the MVPPM code is verified and validated by simulating three instability cases: The first one is a Riemann problem of viscous flow on the shock tube; the second one is the hydrodynamic instability and mixing of gaseous flows under re-shocks; the third one is a half height experiment of interfacial instability, which is conducted on the AWE’s shock tube. By comparing the numerical results with experimental data, good agreement is achieved. Then the MVPPM code is applied to simulate the two cases of the interfacial instabilities of jelly models accelerated by explosion products of a gaseous explosive mixture (GEM), which are adopted in our experiments. The first is implosive dynamic interfacial instability of cylindrical symmetry and mixing. The evolving process of inner and outer interfaces, and the late distribution of mixing mass caused by Rayleigh-Taylor (RT) instability in the center of different radius are given. The second is jelly layer experiment which is initialized with one periodic perturbation with different amplitude and wave length. It reveals the complex processes of evolution of interface, and presents the displacement of front face of jelly layer, bubble head and top of spike relative to initial equilibrium position vs. time. The numerical results are in excellent agreement with that experimental images, and show that the amplitude of initial perturbations affects the evolvement of fluid mixing zone (FMZ) growth rate extremely, especially at late times.

  15. Collisional plateaus. [in earth and Venus lithospheres

    NASA Technical Reports Server (NTRS)

    Morgan, P.; Burke, K.

    1985-01-01

    Aspects of the geology of collisional plateaus formed by the thickening of continental crust are briefly reviewed. The history of studies of collisional plateaus is summarized, and igneous activity in collisional plateaus is discussed. Isostatic considerations pertaining to these plateaus are addressed, developing models of isostatic support of topography which illustrate the importance of compressional tectonics in the creation of high altitude plateaus. Possible analogous environments on Venus are considered. Finally, the paradox of extension associated with compression in the plateaus is discussed.

  16. Proposed Ozark Plateaus Province Hydrologic Observatory

    NASA Astrophysics Data System (ADS)

    Davis, R. K.; Brahana, J. V.; Matlock, M. D.; Chaubey, I.; Pavlowsky, R. T.; Gilzow, F.

    2004-12-01

    The Upper White River, which drains about 40 percent of the Ozarks Ecoregion, is the main drain for the Ozark Plateaus and is characteristic of rivers draining other karst areas within the United States and the world. The proposed Ozark Plateaus Hydrologic Observatory (OPHO) encompasses twelve 8-digit hydrologic units covering about 67,000 km2 in parts of three states (Arkansas, Missouri, and Oklahoma). Six major U.S. Army Corps of Engineers reservoirs are within the OPHO including four on the main stem of the White River and one on the Illinois River. Karst features are prominent in the Salem, Ozark, and Springfield Plateaus of the OPHO, and include numerous solutionally enlarged fractures, caves, sinkholes, and sinking streams. Within the basin are numerous and diverse biological communities, representing influences from 1) eastern deciduous forest, 2) Great Plains prairies, 3) arid southwest, and 4) relicts of northern species from the Pleistocene Ice Age. Also contain in the OPHO is a diverse and unique array of mussels, an imperiled river organism (38 species), and crayfish. In the extensive karst regions of the OPHO are found largely endemic subterranean organisms also dependent on good water quality: for example, the Ozark Cavefish, Bristly Cave Crayfish and the recently federally- listed Tumbling Creek Cave Snail. Mantled karst aquifers characteristic of the Ozark Plateaus Region represent a coupled atmospheric/surface water/groundwater system that is highly susceptible to external forcing. Little attenuation of contaminants occurs as water moves from surface sources into and through the mantled karst aquifer to discharge naturally at springs and streams throughout the Ozark Plateau Region, and to wells. Because of the very open character of the aquifer, extremely dynamic biogeochemical cycling of nutrients occurs. Upper White River Reservoir development, filling and operation historically have altered and continue to alter the hydrologic and ecosystems

  17. Numerical Investigation of Galloping Instabilities in Z-Shaped Profiles

    PubMed Central

    Chavez, Miguel; Valero, Eusebio

    2014-01-01

    Aeroelastic effects are relatively common in the design of modern civil constructions such as office blocks, airport terminal buildings, and factories. Typical flexible structures exposed to the action of wind are shading devices, normally slats or louvers. A typical cross-section for such elements is a Z-shaped profile, made out of a central web and two-side wings. Galloping instabilities are often determined in practice using the Glauert-Den Hartog criterion. This criterion relies on accurate predictions of the dependence of the aerodynamic force coefficients with the angle of attack. The results of a parametric analysis based on a numerical analysis and performed on different Z-shaped louvers to determine translational galloping instability regions are presented in this paper. These numerical analysis results have been validated with a parametric analysis of Z-shaped profiles based on static wind tunnel tests. In order to perform this validation, the DLR TAU Code, which is a standard code within the European aeronautical industry, has been used. This study highlights the focus on the numerical prediction of the effect of galloping, which is shown in a visible way, through stability maps. Comparisons between numerical and experimental data are presented with respect to various meshes and turbulence models. PMID:25054176

  18. Numerical aspects of searching convective/absolute instability transition

    NASA Astrophysics Data System (ADS)

    Suslov, Sergey A.

    2006-02-01

    An overview of various numerical techniques used to determine the spatio-temporal character of instabilities in fluid flows is given. The advantageous features of various previously known individual techniques are discussed and a practical procedure combining them is suggested for a specific task of determining the complete boundary between linearly convectively and absolutely unstable regimes in a multi-parameter space in problems with a fully numerical dispersion relation. Special attention is paid to aspects of automatization of computations as this is a crucial condition for their efficiency. The suggested procedure is successfully used and is shown to provide a high degree of automatism in the physical example of non-Boussinesq mixed convection in a vertical channel. This example comprises most of the major numerical difficulties found in various spatio-temporal instability studies of two-dimensional fluid flows which previously could not be handled without frequent human intervention and visual inspection of intermediate results. This paper focuses on the general numerical aspects of the computations leaving the detailed discussion of the obtained physical results for a separate publication.

  19. Numerical investigation of galloping instabilities in Z-shaped profiles.

    PubMed

    Gomez, Ignacio; Chavez, Miguel; Alonso, Gustavo; Valero, Eusebio

    2014-01-01

    Aeroelastic effects are relatively common in the design of modern civil constructions such as office blocks, airport terminal buildings, and factories. Typical flexible structures exposed to the action of wind are shading devices, normally slats or louvers. A typical cross-section for such elements is a Z-shaped profile, made out of a central web and two-side wings. Galloping instabilities are often determined in practice using the Glauert-Den Hartog criterion. This criterion relies on accurate predictions of the dependence of the aerodynamic force coefficients with the angle of attack. The results of a parametric analysis based on a numerical analysis and performed on different Z-shaped louvers to determine translational galloping instability regions are presented in this paper. These numerical analysis results have been validated with a parametric analysis of Z-shaped profiles based on static wind tunnel tests. In order to perform this validation, the DLR TAU Code, which is a standard code within the European aeronautical industry, has been used. This study highlights the focus on the numerical prediction of the effect of galloping, which is shown in a visible way, through stability maps. Comparisons between numerical and experimental data are presented with respect to various meshes and turbulence models. PMID:25054176

  20. Numerical Simulation of Flow Instability and Heat Transfer

    NASA Astrophysics Data System (ADS)

    Dou, Hua-Shu; Jiang, Gang

    2014-11-01

    This paper numerically investigates the physical mechanism of flow instability and heat transfer of natural convection in a cavity with thin fin(s). The left and the right walls of the cavity are differentially heated. The cavity is given an initial temperature, and the thin fin(s) is fixed on the hot wall in order to control the heat transfer. The finite volume method with the SIMPLE scheme is used to simulate the flow. Distributions of the temperature, the pressure, the velocity and the total pressure are achieved. Then, the energy gradient method is employed to study the physical mechanism of flow instability and the effect of the thin fin(s) on heat transfer. Based on the energy gradient method, the energy gradient function K represents the characteristic of flow instability. It is observed from the simulation results that the positions where instabilities take place in the temperature contours accord well with those of higher K value, which demonstrates that the energy gradient method reveals the physical mechanism of flow instability. Furthermore, the effect of the fin length, the fin position, the fin number, and Ra on heat transfer is also investigated. It is found that the effect of the fin length on heat transfer is negligible when Ra is relatively high. When there is only one fin, the most efficient heat transfer rate is achieved as the fin is fixed at the middle height of the cavity. The fin blocks heat transfer with a relatively small Ra, but the fin enhances heat transfer with a relatively large Ra. The fin(s) enhances heat transfer gradually with the increase of Ra under the influence of the thin fin(s). Finally, it is observed that both Kmax and Ra can reveal the physical mechanism of natural convection from different approaches.

  1. A Numerical Instability in an ADI Algorithm for Gyrokinetics

    SciTech Connect

    E.A. Belli; G.W. Hammett

    2004-12-17

    We explore the implementation of an Alternating Direction Implicit (ADI) algorithm for a gyrokinetic plasma problem and its resulting numerical stability properties. This algorithm, which uses a standard ADI scheme to divide the field solve from the particle distribution function advance, has previously been found to work well for certain plasma kinetic problems involving one spatial and two velocity dimensions, including collisions and an electric field. However, for the gyrokinetic problem we find a severe stability restriction on the time step. Furthermore, we find that this numerical instability limitation also affects some other algorithms, such as a partially implicit Adams-Bashforth algorithm, where the parallel motion operator v{sub {parallel}} {partial_derivative}/{partial_derivative}z is treated implicitly and the field terms are treated with an Adams-Bashforth explicit scheme. Fully explicit algorithms applied to all terms can be better at long wavelengths than these ADI or partially implicit algorithms.

  2. NUMERICAL STUDY OF THE VISHNIAC INSTABILITY IN SUPERNOVA REMNANTS

    SciTech Connect

    Michaut, C.; Cavet, C.; Bouquet, S. E.; Roy, F.; Nguyen, H. C.

    2012-11-10

    The Vishniac instability is thought to explain the complex structure of radiative supernova remnants in their Pressure-Driven Thin Shell (PDTS) phase after a blast wave (BW) has propagated from a central explosion. In this paper, the propagation of the BW and the evolution of the PDTS stage are studied numerically with the two-dimensional (2D) code HYDRO-MUSCL for a finite-thickness shell expanding in the interstellar medium (ISM). Special attention is paid to the adiabatic index, {gamma}, and three distinct values are taken for the cavity ({gamma}{sub 1}), the shell ({gamma}{sub 2}), and the ISM ({gamma}{sub 3}) with the condition {gamma}{sub 2} < {gamma}{sub 1}, {gamma}{sub 3}. This low value of {gamma}{sub 2} accounts for the high density in the shell achieved by a strong radiative cooling. Once the spherical background flow is obtained, the evolution of a 2D-axisymmetric perturbation is computed from the linear to the nonlinear regime. The overstable mechanism, previously demonstrated theoretically by E. T. Vishniac in 1983, is recovered numerically in the linear stage and is expected to produce and enhance anisotropies and clumps on the shock front, leading to the disruption of the shell in the nonlinear phase. The period of the increasing oscillations and the growth rate of the instability are derived from several points of view (the position of the perturbed shock front, mass fluxes along the shell, and density maps), and the most unstable mode differing from the value given by Vishniac is computed. In addition, the influence of several parameters (the Mach number, amplitude and wavelength of the perturbation, and adiabatic index) is examined and for wavelengths that are large enough compared to the shell thickness, the same conclusion arises: in the late stage of the evolution of the radiative supernova remnant, the instability is dampened and the angular initial deformation of the shock front is smoothed while the mass density becomes uniform with the

  3. A numerical and analytical investigation of Rayleigh-Taylor instability in a solid tungsten plate

    SciTech Connect

    Robinson, A.C.; Swegle, J.W.

    1987-07-01

    The Rayleigh-Taylor instability response of an elastic-plastic tungsten plate is investigated by numerical experiments and an approximate modal analysis. The so-called ''minimum amplitude'' instability criteria derived from plasticity analyses is shown to be incomplete as a general indicator of instability or stability at very large driving pressures. Model equations are derived which are able to reproduce the basic qualitative features of the observed instability response given by the numerical calculations. 11 refs., 29 figs.

  4. Plasma self-heating and saturation due to numerical instabilities

    SciTech Connect

    Birdsall, C.K.; Maron, N.

    1980-06-01

    The cold-beam nonphysical instability due to the aliases produced by the numerical spatial grid is presented in detail for momentum-conserving linear weighting codes. Additions to previous work include: linear analysis dispersion diagrams showing large growth rates, ..omega../sub i/< or =0.2 ..omega../sub p/; methods for reducing ..omega../sub i/, effectively broadening the finite-size particle width; simulation results verifying the linear theory, plus plots showing the p=1 alias in phase space (..nu../sub x/, x); the growth of beam, thermal spread (..nu../sub t//sup 2/ in temperature) an loss of energy conservation; end-of-growth (saturation at small thermal spread (lambda/sub D//..delta..x=upsilon/sub t//..omega../sub p/..delta..xapprox. =0.046, for lambda/sub B//..delta..xequivalent upsilon/sub 0//..omega../sub p/..delta..x> or =0.3, i.e., upsilon/sub t/< or =0.14upsilon/sub 0/), with return to near energy conservation (stability); demonstration of no growth for a warm beam, with upsilon/sub t/(initial)> or =upsilon/sub t/(saturation); and the mechanism of stabilization (trapping). A thermal (Maxwellian) plasma, also nonphysically unstable at small lambda/sub D//..delta..x, is also found to approach stabilization by self-heating. The two-stream physical instability is affected by the grid, with the aliasing instability also present; the linear theory for this is presented, with guidelines for minimizing the effects of the grid.

  5. A Cartesian parametrization for the numerical analysis of material instability

    DOE PAGESBeta

    Mota, Alejandro; Chen, Qiushi; Foulk, III, James W.; Ostien, Jakob T.; Lai, Zhengshou

    2016-02-25

    We examine four parametrizations of the unit sphere in the context of material stability analysis by means of the singularity of the acoustic tensor. We then propose a Cartesian parametrization for vectors that lie a cube of side length two and use these vectors in lieu of unit normals to test for the loss of the ellipticity condition. This parametrization is then used to construct a tensor akin to the acoustic tensor. It is shown that both of these tensors become singular at the same time and in the same planes in the presence of a material instability. Furthermore, themore » performance of the Cartesian parametrization is compared against the other parametrizations, with the results of these comparisons showing that in general, the Cartesian parametrization is more robust and more numerically efficient than the others.« less

  6. Convective instability in sedimentation: 3-D numerical study

    NASA Astrophysics Data System (ADS)

    Yu, Xiao; Hsu, Tian-Jian; Balachandar, S.

    2014-11-01

    To provide a probable explanation on the field observed rapid sedimentation process near river mouths, we investigate the convective sedimentation in stably stratified saltwater using 3-D numerical simulations. Guided by the linear stability analysis, this study focuses on the nonlinear interactions of several mechanisms, which lead to various sediment finger patterns, and the effective settling velocity for sediment ranging from clay (single-particle settling velocity V0 = 0.0036 and 0.0144 mm/s, or particle diameter d = 2 and 4 μm) to silt (V0 = 0.36 mm/s, or d = 20 μm). For very fine sediment with V0 = 0.0036 mm/s, the convective instability is dominated by double diffusion, characterized by millimeter-scale fingers. Gravitational settling slightly increases the growth rate; however, it has notable effect on the downward development of vertical mixing shortly after the sediment interface migrates below the salt interface. For sediment with V0 = 0.0144 mm/s, Rayleigh-Taylor instabilities become dominant before double-diffusive modes grow sufficiently large. Centimeter-scale and highly asymmetric sediment fingers are obtained due to nonlinear interactions between different modes. For sediment with V0 = 0.36 mm/s, Rayleigh-Taylor mechanism dominates and the resulting centimeter-scale sediment fingers show a plume-like structure. The flow pattern is similar to that without ambient salt stratification. Rapid sedimentation with effective settling velocity on the order of 1 cm/s is likely driven by convective sedimentation for sediment with V0 greater than 0.1 mm/s at concentration greater than 10-20 g/L.

  7. Dynamics of plume-triple junction interaction: Results from a series of three-dimensional numerical models and implications for the formation of oceanic plateaus

    NASA Astrophysics Data System (ADS)

    Dordevic, Mladen; Georgen, Jennifer

    2016-03-01

    Mantle plumes rising in the vicinity of mid-ocean ridges often generate anomalies in melt production and seafloor depth. This study investigates the dynamical interactions between a mantle plume and a ridge-ridge-ridge triple junction, using a parameter space approach and a suite of steady state, three-dimensional finite element numerical models. The top domain boundary is composed of three diverging plates, with each assigned half-spreading rates with respect to a fixed triple junction point. The bottom boundary is kept at a constant temperature of 1350°C except where a two-dimensional, Gaussian-shaped thermal anomaly simulating a plume is imposed. Models vary plume diameter, plume location, the viscosity contrast between plume and ambient mantle material, and the use of dehydration rheology in calculating viscosity. Importantly, the model results quantify how plume-related anomalies in mantle temperature pattern, seafloor depth, and crustal thickness depend on the specific set of parameters. To provide an example, one way of assessing the effect of conduit position is to calculate normalized area, defined to be the spatial dispersion of a given plume at specific depth (here selected to be 50 km) divided by the area occupied by the same plume when it is located under the triple junction. For one particular case modeled where the plume is centered in an intraplate position 100 km from the triple junction, normalized area is just 55%. Overall, these models provide a framework for better understanding plateau formation at triple junctions in the natural setting and a tool for constraining subsurface geodynamical processes and plume properties.

  8. Numerical Study of a Hydrodynamic Instability Driven by Evaporation

    NASA Astrophysics Data System (ADS)

    Hernandez-Zapata, Sergio; Romo-Cruz, Julio Cesar Ruben; Lopez-Sanchez, Erick Javier; Ruiz-Chavarria, Gerardo

    2013-11-01

    The study of hydrodynamic instabilities in liquid layers produced by evaporation has several applications on industry and technology. In this work we study numerically the conditions under which a liquid layer becomes unstable when evaporation in the vapor-liquid interphase is present. The evaporation process follows the Hertz-Knudsen law (the evaporation rate is proportional to the difference between the saturated vapor pressure at the liquid layer temperature and the vapor partial pressure in the environment). Additionally to the usual boundary conditions on solid walls (for example, the non-slip condition for the velocity), we analyze the boundary conditions in the vapor-liquid interphase where the momentum and energy balances have to be taken into account and where the evaporation plays a crucial role. To solve this problem the linear theory of stability is used; that is, a small perturbation around the basic solution is applied (flow at rest and a temperature stationary field). The equations are solved using the Chebyshev pseudo-spectral method. The results are compared with the more usual Rayleigh-Bénard and Marangoni mechanisms as well as with some experiments carried out by our team. Authors acknowledge DGAPA-UNAM by support under project IN116312, ``Vorticidad y Ondas no lineales en fluidos.''

  9. Numerical investigation of the transverse instability on the radiation-pressure-driven foil

    NASA Astrophysics Data System (ADS)

    Wang, W. Q.; Yin, Y.; Yu, T. P.; Xu, H.; Zou, D. B.; Shao, F. Q.

    2015-12-01

    The development of transverse instability in the radiation-pressure-acceleration dominant laser-foil interaction is numerically examined by two-dimensional particle-in-cell simulations. When a plane laser impinges on a foil with modulated surface, the transverse instability is incited, and periodic perturbations of the proton density develop. The growth rate of the transverse instability is numerically diagnosed. It is found that the linear growth of the transverse instability lasts only a few laser periods, then the instability gets saturated. In order to optimize the modulation wavelength of the target, a method of information entropy is put forward to describe the chaos degree of the transverse instability. With appropriate modulation, the transverse instability shows a low chaos degree, and a quasi-monoenergetic proton beam is produced.

  10. Numerical Prediction of Laminar Instability Noise for NACA 0012 Aerofoil

    NASA Astrophysics Data System (ADS)

    De Gennaro, Michele; Hueppe, Andreas; Kuehnelt, Helmut; Kaltenbacher, Manfred

    2011-09-01

    Aerofoil self-generated noise is recognized to be of fundamental importance in the frame of applied aeroacoustics and the use of computational methods to assess the acoustic behaviour of airframe components challenges an even larger community of engineers and scientists. Several noise generation mechanisms can be found which are mainly related to the physical development of turbulence over the boundary layer. They can be classified in 3 main categories: the Turbulent Boundary Layer—Trailing Edge noise (TBL-TE), the Laminar Boundary Layer—Vortex Shedding (LBL-VS) noise and the Separation Stall (S-S) noise. The TBL-TE is mainly related to the noise generated by turbulent eddies which develop into the boundary layer and usually exhibits a broadband spectrum. The LBL-VS is related to laminar instabilities that can occur within the boundary layer which are responsible for a very late transition and generate a typical peaked tonal noise, while the S-S noise mainly results from the development of large vortices after the separation point. In this paper we propose a numerical analysis targeted to the simulation the LBL-VS noise mechanisms on a NACA 0012 aerofoil, tested at a Reynolds number of 1.1 M and Mach number of 0.2. The aerodynamic simulation is performed with a 2D transient RANS approach using the k-ω transitional turbulence model, while the acoustic computations are performed with the FfowcsWilliams-Hawkings (FW-H) acoustic analogy and with a Finite Element (FE) approach solving Lighthill's wave equation. Computed noise spectra are compared with experimental data published by NASA showing a good agreement both for peak location as well as for the predicted noise level.

  11. Comprehensive numerical methodology for direct numerical simulations of compressible Rayleigh-Taylor instability

    NASA Astrophysics Data System (ADS)

    Reckinger, Scott J.; Livescu, Daniel; Vasilyev, Oleg V.

    2016-05-01

    An investigation of compressible Rayleigh-Taylor instability (RTI) using Direct Numerical Simulations (DNS) requires efficient numerical methods, advanced boundary conditions, and consistent initialization in order to capture the wide range of scales and vortex dynamics present in the system, while reducing the computational impact associated with acoustic wave generation and the subsequent interaction with the flow. An advanced computational framework is presented that handles the challenges introduced by considering the compressive nature of RTI systems, which include sharp interfacial density gradients on strongly stratified background states, acoustic wave generation and removal at computational boundaries, and stratification dependent vorticity production. The foundation of the numerical methodology described here is the wavelet-based grid adaptivity of the Parallel Adaptive Wavelet Collocation Method (PAWCM) that maintains symmetry in single-mode RTI systems to extreme late-times. PAWCM is combined with a consistent initialization, which reduces the generation of acoustic disturbances, and effective boundary treatments, which prevent acoustic reflections. A dynamic time integration scheme that can handle highly nonlinear and potentially stiff systems, such as compressible RTI, completes the computational framework. The numerical methodology is used to simulate two-dimensional single-mode RTI to extreme late-times for a wide range of flow compressibility and variable density effects. The results show that flow compressibility acts to reduce the growth of RTI for low Atwood numbers, as predicted from linear stability analysis.

  12. Lessons Learned from Numerical Simulations of Interfacial Instabilities

    NASA Astrophysics Data System (ADS)

    Cook, Andrew

    2015-11-01

    Rayleigh-Taylor (RT), Richtmyer-Meshkov (RM) and Kelvin-Helmholtz (KH) instabilities serve as efficient mixing mechanisms in a wide variety of flows, from supernovae to jet engines. Over the past decade, we have used the Miranda code to temporally integrate the multi-component Navier-Stokes equations at spatial resolutions up to 29 billion grid points. The code employs 10th-order compact schemes for spatial derivatives, combined with 4th-order Runge-Kutta time advancement. Some of our major findings are as follows: The rate of growth of a mixing layer is equivalent to the net mass flux through the equi-molar plane. RT growth rates can be significantly reduced by adding shear. RT instability can produce shock waves. The growth rate of RM instability can be predicted from known interfacial perturbations. RM vortex projectiles can far outrun the mixing region. Thermal fluctuations in molecular dynamics simulations can seed instabilities along the braids in KH instability. And finally, enthalpy diffusion is essential in preserving the second law of thermodynamics. This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.

  13. On the Numerical Dispersion of Electromagnetic Particle-In-Cell Code : Finite Grid Instability

    SciTech Connect

    Meyers, Michael David; Huang, Chengkun; Zeng, Yong; Yi, Sunghwan; Albright, Brian James

    2014-07-15

    The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the electromagnetic PIC algorithm to analyze the origin of these instabilities. We rigorously derive the faithful 3D numerical dispersion of the PIC algorithm, and then specialize to the Yee FDTD scheme. In particular, we account for the manner in which the PIC algorithm updates and samples the fields and distribution function. Temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme are also explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical 1D modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction.

  14. Numerical Modeling of Extensional Necking Instabilities: Application to Ganymede's Grooved Terrain

    NASA Technical Reports Server (NTRS)

    Bland, M. T.; Showman, A. P.

    2005-01-01

    Ganymede s pervasive 5-10 km-wavelength grooves have been suggested to result from a necking instability during an epoch of lithospheric extension, but to date few quantitative studies of groove formation have been performed. We present two-dimensional numerical models of necking instabilities under conditions that are appropriate to Ganymede at the time of groove formation. Preliminary simulations indicate that extensional necking instabilities can occur under a range of conditions, many of which may be relevant to Ganymede. The form of the surface topography produced by these instabilities varies as a function of the strain rate, amount of extension, initial topographic perturbation, and rheological parameters.

  15. Numerical studies of baroclinic instability at small Richardson number

    NASA Technical Reports Server (NTRS)

    Miller, T. L.

    1984-01-01

    Baroclinic instability at small Richardson number consists of competition between symmetric (Solberg) modes, zonal (Eady) modes, and possibly other modes whose orientation is neither symmetric nor zonal. In FY-84, a study of the fully nonlinear development and energetics of the symmetric modes was completed. The nonlinear effects were quite strong and, in some senses, unexpected. Secondly, a study of the dependency of the energetics of the symmetric modes upon the physical parameters was completed.

  16. Direct Numerical Simulations of Flame Instabilities in Type Ia Supernovae

    NASA Astrophysics Data System (ADS)

    Zingale, M.; Bell, J. B.; Day, M. S.; Rendleman, C. A.; Woosley, S. E.

    2003-12-01

    Instabilities serve an important role in accelerating a thermonuclear flame in a white dwarf to a large fraction of the speed of sound (perhaps to a supersonic detonation), consuming the carbon/oxygen, and producing a Type Ia supernovae. The precise mechanism for this acceleration is not well understood, but large scale simulations show that a deflagration alone can unbind the star. We present fully resolved, multidimensional calculations of Rayleigh-Taylor unstable flames in conditions appropriate to the late stages of Type Ia SNe, using a low Mach number hydrodynamics code. At densities below 1.e7 g/cc, a fundamental change in the burning is observed, as the flame transitions from the wrinkled flame to the distributed burning regime. Significant acceleration is observed for all densities we study, limited only by the size of the domain we can address. We compare with corresponding simulations of the Landau-Darrieus instability. We discuss the physics of these instabilities on the small scales and the implications they have for large scale flame modeling and the possibility for deflagration to detonation transitions. Support for this work was provided by the DOE grant No. DE-FC02-01ER41176 to the Supernova Science Center/UCSC and the Applied Mathematics Program of the DOE Office of Mathematics, Information, and Computational Sciences under the U.S. Department of Energy under contract No. DE-AC03-76SF00098.

  17. The Zombie Instability: Using Numerical Simulation to Design a Laboratory Experiment

    NASA Astrophysics Data System (ADS)

    Wang, Meng; Pei, Suyang; Jiang, Chung-Hsiang; Hassanzadeh, Pedram; Marcus, Philip

    2014-11-01

    A new type of finite amplitude-instability has been found in numerical simulations of stratified, rotating, shear flows. The instability occurs via baroclinic critical layers that create linearly unstable vortex layers, which roll-up into vortices. Under the right conditions, those vortices can form a new generation of vortices, resulting in ``vortex self-replication'' that fills the fluid with vortices. Creating this instability in a laboratory would provide further evidence for the existence of the instability, which we first found in numerical simulations of protoplanetary disks. To design a laboratory experiment we need to know how the flow parameters-- shear, rotation and stratification, etc. affect the instability. To build an experiment economically, we also need to know how the finite-amplitude trigger of the instability scales with viscosity and the size of the domain. In this talk, we summarize our findings. We present a map, in terms of the experimentally controllable parameters, that shows where the instability occurs and whether the instability creates a few isolated transient vortices, a few long-lived vortices, or long-lived, self-replicating vortices that fill the entire flow.

  18. Mountains and Plateaus on Io

    NASA Technical Reports Server (NTRS)

    1997-01-01

    These two views of Io were acquired by NASA's Galileo spacecraft during its seventh orbit (G7) of Jupiter. The images were designed to view large features on Io at low sun angles when the lighting conditions emphasize the topography or relief of the volcanic satellite. Sun angles are low near the terminator which is the day-night boundary near the left side of the images. These images reveal that the topography is very flat near the active volcanic centers such as Loki Patera (the large dark horseshoe-shaped feature near the terminator in the left-hand image) and that a variety of mountains and plateaus exist elsewhere.

    North is to the top of the picture. The resolution is about 6 kilometers per picture element (6.1 for the left hand image and 5.7 for the right). The images were taken on April 4th, 1997 at a ranges of 600,000 kilometers (left image) and 563,000 kilometers (right image) by the solid state imaging (CCD) system on NASA's Galileo spacecraft.

    The Jet Propulsion Laboratory, Pasadena, CA manages the mission for NASA's Office of Space Science, Washington, DC.

    This image and other images and data received from Galileo are posted on the World Wide Web, on the Galileo mission home page at URL http://galileo.jpl.nasa.gov. Background information and educational context for the images can be found at URL http://www.jpl.nasa.gov/galileo/sepo

  19. Numerical analysis of hydrodynamic instability in magnetized laser ablation flow

    NASA Astrophysics Data System (ADS)

    Ohnishi, Naofumi; Ishii, Ayako; Kuramitsu, Yasuhiro; Morita, Taichi; Sakawa, Youichi; Takabe, Hideaki

    2015-12-01

    We have conducted radiation magneto-hydrodynamics (RMHD) simulations of Richtmyer-Meshkov instability (RMI) in a magnetized counter flow produced by intense lasers. A jet-like plasma from a planar plastic target is formed and maintained in several-tens of nanoseconds by expanding plasma from rear side of two separated laser spots, and parallelly located another target is ablated by the radiation from the plasma, reproducing past experimental works. A planar shock driven by the radiation interacts with the jet as a nonuniform density structure, resulting in the RMI. The magnetic field is amplified up to ∼40 times greater than the background value at the interface at which the instability occurs. However, a certain extent of the amplification results from the compression effect induced by the counter flow, and the obtained amplification level is difficult to be measured in the experiments. An experiment for observing a clear amplification must be designed through the RMHD simulations so that the RMI takes place in the low-density area between two targets.

  20. Numerical simulations of multifrequency instability-wave growth and suppression in the Blasius boundary layer

    NASA Astrophysics Data System (ADS)

    Pal, A.; Bower, W. W.; Meyer, G. H.

    1991-02-01

    A mathematical model based on the Orr-Sommerfeld equation is developed to describe the growth and suppression of multifrequency, two-dimensional instability waves in the Blasius boundary layer over a flat place through localized perturbations at the surface caused by time-varying suction/blowing. It is shown for harmonic (single-frequency) perturbations that the instability wave can be decomposed into two components: an idealized Tollmien-Schlichting wave and a second perturbation that approximately cancels the first component upstream of the surface disturbance and becomes small downstream. Because the first component alone fully expresses the instability of the flow, the need to perform numerical Fourier transformation over the wave number is eliminated, permitting easy extension of the analysis to the more general case of arbitrary waveform of the perturbation. Numerical results are presented for examples of instability-wave generation and suppression in the boundary layer.

  1. Direct numerical simulations of the elliptic instability of a vortex pair

    NASA Astrophysics Data System (ADS)

    Laporte, F.; Corjon, A.

    2000-05-01

    The objective of this study is to perform direct numerical simulations (DNS) of the three-dimensional short-wavelength elliptic instability developing in a counter-rotating vortex pair, and to reproduce numerically a water-tank experiment. The main features of the elliptic instability are recovered by the simulations. In particular, the spatial structure and the temporal evolution of the most amplified perturbation mode during the linear regime correspond to both experimental measurements and theoretical predictions. The long-term evolution is also simulated, and the stages leading to transition to turbulence are described. Some elements resulting from simulations related to the interaction between the short-wavelength elliptic instability and the long-wavelength Crow instability are provided.

  2. Scale separation in granular packings: stress plateaus and fluctuations.

    PubMed

    Goldenberg, C; Atman, A P F; Claudin, P; Combe, G; Goldhirsch, I

    2006-04-28

    It is demonstrated, by numerical simulations of a 2D assembly of polydisperse disks, that there exists a range (plateau) of coarse-graining scales for which the stress tensor field in a granular solid is nearly resolution independent, thereby enabling an "objective" definition of this field. Expectedly, it is not the mere size of the system but the (related) magnitudes of the gradients that determine the widths of the plateaus. Ensemble averaging (even over "small" ensembles) extends the widths of the plateaus to subparticle scales. The fluctuations within the ensemble are studied as well. Both the response to homogeneous forcing and to an external compressive localized load (and gravity) are studied. Implications to small solid systems and constitutive relations are briefly discussed. PMID:16712277

  3. Termination of the MRI via parasitic instabilities in core-collapse supernovae: influence of numerical methods

    NASA Astrophysics Data System (ADS)

    Rembiasz, T.; Obergaulinger, M.; Cerdá-Durán, P.; Aloy, M. Á.; Müller, E.

    2016-05-01

    We study the influence of numerical methods and grid resolution on the termination of the magnetorotational instability (MRI) by means of parasitic instabilities in threedimensional shearing-disc simulations reproducing typical conditions found in core-collapse supernovae. Whether or not the MRI is able to amplify weak magnetic fields in this context strongly depends, among other factors, on the amplitude at which its growth terminates. The qualitative results of our study do not depend on the numerical scheme. In all our models, MRI termination is caused by Kelvin-Helmholtz instabilities, consistent with theoretical predictions. Quantitatively, however, there are differences, but numerical convergence can be achieved even at relatively low grid resolutions if high-order reconstruction methods are used.

  4. On physical and numerical instabilities arising in simulations of non-stationary radiatively cooling shocks

    NASA Astrophysics Data System (ADS)

    Badjin, D. A.; Glazyrin, S. I.; Manukovskiy, K. V.; Blinnikov, S. I.

    2016-06-01

    We describe our modelling of the radiatively cooling shocks and their thin shells with various numerical tools in different physical and calculational setups. We inspect structure of the dense shell, its formation and evolution, pointing out physical and numerical factors that sustain its shape and also may lead to instabilities. We have found that under certain physical conditions, the circular shaped shells show a strong bending instability and successive fragmentation on Cartesian grids soon after their formation, while remain almost unperturbed when simulated on polar meshes. We explain this by physical Rayleigh-Taylor-like instabilities triggered by corrugation of the dense shell surfaces by numerical noise. Conditions for these instabilities follow from both the shell structure itself and from episodes of transient acceleration during re-establishing of dynamical pressure balance after sudden radiative cooling onset. They are also easily excited by physical perturbations of the ambient medium. The widely mentioned non-linear thin shell instability, in contrast, in tests with physical perturbations is shown to have only limited chances to develop in real radiative shocks, as it seems to require a special spatial arrangement of fluctuations to be excited efficiently. The described phenomena also set new requirements on further simulations of the radiatively cooling shocks in order to be physically correct and free of numerical artefacts.

  5. On the numerical dispersion of electromagnetic particle-in-cell code: Finite grid instability

    SciTech Connect

    Meyers, M.D.; Huang, C.-K.; Zeng, Y.; Yi, S.A.; Albright, B.J.

    2015-09-15

    The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the Electromagnetic PIC model. We rigorously derive the faithful 3-D numerical dispersion relation of the PIC model, for a simple, direct current deposition scheme, which does not conserve electric charge exactly. We then specialize to the Yee FDTD scheme. In particular, we clarify the presence of alias modes in an eigenmode analysis of the PIC model, which combines both discrete and continuous variables. The manner in which the PIC model updates and samples the fields and distribution function, together with the temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme, is explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1-D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction, which is then verified by simulation. We demonstrate that our analysis is readily extendable to charge conserving models.

  6. On the numerical dispersion of electromagnetic particle-in-cell code: Finite grid instability

    NASA Astrophysics Data System (ADS)

    Meyers, M. D.; Huang, C.-K.; Zeng, Y.; Yi, S. A.; Albright, B. J.

    2015-09-01

    The Particle-In-Cell (PIC) method is widely used in relativistic particle beam and laser plasma modeling. However, the PIC method exhibits numerical instabilities that can render unphysical simulation results or even destroy the simulation. For electromagnetic relativistic beam and plasma modeling, the most relevant numerical instabilities are the finite grid instability and the numerical Cherenkov instability. We review the numerical dispersion relation of the Electromagnetic PIC model. We rigorously derive the faithful 3-D numerical dispersion relation of the PIC model, for a simple, direct current deposition scheme, which does not conserve electric charge exactly. We then specialize to the Yee FDTD scheme. In particular, we clarify the presence of alias modes in an eigenmode analysis of the PIC model, which combines both discrete and continuous variables. The manner in which the PIC model updates and samples the fields and distribution function, together with the temporal and spatial phase factors from solving Maxwell's equations on the Yee grid with the leapfrog scheme, is explicitly accounted for. Numerical solutions to the electrostatic-like modes in the 1-D dispersion relation for a cold drifting plasma are obtained for parameters of interest. In the succeeding analysis, we investigate how the finite grid instability arises from the interaction of the numerical modes admitted in the system and their aliases. The most significant interaction is due critically to the correct representation of the operators in the dispersion relation. We obtain a simple analytic expression for the peak growth rate due to this interaction, which is then verified by simulation. We demonstrate that our analysis is readily extendable to charge conserving models.

  7. Lithospheric Architecture, Heterogenities, Instabilities, Melting - insight form numerical modelling

    NASA Astrophysics Data System (ADS)

    Gorczyk, Weronika; Hobbs, Bruce; Ord, Alison; Gessner, Klaus; Gerya, Taras V.

    2010-05-01

    yield stress of the blocks (representing heterogeneous fused material) nucleates localised deformation and creates conditions for delamination via a Rayleigh-Taylor instability. Above the site of localised delamination of the mantle lithosphere, a series of deep crustal faults develop that may extend into the upper mantle. These deep structures can act as the pathways for mantle derived CO2±H2O fluids and alkaline igneous complexes. Isotherms are commonly elevated throughout the lithosphere in the hanging wall of deep through-going structures and are depressed in the footwalls. This means that some architectures favour devolatilisation and melting in the hanging wall. A large spectrum of behaviour is observed and results from minor changes in the orientation and strength of the blocks.

  8. A numerical simulation of barotropic instability. II Wave-wave interaction

    NASA Technical Reports Server (NTRS)

    Nielsen, J. E.; Schoeberl, M. R.

    1984-01-01

    A fully nonlinear numerical model of the point jet barotropic instability is used to test and confirm the hypothesis that the magnitude of the wave vorticity does not exceed the magnitude of the initial shear. This result arises directly from the local conservation of vorticity following a parcel and the fact that unstable waves are principally confined to the region where the zonal mean vorticity can be smoothed by the wave so as to eliminate the instability. Comparisons are made between fully nonlinear and quasi-linear models of the point jet instability and their tracer transport properties. Differences become particularly evident after wave saturation. The most important effect neglected by the wave-mean flow model appears to be the advection of wave vorticity by the most unstable mode. However, as equilibration of the instability proceeds, the globally averaged properties of both models are found to be similar.

  9. Electron cloud buildup and instability: Numerical simulations for the CERN Proton Synchrotron

    NASA Astrophysics Data System (ADS)

    Giovannozzi, M.; Métral, E.; Métral, G.; Rumolo, G.; Zimmermann, F.

    2003-01-01

    Experimental observations on the electron cloud have been collected at the CERN Proton Synchrotron machine throughout the last two years. At the same time, an intense campaign of simulations has been carried out to understand the observed electron cloud buildup and the related instability. In this paper, the results of the numerical simulations are presented and discussed.

  10. Numerical study of self modulation instability of 1 nC electron bunch at ATF

    SciTech Connect

    Fang Yun; Mori, Warren; Muggli, Patric

    2012-12-21

    The development of self-modulation instability (SMI) is investigated numerically for the 1 nC electron bunch available at Accelerator Test Facility (ATF) of Brookhaven National Laboratory (BNL). Possible experiment based on the simulation results is proposed. All the simulations are performed with the 2D-cylindrically symmetric particle-in-cell code.

  11. Linear and nonlinear evolution of the ion resonance instability in cylindrical traps: A numerical study

    SciTech Connect

    Sengupta, M.; Ganesh, R.

    2015-07-15

    Numerical experiments have been performed to investigate the linear and nonlinear dynamics, and energetics of the ion resonance instability in cylindrically confined nonneutral plasma. The instability is excited on a set of parametrically different unstable equilibria of a cylindrical nonneutral cloud, composed of electrons partially neutralized by a much heavier ion species of single ionization. A particle-in-cell code has been developed and employed to carry out these simulations. The results obtained from the initial exponential growth phase of the instability in these numerical experiments are in agreement with the linearised analytical model of the ion resonance instability. As the simulations delve much further in time beyond the exponential growth phase, very interesting nonlinear phenomena of the ion resonance instability are revealed, such as a process of simultaneous wave breaking of the excited poloidal mode on the ion cloud and pinching of the poloidal perturbations on the electron cloud. This simultaneous nonlinear dynamics of the two components is associated with an energy transfer process from the electrons to the ions. At later stages there is heating induced cross-field transport of the heavier ions and tearing across the pinches on the electron cloud followed by an inverse cascade of the torn sections.

  12. Numerical study on Rayleigh-Taylor instabilities in the lightning return stroke

    SciTech Connect

    Chen, Qiang; Chen, Bin Shi, Lihua; Yi, Yun; Wang, Yangyang

    2015-09-15

    The Rayleigh-Taylor (R-T) instabilities are important hydrodynamics and magnetohydrodynamics (MHD) phenomena that are found in systems in high energy density physics and normal fluids. The formation and evolution of the R-T instability at channel boundary during back-flow of the lightning return stroke are analyzed using the linear perturbation theory and normal mode analysis methods, and the linear growth rate of the R-T instability in typical condition for lightning return stroke channel is obtained. Then, the R-T instability phenomena of lightning return stroke are simulated using a two-dimensional Eulerian finite volumes resistive radiation MHD code. The numerical results show that the evolution characteristics of the R-T instability in the early stage of back-flow are consistent with theoretical predictions obtained by linear analysis. The simulation also yields more evolution characteristics for the R-T instability beyond the linear theory. The results of this work apply to some observed features of the return stroke channel and further advance previous theoretical and experimental work.

  13. High order numerical simulations of the Richtmyer- Meshkov instability in a relativistic fluid

    NASA Astrophysics Data System (ADS)

    Zanotti, O.; Dumbser, M.

    2015-07-01

    We study the Richtmyer-Meshkov (RM) instability of a relativistic perfect fluid by means of high order numerical simulations with adaptive mesh refinement (AMR). The numerical scheme combines a finite volume reconstruction in space, a local space-time discontinuous Galerkin predictor method, a high order one-step time update scheme, and a "cell-by-cell" space-time AMR strategy with time-accurate local time stepping. In this way, third order accurate (both in space and in time) numerical simulations of the RM instability are performed, spanning a wide parameter space. We present results both for the case in which a light fluid penetrates into a higher density one (Atwood number A > 0) and for the case in which a heavy fluid penetrates into a lower density one (Atwood number A < 0). We find that for large Lorentz factors γs of the incident shock wave, the relativistic RM instability is substantially weakened and ultimately suppressed. More specifically, the growth rate of the RM instability in the linear phase has a local maximum which occurs at a critical value of γs ≈ [1.2, 2]. Moreover, we have also revealed a genuinely relativistic effect, absent in Newtonian hydrodynamics, which arises in three dimensional configurations with a non-zero velocity component tangent to the incident shock front. In particular, in A > 0 models, the tangential velocity has a net magnification effect, while in A < 0 models, the tangential velocity has a net suppression effect.

  14. A semi-numerical algorithm for instability of compressible multilayered structures

    NASA Astrophysics Data System (ADS)

    Tang, Shan; Yang, Yang; Peng, Xiang He; Liu, Wing Kam; Huang, Xiao Xu; Elkhodary, Khalil

    2015-07-01

    A computational method is proposed for the analysis and prediction of instability (wrinkling or necking) of multilayered compressible plates and sheets made by metals or polymers under plane strain conditions. In previous works, a basic assumption (or a physical argument) that has been frequently made is that materials are incompressible to simplify mathematical derivations. To account for the compressibility of metals and polymers (the lower Poisson's ratio leads to the more compressible material), we propose a combined semi-numerical algorithm and finite element method for instability analysis. Our proposed algorithm is herein verified by comparing its predictions with published results in literature for thin films with polymer/metal substrates and for polymer/metal systems. The new combined method is then used to predict the effects of compressibility on instability behaviors. Results suggest potential utility for compressibility in the design of multilayered structures.

  15. Numerical investigation of the instability and nonlinear evolution of narrow-band directional ocean waves.

    PubMed

    Eliasson, Bengt; Shukla, P K

    2010-07-01

    The instability and nonlinear evolution of directional ocean waves is investigated numerically by means of simulations of the governing kinetic equation for narrow-band surface waves. Our simulation results reveal the onset of the modulational instability for long-crested wave trains, which agrees well with recent large-scale experiments in wave basins, where it was found that narrower directional spectra lead to self-focusing of ocean waves and an enhanced probability of extreme events. We find that the modulational instability is nonlinearly saturated by a broadening of the wave spectrum, which leads to the stabilization of the water-wave system. Applications of our results to other fields of physics, such as nonlinear optics and plasma physics, are discussed. PMID:20867450

  16. Shear-induced instability and arch filament eruption - A magnetohydrodynamic (MHD) numerical simulation

    NASA Technical Reports Server (NTRS)

    Wu, S. T.; Song, M. T.; Martens, P. C. H.; Dryer, M.

    1991-01-01

    A situation wherein a bipolar magnetic field embedded in a stratified solar atmosphere undergoes symmetrical shear motion at the footpoints is investigated via a 2D (nonplanar) MHD simulation. It was found that the vertical plasma flow velocities grow exponentially, leading to a new type of global MHD instability. The growth rate increases almost linearly until it reaches the same order of magnitude as the Alfven speed. Then a nonlinear MHD instability occurs beyond this point. It was found that the central loops are pinched by opposing Lorentz forces, and the outer closed loops stretch upward with the vertically-rising mass flow. The nonlinear dynamical shearing instability is illustrated by a numerical example that is given for three different values of the plasma beta that span several orders of magnitude.

  17. Experimental demonstration of bow-shock instability and its numerical analysis

    NASA Astrophysics Data System (ADS)

    Kikuchi, Y.; Ohnishi, N.; Ohtani, K.

    2016-07-01

    An experimental demonstration was carried out in a ballistic range at high Mach numbers with the low specific heat ratio gas hydrofluorocarbon HFC-134a to observe the unstable bow-shock wave generated in front of supersonic blunt objects. The shadowgraph images obtained from the experiments showed instability characteristics, in which the disturbances grow and flow downstream and the wake flow appears wavy because of the shock oscillation. Moreover, the influence of the body shape and specific heat ratio on the instability was investigated for various experimental conditions. Furthermore, the observed features, such as wave structure and disturbance amplitude, were captured by numerical simulations, and it was demonstrated that computational fluid dynamics could effectively simulate the physical instability. In addition, it was deduced that the shock instability is induced by sound emissions from the edge of the object. This inference supports the dependence of the instability on the specific heat ratio and Mach number because the shock stand-off distance is affected by these parameters and limits the sound wave propagation.

  18. System dynamic instabilities induced by sliding contact: A numerical analysis with experimental validation

    NASA Astrophysics Data System (ADS)

    Brunetti, J.; Massi, F.; Saulot, A.; Renouf, M.; D`Ambrogio, W.

    2015-06-01

    Mechanical systems present several contact surfaces between deformable bodies. The contact interface can be either static (joints) or in sliding (active interfaces). The sliding interfaces can have several roles and according to their application they can be developed either for maximizing the friction coefficient and the energy dissipation (e.g. brakes) or rather to allow the relative displacement at joints with a maximum efficiency. In both cases the coupling between system and local contact dynamics can bring to system dynamics instabilities (e.g. brake squeal or squeaking of hip prostheses). This results in unstable vibrations of the system, induced by the oscillation of the contact forces. In the literature, a large number of works deal with such kind of instabilities and are mainly focused on applied problems such as brake squeal noise. This paper shows a more general numerical analysis of a simple system constituted by two bodies in sliding contact: a rigid cylinder rotating inside a deformable one. The parametrical Complex Eigenvalue Analysis and the transient numerical simulations show how the friction forces can give rise to in-plane dynamic instabilities due to the interaction between two system modes, even for such a simple system characterized by one deformable body. Results from transient simulations highlight the key role of realistic values of the material damping to have convergence of the model and, consequently, reliable physical results. To this aim an experimental estimation of the material damping has been carried out. Moreover, the simplicity of the system allows for a deeper analysis of the contact instability and a balance of the energy flux among friction, system vibrations and damping. The numerical results have been validated by comparison with experimental ones, obtained by a specific test bench developed to reproduce and analyze the contact friction instabilities.

  19. Experimental and numerical investigations on flashing-induced instabilities in a single channel

    SciTech Connect

    Marcel, Christian P.; Rohde, M.; Van Der Hagen, T.H.J.J.

    2009-11-15

    During the start-up phase, natural circulation BWRs (NC-BWRs) need to be operated at low pressure conditions. Such conditions favor flashing-induced instabilities due to the large hydrostatic pressure drop induced by the tall chimney. Moreover, in novel NC-BWR designs the steam separation is performed in the steam separators which create large pressure drops at the chimney outlet, which effect on stability has not been investigated yet. In this work, flashing-induced oscillations occurring in a tall, bottom heated channel are numerically investigated by using a simple linear model with three regions and an accurate implementation for estimating the water properties. The model is used to investigate flashing-induced instabilities in a channel for different values of the core inlet friction value. The results are compared with experiments obtained by using the CIRCUS facility at the same conditions, showing a good agreement. In addition, the experiments on flashing-induced instabilities are presented in a novel manner allowing visualizing new details of the phenomenon numerical stability investigations on the effect of the friction distribution are also done. It is found that by increasing the total restriction in the channel the system is destabilized. In addition, the chimney outlet restriction has a stronger destabilizing effect than the core inlet restriction. A stable two-phase region is observed prior to the instabilities in the experiments and the numerical simulations which may help to pressurize the vessel of NC-BWRs and thus reducing the effects of flashing instabilities during start-up. (author)

  20. A numerical model for elliptical instability of the Earth's fluid outer core

    NASA Astrophysics Data System (ADS)

    Seyed-Mahmoud, Behnam; Henderson, Gary; Aldridge, Keith

    2000-01-01

    A dynamical model is proposed for the elliptical instability that has been reported by Aldridge et al. [Aldridge, K.D., Seyed-Mahmoud, B., Henderson, G.A., van Wijngaarden, W., 1997. Elliptical instability of the Earth's fluid core. Phys. Earth Planet. Inter., 103, 365-374] in connection with recent experiments on an ellipsoidal shell of rotating fluid. The frequencies and growth rates of the instability are obtained numerically by means of a Galerkin method that is based upon the normal modes of the contained fluid. A finite-element method has been employed to approximately solve the ill-posed Poincaré problem for the normal modes. The numerical results for a special case are compared with their analytical counterparts, and the agreement is to within 0.1% for shells of small ellipticity. Results are presented for other cases, including some where the boundary perturbation is allowed to rotate slowly with respect to the inertial frame. The conclusion is that such investigations are of geophysical interest, since tidal forcing might be sufficient to excite an elliptical instability of the fluid outer core of the Earth and thus contribute to the geomagnetic field.

  1. Instability of surface lenticular vortices: results from laboratory experiments and numerical simulations

    NASA Astrophysics Data System (ADS)

    Lahaye, Noé; Paci, Alexandre; Smith, Stefan Llewellyn

    2016-04-01

    We examine the instability of lenticular vortices -- or lenses -- in a stratified rotating fluid. The simplest configuration is one in which the lenses overlay a deep layer and have a free surface, and this can be studied using a two-layer rotating shallow water model. We report results from laboratory experiments and high-resolution direct numerical simulations of the destabilization of vortices with constant potential vorticity, and compare these to a linear stability analysis. The stability properties of the system are governed by two parameters: the typical upper-layer potential vorticity and the size (depth) of the vortex. Good agreement is found between analytical, numerical and experimental results for the growth rate and wavenumber of the instability. The nonlinear saturation of the instability is associated with conversion from potential to kinetic energy and weak emission of gravity waves, giving rise to the formation of coherent vortex multipoles with trapped waves. The impact of flow in the lower layer is examined. In particular, it is shown that the growth rate can be strongly affected and the instability can be suppressed for certain types of weak co-rotating flow.

  2. Numerical techniques for solving nonlinear instability problems in smokeless tactical solid rocket motors. [finite difference technique

    NASA Technical Reports Server (NTRS)

    Baum, J. D.; Levine, J. N.

    1980-01-01

    The selection of a satisfactory numerical method for calculating the propagation of steep fronted shock life waveforms in a solid rocket motor combustion chamber is discussed. A number of different numerical schemes were evaluated by comparing the results obtained for three problems: the shock tube problems; the linear wave equation, and nonlinear wave propagation in a closed tube. The most promising method--a combination of the Lax-Wendroff, Hybrid and Artificial Compression techniques, was incorporated into an existing nonlinear instability program. The capability of the modified program to treat steep fronted wave instabilities in low smoke tactical motors was verified by solving a number of motor test cases with disturbance amplitudes as high as 80% of the mean pressure.

  3. Numerical simulations of magnetic Kelvin-Helmholtz instability at a twisted solar flux tube

    NASA Astrophysics Data System (ADS)

    Murawski, K.; Chmielewski, P.; Zaqarashvili, T. V.; Khomenko, E.

    2016-04-01

    The paper aims to study the response of a solar small-scale and weak magnetic flux tube to photospheric twisting motions. We numerically solve three-dimensional ideal magnetohydrodynamic equations to describe the evolution of the perturbation within the initially static flux tube, excited by twists in the azimuthal component of the velocity. These twists produce rotation of the magnetic field lines. Perturbation of magnetic field lines propagates upwardly, driving vertical and azimuthal flow as well as plasma compressions and rarefactions in the form of eddies. We conclude that these eddies result from the sheared azimuthal flow which seeds Kelvin-Helmholtz instability (KHI) between the flux tube and the ambient medium. Numerically obtained properties of the KHI confirm the analytical predictions for the occurrence of the instability.

  4. Numerical simulations of magnetic Kelvin-Helmholtz instability at a twisted solar flux tube

    NASA Astrophysics Data System (ADS)

    Murawski, K.; Chmielewski, P.; Zaqarashvili, T. V.; Khomenko, E.

    2016-07-01

    The paper aims to study the response of a solar small-scale and weak magnetic flux tube to photospheric twisting motions. We numerically solve three-dimensional ideal magnetohydrodynamic equations to describe the evolution of the perturbation within the initially static flux tube, excited by twists in the azimuthal component of the velocity. These twists produce rotation of the magnetic field lines. Perturbation of magnetic field lines propagates upwardly, driving vertical and azimuthal flow as well as plasma compressions and rarefactions in the form of eddies. We conclude that these eddies result from the sheared azimuthal flow which seeds Kelvin-Helmholtz instability (KHI) between the flux tube and the ambient medium. Numerically obtained properties of the KHI confirm the analytical predictions for the occurrence of the instability.

  5. Three-dimensional numerical studies of the temperature anisotropy instability in intense charged particle beams

    NASA Astrophysics Data System (ADS)

    Startsev, Edward A.; Davidson, Ronald C.; Qin, Hong

    2005-05-01

    In neutral plasmas with a uniform magnetic field and strongly anisotropic distribution function (T∥/T⊥≪1) an electrostatic Harris-type collective instability may develop if the plasma is sufficiently dense. Such anisotropies develop naturally in accelerators, and a similar instability may lead to a deterioration of the beam quality in a one-component nonneutral charged particle beam. The instability may also lead to an increase in the longitudinal velocity spread, which would make the focusing of the beam difficult and impose a limit on the minimum spot size achievable in heavy ion fusion experiments. This paper reports the results of recent numerical studies of the temperature anisotropy instability using the newly developed Beam Eigenmodes And Spectra (bEASt) code for space-charge-dominated, low-emittance beams with large tune depression (ν/ν0≪1). Such high-intensity beams are relevant to next-step experiments such as the Integrated Beam Experiment (IBX), which would serve as proof-of-principal experiment for heavy-ion fusion.

  6. Numerical Analysis of Standing Accretion Shock Instability with Neutrino Heating in Supernova Cores

    NASA Astrophysics Data System (ADS)

    Ohnishi, Naofumi; Kotake, Kei; Yamada, Shoichi

    2006-04-01

    We have numerically studied the instability of the spherically symmetric standing accretion shock wave against nonspherical perturbations. We have in mind the application to collapse-driven supernovae in the postbounce phase, where the prompt shock wave generated by core bounce is commonly stalled. We take an experimental standpoint in this paper. Using spherically symmetric, completely steady, shocked accretion flows as unperturbed states, we have clearly observed both the linear growth and the subsequent nonlinear saturation of the instability. In so doing, we have employed a realistic equation of state, together with heating and cooling via neutrino reactions with nucleons. We have performed a mode analysis based on the spherical harmonics decomposition and found that the modes with l=1,2 are dominant not only in the linear regime but also after nonlinear couplings generate various modes and saturation occurs. By varying the neutrino luminosity, we have constructed unperturbed states both with and without a negative entropy gradient. We have found that in both cases the growth of the instability is similar, suggesting that convection does not play a dominant role, which also appears to be supported by the recent linear analysis of the convection in accretion flows by Foglizzo et al. The oscillation period of the unstable l=1 mode is found to fit better with the advection time rather than with the sound crossing time. Whatever the cause may be, the instability favors a shock revival.

  7. Direct Numerical Simulations of Small-Scale Gravity Wave Instability Dynamics in Variable Stratification and Shear

    NASA Astrophysics Data System (ADS)

    Mixa, T.; Fritts, D. C.; Laughman, B.; Wang, L.; Kantha, L. H.

    2015-12-01

    Multiple observations provide compelling evidence that gravity wave dissipation events often occur in multi-scale environments having highly-structured wind and stability profiles extending from the stable boundary layer into the mesosphere and lower thermosphere. Such events tend to be highly localized and thus yield local energy and momentum deposition and efficient secondary gravity wave generation expected to have strong influences at higher altitudes [e.g., Fritts et al., 2013; Baumgarten and Fritts, 2014]. Lidars, radars, and airglow imagers typically cannot achieve the spatial resolution needed to fully quantify these small-scale instability dynamics. Hence, we employ high-resolution modeling to explore these dynamics in representative environments. Specifically, we describe numerical studies of gravity wave packets impinging on a sheet of high stratification and shear and the resulting instabilities and impacts on the gravity wave amplitude and momentum flux for various flow and gravity wave parameters. References: Baumgarten, Gerd, and David C. Fritts (2014). Quantifying Kelvin-Helmholtz instability dynamics observed in noctilucent clouds: 1. Methods and observations. Journal of Geophysical Research: Atmospheres, 119.15, 9324-9337. Fritts, D. C., Wang, L., & Werne, J. A. (2013). Gravity wave-fine structure interactions. Part I: Influences of fine structure form and orientation on flow evolution and instability. Journal of the Atmospheric Sciences, 70(12), 3710-3734.

  8. Experimental and numerical study of the shear layer instability between two counter-rotating disks

    NASA Astrophysics Data System (ADS)

    Moisy, F.; Doaré, O.; Pasutto, T.; Daube, O.; Rabaud, M.

    2004-05-01

    The shear layer instability in the flow between two counter-rotating disks enclosed by a cylinder is investigated experimentally and numerically, for radius-to-height ratio Gamma {=} R/h between 2 and 21. For sufficiently large rotation ratio, the internal shear layer that separates two regions of opposite azimuthal velocities is prone to an azimuthal symmetry breaking, which is investigated experimentally by means of visualization and particle image velocimetry. The associated pattern is a combination of a sharp-cornered polygonal pattern, as observed by Lopez et al. (2002) for low aspect ratio, surrounded by a set of spiral arms, first described by Gauthier et al. (2002) for high aspect ratio. The spiral arms result from the interaction of the shear layer instability with the Ekman boundary layer over the faster rotating disk. Stability curves and critical modes are experimentally measured for the whole range of aspect ratios, and are found to compare well with numerical simulations of the three-dimensional time-dependent Navier Stokes equations over an extensive range of parameters. Measurements of a local Reynolds number based on the shear layer thickness confirm that a shear layer instability, with only weak curvature effect, is responsible for the observed patterns. This scenario is supported by the observed onset modes, which scale as the shear layer radius, and by the measured phase velocities.

  9. How common are ALS plateaus and reversals?

    PubMed Central

    Vaughan, Timothy; Wicks, Paul; Heywood, Jamie; Sinani, Ervin; Selsov, Roger; Macklin, Eric A.; Schoenfeld, David; Cudkowicz, Merit; Sherman, Alex

    2016-01-01

    Objective: To determine the frequency of amyotrophic lateral sclerosis (ALS) plateaus and reversals in the Pooled Resource Open-Access ALS Clinical Trials (PRO-ACT) database. Methods: We analyzed Amyotrophic Lateral Sclerosis Functional Rating Scale (ALSFRS) and ALSFRS–revised (ALSFRS-R) data from PRO-ACT participants. The frequencies of participants experiencing plateaus (periods where scores did not change) were calculated over 6-, 12-, and 18-month epochs. The percentage of participants ever experiencing reversals (periods where scores improved) of different lengths were also calculated and plotted. Results: Over 6 months, 25% of 3,132 participants did not decline. Over 12 months, 16% of 2,105 participants did not decline. Over 18 months, 7% of 1,218 participants did not decline. Small ALS reversals were also common, especially over shorter follow-up intervals; 14% of 1,343 participants had a 180-day interval where their ALSFRS-R slope was greater than zero. Fewer than 1% of participants ever experienced improvements of 4 or more ALSFRS-R points lasting at least 12 months. Conclusion: ALS plateaus and small reversals are common, especially over brief intervals. In light of these data, stable disease, especially for a short period of time, should not be interpreted as an ALS treatment effect. Large sustained ALS reversals, on the other hand, are rare, potentially important, and warrant further study. PMID:26658909

  10. Numerical analysis of the Eckhaus instability in travelling-wave convection in binary mixtures.

    PubMed

    Mercader, I; Alonso, A; Batiste, O

    2004-11-01

    The Eckhaus stability boundaries of travelling periodic roll patterns arising in binary fluid convection is analysed using high-resolution numerical methods. We present results corresponding to three different values of the separation ratio used in experiments. Our results show that the subcritical branches of travelling waves bifurcating at the onset of convection suffer sideband instabilities that are restabilised further away in the branch. If this restabilisation is produced after the turning point of the travelling-wave branch, these waves do not become stable in a saddle node bifurcation as would have been the case in a smaller domain. In the regions of instability of the uniform travelling waves we expect to find either transitions between states of different wave number or modulated travelling waves arising in these bifurcations. PMID:15592771

  11. Numerical simulations of turbulent trapping in the weak beam-plasma instability

    NASA Technical Reports Server (NTRS)

    Theilhaber, K.; Laval, G.; Pesme, D.

    1987-01-01

    Numerical simulations of the weak beam-plasma instability were done in the turbulent regime where small-scale trapping is a dominant feature of the instability, a regime with behavior not predicted by quasi-linear theory. The results of the simulations were compared with those of a specific model of the turbulence, the so-called 'turbulent trapping' model, which gives precise formulas for the particle correlation functions, and predicts a growth rate well enhanced over the quasi-linear value. It was found that the model gives accurate predictions for the correlation functions. On the other hand, while growth rates were enhanced over the quasi-linear values, the enhancements observed are smaller than expected from the quantitative predictions of the model.

  12. Numerical study of baroclinic instability associated with thermobaric deep convection at high latitudes: Idealized cases

    NASA Astrophysics Data System (ADS)

    Akitomo, Kazunori

    2005-06-01

    Numerical experiments with a three-dimensional nonhydrostatic model in a rotating frame have been executed to investigate baroclinic instability associated with thermobaric deep convection in weakly stratified polar oceans and its role in the transport processes. The model ocean has a two-layered structure with the cold, fresh mixed layer overlying the warm, saline deep water cell, as in the Weddell Sea. In contrast with a scenario based on the linear equation of state, thermobaric overturning of the water column enhances the horizontal density gradient (baroclinicity) through nonlinearity of the equation of state. If temperature controls water density (TEM cases), baroclinicity is intensified at the bottom of the overturned layer while at the surface if salinity does (SAL cases). Such intensification causes further development of baroclinic instability or baroclinic destabilization and more effective vertical heat transport. In the post-overturning stage, on the other hand, surface cooling (convective motion) has two oppositely operating effects on baroclinic instability and the associated heat transport. One is that horizontal convergence due to convective motion enhances baroclinic instability in the surface layer, as in previous studies focusing on strongly stratified oceans. This is observed in SAL cases with weak cooling, but not in TEM cases. The other is that strong cooling suppresses baroclinic instability by homogenizing the overturned layer vertically. This effect has not been found in the strongly stratified oceans. As a result, the vertical heat transport is most effective at low cooling rates (˜125Wm-2) in SAL cases while it monotonically decreases with cooling rate in TEM cases. When baroclinicity is initially weak as in the Weddell Sea, the most effective transport occurs with the cooling rate of 25Wm-2 which is a possible value under sea-ice cover in the actual situation.

  13. Elimination of the numerical Cerenkov instability for spectral EM-PIC codes

    NASA Astrophysics Data System (ADS)

    Yu, Peicheng; Xu, Xinlu; Decyk, Viktor K.; Fiuza, Frederico; Vieira, Jorge; Tsung, Frank S.; Fonseca, Ricardo A.; Lu, Wei; Silva, Luis O.; Mori, Warren B.

    2015-07-01

    When using an electromagnetic particle-in-cell (EM-PIC) code to simulate a relativistically drifting plasma, a violent numerical instability known as the numerical Cerenkov instability (NCI) occurs. The NCI is due to the unphysical coupling of electromagnetic waves on a grid to wave-particle resonances, including aliased resonances, i.e., ω + 2 πμ / Δt =(k1 + 2 πν1 / Δx1) v0, where μ and ν1 refer to the time and space aliases and the plasma is drifting relativistically at velocity v0 in the 1 ˆ -direction. We extend our previous work Xu et al. (2013) by recasting the numerical dispersion relation of a relativistically drifting plasma into a form which shows explicitly how the instability results from the coupling modes which are purely transverse electromagnetic (EM) modes and purely longitudinal modes in the rest frame of the plasma for each time and space aliasing. The dispersion relation for each μ and ν1 is the product of the dispersion relation of these two modes set equal to a coupling term that vanishes in the continuous limit. The new form of the numerical dispersion relation provides an accurate method of systematically calculating the growth rate and location of the mode in the fundamental Brillouin zone for any Maxwell solver for each μ and ν1. We then focus on the spectral Maxwell solver and systematically discuss its NCI modes. We show that the second fastest growing NCI mode for the spectral solver corresponds to μ =ν1 = 0, that it has a growth rate approximately one order of magnitude smaller than the fastest growing μ = 0 and ν1 = 1 mode, and that its location in the k space fundamental Brillouin zone is sensitive to the grid size and time step. Based on these studies, strategies to systematically eliminate the NCI modes for a spectral solver are developed. We apply these strategies to both relativistic collisionless shock and LWFA simulations, and demonstrate that high-fidelity multi-dimensional simulations of drifting plasmas

  14. Fluid Instabilities in the Crab Nebula Jet: Results from Numerical Simulations

    NASA Astrophysics Data System (ADS)

    Mignone, A.; Striani, E.; Bodo, G.; Anjiri, M.

    2014-09-01

    We present an overview of high-resolution relativistic MHD numerical simulations of the Crab Nebula South-East jet. The models are based on hot and relativistic hollow outflows initially carrying a purely toroidal magnetic field. Our results indicate that weakly relativistic (γ˜ 2) and strongly magnetized jets are prone to kink instabilities leading to a noticeable deflection of the jet. These conclusions are in good agreement with the recent X-ray (Chandra) data of Crab Nebula South-East jet indicating a change in the direction of propagation on a time scale of the order of few years.

  15. Numerical Simulations Studies of the Convective Instability Onset in a Supercritical Fluid

    NASA Technical Reports Server (NTRS)

    Furukawa, A.; Meyer, H.; Onuki, A.

    2004-01-01

    Numerical simulation studies are reported for the convection of a supercritical fluid, He-3, in a Rayleigh-Benard cell. The calculations provide the temporal profile DeltaT(t) of the temperature drop across the fluid layer. In a previous article, systematic delays in the onset of the convective instability in simulations relative to experiments were reported, as seen from the DeltaT(t) profiles. They were attributed to the smallness of the noise which is needed to start the instability. Therefore i) homogeneous temperature noise and ii) spatial lateral periodic temperature variations in the top plate were programmed into the simulations, and DeltaT(t) compared with that of an experiment with the same fluid parameters. An effective speed-up in the instability onset was obtained, with the best results obtained through the spatial temperature variations with a period of 2L, close to the wavelength of a pair of convections rolls. For a small amplitude of 0.5 micro-K, this perturbation gave a semiquantitative agreement with experimental observations. Results for various noise amplitudes are presented and discussed in relation to predictions by El Khouri and Carl es.

  16. Numerical simulations of the two-dimensional multimode Richtmyer-Meshkov instability

    SciTech Connect

    Thornber, B.; Zhou, Y.

    2015-03-15

    The two-dimensional Richtmyer-Meshkov instability occurs as shock waves pass through a perturbed material interface, triggering transition to an inhomogeneous turbulence variable density flow. This paper presents a series of large-eddy-simulations of the two dimensional turbulent RM instability and compares the results to the fully three dimensional simulations. There are two aims for this paper, the first is to explore what numerical resolution is required for a statistically converged solution for a two dimensional inhomogeneous flow field. The second aim is to elucidate the key differences in flow physics between the two dimensional and three dimensional Richtmyer-Meshkov instabilities, particularly their asymptotic self-similar regime. Convergence is achieved using 64 independent realisations and grid resolutions up to 4096{sup 2} in the plane. It is shown that for narrowband cases the growth rate θ = 0.48 which is substantially higher than the three-dimensional equivalent. Mix measures are consistently lower compared to three-dimensional, and the kinetic energy distribution is homogeneous at late time. The broadband case has a similar initial growth rate as the three-dimensional case, with a marginally lower θ = 0.63. Mix is similar in magnitude, but is reducing at late time. The spectra in both cases exhibit the dual-cascade expected from two-dimensional turbulence.

  17. Analytical and numerical treatment of resistive drift instability in a plasma slab

    NASA Astrophysics Data System (ADS)

    Mirnov, V. V.; Sauppe, J. P.; Hegna, C. C.; Sovinec, C. R.

    2016-05-01

    An analytic approach combining the effect of equilibrium diamagnetic flows and the finite ionsound gyroradius associated with electron-ion decoupling and kinetic Alfvén wave dispersion is derived to study resistive drift instabilities in a plasma slab. Linear numerical computations using the NIMROD code are performed with cold ions and hot electrons in a plasma slab with a doubly periodic box bounded by two perfectly conducting walls. A linearly unstable resistive drift mode is observed in computations with a growth rate that is consistent with the analytic dispersion relation. The resistive drift mode is expected to be suppressed by magnetic shear in unbounded domains, but the mode is observed in numerical computations with and without magnetic shear. In the slab model, the finite slab thickness and the perfectly conducting boundary conditions are likely to account for the lack of suppression.

  18. Numerical analysis of instability processes in underground cavities and of the related effects at the surface

    NASA Astrophysics Data System (ADS)

    Lollino, Piernicola; Parise, Mario

    2010-05-01

    Natural and anthropogenic caves may represent a potential hazard for the built-up environment, due to the occurrence of underground instability processes, that may propagate upward and eventually reach the ground surface, thus inducing the occurrence of sinkholes. Especially when the caves are at shallow depth, the effects at the ground surface may result extremely severe. In the Apulia region of southern Italy, there are many sites where underground quarrying developed in the past, due to presence at a certain depth of rock of good quality for building purposes. Development of underground quarries, rather than open pit mines, was also favoured by the preservation of the terrains on the ground surface for agricultural practices. The Pliocene-Pleistocene calcarenite (a typical soft rock) was therefore quarried underground, by digging extensive networks of galleries in those levels within the local geological succession most suitable for the quarrying activity. With time, these underground activities have progressively been abandoned, and later on many quarries were used for other purposes, including illegal discharge of solid and liquid wastes. Many Apulian towns are nowadays located just above these caves, due to urban expansion in the last decades and loss of memory of the presence of the underground quarries. Thus, a serious risk exists for civil society, which should not be left uninvestigated. The present contribution deals with the analysis of the main factors at the origin of the instability processes described, also including those causing weathering of the soft rock wihich induces gradual decay of the physical and mechanical properties of the rock mass. Aimed at exploring the evolution with time of the stability conditions within the cavities, numerical analysis have been implemented by using finite element methods with respect to ideal situations which are representative of typical case studies in Apulia. Both the effects of local instability processes

  19. Numerical investigation of a single-mode chemically reacting Richtmyer-Meshkov instability

    NASA Astrophysics Data System (ADS)

    Attal, N.; Ramaprabhu, P.

    2015-07-01

    We report on high-resolution, numerical simulations of a single-mode, chemically reacting, Richtmyer-Meshkov (RM) instability, at different interface thicknesses. The gases on either side of the diffuse interface were Hydrogen (H and Oxygen (O, with a pre-shock Atwood number of 0.5. An incident shock with a Mach number of 1.2 is allowed to traverse from the light (H to the heavy (O medium in the 2D numerical shock tube. The simulations were performed using the astrophysical FLASH code developed at the University of Chicago, with extensive modifications implemented by the authors to describe detailed H-O chemistry, temperature-dependent specific heats, and multi-species equation of state. The interface thickness was systematically varied in the simulations to study the effect of the total mass of fuel burnt and heat added on the hydrodynamic instability growth rates. In the absence of an incident shock, burning results in the formation of so-called combustion waves, which spontaneously trigger RM and Rayleigh-Taylor like instability growth of the interface. We are able to obtain the resulting growth rates of an imposed sinusoidal perturbation, and compare them with the predictions of an impulsive model, with simple modifications to account for the finite thickness of the interface, density changes due to heat addition, and compression of the material line due to the combustion wave. When additionally an incident shock is present, we observe complex interactions between the shock and the aforementioned combustion waves, resulting in significant non-planar distortions of each. When the unstable interface is subjected to a reshock, significant mixing enhancement is observed, accompanied by a dramatic increase in combustion product formation, and combustion efficiency.

  20. Numerical Simulation of Liquid Sheet Instability in a Multiphase Flow Domain

    NASA Astrophysics Data System (ADS)

    Souvick, Chatterjee; Mahapatra, Soumik; Mukhopadhyay, Achintya; Sen, Swarnendu

    2013-11-01

    Instability of a liquid sheet leading to the formation of droplets is a classical problem finding a wide range of multi-scale applications like gas turbine engines and inkjet printers. Numerical simulation of such a phenomenon is crucial because of its cost and time effective nature. In this work, the hydrodynamics in a custom designed nozzle is analyzed using Volume of Fluid method in Ansys Fluent. This innovative nozzle design includes an annular liquid sheet sandwiched between two air streams such that the inner air channel is recessed to a certain length. Such a recession leads to interaction between the two multiphase streams inside the atomizer resulting to an increased shear layer instability which augments the disintegration process. The numerical technique employed in this work couples Navier Stokes equation with VoF surface tracking technique. A parametric study with the hydrodynamic parameters involved in the problem, as well as the recession length, is performed while monitoring the axial and tangential exit velocities along with the spray cone angle. Comparison between the full 3D model and two different equivalent 2D axisymmetric models have been shown. The two axisymmetric models vary based on conserving different physical parameters between the 2D and 3D cases.

  1. Compressible magnetic Rayleigh-Taylor instability in stratified plasmas: Comparison of analytical and numerical results in the linear regime

    SciTech Connect

    Liberatore, S.; Jaouen, S.; Tabakhoff, E.; Canaud, B.

    2009-04-15

    Magnetic Rayleigh-Taylor instability is addressed in compressible hydrostatic media. A full model is presented and compared to numerical results from a linear perturbation code. A perfect agreement between both approaches is obtained in a wide range of parameters. Compressibility effects are examined and substantial deviations from classical Chandrasekhar growth rates are obtained and confirmed by the model and the numerical calculations.

  2. Direct numerical simulation of instabilities in parallel flow with spherical roughness elements

    NASA Technical Reports Server (NTRS)

    Deanna, R. G.

    1992-01-01

    Results from a direct numerical simulation of laminar flow over a flat surface with spherical roughness elements using a spectral-element method are given. The numerical simulation approximates roughness as a cellular pattern of identical spheres protruding from a smooth wall. Periodic boundary conditions on the domain's horizontal faces simulate an infinite array of roughness elements extending in the streamwise and spanwise directions, which implies the parallel-flow assumption, and results in a closed domain. A body force, designed to yield the horizontal Blasius velocity in the absence of roughness, sustains the flow. Instabilities above a critical Reynolds number reveal negligible oscillations in the recirculation regions behind each sphere and in the free stream, high-amplitude oscillations in the layer directly above the spheres, and a mean profile with an inflection point near the sphere's crest. The inflection point yields an unstable layer above the roughness (where U''(y) is less than 0) and a stable region within the roughness (where U''(y) is greater than 0). Evidently, the instability begins when the low-momentum or wake region behind an element, being the region most affected by disturbances (purely numerical in this case), goes unstable and moves. In compressible flow with periodic boundaries, this motion sends disturbances to all regions of the domain. In the unstable layer just above the inflection point, the disturbances grow while being carried downstream with a propagation speed equal to the local mean velocity; they do not grow amid the low energy region near the roughness patch. The most amplified disturbance eventually arrives at the next roughness element downstream, perturbing its wake and inducing a global response at a frequency governed by the streamwise spacing between spheres and the mean velocity of the most amplified layer.

  3. Direct numerical simulation of electrokinetic instability and transition to chaotic motion

    SciTech Connect

    Demekhin, E. A.; Nikitin, N. V.; Shelistov, V. S.

    2013-12-15

    A new type of instability—electrokinetic instability—and an unusual transition to chaotic motion near a charge-selective surface (semiselective electric membrane, electrode, or system of micro-/nanochannels) was studied by the numerical integration of the Nernst-Planck-Poisson-Stokes system and a weakly nonlinear analysis near the threshold of instability. A special finite-difference method was used for the space discretization along with a semi-implicit 31/3 -step Runge-Kutta scheme for the integration in time. Two kinds of initial conditions were considered: (a) white-noise initial conditions to mimic “room disturbances” and subsequent natural evolution of the solution, and (b) an artificial monochromatic ion distribution with a fixed wave number to simulate regular wave patterns. The results were studied from the viewpoint of hydrodynamic stability and bifurcation theory. The threshold of electroconvective movement was found by the linear spectral stability theory, the results of which were confirmed by numerical simulation of the entire system. Our weakly nonlinear analysis and numerical integration of the entire system predict possibility of both kinds of bifurcations at the critical point, supercritical and subcritical, depending on the system parameters. The following regimes, which replace each other as the potential drop between the selective surfaces increases, were obtained: one-dimensional steady solution, two-dimensional steady electroconvective vortices (stationary point in a proper phase space), unsteady vortices aperiodically changing their parameters (homoclinic contour), periodic motion (limit cycle), and chaotic motion. The transition to chaotic motion does not include Hopf bifurcation. The numerical resolution of the thin concentration polarization layer showed spike-like charge profiles along the surface, which could be, depending on the regime, either steady or aperiodically coalescent. The numerical investigation confirmed the

  4. Analytical and numerical study of the transverse Kelvin-Helmholtz instability in tokamak edge plasmas

    DOE PAGESBeta

    Myra, James R.; D'Ippolito, Daniel A.; Russell, David A.; Umansky, Maxim V.; Baver, Derek A.

    2016-04-11

    Sheared flows perpendicular to the magnetic field can be driven by the Reynolds stress or ion pressure gradient effects and can potentially influence the stability and turbulent saturation level of edge plasma modes. On the other hand, such flows are subject to the transverse Kelvin- Helmholtz (KH) instability. Here, the linear theory of KH instabilities is first addressed with an analytic model in the asymptotic limit of long wavelengths compared with the flow scale length. The analytic model treats sheared ExB flows, ion diamagnetism (including gyro-viscous terms), density gradients and parallel currents in a slab geometry, enabling a unified summarymore » that encompasses and extends previous results. In particular, while ion diamagnetism, density gradients and parallel currents each individually reduce KH growth rates, the combined effect of density and ion pressure gradients is more complicated and partially counteracting. Secondly, the important role of realistic toroidal geometry is explored numerically using an invariant scaling analysis together with the 2DX eigenvalue code to examine KH modes in both closed and open field line regions. For a typical spherical torus magnetic geometry, it is found that KH modes are more unstable at and just outside the separatrix as a result of the distribution of magnetic shear. Lastly implications for reduced edge turbulence modeling codes are discussed.« less

  5. Analytical and numerical study of the transverse Kelvin-Helmholtz instability in tokamak edge plasmas

    NASA Astrophysics Data System (ADS)

    Myra, J. R.; D’Ippolito, D. A.; Russell, D. A.; Umansky, M. V.; Baver, D. A.

    2016-04-01

    > Sheared flows perpendicular to the magnetic field can be driven by the Reynolds stress or ion pressure gradient effects and can potentially influence the stability and turbulent saturation level of edge plasma modes. On the other hand, such flows are subject to the transverse Kelvin-Helmholtz (KH) instability. Here, the linear theory of KH instabilities is first addressed with an analytic model in the asymptotic limit of long wavelengths compared with the flow scale length. The analytic model treats sheared flows, ion diamagnetism (including gyro-viscous terms), density gradients and parallel currents in a slab geometry, enabling a unified summary that encompasses and extends previous results. In particular, while ion diamagnetism, density gradients and parallel currents each individually reduce KH growth rates, the combined effect of density and ion pressure gradients is more complicated and partially counteracting. Secondly, the important role of realistic toroidal geometry is explored numerically using an invariant scaling analysis together with the 2DX eigenvalue code to examine KH modes in both closed and open field line regions. For a typical spherical torus magnetic geometry, it is found that KH modes are more unstable at, and just outside of, the separatrix as a result of the distribution of magnetic shear. Finally implications for reduced edge turbulence modelling codes are discussed.

  6. Analytical and numerical treatment of drift-tearing and resistive drift instabilities in plasma slab

    NASA Astrophysics Data System (ADS)

    Mirnov, V. V.; Hegna, C. C.; Sauppe, J. P.; Sovinec, C. R.

    2015-11-01

    We consider modification to linear resistive MHD instability theory in a slab due to two categories of non-MHD effects: (1) electron and ion diamagnetic flows caused by equilibrium pressure gradients and (2) electron and ion decoupling on short scales associated with kinetic Alfven and whistler waves. The relationship between the expected stabilizing response due to the effects (1) and the destabilizing contribution caused by the dispersive waves (2) is investigated. An analytic solution combining the effect of diamagnetic flows and the ion-sound gyroradius contribution is derived using a perturbative approach. Linear numerical simulations using the NIMROD code are performed with cold ions and hot electrons in plasma slab with a doubly periodic box bounded by two perfectly conducting walls. Configurations with magnetic shear are unstable to current-driven drift-tearing instability. A second linearly unstable resistive drift type mode with largely electrostatic perturbations is also observed in simulations. The resistive-drift mode is suppressed by magnetic shear in unbounded domains but can remain unstable in the simulations with finite slab thickness and perfectly conducting wall. Additionally, the growth rate is sensitive to the magnetic shear length. We analyze whether these modes can be unstable in cylindrical configurations with magnetic shear typical for reversed field pinches. The material is based on work supported by the U.S. DOE and NSF.

  7. Numerical Study of Electrolytic Flow Instabilities Driven by an Azimuthal Lorentz Force in a Cylindrical Geometry

    NASA Astrophysics Data System (ADS)

    Pérez-Barrera, James; Pérez-Espinoza, José Enrique; Ortíz, Alejandro; Cuevas, Sergio; Ramos, Eduardo

    2014-11-01

    We present numerical simulations of the flow produced by an azimuthal Lorentz force in an electromagnetic stirrer. The stirrer consists of a cylindrical cavity with two copper concentric cylindrical electrodes, filled with an electrolytic solution. Underneath the cavity, a permanent magnet creates an almost uniform magnetic field, perpendicular to the circular section of the stirrer. An electric potential difference between the electrodes produces a radial D.C. current that passes through the fluid and interacts with the axial magnetic field, generating an azimuthal Lorentz force that drives the fluid. Experiments have shown the appearance of a flow instability that gives rise to a varying number of anticyclonic vortices for given values of the current intensity and fluid layer thickness. The MHD governing equations are expressed in terms of the velocity, pressure and electric potential. Numerical simulations are carried out using a hybrid Finite volume-Fourier method to ensure periodicity in the azimuthal direction. Numerical results show the formation of different modes of perturbation in the velocity field, which give rise to a varying number of traveling vortical structures. Work supported by CONACYT, Mexico under Project 131399. JPB acknowledges a Grant from CONACYT.

  8. Numerical Experiments with a Turbulent Single-Mode Rayleigh-Taylor Instability

    SciTech Connect

    Cloutman, L.D.

    2000-04-01

    Direct numerical simulation is a powerful tool for studying turbulent flows. Unfortunately, it is also computationally expensive and often beyond the reach of the largest, fastest computers. Consequently, a variety of turbulence models have been devised to allow tractable and affordable simulations of averaged flow fields. Unfortunately, these present a variety of practical difficulties, including the incorporation of varying degrees of empiricism and phenomenology, which leads to a lack of universality. This unsatisfactory state of affairs has led to the speculation that one can avoid the expense and bother of using a turbulence model by relying on the grid and numerical diffusion of the computational fluid dynamics algorithm to introduce a spectral cutoff on the flow field and to provide dissipation at the grid scale, thereby mimicking two main effects of a large eddy simulation model. This paper shows numerical examples of a single-mode Rayleigh-Taylor instability in which this procedure produces questionable results. We then show a dramatic improvement when two simple subgrid-scale models are employed. This study also illustrates the extreme sensitivity to initial conditions that is a common feature of turbulent flows.

  9. Numerical analysis of gas-dynamic instabilities during the laser drilling process

    NASA Astrophysics Data System (ADS)

    Khan, A. H.; O'Neill, W.; Tunna, L.; Sutcliffe, C. J.

    2006-08-01

    The use of high-pressure gas jets in the laser-drilling process has significant influence on the melt ejection mechanism. These jets are highly unstable and this directly relates to the gas pressure and the geometry of the hole being drilled. The evolution of gas-dynamic instabilities during the laser-drilling process was investigated numerically. A minimum length nozzle (MLN) with a 300 μm throat diameter was modelled at various gas pressures, with the gas jet impinging on a range of simulated holes with different aspect ratios. The simulations predict the formation of surface pressure fluctuations that have a broad spectrum due to both the turbulent nature of the jet and the blunt shock oscillation on the surface. The surface pressure variations and the blunt shock oscillation govern the gas dynamic conditions inside the hole, which strongly influence the melt ejection phenomena during the laser-drilling process.

  10. Experimental and numerical study of plastic shear instability under high-speed loading conditions

    SciTech Connect

    Sokovikov, Mikhail E-mail: naimark@icmm.ru; Chudinov, Vasiliy E-mail: naimark@icmm.ru; Bilalov, Dmitry E-mail: naimark@icmm.ru; Oborin, Vladimir E-mail: naimark@icmm.ru; Uvarov, Sergey E-mail: naimark@icmm.ru; Plekhov, Oleg E-mail: naimark@icmm.ru; Terekhina, Alena E-mail: naimark@icmm.ru; Naimark, Oleg E-mail: naimark@icmm.ru

    2014-11-14

    The behavior of specimens dynamically loaded during the split Hopkinson (Kolsky) bar tests in a regime close to simple shear conditions was studied. The lateral surface of the specimens was investigated in a real-time mode with the aid of a high-speed infra-red camera CEDIP Silver 450M. The temperature field distribution obtained at different time made it possible to trace the evolution of plastic strain localization. The process of target perforation involving plug formation and ejection was examined using a high-speed infra-red camera and a VISAR velocity measurement system. The microstructure of tested specimens was analyzed using an optical interferometer-profilometer and a scanning electron microscope. The development of plastic shear instability regions has been simulated numerically.

  11. Direct Numerical Simulation of Nanofilm Instability Driven by Liquid/Solid Interactions

    NASA Astrophysics Data System (ADS)

    Mahady, Kyle; Afkhami, Shahriar; Kondic, Lou

    2015-11-01

    The nanoscale interaction between liquid and solid molecules underlies fundamental phenomena for systems involving liquids on surfaces. In addition to giving rise to the contact angle of drops, this interaction drives the spontaneous rupture of nanofilms. We study this process by means of direct simulation of the Navier-Stokes equations using the Volume of Fluid interface tracking method. Our numerical method simulates the liquid/solid interaction, and permits the study of the film rupture process with inertial effects and arbitrarily large contact angles, in both two and three dimensions. We focus in particular on the evolution of length scales in a perturbed film as it breaks up, and the spatial organization of the resulting drops. We compare our results to recent experiments, where this instability mechanism has been harnessed for the self-assembly of ordered arrays of metallic nanoparticles (ACS App. Mat. and Int., 2014, 6, 5835).

  12. Numerical simulations of the Single-mode, Doubly-shocked Richtmyer-Meshkov (RM) Instability

    NASA Astrophysics Data System (ADS)

    Karkhanis, Varad; Ramaprabhu, Praveen

    2014-11-01

    We describe results from numerical simulations of a single-mode, doubly-shocked material interface between gases of different densities. The time interval between the shocks was varied to observe interfacial growth due to Richtmyer-Meshkov Instability initialized with different amplitudes. The simulations were performed with low and high density ratio fluids (A = 0.15 and A = -0.99), where the latter case is relevant to ejecta formation. We compare the growth rates from our simulations after the first and second shocks with linear, nonlinear and ejecta models. In the heavy to light configuration (A = -0.99), we observe two consecutive phase inversions following each shock. We have also investigated the effect of variations in the initial interface perturbation to include sine, chevron, sawtooth, and square-wave form, and find our results to be of relevance to machined target experiments.

  13. Full-Scale Direct Numerical Simulation of Two- and Three-Dimensional Instabilities and Rivulet Formulation in Heated Falling Films

    NASA Technical Reports Server (NTRS)

    Krishnamoorthy, S.; Ramaswamy, B.; Joo, S. W.

    1995-01-01

    A thin film draining on an inclined plate has been studied numerically using finite element method. Three-dimensional governing equations of continuity, momentum and energy with a moving boundary are integrated in an arbitrary Lagrangian Eulerian frame of reference. Kinematic equation is solved to precisely update interface location. Rivulet formation based on instability mechanism has been simulated using full-scale computation. Comparisons with long-wave theory are made to validate the numerical scheme. Detailed analysis of two- and three-dimensional nonlinear wave formation and spontaneous rupture forming rivulets under the influence of combined thermocapillary and surface-wave instabilities is performed.

  14. Numerical simulations of combustion instabilities in gas turbine combustors, with applications

    NASA Astrophysics Data System (ADS)

    Swenson, Grant Douglas

    Recent advances in technology have opened up a potential market for small gas turbine power systems in the 50--100 MW range. In an effort to improve their systems, the gas-turbine industry is interested in understanding and controlling combustion instabilities as well as reducing pollutant production. To understand the dynamics inherent in a combustion system, information about the flow field behavior is required. Because of a scarcity of available experimental or numerical results for full-scale gas-turbine combustors, we decided to use numerical simulations to provide the required information about the flow field dynamics. The ability of the numerical simulations to reproduce unstable behavior in combustion environments will be presented. The investigation of the flow field dynamics has been conducted for three test cases; a planar heat source in a tube, premixed flow in a dump combustor, and premixed and diffusion flames in a full-scale gas turbine combustor. The numerically determined unsteady acoustic modes will be shown to compare well with theory and experiments. An investigation of the local heat release response to an unsteady flow field is conducted for incorporation into an approximate analysis method. The results of including a Helmholtz resonator in a dump combustor as a passive control mechanism will be presented. The production of NOx and CO will be compared between stable and unstable flow configurations. The pollutant results indicate that for the planar flame in a tube and the dump combustor, the NOx levels at the exit plane are reduced when the system is unstable.

  15. Numerical simulation of multi-material mixing in an inclined interface Richtmyer-Meshkov instability

    NASA Astrophysics Data System (ADS)

    Subramaniam, Akshay; Lele, Sanjiva

    2015-06-01

    The Richtmyer-Meshkov instability arises when a shock wave interacts with an interface separating two fluids. In this work, high fidelity simulations of shock induced multi-material mixing between air and SF6 in a shock tube are performed for a Mach 1.5 shock interacting with a planar material interface that is inclined with respect to the shock propagating direction. In the current configuration, unlike in the classical sinusoidal interface case, the evolution of the interface is fully non-linear from early time onwards. The simulations attempt to replicate an experiment conducted at the Texas A&M fluid mixing shock tube facility. Simulations of this problem at multiple spatial resolutions (upto 270 million grid points) have shown that even low order statistics like the net circulation are hard to capture at resolutions where the classical RM cases yield good results. Tight coupling between numerics and flow physics and large range of spatial scales make this a challenging problem to simulate numerically. Simulations shown are conducted with an extended version of the MIRANDA solver developed by Cook et al. (2007) which combines high-order compact finite differences with localized non-linear artificial properties for shock and interface capturing.

  16. Numerical simulation on macro-instability of coupling flow field structure in jet-stirred tank

    NASA Astrophysics Data System (ADS)

    Luan, D. Y.; Lu, J. P.; Bu, Q. X.; Zhang, S. F.; Zheng, S. X.

    2016-05-01

    The velocity field macro-instability (MI) can help to improve the mixing efficiency. In this work, the MI features of flow field induced by jet-stirred coupling action is studied by using computational fluid dynamics (CFD) simulations. The numerical simulation method of jet-stirred model was established based on standard turbulent equations, and the impeller rotation was modeled by means of the Sliding Mesh (SM) technology. The numerical results of test fluid (water) power consumption were compared with the data obtained by power test experiments. The effects of jet flow velocity and impeller speed on MI frequency were analyzed thoroughly. The results show that the calculated values of power consumption agree well with the experiment measured data, which validates the turbulent model, and the flow structure and MI frequency distribution are affected by both impeller speed and jet flow rate. The amplitude of MI frequency increases obviously with the increasing rotation speed of impeller and the eccentric jet rate, and it can be enhanced observably by eccentric jet rate, in condition of comparatively high impeller speed. At this time, the MI phenomenon disappears with the overall chaotic mixing.

  17. Perturbation theory and numerical modelling of weakly and moderately nonlinear incompressible Richtmyer-Meshkov instability

    NASA Astrophysics Data System (ADS)

    Herrmann, M.; Velikovich, A. L.; Abarzhi, S. I.

    2014-10-01

    A study of incompressible two-dimensional Richtmyer-Meshkov instability by means of high-order Eulerian perturbation theory and numerical simulations is reported. Nonlinear corrections to Richtmyer's impulsive formula for the bubble and spike growth rates have been calculated analytically for arbitrary Atwood number and an explicit formula has been obtained for it in the Boussinesq limit. Conditions for early-time acceleration and deceleration of the bubble and the spike have been derived. In our simulations we have solved 2D unsteady Navier-Stokes equations for immiscible incompressible fluids using the finite volume fractional step flow solver NGA developed by, coupled to the level set based interface solver LIT,. The impact of small amounts of viscosity and surface tension on the RMI flow dynamics is studied numerically. Simulation results are compared to the theory to demonstrate successful code verification and highlight the influence of the theory's ideal inviscid flow assumption. Theoretical time histories of the interface curvature at the bubble and spike tip and the profiles of vertical and horizontal velocities have been favorably compared to simulation results, which converge to the theoretical predictions as the Reynolds and Weber numbers are increased. Work supported by the US DOE/NNSA.

  18. Numerical simulation of three-dimensional reconnection due to the instability of collisionless current sheets

    NASA Astrophysics Data System (ADS)

    Büchner, J.; Kuska, J.-P.

    1997-01-01

    Based on analytical calculations we have currently argued that spontaneous reconnection through thin collisionless current sheets is an essentially three-dimensional (3 D) process (Büchner, 1996 a, b). Since 3 D kinetic PIC codes have become available, the three dimensional nature of the collisionless current sheet decay are now illustrated by numerical simulations (Büchner and Kuska, 1996; Pritchett and Coroniti, 1996; Zhu and Winglee, 1996). While the latter two claim a coupling to a longer wavelength kink mode as a main factor, destabilizing thin current sheets in 3 D, our simulations have revealed that even shorter scale perturbations in the current direction suffice to destabilize thin sheets very quickly. Since past simulation runs, however, were limited to mass ratios near unity, the influence of the electrons was not treated adequately. We have now investigated the stability of thin collisionless current sheets including 64 times lighter negatively charged particles. We can now show that while the two-dimensional tearing instability slows down for M = M_p/m_e = 64, the three-dimensional current sheet decay is a much faster process - practically as fast as the mass ratio M = 1 3 D sheet decay, even without kinking the sheet. We further conclude that, unlike the two-dimensional tearing instability, the three-dimensional decay of thin current sheets is not controlled by the electrons. For a sheet width comparable with the ion inertial length, we also recovered signatures of the Hall effect as predicted by Vasyliunas (1975) in the mass ratio M = 64 case. The ion inertial length seems to be the critical scale at which the sheet starts to decay.

  19. Adaptive Wavelet-Based Direct Numerical Simulations of Rayleigh-Taylor Instability

    NASA Astrophysics Data System (ADS)

    Reckinger, Scott J.

    The compressible Rayleigh-Taylor instability (RTI) occurs when a fluid of low molar mass supports a fluid of higher molar mass against a gravity-like body force or in the presence of an accelerating front. Intrinsic to the problem are highly stratified background states, acoustic waves, and a wide range of physical scales. The objective of this thesis is to develop a specialized computational framework that addresses these challenges and to apply the advanced methodologies for direct numerical simulations of compressible RTI. Simulations are performed using the Parallel Adaptive Wavelet Collocation Method (PAWCM). Due to the physics-based adaptivity and direct error control of the method, PAWCM is ideal for resolving the wide range of scales present in RTI growth. Characteristics-based non-reflecting boundary conditions are developed for highly stratified systems to be used in conjunction with PAWCM. This combination allows for extremely long domains, which is necessary for observing the late time growth of compressible RTI. Initial conditions that minimize acoustic disturbances are also developed. The initialization is consistent with linear stability theory, where the background state consists of two diffusively mixed stratified fluids of differing molar masses. The compressibility effects on the departure from the linear growth, the onset of strong non-linear interactions, and the late-time behavior of the fluid structures are investigated. It is discovered that, for the thermal equilibrium case, the background stratification acts to suppress the instability growth when the molar mass difference is small. A reversal in this monotonic behavior is observed for large molar mass differences, where stratification enhances the bubble growth. Stratification also affects the vortex creation and the associated induced velocities. The enhancement and suppression of the RTI growth has important consequences for a detailed understanding of supernovae flame front

  20. Numerical simulations of resistive magnetohydrodynamic instabilities in a poloidal divertor tokamak

    NASA Astrophysics Data System (ADS)

    Uchimoto, E.

    1988-03-01

    A new 3-D resistive MHD initial value code RPD has been successfully developed from scratch to study the linear and nonlinear evolution of long wavelength resistive MHD instabilities in a square cross-section tokamak with or without a poloidal divertor. The code numerically advances the full set of compressible resistive MHD equations in a toroidal geometry, with an important option of permitting the divertor separatrix and the region outside it to be in the computational domain. A severe temporal step size restriction for numerical stability imposed by the fast compressional waves was removed by developing and implementing a new, efficient semi-implicit scheme extending one first proposed by Harned and Kerner. As a result, the code typically runs faster than that with a mostly explicit scheme by a factor of about the aspect ratio. The equilibrium input for RPD is generated by a new 2-D code EQPD that is based on the Chodura-Schluter method. The RPD code, as well as the new semi-implicit scheme, has passed very extensive numerical tests in both divertor and divertorless geometries. Linear and nonlinear simulations in a divertorless geometry have reproduced the standard, previously known results. In a geometry with a four-node divertor the m = 2, n = 1 (2/1) tearing mode tends to be linearly stabilized as the q = 2 surface approaches the divertor separatrix. However, the m = 1, n = 1 (1/1) resistive kink mode remains relatively unaffected by the nearness of the q = 1 surface to the divertor separatrix. When plasma current is added to the region outside the divertor separatrix, the 2/1 tearing mode is linearly stabilized not by this current, but by the profile modifications induced near the q = 2 surface and the divertor separatrix. A similar stabilization effect is seen for the 1/1 resistive kink mode, but to a lesser extent.

  1. Are oceanic plateaus sites of komatiite formation?

    NASA Astrophysics Data System (ADS)

    Storey, M.; Mahoney, J. J.; Kroenke, L. W.; Saunders, A. D.

    1991-04-01

    During Cretaceous and Tertiary time a series of oceanic terranes were accreted onto the Pacific continental margin of Colombia. The island of Gorgona is thought to represent part of the most recent, early Eocene, terrane-forming event. Gorgona is remarkable for the occurrence of komatiites of middle Cretaceous age, having MgO contents up to 24%. The geochemistry of spatially and temporally associated tholeiites suggests that Gorgona is an obducted fragment of the oceanic Caribbean Plateau, postulated by Duncan and Hargraves (1984) to have formed at 100 to 75 Ma over the Galapagos hotspot. Further examples of high-MgO oceanic lavas that may represent fragments of the Caribbean Plateau occur in allochthonous terranes on the island of Curaçao in the Netherlands Antilles and in the Romeral zone ophiolites in the southwestern Colombian Andes. These and other examples suggest that the formation of high-MgO liquids may be a feature of oceanic-plateau settings. The association of Phanerozoic komatiites with oceanic plateaus, coupled with thermal considerations, provides a plausible analogue for the origin of some komatiite-tholeiite sequences in Archean greenstone belts.

  2. Numerical analysis of flow instability in the water wall of a supercritical CFB boiler with annular furnace

    NASA Astrophysics Data System (ADS)

    Xie, Beibei; Yang, Dong; Xie, Haiyan; Nie, Xin; Liu, Wanyu

    2016-08-01

    In order to expand the study on flow instability of supercritical circulating fluidized bed (CFB) boiler, a new numerical computational model considering the heat storage of the tube wall metal was presented in this paper. The lumped parameter method was proposed for wall temperature calculation and the single channel model was adopted for the analysis of flow instability. Based on the time-domain method, a new numerical computational program suitable for the analysis of flow instability in the water wall of supercritical CFB boiler with annular furnace was established. To verify the code, calculation results were respectively compared with data of commercial software. According to the comparisons, the new code was proved to be reasonable and accurate for practical engineering application in analysis of flow instability. Based on the new program, the flow instability of supercritical CFB boiler with annular furnace was simulated by time-domain method. When 1.2 times heat load disturbance was applied on the loop, results showed that the inlet flow rate, outlet flow rate and wall temperature fluctuated with time eventually remained at constant values, suggesting that the hydrodynamic flow was stable. The results also showed that in the case of considering the heat storage, the flow in the water wall is easier to return to stable state than without considering heat storage.

  3. Laboratory and numerical simulation of internal wave attractors and their instability.

    NASA Astrophysics Data System (ADS)

    Brouzet, Christophe; Dauxois, Thierry; Ermanyuk, Evgeny; Joubaud, Sylvain; Sibgatullin, Ilias

    2015-04-01

    Internal wave attractors are formed as result of focusing of internal gravity waves in a confined domain of stably stratified fluid due to peculiarities of reflections properties [1]. The energy injected into domain due to external perturbation, is concentrated along the path formed by the attractor. The existence of attractors was predicted theoretically and proved both experimentally and numerically [1-4]. Dynamics of attractors is greatly influenced by geometrical focusing, viscous dissipation and nonlinearity. The experimental setup features Schmidt number equal to 700 which impose constraints on resolution in numerical schemes. Also for investigation of stability on large time intervals (about 1000 periods of external forcing) numerical viscosity may have significant impact. For these reasons, we have chosen spectral element method for investigation of this problem, what allows to carefully follow the nonlinear dynamics. We present cross-comparison of experimental observations and numerical simulations of long-term behavior of wave attractors. Fourier analysis and subsequent application of Hilbert transform are used for filtering of spatial components of internal-wave field [5]. The observed dynamics shows a complicated coupling between the effects of local instability and global confinement of the fluid domain. The unstable attractor is shown to act as highly efficient mixing box providing the efficient energy pathway from global-scale excitation to small-scale wave motions and mixing. Acknowledgement, IS has been partially supported by Russian Ministry of Education and Science (agreement id RFMEFI60714X0090) and Russian Foundation for Basic Research, grant N 15-01-06363. EVE gratefully acknowledges his appointment as a Marie Curie incoming fellow at Laboratoire de physique ENS de Lyon. This work has been partially supported by the ONLITUR grant (ANR-2011-BS04-006-01) and achieved thanks to the resources of PSMN from ENS de Lyon 1. Maas, L. R. M. & Lam, F

  4. Three-dimensional numerical simulation of Marangoni instabilities in liquid bridges: influence of geometrical aspect ratio

    NASA Astrophysics Data System (ADS)

    Lappa, M.; Savino, R.; Monti, R.

    2001-05-01

    Oscillatory Marangoni convection in silicone oil-liquid bridges with different geometrical aspect ratios is investigated by three-dimensional and time-dependent numerical simulations, based on control volume methods in staggered cylindrical non-uniform grids. The three-dimensional oscillatory flow regimes are studied and compared with previous experimental and theoretical results. The results show that the critical wavenumber (m), related to the azimuthal spatio-temporal flow structure, is a monotonically decreasing function of the geometrical aspect ratio of the liquid bridge (defined as the ratio of length to diameter). For this function, a general correlation formula is found, which is in agreement with the previous experimental findings. The critical Marangoni number and the oscillation frequency are decreasing functions of the aspect ratio; however, the critical Marangoni number, based on the axial length of the bridge, does not change much with the aspect ratio. For each aspect ratio investigated, the onset of the instability from the axisymmetric steady state to the three-dimensional oscillatory one is characterized by the appearance of a standing wave regime that exhibits, after a certain time, a second transition to a travelling wave regime. The standing wave regime is more stable for lower aspect ratios since it lasts for a long time. This behaviour is explained on the basis of the propagation velocity of the disturbances in the liquid phase. For this velocity, a general correlation law is found as a function of the aspect ratio and of the Marangoni number. Copyright

  5. Evolution of length scales and statistics of Richtmyer-Meshkov instability from direct numerical simulations.

    PubMed

    Tritschler, V K; Zubel, M; Hickel, S; Adams, N A

    2014-12-01

    In this study we present direct numerical simulation results of the Richtmyer-Meshkov instability (RMI) initiated by Ma=1.05,Ma=1.2, and Ma=1.5 shock waves interacting with a perturbed planar interface between air and SF(6). At the lowest shock Mach number the fluids slowly mix due to viscous diffusion, whereas at the highest shock Mach number the mixing zone becomes turbulent. When a minimum critical Taylor microscale Reynolds number is exceeded, an inertial range spectrum emerges, providing further evidence of transition to turbulence. The scales of turbulent motion, i.e., the Kolmogorov length scale, the Taylor microscale, and the integral length, scale are presented. The separation of these scales is found to increase as the Reynolds number is increased. Turbulence statistics, i.e., the probability density functions of the velocity and its longitudinal and transverse derivatives, show a self-similar decay and thus that turbulence evolving from RMI is not fundamentally different from isotropic turbulence, though nominally being only isotropic and homogeneous in the transverse directions. PMID:25615181

  6. Numerical simulations of the magnetorotational instability in protoneutron stars - I. Influence of buoyancy

    NASA Astrophysics Data System (ADS)

    Guilet, Jérôme; Müller, Ewald

    2015-06-01

    The magnetorotational instability (MRI) is considered to be a promising mechanism to amplify the magnetic field in fast-rotating protoneutron stars. In contrast to accretion discs, radial buoyancy driven by entropy and lepton fraction gradients is expected to have a dynamical role as important as rotation and shear. We investigate the poorly known impact of buoyancy on the non-linear phase of the MRI, by means of three-dimensional numerical simulations of a local model in the equatorial plane of a protoneutron star. The use of the Boussinesq approximation allows us to utilize a shearing box model with clean shearing periodic boundary conditions, while taking into account the buoyancy driven by radial entropy and composition gradients. We find significantly stronger turbulence and magnetic fields in buoyantly unstable flows. On the other hand, buoyancy has only a limited impact on the strength of turbulence and magnetic field amplification for buoyantly stable flows in the presence of a realistic thermal diffusion. The properties of the turbulence are, however, significantly affected in the latter case. In particular, the toroidal components of the magnetic field and of the velocity become even more dominant with respect to the poloidal ones. Furthermore, we observed in the regime of stable buoyancy the formation of long-lived coherent structures such as channel flows and zonal flows. Overall, our results support the ability of the MRI to amplify the magnetic field significantly even in stably stratified regions of protoneutron stars.

  7. Direct Numerical Simulation of Richtmeyer-Meshkov Instability Using pWAMR

    NASA Astrophysics Data System (ADS)

    Grenga, Temistocle; Paolucci, Samuel

    2015-11-01

    The parallel Wavelet Adaptive Multiresolution Representation (pWAMR) method is used to simulate the Richtmyer-Meshkov instability caused by a shock interacting with a density-stratified interface. The physical problem is studied in several configurations. We present results of numerical studies that investigate the influence of initial condition parameters (amplitude and wavelength of perturbations) on mixing and transition. In addition, the evaluation of turbulence statistics provides a measure of the mixing across the scales and the correlation with the initial condition parameters. The problem is modeled using the compressible reactive Navier-Stokes equations for a gas mixture, including multi-component diffusion, Soret and Dufour effects, and state dependent thermodynamic and transport properties. Since the amplitudes of wavelets provide a direct measure of the local error, the method is able to efficiently capture to any desired accuracy a wide range of spatial scales using a relatively small number of degrees of freedom by evolving the dynamically adaptive grid. In an effective fashion, the multilevel structure of the algorithm provides a simple way to adapt computational refinements to local demands of the solution, thus automatically producing verified solutions. Supported by C-SWARM through the Department of Energy, National Nuclear Security Administration, under Award Number DE-NA0002377.

  8. Numerical study of acoustic instability in a partly lined flow duct using the full linearized Navier-Stokes equations

    NASA Astrophysics Data System (ADS)

    Xin, Bo; Sun, Dakun; Jing, Xiaodong; Sun, Xiaofeng

    2016-07-01

    Lined ducts are extensively applied to suppress noise emission from aero-engines and other turbomachines. The complex noise/flow interaction in a lined duct possibly leads to acoustic instability in certain conditions. To investigate the instability, the full linearized Navier-Stokes equations with eddy viscosity considered are solved in frequency domain using a Galerkin finite element method to compute the sound transmission in shear flow in the lined duct as well as the flow perturbation over the impedance wall. A good agreement between the numerical predictions and the published experimental results is obtained for the sound transmission, showing that a transmission peak occurs around the resonant frequency of the acoustic liner in the presence of shear flow. The eddy viscosity is an important influential factor that plays the roles of both providing destabilizing and making coupling between the acoustic and flow motions over the acoustic liner. Moreover, it is shown from the numerical investigation that the occurrence of the sound amplification and the magnitude of transmission coefficient are closely related to the realistic velocity profile, and we find it essential that the actual variation of the velocity profile in the axial direction over the liner surface be included in the computation. The simulation results of the periodic flow patterns possess the proper features of the convective instability over the liner, as observed in Marx et al.'s experiment. A quantitative comparison between numerical and experimental results of amplitude and phase of the instability is performed. The corresponding eigenvalues achieve great agreement.

  9. Numerical Experiments Using a Convective Flux Limiter on a Turbulent Single-Mode Rayleigh-Taylor Instability

    SciTech Connect

    Cloutman, L.D.

    2000-07-10

    Direct numerical simulation and large eddy simulations are powerful tools for studying turbulent flows. Unfortunately, they are computationally demanding in terms of run times, storage, and accuracy of the numerical method used. In particular, high order methods promise high accuracy on a given grid, but they often fail to deliver the expected accuracy due to dispersive truncation errors that appear as unphysical oscillations in the numerical solutions. This report describes a nonlinear flux limiter that has been applied to the second-order tensor viscosity method and markedly reduces the dispersive truncation errors. A Rayleigh-Taylor instability is simulated to show how well the flux limiter works.

  10. Direct numerical simulation of transitional mixed convection flows: Viscous and inviscid instability mechanisms

    NASA Astrophysics Data System (ADS)

    Sengupta, T. K.; Bhaumik, S.; Bose, Rikhi

    2013-09-01

    Receptivity studies using direct numerical simulation require computations of equilibrium flow and its response to deterministic excitation. Equivalent flow problem, without heat interaction for zero-pressure gradient boundary layer, has been studied with respect to wall-excitation by a finite difference high accuracy method based on the solution of Navier-Stokes equation in Sengupta and Bhaumik [Phys. Rev. Lett. 107, 154501 (2011)] and Sengupta et al. [Phys. Rev. E 85, 026308 (2012)]. One of the key features of this study has been that the same methodology is used for computing the equilibrium flow and the disturbance field. Computation of equilibrium flow was performed by solving Navier-Stokes equation to include the leading edge of the plate, so that the effects of leading edge singularity and the growth of the boundary layer is included in the nonlinear framework. When the same methodology is attempted for mixed convection flows past horizontal plate (with Boussinesq approximation to model heat transfer effects) some of the equilibrium flow features could not be explained with linear viscous instability theory results. For horizontal hot flat plate with adiabatic wall conditions, the equilibrium flow could be computed and its receptivity could be correlated with linear spatial theory for lower buoyancy parameter. Here, we focus on receptivity of mixed convection flows to wall excitation for the following cases which do not allow computing the equilibrium flows from the solution of Navier-Stokes equation: (i) Aadiabatic horizontal flat plate cooled significantly at the leading edge only and (ii) strongly heated isothermal wedge flow for a wedge angle of 60°. The cold plate case is particularly interesting as the linear spatial theory indicates enhanced stabilization for higher magnitude of the buoyancy parameter. Results presented for the cold plate case indicates disturbance growth outside the shear layer. This prompted us to re-investigate various mechanisms

  11. Experimental and numerical study on the wavy instability in a Rayleigh-Bénard-Poiseuille flow: non linear effects

    NASA Astrophysics Data System (ADS)

    Seychelles, F.; Mergui, S.; Nicolas, X.

    2012-11-01

    A combined experimental and numerical study of a Rayleigh-Bénard-Poiseuille air flow in a rectangular channel is presented. The aim of the present paper is to characterize a secondary instability, referred to as wavy instability and known to be a convective instability, with the objective to identify the best situation for an optimal homogenization of heat transfers in the system. A periodic mechanical excitation is introduced at channel inlet and the spatial and temporal evolution of the temperature fluctuations are analyzed, depending on the Rayleigh number, the frequency and the amplitude of the perturbation. The Reynolds number is fixed. As the saturated state is a priori the best situation to homogenize the transfers, the objective is to expand the saturation area and to generate a maximum saturation amplitude value. It is shown that the best choice is a high Rayleigh number or/and a large magnitude of perturbation associated with a specific low value of the forcing frequency.

  12. Numerical simulation of unsteady heat release of low frequency instabilities in a dump combustor

    NASA Astrophysics Data System (ADS)

    Laverdant, A.

    The influence of combustion instabilities on heat transfer is investigated using an adaptation of KIVA code. A simulation of low-frequency instabilities observed on a small burner is described. It is shown that the turbulence is distributed in the flame zone, and the heat transfer increases by acoustic pulsation emitted from the entrance plane of the cavity.

  13. Material property discontinuities in intervertebral disc porohyperelastic finite element models generate numerical instabilities due to volumetric strain variations.

    PubMed

    Ruiz, C; Noailly, J; Lacroix, D

    2013-10-01

    Numerical studies of the intervertebral disc (IVD) are important to better understand the load transfer and the mechanobiological processes within the disc. Among the relevant calculations, fluid-related outputs are critical to describe and explore accurately the tissue properties. Porohyperelastic finite element models of IVD can describe accurately the disc behaviour at the organ level and allow the inclusion of fluid effects. However, results may be affected by numerical instabilities when fast load rates are applied. We hypothesized that such instabilities would appear preferentially at material discontinuities such as the annulus-nucleus boundary and should be considered when testing mesh convergence. A L4-L5 IVD model including the nucleus, annulus and cartilage endplates were tested under pure rotational loads, with different levels of mesh refinement. The effect of load relaxation and swelling were also studied. Simulations indicated that fluid velocity oscillations appeared due to numerical instability of the pore pressure spatial derivative at material discontinuities. Applying local refinement only was not enough to eliminate these oscillations. In fact, mesh refinements had to be local, material-dependent, and supplemented by the creation of a material transition zone, including interpolated material properties. Results also indicated that oscillations vanished along load relaxation, and faster attenuation occurred with the incorporation of the osmotic pressure. We concluded that material discontinuities are a major cause of instability for poromechanical calculations in multi-tissue models when load velocities are simulated. A strategy was presented to address these instabilities and recommendations on the use of IVD porohyperelastic models were given. PMID:23796430

  14. Analysis of Stationary Crossflow Instability on HIFiRE-5 Using Direct Numerical Simulation

    NASA Astrophysics Data System (ADS)

    Dinzl, Derek John

    Direct numerical simulation is performed on a 38.1% scale HIFiRE-5 forebody to study stationary crossflow instability. Computations use the US3D Navier-Stokes solver to simulate Mach 6 flow at Reynolds numbers of 8.1 x 10 6 /m and 11.8 x 106 /m, which are conditions used by quiet tunnel experiments at Purdue University. Distributed roughness with point-to-point height variation on the computational grid and maximum heights of 0.5-4.0 mum is used with the intent to emulate smooth-body transition and excite the naturally-occurring most unstable disturbance wavenumber. Cases at the low Reynolds number condition use three grid sizes, and hence three different roughness patterns of varying wavelength, and demonstrate that the final flow solution is extremely dependent on the particular roughness pattern. The same roughness pattern is interpolated onto each grid which yields similar solutions, indicating grid convergence. At the high Reynolds number condition, a steady physical mechanism is introduced which explains sharp increases seen in the wall heat flux for both computations and experiment. Namely, the sharp increase is caused by large streamwise velocity disturbances impinging on the wall. Evolution of disturbance spanwise wavelength is computed, and it is found that this wavelength is more sensitive to Reynolds number than roughness, indicating that the disturbance wavelength is primarily flow--selected for these cases. The calculation of disturbance growth rates shows the region over which crossflow disturbances behave linearly and where nonlinear effects become important. The effect of roughness height and nose sharpness are considered, and both were found to have a large effect on the resulting heating pattern. Crossflow vortex coalescence is observed and a possible cause is discussed.

  15. Thermal and Magnetorotational Instability in the Interstellar Medium: Two-dimensional Numerical Simulations

    NASA Astrophysics Data System (ADS)

    Piontek, Robert A.; Ostriker, Eve C.

    2004-02-01

    The structure and dynamics of diffuse gas in the Milky Way and other disk galaxies may be strongly influenced by thermal and magnetorotational instabilities (TI and MRI, respectively) on scales ~1-100 pc. We initiate a study of these processes, using two-dimensional numerical hydrodynamic and magnetohydrodynamic simulations with conditions appropriate for the atomic interstellar medium (ISM). Our simulations incorporate thermal conduction and adopt local ``shearing-periodic'' equations of motion and boundary conditions to study dynamics of a (100 pc) 2 radial-vertical section of the disk. We demonstrate, consistent with previous work, that nonlinear development of ``pure TI'' produces a network of filaments that condense into cold clouds at their intersections, yielding a distinct two-phase warm/cold medium within ~20 Myr. TI-driven turbulent motions of the clouds and warm intercloud medium are present but saturate at quite subsonic amplitudes for uniform initial P/k=2000 K cm -3. MRI has previously been studied in near-uniform media; our simulations include both TI+MRI models, which begin from uniform-density conditions, and cloud+MRI models, which begin with a two-phase cloudy medium. Both the TI+MRI and cloud+MRI models show that MRI develops within a few galactic orbital times, just as for a uniform medium. The mean separation between clouds can affect which MRI mode dominates the evolution. Provided intercloud separations do not exceed half the MRI wavelength, we find the MRI growth rates are similar to those for the corresponding uniform medium. This opens the possibility that if low cloud volume filling factors increase MRI dissipation times compared to those in a uniform medium, then MRI-driven motions in the ISM could reach amplitudes comparable to observed H I turbulent line widths.

  16. Numerical study of particle-induced Rayleigh-Taylor instability: Effects of particle settling and entrainment

    NASA Astrophysics Data System (ADS)

    Chou, Yi-Ju; Shao, Yun-Chuan

    2016-04-01

    In this study, we investigate Rayleigh-Taylor instability in which the density stratification is caused by the suspension of particles in liquid flows using the conventional single-phase model and Euler-Lagrange (EL) two-phase model. The single-phase model is valid only when the particles are small and number densities are large, such that the continuum approximation applies. The present single-phase results show that the constant settling of the particle concentration restricts the lateral development of the vortex ring, which results in a decrease of the rising speed of the Rayleigh-Taylor bubbles. The EL model enables the investigation of particle-flow interaction and the influence of particle entrainment, resulting from local non-uniformity in the particle distribution. We compare bubble dynamics in the single-phase and EL cases, and our results show that the deviation between the two cases becomes more pronounced when the particle size increases. The main mechanism responsible for the deviation is particle entrainment, which can only be resolved in the EL model. We provide a theoretical argument for the small-scale local entrainment resulting from the local velocity shear and non-uniformity of the particle concentration. The theoretical argument is supported by numerical evidence. Energy budget analysis is also performed and shows that potential energy is released due to the interphase drag and buoyant effect. The buoyant effect, which results in the transformation of potential energy into kinetic energy and shear dissipation, plays a key role in settling enhancement. We also find that particle entrainment increases the shear dissipation, which in turn enhances the release of potential energy.

  17. Direct numerical simulation of interfacial instabilities: A consistent, conservative, all-speed, sharp-interface method

    NASA Astrophysics Data System (ADS)

    Chang, Chih-Hao; Deng, Xiaolong; Theofanous, Theo G.

    2013-06-01

    We present a conservative and consistent numerical method for solving the Navier-Stokes equations in flow domains that may be separated by any number of material interfaces, at arbitrarily-high density/viscosity ratios and acoustic-impedance mismatches, subjected to strong shock waves and flow speeds that can range from highly supersonic to near-zero Mach numbers. A principal aim is prediction of interfacial instabilities under superposition of multiple potentially-active modes (Rayleigh-Taylor, Kelvin-Helmholtz, Richtmyer-Meshkov) as found for example with shock-driven, immersed fluid bodies (locally oblique shocks)—accordingly we emphasize fidelity supported by physics-based validation, including experiments. Consistency is achieved by satisfying the jump discontinuities at the interface within a conservative 2nd-order scheme that is coupled, in a conservative manner, to the bulk-fluid motions. The jump conditions are embedded into a Riemann problem, solved exactly to provide the pressures and velocities along the interface, which is tracked by a level set function to accuracy of O(Δx5, Δt4). Subgrid representation of the interface is achieved by allowing curvature of its constituent interfacial elements to obtain O(Δx3) accuracy in cut-cell volume, with attendant benefits in calculating cell- geometric features and interface curvature (O(Δx3)). Overall the computation converges at near-theoretical O(Δx2). Spurious-currents are down to machine error and there is no time-step restriction due to surface tension. Our method is built upon a quadtree-like adaptive mesh refinement infrastructure. When necessary, this is supplemented by body-fitted grids to enhance resolution of the gas dynamics, including flow separation, shear layers, slip lines, and critical layers. Comprehensive comparisons with exact solutions for the linearized Rayleigh-Taylor and Kelvin-Helmholtz problems demonstrate excellent performance. Sample simulations of liquid drops subjected to

  18. Numerical methods for instability mitigation in the modeling of laser wakefield accelerators in a Lorentz-boosted frame

    SciTech Connect

    Vay, J.-L.; Geddes, C.G.R.; Cormier-Michel, E.; Grote, D.P.

    2011-07-01

    Modeling of laser-plasma wakefield accelerators in an optimal frame of reference has been shown to produce orders of magnitude speed-up of calculations from first principles. Obtaining these speedups required mitigation of a high-frequency instability that otherwise limits effectiveness. In this paper, methods are presented which mitigated the observed instability, including an electromagnetic solver with tunable coefficients, its extension to accommodate Perfectly Matched Layers and Friedman's damping algorithms, as well as an efficient large bandwidth digital filter. It is observed that choosing the frame of the wake as the frame of reference allows for higher levels of filtering or damping than is possible in other frames for the same accuracy. Detailed testing also revealed the existence of a singular time step at which the instability level is minimized, independently of numerical dispersion. A combination of the techniques presented in this paper prove to be very efficient at controlling the instability, allowing for efficient direct modeling of 10 GeV class laser plasma accelerator stages. The methods developed in this paper may have broader application, to other Lorentz-boosted simulations and Particle-In-Cell simulations in general.

  19. Fundamental form of the electrostatic δf-PIC algorithm and discovery of a converged numerical instability

    NASA Astrophysics Data System (ADS)

    Wilkie, George J.; Dorland, William

    2016-05-01

    The δf particle-in-cell algorithm has been a useful tool in studying the physics of plasmas, particularly turbulent magnetized plasmas in the context of gyrokinetics. The reduction in noise due to not having to resolve the full distribution function indicates an efficiency advantage over the standard ("full-f") particle-in-cell. Despite its successes, the algorithm behaves strangely in some circumstances. In this work, we document a fully resolved numerical instability that occurs in the simplest of multiple-species test cases: the electrostatic ΩH mode. There is also a poorly understood numerical instability that occurs when one is under-resolved in particle number, which may require a prohibitively large number of particles to stabilize. Both of these are independent of the time-stepping scheme, and we conclude that they exist if the time advancement were exact. The exact analytic form of the algorithm is presented, and several schemes for mitigating these instabilities are also presented.

  20. High resolution numerical investigation on the effect of convective instability on long term CO2 storage in saline aquifers

    NASA Astrophysics Data System (ADS)

    Lu, C.; Lichtner, P. C.

    2007-07-01

    CO2 sequestration (capture, separation, and long term storage) in various geologic media including depleted oil reservoirs, saline aquifers, and oceanic sediments is being considered as a possible solution to reduce green house gas emissions. Dissolution of supercritical CO2 in formation brines is considered an important storage mechanism to prevent possible leakage. Accurate prediction of the plume dissolution rate and migration is essential. Analytical analysis and numerical experiments have demonstrated that convective instability (Rayleigh instability) has a crucial effect on the dissolution behavior and subsequent mineralization reactions. Global stability analysis indicates that a certain grid resolution is needed to capture the features of density-driven fingering phenomena. For 3-D field scale simulations, high resolution leads to large numbers of grid nodes, unfeasible for a single workstation. In this study, we investigate the effects of convective instability on geologic sequestration of CO2 by taking advantage of parallel computing using the code PFLOTRAN, a massively parallel 3-D reservoir simulator for modeling subsurface multiphase, multicomponent reactive flow and transport based on continuum scale mass and energy conservation equations. The onset, development and long-term fate of a supercritical CO2 plume will be resolved with high resolution numerical simulations to investigate the rate of plume dissolution caused by fingering phenomena.

  1. An Experimental and Numerical Study of Roughness-Induced Instabilities in a Mach 3.5 Boundary Layer

    NASA Technical Reports Server (NTRS)

    Kegerise, Michael A.; King, Rudolph A.; Owens, Lewis R.; Choudhari, Meelan M.; Norris, Andrew T.; Li, Fei; Chang, Chau-Layn

    2012-01-01

    Progress on a joint experimental and numerical study of laminar-to-turbulent transition induced by an isolated roughness element in a high-speed laminar boundary layer is reported in this paper. The numerical analysis suggests that transition is driven by the instability of high- and low-speed streaks embedded in the wake of the isolated roughness element. In addition, spatial stability analysis revealed that the wake flow supports multiple modes (even and odd) of convective instabilities that experience strong enough growth to cause transition. The experimental measurements, which included hot-wire and pitot-probe surveys, confirmed the existence of embedded high- and low-speed streaks in the roughness wake. Furthermore, the measurements indicate the presence of both even and odd modes of instability, although their relative magnitude depends on the specifics of the roughness geometry and flow conditions (e.g., the value of Re(sub kk) or k/delta. For the two test cases considered in the measurements (Re(sub kk) values of 462 and 319), the even mode and the odd mode were respectively dominant and appear to play a primary role in the transition process.

  2. Numerical studies of the Weibel Instability in Intense Charged Particle Beams with Large Energy Anisotropy

    NASA Astrophysics Data System (ADS)

    Lee, Wei-Li; Startsev, Edward A.; Davidson, Ronald C.

    2004-11-01

    In intense charged particle beams with large temperature anisotropy free energy is available to drive a transverse electromagnetic Weibel-type instability. The finite transverse geometry of the confined beam makes a detailed theoretical investigation difficult. In this paper the newly developed bEASt (beam eigenmode and spectra) code which solves the linearized Vlasov-Maxwell equations is used to investigate the detailed properties of the Weibel instability for a long charge bunch propagating through a cylindrical pipe of radius r_w. The stability analysis is carried out for azimuthally symmetric perturbations about a two-temperature thermal equilibrium distribution in the smooth-focusing approximation. To study the nonlinear stage of the instability, the Darwin model is being developed and incorporated into the Beam Equilibrium Stability and Transport(BEST) code.

  3. Release of potential instability by mesoscale triggering - An objective model simulation. [in precipitation numerical weather forecasting

    NASA Technical Reports Server (NTRS)

    Matthews, D. A.

    1978-01-01

    The effects of mesoscale triggering on organized nonsevere convective cloud systems in the High Plains are considered. Two experiments were conducted to determine if a one-dimensional quasi-time dependent model could (1) detect soundings which were sensitive to mesoscale triggering, and (2) discriminate between cases which had mesoscale organized convection and those with no organized convection. The MESOCU model was used to analyze the available potential instability and thermodynamic potential for cloud growth. It is noted that lifting is a key factor in the release of available potential instability on the High Plains.

  4. Axisymmetric numerical and analytical studies of the magnetorotational instability in a magnetized Taylor-Couette flow

    NASA Astrophysics Data System (ADS)

    Liu, Wei

    2007-08-01

    The magnetorotational instability (MRI) is probably the main cause of turbulence and accretion in sufficiently ionized astrophysical disks. However, despite much theoretical and computational work, the nonlinear saturation of MRI is imperfectly understood. In Chap. 2 and Chap. 3 of this thesis we present non-ideal magnetohydrodynamic simulations of the Princeton MRI experiment. In vertically infinite or periodic cylinders, MRI saturates in a resistive current-sheet with a significant reduction of the mean shear, and with poloidal circulation scaling as the square root of resistivity. Angular momentum transport scales as the reciprocal square root of viscosity but only weakly depends on resistivity. For finite cylinders with insulating end caps, a method for implementing the fully insulating boundary condition is introduced. MRI grows with a clear linear phase from small amplitudes at rates in good agreement with linear analysis. In the final state one inflowing "jet" opposite to the usual Ekman "jet" is found near the inner cylinder. The MRI enhances the angular momentum transport at saturation. Under proper conditions, our experimental facility is a good platform to show that MRI could be suppressed by a strong magnetic field. Recently, Hollerbach and Rüdiger have reported that MRI modes may grow at much reduced magnetic Reynolds number ( Re m ) and Lundquist number S in the presence of a helical background field, a current-free combination of axial and toroidal field. We have investigated these helical MRI modes in Chap. 4 and Chap. 5. In vertically infinite or periodic cylinders, resistive HMRI is a weakly destabilized hydrodynamic inertial oscillation propagating axially along the background Poynting flux. Growth rates are small, however, and require large axial currents. Furthermore, finite cylinders with insulating endcaps were shown to reduce the growth rate and to stabilize highly resistive, inviscid flows entirely, and the new mode is stable in Keplerian

  5. Comprehensive experimental and numerical analysis of instability phenomena in pump turbines

    NASA Astrophysics Data System (ADS)

    Gentner, Ch; Sallaberger, M.; Widmer, Ch; Bobach, B.-J.; Jaberg, H.; Schiffer, J.; Senn, F.; Guggenberger, M.

    2014-03-01

    The changes in the electricity market have led to changed requirements for the operation of pump turbines. Utilities need to change fast and frequently between pumping and generating modes and increasingly want to operate at off-design conditions for extended periods. Operation of the units in instable areas of the machine characteristic is not acceptable and may lead to self-excited vibration of the hydraulic system. In turbine operation of pump turbines unstable behaviour can occur at low load off-design operation close to runaway conditions (S-shape of the turbine characteristic). This type of instability may impede the synchronization of the machine in turbine mode and thus increase start-up and switch over times. A pronounced S-shaped instability can also lead to significant drop of discharge in the event of load rejection. Low pressure on the suction side and in the tail-race tunnel could cause dangerous separation of the water column. Understanding the flow features that lead to the instable behaviour of pump turbines is a prerequisite to the design of machines that can fulfil the growing requirements relating to operational flexibility. Flow simulation in these instability zones is demanding due to the complex and highly unsteady flow patterns. Only unsteady simulation methods are able to reproduce the governing physical effects in these operating regions. ANDRITZ HYDRO has been investigating the stability behaviour of pump turbines in turbine operation in cooperation with several universities using simulation and measurements. In order to validate the results of flow simulation of unstable operating points, the Graz University of Technology (Austria) performed detailed experimental investigations. Within the scope of a long term research project, the operating characteristics of several pump turbine runners have been measured and flow patterns in the pump turbine at speed no load and runaway have been examined by 2D Laser particle image velocimetry (PIV

  6. Conceptual model of hydrogeology in the Ozark Plateaus region during Pennsylvanian time

    SciTech Connect

    Brahana, J.V. )

    1993-03-01

    Recently completed studies of the Ozark Plateaus region of southern Missouri and northern Arkansas provide a conceptual framework for understanding current hydrogeology, and form the basis for numerical models that can be used to quantitatively assess flow and solute transport in the aquifers of this area. Three separate investigations were completed as part of the Regional Aquifer-Systems Analysis (RASA) program of the US Geological Survey during 1985--1993. Although the objectives of these RASA studies [Northern Midwest (NM) RASA, Gulf Coast (GC) RASA, and Central Midwest (CM) RASA] focused on recent hydrologic conditions, each study has contributed o increased understanding of the evolution of the hydrogeology of the region.

  7. Effects of chemical reactions on density-dependent fluid flow: On the numerical formulation and the development of instabilities

    SciTech Connect

    Freedman, Vicky L.; Ibaraki, Motomu

    2002-04-01

    A three-dimensional, reactive numerical flow model is developed that couples chemical reactions with density-dependent mass transport and fluid flow. The model includes equilibrium reactions for the aqueous species, kinetic reactions between the solid and aqueous phases, and full coupling of porosity and permeability changes that result from precipitation and dissolution reactions in the porous media. The coupled effects of the processes driving flow and the chemical reactions occurring during solute transport is studied using a carbonate system in fully saturated porous media. Results demonstrate that instability development is sensitive to the initial perturbation caused by density differences between the solute plume and the ambient ground water. If the initial perturbation is large, then it acts as a''trigger'' in the flow system that causes instabilities to develop in a planar reaction front. When permeability changes occur due to dissolution reactions occurring in the porous media, a reactive feedback loop is created by calcite dissolution and the mixed convective transport of the system. Although the feedback loop does not have a significant impact on plume shape, complex concentration distributions develop as a result of the instabilities generated in the flow system.

  8. Numerical and experimental studies of magnetic Rayleigh-Taylor instabilities in solid liners

    SciTech Connect

    Faehl, R.J.; Atchison, W.L.; Reinovsky, R.E.; Morgan, D.V.

    1997-12-31

    The authors have studied the nonlinear evolution of Magnetic Rayleigh-Taylor (MRT) instability in solid aluminum liners. Two-dimensional MHD modeling of the configuration have been performed using the measured current pulseform of the PEGASUS II capacitor bank as a magnetic field source. The liner configuration consisted of a solid aluminum annulus 20 mm long and 0.4 mm thick. The outer surface of the liner was machined with a sinusoidal ripple. Half of the liner was scribed with one wavelength, the other half with another wavelength. A series of experiments was conducted with a smooth surface and ripple wavelengths of 0.5, 0.75, and 2.0 mm. The peak-to-peak ripple amplitude was varied between 0.025 and 0.10 mm, but the ratio of the initial amplitude to wavelength never exceeded 7%. Diagnostics consisted of radiography, a center-conductor B-dot probe, and transmission line pulsed power diagnostics. The inner surfaces of the liners were coated with a thin film of gold (Au) to permit measurement of material deformation at various stages of the instability growth. Three radiographs, at different times, were taken for each experiment. The center-conductor magnetic probe was used to measure the magnetic field penetration through the liner. Initial comparisons between the data and the MHD calculations showed excellent agreement with respect to nonlinear instability growth and with the time of field rupture of the liner.

  9. Direct numerical simulations of type Ia supernovae flames II: The Rayleigh-Taylor instability

    SciTech Connect

    Bell, J.B.; Day, M.S.; Rendleman, C.A.; Woosley, S.E.; Zingale, M.

    2004-01-12

    A Type Ia supernova explosion likely begins as a nuclear runaway near the center of a carbon-oxygen white dwarf. The outward propagating flame is unstable to the Landau-Darrieus, Rayleigh-Taylor, and Kelvin-Helmholtz instabilities, which serve to accelerate it to a large fraction of the speed of sound. We investigate the Rayleigh-Taylor unstable flame at the transition from the flamelet regime to the distributed-burning regime, around densities of 10e7 gm/cc, through detailed, fully resolved simulations. A low Mach number, adaptive mesh hydrodynamics code is used to achieve the necessary resolution and long time scales. As the density is varied, we see a fundamental change in the character of the burning--at the low end of the density range the Rayleigh-Taylor instability dominates the burning, whereas at the high end the burning suppresses the instability. In all cases, significant acceleration of the flame is observed, limited only by the size of the domain we are able to study. We discuss the implications of these results on the potential for a deflagration to detonation transition.

  10. Numerical study of the ablative Richtmyer-Meshkov instability of laser-irradiated deuterium and deuterium-tritium targets

    NASA Astrophysics Data System (ADS)

    Marocchino, Alberto; Atzeni, Stefano; Schiavi, Angelo

    2010-11-01

    The Richtmyer-Meshkov instability (RMI) at the ablation front of laser-irradiated planar targets is investigated by two-dimensional numerical hydrodynamics simulations. The linear evolution of perturbations seeded either by surface roughness or target inhomogeneity is studied for perturbation wavelengths in the range 10≤λ≤400 μm and laser intensity 4×1012≤I≤4×1014 W/cm2 (with laser wavelength λlaser=0.35 μm). Thin and thick cryogenic deuterium or deuterium-tritium (DT) planar targets are considered. For targets irradiated at constant intensity, it is found that perturbations with wavelength below a given threshold perform damped oscillations, while perturbations above such a threshold are unstable and oscillate with growing amplitude. This is qualitatively in agreement with theoretical predictions by Goncharov et al. [Phys. Plasmas 13, 012702 (2006)], according to which ablation related processes stabilize perturbations with kDc≫1, where Dc is the distance between the ablation front and critical density for laser propagation. For kDc<1 a weakly growing Landau-Darrieus instability (LDI) is instead excited. The stability threshold increases substantially with laser intensity, given the dependence of Dc on laser intensity I (roughly Dc∝I, according to the present simulations). Direct-drive laser fusion targets are irradiated by time-shaped pulses, with a low intensity initial foot. In this case, perturbations with wavelengths below some threshold (about 10 μm, for typical ignition-class all-DT targets) are damped after an initial growth. In a thin target, initial perturbations, either damped or amplified by RMI and LDI, seed the subsequent Rayleigh-Taylor instability. Finally, it is shown that RMI growth of fusion targets can be reduced by using laser pulses including an initial adiabat-shaping picket (originally proposed to reduce the growth of Rayleigh-Taylor instability).

  11. A method for the direct numerical simulation of hypersonic boundary-layer instability with finite-rate chemistry

    SciTech Connect

    Marxen, Olaf; Magin, Thierry E.; Shaqfeh, Eric S.G.; Iaccarino, Gianluca

    2013-12-15

    A new numerical method is presented here that allows to consider chemically reacting gases during the direct numerical simulation of a hypersonic fluid flow. The method comprises the direct coupling of a solver for the fluid mechanical model and a library providing the physio-chemical model. The numerical method for the fluid mechanical model integrates the compressible Navier–Stokes equations using an explicit time advancement scheme and high-order finite differences. This Navier–Stokes code can be applied to the investigation of laminar-turbulent transition and boundary-layer instability. The numerical method for the physio-chemical model provides thermodynamic and transport properties for different gases as well as chemical production rates, while here we exclusively consider a five species air mixture. The new method is verified for a number of test cases at Mach 10, including the one-dimensional high-temperature flow downstream of a normal shock, a hypersonic chemical reacting boundary layer in local thermodynamic equilibrium and a hypersonic reacting boundary layer with finite-rate chemistry. We are able to confirm that the diffusion flux plays an important role for a high-temperature boundary layer in local thermodynamic equilibrium. Moreover, we demonstrate that the flow for a case previously considered as a benchmark for the investigation of non-equilibrium chemistry can be regarded as frozen. Finally, the new method is applied to investigate the effect of finite-rate chemistry on boundary layer instability by considering the downstream evolution of a small-amplitude wave and comparing results with those obtained for a frozen gas as well as a gas in local thermodynamic equilibrium.

  12. Numerical aspects in modeling high Deborah number flow and elastic instability

    NASA Astrophysics Data System (ADS)

    Kwon, Youngdon

    2014-05-01

    Investigating highly nonlinear viscoelastic flow in 2D domain, we explore problem as well as property possibly inherent in the streamline upwinding technique (SUPG) and then present various results of elastic instability. The mathematically stable Leonov model written in tensor-logarithmic formulation is employed in the framework of finite element method for spatial discretization of several representative problem domains. For enhancement of computation speed, decoupled integration scheme is applied for shear thinning and Boger-type fluids. From the analysis of 4:1 contraction flow at low and moderate values of the Deborah number (De) the solution with SUPG method does not show noticeable difference from the one by the computation without upwinding. On the other hand, in the flow regime of high De, especially in the state of elastic instability the SUPG significantly distorts the flow field and the result differs considerably from the solution acquired straightforwardly. When the strength of elastic flow and thus the nonlinearity further increase, the computational scheme with upwinding fails to converge and evolutionary solution does not become available any more. All this result suggests that extreme care has to be taken on occasions where upwinding is applied, and one has to first of all prove validity of this algorithm in the case of high nonlinearity. On the contrary, the straightforward computation with no upwinding can efficiently model representative phenomena of elastic instability in such benchmark problems as 4:1 contraction flow, flow over a circular cylinder and flow over asymmetric array of cylinders. Asymmetry of the flow field occurring in the symmetric domain, enhanced spatial and temporal fluctuation of dynamic variables and flow effects caused by extension hardening are properly described in this study.

  13. Gas Evolution Dynamics in Godunov-Type Schemes and Analysis of Numerical Shock Instability

    NASA Technical Reports Server (NTRS)

    Xu, Kun

    1999-01-01

    In this paper we are going to study the gas evolution dynamics of the exact and approximate Riemann solvers, e.g., the Flux Vector Splitting (FVS) and the Flux Difference Splitting (FDS) schemes. Since the FVS scheme and the Kinetic Flux Vector Splitting (KFVS) scheme have the same physical mechanism and similar flux function, based on the analysis of the discretized KFVS scheme the weakness and advantage of the FVS scheme are closely observed. The subtle dissipative mechanism of the Godunov method in the 2D case is also analyzed, and the physical reason for shock instability, i.e., carbuncle phenomena and odd-even decoupling, is presented.

  14. Collisional interchange instability 2. A comparison of the numerical simulations with the in situ experimental data

    SciTech Connect

    Kelley, M.C.; Seyler, C.E.; Zargham, S. )

    1987-09-01

    A two-dimensional model applicable to F region ionosphere plasma instabilities has been developed and described in a comparison paper. Here the authors apply the model to equatorial F region irregularities and in particular test the model against rocket and satellite data. As a diagnostic they create simulated data sets similar to the one-dimensional measurements of plasma density performed by space probes and take the Fourier transform of these data in the same manner and used by the space experimentors. To their knowledge this is the first attempt to directly relate a simulation to in situ data in this way, diagnostic technique which resolves the ambiguity inherent to one-dimensional measurements. The comparison is in excellent agreement with the in situ data. In particular, they have been able to resolve the apparent differences between satellite quasi-sinusoidal observations with the shocklike observations made on rockets. Unlike previous simulations of this phenomenon they find an inherent anisotropy in the instability development which is mirrored in the in situ data. They also present evidence that the shallow spectral slopes which often characterize spread F rocket spectra near the F peak may be due to a change in the angle between the rocket velocity vector and the characteristic directions in the medium.

  15. Continental accretion: From oceanic plateaus to allochthonous terranes

    USGS Publications Warehouse

    Ben-Avraham, Z.; Nur, A.; Jones, D.; Cox, A.

    1981-01-01

    Some of the regions of the anomalously high sea-floor topography in today's oceans may be modern allochthonous terranes moving with their oceanic plates. Fated to collide with and be accreted to adjacent continents, they may create complex volcanism, cut off and trap oceanic crust, and cause orogenic deformation. The accretion of plateaus during subduction of oceanic plates may be responsible for mountain building comparable to that produced by the collision of continents. Copyright ?? 1981 AAAS.

  16. Shock tube experiments and numerical simulation of the single-mode, three-dimensional Richtmyer-Meshkov instability

    NASA Astrophysics Data System (ADS)

    Long, C. C.; Krivets, V. V.; Greenough, J. A.; Jacobs, J. W.

    2009-11-01

    A vertical shock tube is used to perform experiments in which an interface is formed using opposed flows of air and SF6. A three-dimensional single-mode perturbation is created by the periodic vertical motion of the gases within the shock tube. Richtmyer-Meshkov instability is produced by an impulsive acceleration by a weak shock wave (Ms=1.2). Planar laser induced fluorescence produces still images, and planar Mie scattering produces movies of the experiment. A three-dimensional numerical simulation of this experiment utilizing the Eulerian adaptive mesh refinement code, RAPTOR, was also conducted. Good agreement is obtained between experiments and the simulations. However, existing late time models, which have a 1/t dependence, disagree with measurements of the late time instability development. In contrast, both the experiments and simulation suggest a t-0.54 late time dependence for the overall growth rate. Comparisons with individual bubble and spike velocities show the bubbles appear to decay approximately at 1/t and the spikes to decay at a much slower rate of t-0.38.

  17. Role of Magnetic Field Strength and Numerical Resolution in Simulations of the Heat-flux-driven Buoyancy Instability

    NASA Astrophysics Data System (ADS)

    Avara, Mark J.; Reynolds, Christopher S.; Bogdanović, Tamara

    2013-08-01

    The role played by magnetic fields in the intracluster medium (ICM) of galaxy clusters is complex. The weakly collisional nature of the ICM leads to thermal conduction that is channeled along field lines. This anisotropic heat conduction profoundly changes the instabilities of the ICM atmosphere, with convective stabilities being driven by temperature gradients of either sign. Here, we employ the Athena magnetohydrodynamic code to investigate the local non-linear behavior of the heat-flux-driven buoyancy instability (HBI) relevant in the cores of cooling-core clusters where the temperature increases with radius. We study a grid of two-dimensional simulations that span a large range of initial magnetic field strengths and numerical resolutions. For very weak initial fields, we recover the previously known result that the HBI wraps the field in the horizontal direction, thereby shutting off the heat flux. However, we find that simulations that begin with intermediate initial field strengths have a qualitatively different behavior, forming HBI-stable filaments that resist field-line wrapping and enable sustained vertical conductive heat flux at a level of 10%-25% of the Spitzer value. While astrophysical conclusions regarding the role of conduction in cooling cores require detailed global models, our local study proves that systems dominated by the HBI do not necessarily quench the conductive heat flux.

  18. Numerical investigations and analysis on the effects of geometrical parameters on the group velocity of transverse magnetic pump mode and free electron laser instability

    SciTech Connect

    Sharma, B. S.; Jaiman, N. K.

    2009-02-15

    In this paper, we have numerically investigated the effects of various geometrical parameters of a backward wave oscillator, filled with a magnetized plasma of uniform density and driven by a mild relativistic solid electron beam, on the instability growth rate R{sub 0} of a seeded free electron laser. On changing mean radius corrugation amplitude h and corrugation period z{sub 0} of backward wave oscillator; the ponderomotive potential of space charge wave changes. This in turn, changes the coupling strength of TM mode with negative beam space charge mode and hence the growth rate of parametric instability of free electron laser. A dispersion relation is derived and numerically solved for various geometrical parameters of backward wave oscillator and beam profile. A relation for {Gamma} is also derived and computed numerically. The instability growth scales directly to the square root of beam density and inversely as seven power of relativistic gamma factor {gamma}{sub 0}.

  19. Numerical Simulation of Multi-Material Mixing in an Inclined Interface Richtmyer-Meshkov Instability

    NASA Astrophysics Data System (ADS)

    Subramaniam, Akshay; Lele, Sanjiva

    2015-11-01

    The Richtmyer-Meshkov instability arises when a shock wave interacts with an interface separating two fluids. In this work, high fidelity simulations of shock induced multi-material mixing between N2 and CO2 in a shock tube are performed for a Mach 1.55 shock interacting with a planar material interface that is inclined with respect to the shock propagation direction. In the current configuration, unlike the classical perturbed flat interface case, the evolution of the interface is non-linear from early time onwards. Our previous simulations of this problem at multiple spatial resolutions have shown that very small 3D perturbations have a large effect on vortex breakdown mechanisms and hence fine scale turbulence. We propose a comparison of our simulations to the experiments performed at the Georgia Tech Shock Tube and Advanced Mixing Laboratory (STAML). Results before and after reshock of the interface will be shown. Results from simulations of a second case with a more complex initial interface will also be presented. Simulations shown are conducted with an extended version of the Miranda solver developed by Cook et al. (2007) which combines high-order compact finite differences with localized non-linear artificial properties for shock and interface capturing. This research is part of the Blue Waters sustained-petascale computing project, which is supported by the National Science Foundation (awards OCI-0725070 and ACI-1238993) and the state of Illinois.

  20. Numerical 3D Hydrodynamics Study of Gravitational Instabilities in a Circumbinary Disk

    NASA Astrophysics Data System (ADS)

    Desai, Karna Mahadev; Steiman-Cameron, Thomas Y.; Michael, Scott; Cai, Kai; Durisen, Richard H.

    2016-01-01

    We present a 3D hydrodynamical study of gravitational instabilities (GIs) in a circumbinary protoplanetary disk around a Solar mass star and a brown dwarf companion (0.02 M⊙). GIs can play an important, and at times dominant, role in driving the structural evolution of protoplanetary disks. The reported simulations were performed employing CHYMERA, a radiative 3D hydrodynamics code developed by the Indiana University Hydrodynamics Group. The simulations include disk self-gravity and radiative cooling governed by realistic dust opacities. We examine the role of GIs in modulating the thermodynamic state of the disks, and determine the strengths of GI-induced density waves, non-axisymmetric density structures, radial mass transport, and gravitational torques. The principal goal of this study is to determine how the presence of the companion affects the nature and strength of GIs. Results are compared with a parallel simulation of a protoplanetary disk without the presence of the brown dwarf binary companion. We detect no fragmentation in either disk. A persistent vortex forms in the inner region of both disks. The vortex seems to be stabilized by the presence of the binary companion.

  1. Geodynamic models of terrane accretion: Testing the fate of island arcs, oceanic plateaus, and continental fragments in subduction zones

    NASA Astrophysics Data System (ADS)

    Tetreault, J. L.; Buiter, S. J. H.

    2012-08-01

    Crustal growth at convergent margins can occur by the accretion of future allochthonous terranes (FATs), such as island arcs, oceanic plateaus, submarine ridges, and continental fragments. Using geodynamic numerical experiments, we demonstrate how crustal properties of FATs impact the amount of FAT crust that is accreted or subducted, the type of accretionary process, and the style of deformation on the overriding plate. Our results show that (1) accretion of crustal units occurs when there is a weak detachment layer within the FAT, (2) the depth of detachment controls the amount of crust accreted onto the overriding plate, and (3) lithospheric buoyancy does not prevent FAT subduction during constant convergence. Island arcs, oceanic plateaus, and continental fragments will completely subduct, despite having buoyant lithospheric densities, if they have rheologically strong crusts. Weak basal layers, representing pre-existing weaknesses or detachment layers, will either lead to underplating of faulted blocks of FAT crust to the overriding plate or collision and suturing of an unbroken FAT crust. Our experiments show that the weak, ultramafic layer found at the base of island arcs and oceanic plateaus plays a significant role in terrane accretion. The different types of accretionary processes also affect deformation and uplift patterns in the overriding plate, trench migration and jumping, and the dip of the plate interface. The resulting accreted terranes produced from our numerical experiments resemble observed accreted terranes, such as the Wrangellia Terrane and Klamath Mountain terranes in the North American Cordilleran Belt.

  2. Numerical analysis of MHD instability suppression in a double gas puff

    SciTech Connect

    Glazyrin, Igor V.; Diyankov, Oleg V.; Karlykhanov, Nikolai G.; Koshelev, Sergei V.

    1997-05-05

    The process of double gas-puff implosion is numerically analyzed. Various liner configurations are considered: 1) the outer cascade is ten times lighter then the inner one, 2) the masses of cascades are equal, 3) the outer cascade is ten times heavier then the inner one. The mechanism of liner implosions stabilization is discussed. The comparison of two-cascade and single-cascade schemes is performed. The role of magnetic field penetration at the initial time into liner's body is also discussed.

  3. Numerical simulations of instabilities in the implosion process of inertial confined fusion in 2D cylindrical coordinates

    NASA Astrophysics Data System (ADS)

    Yong, Heng; Zhai, ChuanLei; Jiang, Song; Song, Peng; Dai, ZhenSheng; Gu, JianFa

    2016-01-01

    In this paper, we introduce a multi-material arbitrary Lagrangian and Eulerian method for the hydrodynamic radiative multi-group diffusion model in 2D cylindrical coordinates. The basic idea in the construction of the method is the following: In the Lagrangian step, a closure model of radiation-hydrodynamics is used to give the states of equations for materials in mixed cells. In the mesh rezoning step, we couple the rezoning principle with the Lagrangian interface tracking method and an Eulerian interface capturing scheme to compute interfaces sharply according to their deformation and to keep cells in good geometric quality. In the interface reconstruction step, a dual-material Moment-of-Fluid method is introduced to obtain the unique interface in mixed cells. In the remapping step, a conservative remapping algorithm of conserved quantities is presented. A number of numerical tests are carried out and the numerical results show that the new method can simulate instabilities in complex fluid field under large deformation, and are accurate and robust.

  4. Geological and numerical modelling of sinkholes induced by instability of man-made underground caves

    NASA Astrophysics Data System (ADS)

    Lollino, P.; Margiotta, S.; Parise, M.

    2012-04-01

    An increasing number of areas in southern Italy are being interested by sinkholes related to the presence of artificial cavities, and in particular of underground quarries. Aimed at defining the most proper criteria for the study of such phenomena, the case of Cutrofiano (Apulia, SE Italy) is here presented as an example. This territory is well known for the underground quarrying activity of soft calcarenitic rocks that caused in the last decades widespread phenomena of subsidence at the surface, with extensive damages and problems to the main communication routes in the area. A very complex and intricate network of subterranean galleries is present underground, for a total development on the order of several tens of kilometres. The Gravina Calcarenite Formation, that is the object of the underground quarries, consists of whitish calcarenite with an upper greenish clayey-sandy interval rich in fossils. This formation is overlain by grey sandy clays (Subapennine Clays) grading upward to fossil-rich sands (Brindisi Sands). In recent years, clay mining has been resumed at the ground surface, following to the opening of a cement factory, and resulting in the realization of wide open quarries. The local sedimentary sequence is closed by terraced calcarenite deposits rich in ostreids, holding a phreatic groundwater body and overlain by a shallow soil cover. The results of detailed geological, geomorphological and geomechanical surveys, supported by laboratory tests, show that in the Cutrofiano area the thickness of the stratigraphical succession, the depth of the underground galleries, the structural conditions of the rock mass, and the failure mechanisms observed within the quarries are variable. The geological model reconstructed represents the base for the numerical simulations, which are aimed at defining the eventual mechanisms of rock failure, up to the formation of the sinkhole. In particular, two different geological settings have been considered: the first is

  5. Numerical and physical instabilities in massively parallel LES of reacting flows

    NASA Astrophysics Data System (ADS)

    Poinsot, Thierry

    LES of reacting flows is rapidly becoming mature and providing levels of precision which can not be reached with any RANS (Reynolds Averaged) technique. In addition to the multiple subgrid scale models required for such LES and to the questions raised by the required numerical accurcay of LES solvers, various issues related the reliability, mesh independence and repetitivity of LES must still be addressed, especially when LES is used on massively parallel machines. This talk discusses some of these issues: (1) the existence of non physical waves (known as `wiggles' by most LES practitioners) in LES, (2) the effects of mesh size on LES of reacting flows, (3) the growth of rounding errors in LES on massively parallel machines and more generally (4) the ability to qualify a LES code as `bug free' and `accurate'. Examples range from academic cases (minimum non-reacting turbulent channel) to applied configurations (a sector of an helicopter combustion chamber).

  6. Numerical 3D study of FZ growth: dependence on growth parameters and melt instability

    NASA Astrophysics Data System (ADS)

    Ratnieks, G.; Muižnieks, A.; Mühlbauer, A.; Raming, G.

    2001-08-01

    Three-dimensional modelling of the floating zone (needle-eye) crystal growth process is carried out to analyse numerically the stability of the melt flow and the influence of the crystal rotation rate and inductor slit width on the 3D flow field and on the grown crystal resistivity. The unsteadiness of the melt is simulated and it is found that for the considered growth parameters a steady-state flow can be a reasonable approximation to the unsteady melt motion. The parametric studies have shown that increasing the rotation rate essentially changes the flow pattern and weakens the rotational striations, while the inductor slit width has a more local influence on these characteristics.

  7. Numerical Study of Three-dimensional Spatial Instability of a Supersonic Flat Plate Boundary Layer

    NASA Technical Reports Server (NTRS)

    Maestrello, Lucio; Bayliss, A.; Krishnan, R.

    1989-01-01

    The behavior of spatially growing three-dimensional waves in a supersonic boundary layer was studied numerically by solving the complete Navier-Stokes equations. Satisfactory comparison with linear parallel and non-parallel stability theories, and experiment are obtained when a small amplitude inflow disturbance is used. The three-dimensional unsteady Navier-Stokes equations are solved by a finite difference method which is fourth-order and second-order accurate in the convection and viscous terms respectively, and second-order accurate in time. Spanwise periodicity is assumed. The inflow disturbance is composed of eigenfunctions from linear stability theory. By increasing the amplitude of the inflow disturbance, nonlinear effects in the form of a relaxation type oscillation of the time signal of rho(u) are observed.

  8. Three-dimensional instabilities of pantographic sheets with parabolic lattices: numerical investigations

    NASA Astrophysics Data System (ADS)

    Scerrato, Daria; Giorgio, Ivan; Rizzi, Nicola Luigi

    2016-06-01

    In this paper, we determine numerically a large class of equilibrium configurations of an elastic two-dimensional continuous pantographic sheet in three-dimensional deformation consisting of two families of fibers which are parabolic prior to deformation. The fibers are assumed (1) to be continuously distributed over the sample, (2) to be endowed of bending and torsional stiffnesses, and (3) tied together at their points of intersection to avoid relative slipping by means of internal (elastic) pivots. This last condition characterizes the system as a pantographic lattice (Alibert and Della Corte in Zeitschrift für angewandte Mathematik und Physik 66(5):2855-2870, 2015; Alibert et al. in Math Mech Solids 8(1):51-73, 2003; dell'Isola et al. in Int J Non-Linear Mech 80:200-208, 2016; Int J Solids Struct 81:1-12, 2016). The model that we employ here, developed by Steigmann and dell'Isola (Acta Mech Sin 31(3):373-382, 2015) and first investigated in Giorgio et al. (Comptes rendus Mecanique 2016, doi: 10.1016/j.crme.2016.02.009), is applicable to fiber lattices in which three-dimensional bending, twisting, and stretching are significant as well as a resistance to shear distortion, i.e., to the angle change between the fibers. Some relevant numerical examples are exhibited in order to highlight the main features of the model adopted: In particular, buckling and post-buckling behaviors of pantographic parabolic lattices are investigated. The fabric of the metamaterial presented in this paper has been conceived to resist more effectively in the extensional bias tests by storing more elastic bending energy and less energy in the deformation of elastic pivots: A comparison with a fabric constituted by beams which are straight in the reference configuration shows that the proposed concept is promising.

  9. A joint experimental and numerical study of mechanisms associated to instability of partial cavitation on two-dimensional hydrofoil

    NASA Astrophysics Data System (ADS)

    Leroux, Jean-Baptiste; Coutier-Delgosha, Olivier; Astolfi, Jacques André

    2005-05-01

    The present work was carried out in the scope of a numerical-experimental collaborative research program, whose main objective is to understand the mechanisms of instabilities in partial cavitating flow. Experiments were conducted in the configuration of a rectangular foil located in a cavitation tunnel. Partial cavitation was investigated by multipoint wall-pressure measurements together with lift and drag measurements and numerical videos. The computations were conducted on two-dimensional hydrofoil section and are based on a single fluid model of cavitation: the liquid/vapor mixture is considered as a homogeneous fluid whose composition is regulated by a barotropic state law. The algorithm of resolution is derived from the SIMPLE approach, modified to take into account the high compressibility of the medium. Several physical features were pointed out by this joint approach. Particularly two distinct cavity self-oscillation dynamics characterized by two different frequencies (dynamics 1 and dynamics 2) were obtained experimentally and numerically at the angles of incidence of 6° and 8°. In both cases, the reentrant jet was found to be mainly responsible for the cavity breakdown. Dynamics 2 corresponds to the "classical" cavity breakdown and resulting cloud cavitation. A more complex flow pattern was evidenced for dynamics 1. In this case the growth/breakdown cycle of the cavity was observed at a lower Strouhal number (˜0.07/0.09) than dynamics 2 (˜0.3). Moreover, the mechanism is composed of two successive steps: (i) an interaction between the reentrant jet and the cavity interface in the closure region leading to the periodic shedding of secondary cavitation clouds before the main cloud detachment occurs, and (ii) a shock wave induced by the collapse of the main cloud, which influences the growth of the residual cavity.

  10. Direct numerical simulations of trailing-edge noise generated by boundary-layer instabilities

    NASA Astrophysics Data System (ADS)

    Sandberg, R. D.; Sandham, N. D.; Joseph, P. F.

    2007-07-01

    Direct numerical simulations (DNS) are conducted of noise generated at an infinitely thin trailing edge (TE). The aim is to predict the far-field sound and the near-field hydrodynamics, thereby providing an insight into the physical mechanisms of sound generation at airfoil TEs and potentially helping to validate acoustic theories. One of the theories widely used is the classical inviscid theory of Amiet, where the far-field sound can be evaluated in closed form if the convecting surface pressure spectrum upstream of the TE is known. For the first time, data from DNS including viscous effects are compared to the classical inviscid TE noise theory. In the present investigation, Tollmien-Schlichting waves are introduced close to the inflow boundary. The disturbances propagate downstream producing pressure fluctuations at the TE. In conducting two-dimensional DNS the theoretical method requires modification to account for the radiation of the total pressure difference in two dimensions only, as opposed to the three-dimensional sound radiation originally considered by Amiet. The modified theoretical analysis and a comparison between DNS and theoretical results are presented, scrutinizing the assumptions made in the derivation. Amiet's surface pressure jump transfer function is found to predict the scattered pressure field accurately. Directivity plots of DNS data show that viscous effects appear to smear individual lobes and that a downstream pointing lobe is present at higher Mach number which is attributed to an additional wake source.

  11. Numerical study of chiral plasma instability within the classical statistical field theory approach

    NASA Astrophysics Data System (ADS)

    Buividovich, P. V.; Ulybyshev, M. V.

    2016-07-01

    We report on a numerical study of real-time dynamics of electromagnetically interacting chirally imbalanced lattice Dirac fermions within the classical statistical field theory approach. Namely, we perform exact simulations of the real-time quantum evolution of fermionic fields coupled to classical electromagnetic fields, which are in turn coupled to the vacuum expectation value of the fermionic electric current. We use Wilson-Dirac Hamiltonian for fermions, and noncompact action for the gauge field. In general, we observe that the backreaction of fermions on the electromagnetic field prevents the system from acquiring chirality imbalance. In the case of chirality pumping in parallel electric and magnetic fields, the electric field is screened by the produced on-shell fermions and the accumulation of chirality is hence stopped. In the case of evolution with initially present chirality imbalance, axial charge tends to transform to helicity of the electromagnetic field. By performing simulations on large lattices we show that in most cases this decay process is accompanied by the inverse cascade phenomenon, which transfers energy from short-wavelength to long-wavelength electromagnetic fields. In some simulations, however, we observe a very clear signature of inverse cascade for the helical magnetic fields that is not accompanied by the axial charge decay. This suggests that the relation between the inverse cascade and axial charge decay is not as straightforward as predicted by the simplest form of anomalous Maxwell equations.

  12. Partition function zeros and magnetization plateaus of the spin-1 Ising-Heisenberg diamond chain

    NASA Astrophysics Data System (ADS)

    Hovhannisyan, V. V.; Ananikian, N. S.; Kenna, R.

    2016-07-01

    We study the properties of the generalized spin-1 Ising-Heisenberg model on a diamond chain, which can be considered as a theoretical model for the homometallic magnetic complex [Ni3(C4H2O4)2 -(μ3 - OH) 2(H2O)4 ] n ṡ(2H2 O) n. The model possesses a large variety of ground-state phases due to the presence of biquadratic and single-ion anisotropy parameters. Magnetization and quadrupole moment plateaus are observed at one- and two-thirds of the saturation value. The distributions of Yang-Lee and Fisher zeros are studied numerically for a variety of values of the model parameters. The usual value σ = -1/2 alongside an unusual value σ = -2/3 ​is determined for the Yang-Lee edge singularity exponents.

  13. Numerical simulations of two-fluid turbulent mixing at large density ratios and applications to the Rayleigh-Taylor instability.

    PubMed

    Livescu, D

    2013-11-28

    A tentative review is presented of various approaches for numerical simulations of two-fluid gaseous mixtures at high density ratios, as they have been applied to the Rayleigh-Taylor instability (RTI). Systems exhibiting such RTI behaviour extend from atomistic sizes to scales where the continuum approximation becomes valid. Each level of description can fit into a hierarchy of theoretical models and the governing equations appropriate for each model, with their assumptions, are presented. In particular, because the compressible to incompressible limit of the Navier-Stokes equations is not unique and understanding compressibility effects in the RTI critically depends on having the appropriate basis for comparison, two relevant incompressible limits are presented. One of these limits has not been considered before. Recent results from RTI simulations, spanning the levels of description presented, are reviewed in connection to the material mixing problem. Owing to the computational limitations, most in-depth RTI results have been obtained for the incompressible case. Two such results, concerning the asymmetry of the mixing and small-scale anisotropy anomaly, as well as the possibility of a mixing transition in the RTI, are surveyed. New lines for further investigation are suggested and it is hoped that bringing together such diverse levels of description may provide new ideas and increased motivation for studying such flows. PMID:24146007

  14. Sn Attenuation in the Iranian and Tibetan Plateaus

    NASA Astrophysics Data System (ADS)

    Ku, W.; Kaviani, A.; Bao, X.; Christopher, J.; Sandvol, E. A.

    2014-12-01

    The Iranian and Tibetan plateaus are both part of the Alpine-Himalayan mountain belt and both formed as a result of the continental collision (between Arabian plate, Indian plate and Eurasia plates). The two plateaus illustrate many of the same processes but at different stages of development. The Iranian plateau is a much less mature continental plateau than the Tibetan plateau with about half the elevation and a much thinner crust. In order to better understand the nature of the lithosphere mantle and origin of the measure seismic velocity anomalies we have made detailed measurements of the uppermost mantle attenuation using the high frequency regional phase Sn. In order to measure Sn attenuation. We have collected a large data set consisting of 18 years (1995-2012) of waveforms recorded by 305 permanent and temporary stations in Iran and 3 years (2007-2009) of waveforms recorded by 136 permanent and temporary stations in Tibet separately. We used a bandpass filter (0.1-0.5Hz) to identify efficient longer period Sn phases. In order to determine Sn Q we applied a Two Station Method (TSM) and Reverse Two Station Method (RTM) to eliminate the source effects. We have used the LSQR algorithm to tomographically map Sn attenuation tomography across the Iranian plateau. Our results show moderately low Q values beneath the Iranian plateau (~250) and high Q values beneath the south Caspian sea (~400) and Arabian shield (~400). The pattern we observe is roughly consistent with previous studies. We also observe high Q values beneath the Zagros mountains (~450) that is consistent with the Arabian plate underthrusting beneath the Eurasia plate. In Tibet we find high Q values beneath the Qaidam basin (~500) and low Q values beneath much of the Kunlun fault (~200). We also observe high Q values beneath the Bangong-Nujiang suture (~400) which can be explained that the Indian plate is underthrusting beneath the Eurasia plate. Comparing the Sn Q values beneath the Iranian and

  15. Geoid height versus topography for oceanic plateaus and swells

    NASA Technical Reports Server (NTRS)

    Sandwell, David T.; Mackenzie, Kevin R.

    1989-01-01

    Gridded geoid height data (Marsh et al.l, 1986) and gridded bathymetry data (Van Wykhouse, 1973) are used to estimate the average compensation depths of 53 oceanic swells and plateaus. The relationship between geoid height and topography is examined using Airy and thermal compensation models. It is shown that geoid height is linearly related to topography between wavelengths of 400 and 4000 m as predicted by isostatic compensation models. The geoid/topography ratio is dependent on the average depth of compensation. The intermediate geoid/topography ratios of most thermal swells are interpreted as a linear combination of the decaying thermal swell signature and that of the persisting Airy-compensated volcanic edifice.

  16. On the numerical simulation of the ablative Rayleigh-Taylor instability in laser-driven ICF targets using the FastRad3D code

    NASA Astrophysics Data System (ADS)

    Bates, Jason; Schmitt, Andrew; Zalesak, Steve

    2015-11-01

    The ablative Rayleigh-Taylor (RT) instability is a key factor in the performance of directly-drive inertial-confinement-fusion (ICF) targets. Although this subject has been studied for quite some time, the accurate simulation of the ablative RT instability has proven to be a challenging task for many radiation hydrodynamics codes, particularly when it comes to capturing the ablatively-stabilized region of the linear dispersion spectrum and modeling ab initio perturbations. In this poster, we present results from recent two-dimensional numerical simulations of the ablative RT instability that were performed using the Eulerian code FastRad3D at the U.S. Naval Research Laboratory. We consider both planar and spherical geometries, low and moderate-Z target materials, different laser wavelengths and where possible, compare our findings with experiment data, linearized theory and/or results from other radiation hydrodynamics codes. Overall, we find that FastRad3D is capable of simulating the ablative RT instability quite accurately, although some uncertainties/discrepancies persist. We discuss these issues, as well as some of the numerical challenges associated with modeling this class of problems. Work supported by U.S. DOE/NNSA.

  17. Magnetization plateaus and phase diagrams of the Ising model on the Shastry-Sutherland lattice

    NASA Astrophysics Data System (ADS)

    Deviren, Seyma Akkaya

    2015-11-01

    The magnetization properties of a two-dimensional spin-1/2 Ising model on the Shastry-Sutherland lattice are studied within the effective-field theory (EFT) with correlations. The thermal behavior of the magnetizations is investigated in order to characterize the nature (the first- or second-order) of the phase transitions as well as to obtain the phase diagrams of the model. The internal energy, specific heat, entropy and free energy of the system are also examined numerically as a function of the temperature in order to confirm the stability of the phase transitions. The applied field dependence of the magnetizations is also examined to find the existence of the magnetization plateaus. For strong enough magnetic fields, several magnetization plateaus are observed, e.g., at 1/9, 1/8, 1/3 and 1/2 of the saturation. The phase diagrams of the model are constructed in two different planes, namely (h/|J|, |J‧|/|J|) and (h/|J|, T/|J|) planes. It was found that the model exhibits first- and second-order phase transitions; hence tricitical point is also observed in additional to the zero-temperature critical point. Moreover the Néel order (N), collinear order (C) and ferromagnetic (F) phases are also found with appropriate values of the system parameters. The reentrant behavior is also obtained whenever model displays two Néel temperatures. These results are compared with some theoretical and experimental works and a good overall agreement has been obtained.

  18. Role of Magnetic Field Strength and Numerical Resolution in Simulations of the Heat-flux Driven Buoyancy Instability

    NASA Astrophysics Data System (ADS)

    Avara, Mark J.; Reynolds, C. S.; Bogdanovic, T.

    2013-04-01

    The role played by magnetic fields in the intracluster medium (ICM) of galaxy clusters is complex. The weakly collisional nature of the ICM leads to thermal conduction that is channeled along field lines. This anisotropic heat conduction profoundly changes the stability of the ICM atmosphere, with convective stabilities being driven by temperature gradients of either sign. We employ the Athena magnetohydrodynamic code to investigate the local non-linear behavior of the heat-flux driven buoyancy instability (HBI), relevant in the cores of cooling-core clusters where the temperature increases with radius. We study a grid of 2-d simulations that span a large range of initial magnetic field strengths and numerical resolutions. For very weak initial fields, we recover the previously known result that the HBI wraps the field in the horizontal direction thereby shutting off the heat flux. However, we find that simulations which begin with intermediate initial field strengths have a qualitatively different behavior, forming HBI-stable filaments that resist field-line wrapping and enable sustained vertical conductive heat flux at a level of 10-25% of the Spitzer value. We explain the presence and persistence of these filaments in terms of the linear stability of the HBI and the total energetics of the plasma. A complimentary 3-d simulation of high resolution confirms the presence of sustained filaments and shows they can be formed in the ideal MHD regime, even without anisotropic viscosity, previously thought to be necessary. While astrophysical conclusions regarding the role of conduction in cooling cores require detailed global models and a better understanding of conduction in the ICM, our local study proves that systems dominated by HBI do not necessarily quench the conductive heat flux.

  19. Modelling Rayleigh Taylor instability of a sedimenting suspension of several thousand circular particles in a direct numerical simulation

    NASA Astrophysics Data System (ADS)

    Pan, T. W.; Joseph, D. D.; Glowinski, R.

    2001-05-01

    In this paper we study the sedimentation of several thousand circular particles in two dimensions using the method of distributed Lagrange multipliers for solid liquid flow. The simulation gives rise to fingering which resembles Rayleigh Taylor instabilities. The waves have a well-defined wavelength and growth rate which can be modelled as a conventional Rayleigh Taylor instability of heavy fluid above light. The heavy fluid is modelled as a composite solid liquid fluid with an effective composite density and viscosity. Surface tension cannot enter this problem and the characteristic shortwave instability is regularized by the viscosity of the solid liquid dispersion. The dynamics of the Rayleigh Taylor instability are studied using viscous potential flow, generalizing work of Joseph, Belanger & Beavers (1999) to a rectangular domain bounded by solid walls; an exact solution is obtained.

  20. Development of a 3D numerical model to evaluate the Stromboli NW flank instability in relation to magma intrusion

    NASA Astrophysics Data System (ADS)

    Apuani, T.; Merri, A.

    2009-04-01

    A stress-strain analysis of the Stromboli volcano was performed using a three-dimensional explicit finite difference numerical code (FLAC 3D, ITASCA, 2005), to evaluate the effects associated to the presence of magma pressure in magmatic conduit and to foresee the evolution of the magmatic feeding complex. The simulations considered both the ordinary state for the Stromboli, characterized by a partial fill of the active dyke with regular emission of gas and lava fountains and the paroxysmal conditions observed during the March 2007's eruptive crisis, with the magma level in the active dyke reaching the topographic surface along the Sciara del Fuoco depression. The modeling contributes to identify the most probable directions of propagation of new dikes, and the effects of their propagation on the stability of the volcano edifice. The numerical model extends 6 x 6 x 2.6 km3, with a mesh resolution of 100 m, adjusting the grid to fit the shape of the object to be modeled. An elasto-plastic constitutive law was adopted and an homogeneous Mohr-Coulomb strength criterion was chosen for the volcanic cone, assuming one lithotechnical unit (alternation of lava and breccia layers "lava-breccia unit"- Apuani et al 2005). The dykes are represented as discontinuities of the grid, and are modeled by means of interfaces. The magmatic pressure is imposed to the model as normal pressure applied on both sides of the interfaces. The magmastatic pressure was calculated as Pm=d•h, where d is the magma unit weight assumed equal to 25 KN/m3, and h (m) is the height of the magma column. Values of overpressure between 0 and 1 MPa were added to simulate the paroxysmal eruption. The simulation was implemented in successive stages, assuming the results of the previous stages as condition for the next one. A progressive propagation of the dike was simulated, in accordance with the stress conditions identified step by step, and in accordance with the evidences detected by in situ survey, and

  1. Line defects and quantum Hall plateaus in graphene.

    PubMed

    Dal Lago, V; Torres, L E F Foa

    2015-04-15

    Line defects in graphene can be either tailored-growth or arise naturally and are at the center of many discussions. Here we study the multiterminal conductance of graphene with an extended line defect in the quantum Hall regime analyzing the effects of the geometry of the setup, disorder and strain on the quantum Hall plateaus. We show that the defect turns out to affect the local and non-local conductance in very different ways depending on the geometrical configuration. When the defect is parallel to the sample edges one gets an equivalent circuit formed by parallel resistors. In contrast, when the defect bridges opposite edges, the Hall conductance may remain unaltered depending on the geometry of the voltage/current probes. The role of disorder, strain and the microscopic details of the defect in our results is also discussed. We show that the defect provides a realization of the electrical analog of an optical beam splitter. Its peculiar energy dependent inter-edge transmission allows it to be turned on or off at will and it may be used for routing the chiral edge states. PMID:25786352

  2. Report on the geology of the high plateaus of Utah

    USGS Publications Warehouse

    Dutton, C.E.; Powell, John Wesley

    1880-01-01

    In the year 1874 my kind friend Prof. J.W. Powell proposed to me that I should undertake, under his direction, the study of a large volcanic tract in the Territory of Utah, provided the consent of proper authority could be entertained. Distrusting my own fitness for the work, I felt that it would be better for him if his proposals were thankfully declined. In 1875, however, he renewed the proposition in such a friendly and complimentary manner that a refusal seemed ungracious. He therefore laid the matter before the Secretary of War, the General of the Army, and the Chief of Ordnance, all of whom gave their cordial approbation; and by order of the Ward Department I was detailed for duty in connection with the survey of the Rocky Mountain Region in charge of Professor Powell. The field which he assigned me to study was the District of the High Plateaus, and the investigations were made during the summers of 1875, 1876, and 1877. The preparation of a report or monograph upon the district has several times between interrupted by the pressure of other official duties to which the writer has been assigned during the last three years.

  3. Karst morphology and groundwater vulnerability of high alpine karst plateaus

    NASA Astrophysics Data System (ADS)

    Plan, Lukas; Decker, Kurt; Faber, Robert; Wagreich, Michael; Grasemann, Bernhard

    2009-07-01

    High alpine karst plateaus are recharge areas for major drinking water resources in the Alps and many other regions. Well-established methods for the vulnerability mapping of groundwater to contamination have not been applied to such areas yet. The paper characterises this karst type and shows that two common vulnerability assessment methods (COP and PI) classify most of the areas with high vulnerability classes. In the test site on the Hochschwab plateau (Northern Calcareous Alps, Austria), overlying layers are mostly absent, not protective or even enhance point recharge, where they have aquiclude character. The COP method classifies 82% of the area as highly or extremely vulnerable. The resulting maps are reasonable, but do not differentiate vulnerabilities to the extent that the results can be used for protective measures. An extension for the upper end of the vulnerability scale is presented that allows identifying ultra vulnerable areas. The proposed enhancement of the conventional approach points out that infiltration conditions are of key importance for vulnerability. The method accounts for karst genetical and hydrologic processes using qualitative and quantitative properties of karst depressions and sinking streams including parameters calculated from digital elevations models. The method is tested on the Hochschwab plateau where 1.7% of the area is delineated as ultra vulnerable. This differentiation could not be reached by the COP and PI methods. The resulting vulnerability map highlights spots of maximum vulnerability and the combination with a hazard map enables protective measures for a manageable area and number of sites.

  4. Clustering of cochlear oscillations in frequency plateaus as a tool to investigate SOAE generation

    NASA Astrophysics Data System (ADS)

    Epp, Bastian; Wit, Hero; van Dijk, Pim

    2015-12-01

    Spontonaeous otoacoustic emissions (SOAE) reflect the net effect of self-sustained activity in the cochlea, but do not directly provide information about the underlying mechanism and place of origin within the cochlea. The present study investigates if frequency plateaus as found in a linear array of coupled oscillators (OAM) [7] are also found in a transmission line model (TLM) which is able to generate realistic SOAEs [2] and if these frequency plateaus can be used to explain the formation of SOAEs. The simulations showed a clustering of oscillators along the simulated basilar membrane Both, the OAM and the TLM show traveling-wave like behavior along the oscillators coupled into one frequency plateau. While in the TLM roughness is required in order to produce SOAEs, no roughness is required to trigger frequency plateaus in the linear array of oscillators. The formation of frequency plateaus as a consequence of coupling between neighbored active oscillators might be the mechanism underlying SOAEs.

  5. Crustal volumes of the continents and of oceanic and continental submarine plateaus

    NASA Technical Reports Server (NTRS)

    Schubert, G.; Sandwell, D.

    1989-01-01

    Using global topographic data and the assumption of Airy isostasy, it is estimated that the crustal volume of the continents is 7182 X 10 to the 6th cu km. The crustal volumes of the oceanic and continental submarine plateaus are calculated at 369 X 10 to the 6th cu km and 242 X 10 to the 6th cu km, respectively. The total continental crustal volume is found to be 7581 X 10 to the 6th cu km, 3.2 percent of which is comprised of continental submarine plateaus on the seafloor. An upper bound on the contintental crust addition rate by the accretion of oceanic plateaus is set at 3.7 cu km/yr. Subduction of continental submarine plateaus with the oceanic lithosphere on a 100 Myr time scale yields an upper bound to the continental crustal subtraction rate of 2.4 cu km/yr.

  6. Numerical Simulations of Hydrodynamic Instability Growth in Polar-Direct-Drive Implosions at the National Ignition Facility

    NASA Astrophysics Data System (ADS)

    Shvydky, A.; Hohenberger, M.; Radha, P. B.; Rosenberg, M. J.; Craxton, R. S.; Goncharov, V. N.; Marozas, J. A.; Marshall, F. J.; McKenty, P. W.; Regan, S. P.; Sangster, T. C.

    2015-11-01

    Control of shell nonuniformities imprinted by the laser and amplified by hydrodynamic instabilities in the imploding target is critical to the success of polar-direct-drive ignition at the National Ignition Facility (NIF). To develop a platform for laser-imprint studies, hydrodynamic instability growth experiments in laser-driven implosions were performed on the NIF. The experiments used cone-in-shell targets with sinusoidal modulations of various wavelengths and amplitudes machined on the surface. Throughshell x-ray radiography was used to measure optical depth variations, from which the amplitudes of the shell areal-density modulations were extracted. Results of DRACO simulations of the growth of preimposed modulations and imprint-seeded perturbations will be presented and compared with the experimental data. This material is based upon work supported by the Department of Energy National Nuclear Security Administration under Award Number DE-NA0001944.

  7. TWO-DIMENSIONAL NUMERICAL STUDY FOR RAYLEIGH-TAYLOR AND RICHTMYER-MESHKOV INSTABILITIES IN RELATIVISTIC JETS

    SciTech Connect

    Matsumoto, Jin; Masada, Youhei

    2013-07-20

    We study the stability of a non-rotating single-component jet using two-dimensional special relativistic hydrodynamic simulations. By assuming translational invariance along the jet axis, we exclude the destabilization effect by Kelvin-Helmholtz mode. The nonlinear evolution of the transverse structure of the jet with a normal jet velocity is highlighted. An intriguing finding in our study is that Rayleigh-Taylor and Richtmyer-Meshkov type instabilities can destroy cylindrical jet configuration as a result of spontaneously induced radial oscillating motion. This is powered by in situ energy conversion between the thermal and bulk kinetic energies. The effective inertia ratio of the jet to the surrounding medium {eta} determines a threshold for the onset of instabilities. The condition {eta} < 1 should be satisfied for the transverse structure of the jet being persisted.

  8. Self-Consistant Numerical Modeling of E-Cloud Driven Instability of a Bunch Train in the CERN SPS

    SciTech Connect

    Vay, J-L.; Furman, M.A.; Secondo, R.; Venturini, M.; Fox, J.D.; Rivetta, C.H,

    2010-09-01

    The simulation package WARP-POSINST was recently upgraded for handling multiple bunches and modeling concurrently the electron cloud buildup and its effect on the beam, allowing for direct self-consistent simulation of bunch trains generating, and interacting with, electron clouds. We have used the WARP-POSINST package on massively parallel supercomputers to study the growth rate and frequency patterns in space-time of the electron cloud driven transverse instability for a proton bunch train in the CERN SPS accelerator. Results suggest that a positive feedback mechanism exists between the electron buildup and the e-cloud driven transverse instability, leading to a net increase in predicted electron density. Comparisons to selected experimental data are also given. Electron clouds have been shown to trigger fast growing instabilities on proton beams circulating in the SPS and other accelerators. So far, simulations of electron cloud buildup and their effects on beam dynamics have been performed separately. This is a consequence of the large computational cost of the combined calculation due to large space and time scale disparities between the two processes. We have presented the latest improvements of the simulation package WARP-POSINST for the simulation of self-consistent ecloud effects, including mesh refinement, and generation of electrons from gas ionization and impact at the pipe walls. We also presented simulations of two consecutive bunches interacting with electrons clouds in the SPS, which included generation of secondary electrons. The distribution of electrons in front of the first beam was initialized from a dump taken from a preceding buildup calculation using the POSINST code. In this paper, we present an extension of this work where one full batch of 72 bunches is simulated in the SPS, including the entire buildup calculation and the self-consistent interaction between the bunches and the electrons. Comparisons to experimental data are also given.

  9. Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit

    NASA Astrophysics Data System (ADS)

    Yu, Peicheng; Xu, Xinlu; Tableman, Adam; Decyk, Viktor K.; Tsung, Frank S.; Fiuza, Frederico; Davidson, Asher; Vieira, Jorge; Fonseca, Ricardo A.; Lu, Wei; Silva, Luis O.; Mori, Warren B.

    2015-12-01

    A hybrid Maxwell solver for fully relativistic and electromagnetic (EM) particle-in-cell (PIC) codes is described. In this solver, the EM fields are solved in k space by performing an FFT in one direction, while using finite difference operators in the other direction(s). This solver eliminates the numerical Cerenkov radiation for particles moving in the preferred direction. Moreover, the numerical Cerenkov instability (NCI) induced by the relativistically drifting plasma and beam can be eliminated using this hybrid solver by applying strategies that are similar to those recently developed for pure FFT solvers. A current correction is applied for the charge conserving current deposit to ensure that Gauss's Law is satisfied. A theoretical analysis of the dispersion properties in vacuum and in a drifting plasma for the hybrid solver is presented, and compared with PIC simulations with good agreement obtained. This hybrid solver is applied to both 2D and 3D Cartesian and quasi-3D (in which the fields and current are decomposed into azimuthal harmonics) geometries. Illustrative results for laser wakefield accelerator simulation in a Lorentz boosted frame using the hybrid solver in the 2D Cartesian geometry are presented, and compared against results from 2D UPIC-EMMA simulation which uses a pure spectral Maxwell solver, and from OSIRIS 2D lab frame simulation using the standard Yee solver. Very good agreement is obtained which demonstrates the feasibility of using the hybrid solver for high fidelity simulation of relativistically drifting plasma with no evidence of the numerical Cerenkov instability.

  10. Direct Numerical Modeling of E-Cloud Driven Instability of a Bunch Train in the CERN SPS

    SciTech Connect

    Vay, J-L.; Furman, M.A.; Venturini, M.

    2011-03-01

    The simulation package WARP-POSINST was recently upgraded for handling multiple bunches and modeling concurrently the electron cloud buildup and its effect on the beam, allowing for direct self-consistent simulation of bunch trains generating, and interacting with, electron clouds. We have used the WARP-POSINST package on massively parallel supercomputers to study the buildup and interaction of electron clouds with a proton bunch train in the CERN SPS accelerator. Results suggest that a positive feedback mechanism exists between the electron buildup and the e-cloud driven transverse instability, leading to a net increase in predicted electron density. Electron clouds have been shown to trigger fast growing instabilities on proton beams circulating in the SPS and other accelerators. So far, simulations of electron cloud buildup and their effects on beam dynamics have been performed separately. This is a consequence of the large computational cost of the combined calculation due to large space and time scale disparities between the two processes. We have presented the latest improvements of the simulation package WARP-POSINST for the simulation of self-consistent ecloud effects, including mesh refinement, and generation of electrons from gas ionization and impact at the pipe walls. We also presented simulations of two consecutive bunches interacting with electrons clouds in the SPS, which included generation of secondary electrons. The distribution of electrons in front of the first beam was initialized from a dump taken from a preceding buildup calculation using the POSINST code. In this paper, we present an extension of this work where one full batch of 72 bunches is simulated in the SPS, including the entire buildup calculation and the self-consistent interaction between the bunches and the electrons.

  11. Magnetization plateaus in the antiferromagnetic Ising chain with single-ion anisotropy and quenched disorder.

    PubMed

    Neto, Minos A; de Sousa, J Ricardo; Branco, N S

    2015-05-01

    We have studied the presence of plateaus on the low-temperature magnetization of an antiferromagnetic spin-1 chain, as an external uniform magnetic field is varied. A crystal-field interaction is present in the model and the exchange constants follow a random quenched (Bernoulli or Gaussian) distribution. Using a transfer-matrix technique we calculate the largest Lyapunov exponent and, from it, the magnetization at low temperatures as a function of the magnetic field, for different values of the crystal field and the width of the distributions. For the Bernoulli distribution, the number of plateaus increases, with respect to the uniform case [Litaiff et al., Solid State Commun. 147, 494 (2008)] and their presence can be linked to different ground states, when the magnetic field is varied. For the Gaussian distributions, the uniform scenario is maintained, for small widths, but the plateaus structure disappears as the width increases. PMID:26066165

  12. Franciscan complex calera limestones: Accreted remnants of farallon plate oceanic plateaus

    USGS Publications Warehouse

    Tarduno, J.A.; McWilliams, M.; Debiche, M.G.; Sliter, W.V.; Blake, M.C.

    1985-01-01

    The Calera Limestone, part of the Franciscan Complex of northern California, may have formed in a palaeoenvironment similar to Hess and Shatsky Rises of the present north-west Pacific1. We report here new palaeomagnetic results, palaeontological data and recent plate-motion models that reinforce this assertion. The Calera Limestone may have formed on Farallon Plate plateaus, north of the Pacific-Farallon spreading centre as a counterpart to Hess or Shatsky Rises. In one model2, the plateaus were formed by hotspots close to the Farallon_Pacific ridge axis. On accretion to North America, plateau dissection in the late Cretaceous to Eocene (50-70 Myr) could explain the occurrence of large volumes of pillow basalt and exotic blocks of limestone in the Franciscan Complex. Partial subduction of the plateaus could have contributed to Laramide (70-40 Myr) compressional events3. ?? 1985 Nature Publishing Group.

  13. Pressure dependence of the magnetization plateaus of SrCu2(BO3)2

    NASA Astrophysics Data System (ADS)

    Schneider, David A.; Coester, Kris; Mila, Frédéric; Schmidt, Kai Phillip

    2016-06-01

    We show that the critical fields of the magnetization plateaus of the Shastry-Sutherland model decrease significantly upon increasing the ratio of inter- to intradimer coupling and accordingly that the magnetization plateaus of SrCu2(BO3)2 shift to lower field under pressure, making the first two plateaus at 1 /8 and 2 /15 potentially accessible to neutron scattering experiments. These conclusions are based on the derivation of an effective classical model of interacting pinwheel-shaped spin-2 bound states using a combination of perturbative and graph-based continuous unitary transformations, showing that pinwheel crystals are indeed the lowest-energy plateau structures at low magnetization and that a simple model of intermediate-range two-body repulsion between pinwheels is able to account quantitatively for the plateau sequence.

  14. Theoretical and numerical treatment of modal instability in high-power core and cladding-pumped Raman fiber amplifiers.

    PubMed

    Naderi, Shadi; Dajani, Iyad; Grosek, Jacob; Madden, Timothy

    2016-07-25

    Raman fiber lasers have been proposed as potential candidates for scaling beyond the power limitations imposed on near diffraction-limited rare-earth doped fiber lasers. One limitation is the modal instability (MI) and we explore the physics of this phenomenon in Raman fiber amplifiers (RFAs). By utilizing the conservation of number of photons and conservation of energy in the absence of loss, the 3 × 3 governing system of nonlinear equations describing the pump and the signal modal content are decoupled and solved analytically for cladding-pumped RFAs. By comparing the extracted signal at MI threshold for the same step index-fiber, it is found that the MI threshold is independent of the length of the amplifier or whether the amplifier is co-pumped or counter-pumped; dictated by the integrated heat load along the length of fiber. We extend our treatment to gain-tailored RFAs and show that this approach is of limited utility in suppressing MI. Finally, we formulate the physics of MI in core-pumped RFAs where both pump and signal interferences participate in writing the time-dependent index of refraction grating. PMID:27464110

  15. Crustal structure in the eastern Colorado Plateaus Provence from seismic-refraction measurements

    USGS Publications Warehouse

    Roller, John C.

    1964-01-01

    A reversed seismic-refraction profile was recorded in the Colorado Plateaus Province from Hanksville, Utah, to Chinle, Arizona, The velocity of Pg is 6.2 km/sec, and the true velocity of Pn is 7.8 km/sec, Waves identified as reflections indicate that an intermediate layer in the crust has a velocity of approximately 6.8 km/sec. Thickness of the crust is 43 km at Chinle and 40 km at Hanksville. The Pn velocity in the Colorado Plateaus Province is the same as that in the Basin and Range Province, but is significantly lower than Pn in the High Plains of Colorado.

  16. Erosion and tectonics at the margins of continental plateaus

    NASA Technical Reports Server (NTRS)

    Masek, Jeffrey G.; Isacks, Bryan L.; Gubbels, Timothy L.; Fielding, Eric J.

    1994-01-01

    We hypothesize that the steep frontal slope and high peaks of the Beni region and Himalayan front largely reflect the high orographic precipitation and high erosion rates occurring in these regions and that the more gentle topography of the semiarid Pilcomayo region reflects a tectonic landform only slightly modified by erosion. We propose that orographic precipitation impinging on a plateau margin will generally tend to drop moisture low on the slope, eroding back the plateau while enhancing or maintaining the steep long-wavelength slope. A numerical model coupling orographic precipitation, erosion, and tectonic uplift demonstrates the plausibility of this hypothesis. The erosional efflux in both the Beni and Nepal Himalaya have been considerable, and simple mass balance calculations for the Himalaya suggest that during the Neogene, the erosional mass efflux has generally outpaced the tectonic mass influx. This contrasts with the apparent prior domination of tectonic influx and may reflect a decrease in the rate of tectonic addition during the same period, and/or increased late Cenozoic erosion rates.

  17. Crustal thinning between the Ethiopian and East African Plateaus from modeling Rayleigh wave dispersion

    SciTech Connect

    Benoit, M H; Nyblade, A A; Pasyanos, M E

    2006-01-17

    The East African and Ethiopian Plateaus have long been recognized to be part of a much larger topographic anomaly on the African Plate called the African Superswell. One of the few places within the African Superswell that exhibit elevations of less than 1 km is southeastern Sudan and northern Kenya, an area containing both Mesozoic and Cenozoic rift basins. Crustal structure and uppermost mantle velocities are investigated in this area by modeling Rayleigh wave dispersion. Modeling results indicate an average crustal thickness of 25 {+-} 5 km, some 10-15 km thinner than the crust beneath the adjacent East African and Ethiopian Plateaus. The low elevations can therefore be readily attributed to an isostatic response from crustal thinning. Low Sn velocities of 4.1-4.3 km/s also characterize this region.

  18. The tectonic setting of the Seychelles, Mascarene and Amirante Plateaus in the Western Equatorial Indian Ocean

    NASA Technical Reports Server (NTRS)

    Mart, Y.

    1988-01-01

    A system of marine plateaus occurs in the western equatorial Indian Ocean, forming an arcuate series of wide and shallow banks with small islands in places. The oceanic basins that surround the Seychelles - Amirante region are of various ages and reflect a complex seafloor spreading pattern. The structural analysis of the Seychelle - Amirante - Mascarene region reflects the tectonic evolution of the western equatorial Indian Ocean. It is suggested that due to the seafloor spreading during a tectonic stage, the Seychelles continental block drifted southwestwards to collide with the oceanic crust of the Mascarene Basin, forming an elongated folded structure at first, and then a subduction zone. The morphological similarity, the lithological variability and the different origin of the Seychelles Bank, the Mascarene Plateau and the Amirante Arc emphasizes the significant convergent effects of various plate tectonic processes on the development of marine plateaus.

  19. Digital surfaces and thicknesses of selected hydrogeologic units within the Ozark Plateaus aquifer system, northwestern Arkansas

    USGS Publications Warehouse

    Czarnecki, John B.; Bolyard, Susan E.; Hart, Rheannon M.; Clark, Jimmy M.

    2014-01-01

    Digital surfaces and thicknesses of nine hydrogeologic units of the Ozark Plateaus aquifer system from land surface to the top of the Gunter Sandstone in northwestern Arkansas were created using geophysical logs, drillers’ logs, geologist-interpreted formation tops, and previously published maps. The 6,040 square mile study area in the Ozark Plateaus Province includes Benton, Washington, Carroll, Madison, Boone, Newton, Marion, and Searcy Counties. The top of each hydrogeologic unit delineated on geophysical logs was based partly on previously published reports and maps and also from drillers’ logs. These logs were then used as a basis to contour digital surfaces showing the top and thickness of the Fayetteville Shale, the Boone Formation, the Chattanooga Shale, the Everton Formation, the Powell Dolomite, the Cotter Dolomite, the Roubidoux Formation, the Gasconade Dolomite, and the Gunter Sandstone.

  20. Identification and occurrence of uranium and vanadium minerals from the Colorado Plateaus

    USGS Publications Warehouse

    Weeks, A.D.; Thompson, M.E.

    1954-01-01

    This report, designed to make available to field geologists and others informa- tion obtained in recent investigations by the Geological Survey on identification and occurrence of uranium minerals of the Colorado Plateaus, contains descriptions of the physical properties, X-ray data, and in some instances results of chemical and spectrographic analysis of 48 uranium and vanadium minerals. Also included is a list of locations of mines from which the minerals have been identified.

  1. Retired flies, hidden plateaus, and the evolution of senescence in Drosophila melanogaster.

    PubMed

    Curtsinger, James W

    2016-06-01

    Late-life plateaus in age-specific mortality have been an evolutionary and biodemographic puzzle for decades. Although classic theory on the evolution of senescence predicts late-life walls of death, observations in experimental organisms document the opposite trend: a slowing in the rate of increase of mortality at advanced ages. Here, I analyze published life-history data on individual Drosophila melanogaster females and argue for a fundamental change in our understanding of mortality in this important model system. Mortality plateaus are not, as widely assumed, exclusive to late life, and are not explained by population heterogeneity-they are intimately connected to individual fecundity. Female flies begin adult life in the working stage, a period of active oviposition and low but accelerating mortality. Later they transition to the retired stage, a terminal period characterized by limited fecundity and relatively constant mortality. Because ages of transition differ between flies, age-synchronized cohorts contain a mix of working and retired flies. Early- and mid-life plateaus are obscured by the presence of working flies, but can be detected when cohorts are stratified by retirement status. Stage-specificity may be an important component of Drosophila life-history evolution. PMID:27166620

  2. Superplumes and ridge interactions between Ontong Java and Manihiki Plateaus and the Nova-Canton Trough

    NASA Astrophysics Data System (ADS)

    Larson, Roger L.

    1997-09-01

    The initial pulse of volcanism on the Ontong Java and Manihiki Plateaus began prior to 123 124 Ma and had largely ceased by about 122 Ma, whereas the intervening Pacific-Phoenix spreading ridge probably was disrupted between 120 and 115 Ma by formation of the Nova-Canton Trough rift system. Plateau formation may have resulted from a large plume head that rose beneath each of the two plateau sites. During the plateau eruption phase, these two plumes spread beyond the plateaus and within the asthenosphere. They eventually coalesced beneath the intervening spreading ridge, where the Nova-Canton Trough rift system is now located, at about the time volcanism ceased on the plateaus. The two plumes combined to form the Nova-Canton Trough rift system by reheating the nascent lithosphere near the ridge crest and widening the plate boundary from the usual 1 10 km to a much broader region. The present rift system is a maximum of 500 km wide, and much of this area may have formed by normal faulting and crustal accretion across the broadened plate boundary as plate divergence continued. Eventually, coherent spreading resumed somewhere south of the Nova-Canton Trough rift system.

  3. Detection of vertebral plateaus in lateral lumbar spinal X-ray images with Gabor filters.

    PubMed

    Alvarez Ribeiro, Eduardo; Nogueira-Barbosa, Marcello Henrique; Rangayyan, Rangaraj M; Azevedo-Marques, Paulo M

    2010-01-01

    A few recent studies have proposed computed-aided methods for the detection and analysis of vertebral bodies in radiographic images. This paper presents a method based on Gabor filters. Forty-one lateral lumbar spinal X-ray images from different patients were included in the study. For each image, a radiologist manually delineated the vertebral plateaus of L1, L2, L3, and L4 using a software tool for image display and mark-up. Each original image was filtered with a bank of 180 Gabor filters. The angle of the Gabor filter with the highest response at each pixel was used to derive a measure of the strength of orientation or alignment. In order to limit the spatial extent of the image data and the derived features in further analysis, a semi-automated procedure was applied to the original image. A neural network utilizing the logistic sigmoid function was trained with pixel intensity from the original image, the result of manual delineation of the plateaus, the Gabor magnitude response, and the alignment image. The average overlap between the results of detection by image processing and manual delineation of the plateaus of L1-L4 in the 41 images tested was 0.917. The results are expected to be useful in the analysis of vertebral deformities and fractures. PMID:21097095

  4. Two-fluid instability

    SciTech Connect

    Stewart, H.B.

    1984-01-01

    The ability of interpenetrating flow models to represent multidimensional instabilities is probed by numerical experiments with an L-shaped two-fluid jet. Periodic and nonperiodic oscillations of various types are observed, and a partial phase portrait is constructed. The numerical experiments suggest new approaches to verifying transient interpenetrating flow models. 18 references.

  5. Numerical study of friction-induced instability and acoustic radiation - Effect of ramp loading on the squeal propensity for a simplified brake model

    NASA Astrophysics Data System (ADS)

    Soobbarayen, K.; Sinou, J.-J.; Besset, S.

    2014-10-01

    This paper presents a numerical study of the influence of loading conditions on the vibrational and acoustic responses of a disc brake system subjected to squeal. A simplified model composed of a circular disc and a pad is proposed. Nonlinear effects of contact and friction over the frictional interface are modelled with a cubic law and a classical Coulomb's law with a constant friction coefficient. The stability analysis of this system shows the presence of two instabilities with one and two unstable modes that lead to friction-induced nonlinear vibrations and squeal noise. Nonlinear time analysis by temporal integration is conducted for two cases of loadings and initial conditions: a static load near the associated sliding equilibrium and a slow and a fast ramp loading. The analysis of the time responses shows that a sufficiently fast ramp loading can destabilize a stable configuration and generate nonlinear vibrations. Moreover, the fast ramp loading applied for the two unstable cases generates higher amplitudes of velocity than for the static load cases. The frequency analysis shows that the fast ramp loading generates a more complex spectrum than for the static load with the appearance of new resonance peaks. The acoustic responses for these cases are estimated by applying the multi-frequency acoustic calculation method based on the Fourier series decomposition of the velocity and the Boundary Element Method. Squeal noise emissions for the fast ramp loading present lower or higher levels than for the static load due to the different amplitudes of velocities. Moreover, the directivity is more complex for the fast ramp loading due to the appearance of new harmonic components in the velocity spectrum. Finally, the sound pressure convergence study shows that only the first harmonic components are sufficient to well describe the acoustic response.

  6. Correlating Cordilleran Ice Sheet Collapse with North Atlantic Heinrich Events using Global Radiocarbon Plateaus.

    NASA Astrophysics Data System (ADS)

    Hendy, I. L.; Cosma, T.

    2006-12-01

    The small, ephemeral Cordilleran Ice Sheet, present in Alaska, British Columbia and northern Washington during the last glacial cycle is believed to have behaved differently than the larger Laurentice Ice Sheet to climate forcing. High quality chronology is required to understand the relationship between this ice sheet and global climate change. Presently MD02-2496 (48°58.47N: 127°02.14W; 1190m water depth) is the highest resolution paleoclimate record available for the last 50 Ka in the Pacific northwest and contains intervals of glacial-marine sedimentation. High resolution dating based on 36 radiocarbon dates provide a chronology that includes radiocarbon age plateaus, while planktonic foraminiferal stable isotopes offer a continuous record of climate change. Glacial-marine sediments collected from this site on the continental slope west of Vancouver Island, British Columbia, document three intervals of iceberg discharge during the last ~50 Ka. Gradually increasing quantities of ice rafted detritus (grains >250μm, g-1; IRD) followed by abrupt cessation within ~500 years is strongly suggestive of catastrophic iceberg discharge. The penultimate event is correlated to marine invasion of the Juan de Fuca Strait, and Puget Sound, while the final IRD event with that of Georgia Strait. We posit that these previously unknown IRD events represent repeated rapid iceberg discharge related to Cordilleran Ice Sheet collapse. The events occur near the end of radiocarbon plateaus at 13.35 ±90 and 14.05 ±70 14C Kyr BP (not reservoir corrected). If these plateaus correlate with the 12.2 and 13.3 14C Kyr BP plateaus recorded in Carriaco Basin and elsewhere, local reservoir ages can be calculated and vary between 1,150 and 1,550 years similar to those derived locally from glacial wood-shell pairs. Furthermore, if the plateaus result from reduced North Atlantic Deep Water export and consequently Heinrich Events, the Cordilleran IRD events are related to North Atlantic iceberg

  7. Effects of Coarse Legacy Sediment on Rivers of the Ozark Plateaus and Implications for Native Mussel Fauna

    NASA Astrophysics Data System (ADS)

    Erwin, S. O.; Jacobson, R. B.; Eric, A. B.; Jones, J. C.; Anderson, B. W.

    2015-12-01

    Perturbations to sediment regimes due to anthropogenic activities may have long lasting effects, especially in systems dominated by coarse sediment where travel times are relatively long. Effectively evaluating management alternatives requires understanding the future trajectory of river response at both the river network and reach scales. The Ozark Plateaus physiographic province is a montane region in the interior US composed primarily of Paleozoic sedimentary rock. Historic land-use practices around the turn of the last century accelerated delivery of coarse sediment to river channels. Effects of this legacy sediment persist in two national parks, Ozark National Scenic Riverways, MO and Buffalo National River, AR, and are of special concern for management of native mussel fauna. These species require stable habitat, yet they occupy inherently dynamic environments: alluvial rivers. At the river-network scale, analysis of historical data reveals the signature of sediment waves moving through river networks in the Ozarks. Channel planform alternates between relatively stable, straight reaches, and wider, multithread reaches which have been more dynamic over the past several decades. These alternate planform configurations route and store sediment differently, and translate into different patterns of bed stability at the reach scale, which in turn affects the distribution and availability of habitat for native biota. Geomorphic mapping and hydrodynamic modeling reveal the complex relations between planform (in)stability, flow dynamics, bed mobility, and aquatic habitat in systems responding to increased sediment supply. Reaches that have a more dynamic planform may provide more hydraulic refugia and habitat heterogeneity compared to stable, homogeneous reaches. This research provides new insights that may inform management of sediment and mussel habitat in rivers subject to coarse legacy sediment.

  8. Shoulder Instability

    MedlinePlus

    ... Risk Factors Is shoulder instability the same as shoulder dislocation? No. The signs of dislocation and instability might ... the same to you--weakness and pain. However, dislocation occurs when your shoulder goes completely out of place. The shoulder ligaments ...

  9. Water quality in the Ozark Plateaus, Arkansas, Kansas, Missouri, and Oklahoma, 1992-95

    USGS Publications Warehouse

    Petersen, James C.; Adamski, James C.; Bell, Richard W.; Davis, Jerri V.; Femmer, Suzanne R.; Freiwald, David A.; Joseph, Robert L.

    1998-01-01

    This report is intended to summarize major findings that emerged between 1992 and 1995 from the water-quality assessment of the Ozark Plateaus Study Unit and to relate these findings to water-quality issues of regional and national concern. The information is primarily intended for those who are involved in water-resource management. Yet, the information contained here may also interest those who simply wish to know more about the quality of water in the rivers and aquifers in the area where they live.

  10. Key wintertime meteorological features of the Grand Canyon and the Colorado Plateaus Basin

    SciTech Connect

    Whiteman, C.D.; Allwine, K.J.

    1992-06-01

    In the winter of 1989--1990 a major meteorological and air pollution experiment was conducted in the Colorado Plateaus Basin (Richards et al., 1991). The focus of the experiment, conducted by Arizona`s Soft River Project, was to investigate the influence of three 750-MW coal-fired power plant units at the Navajo Generating Station near Page, Arizona, on visibility at Grand Canyon National Park. As part of the meteorological experiment, surface and upper air data were collected from multiple sites within the basin. This data set is the most comprehensive meteorological data set ever collected within the region, and the purpose of this paper is to briefly summarize the key wintertime meteorological features of the Colorado Plateaus Basin and the Grand Canyon, through which the basin drains, using analyses of the Winter Visibility Study data. Our analyses focused primarily on thermally driven circulations within the basin and the Grand Canyon, but we also investigated the surface energy budget that drives these circulations and the interactions between the thermal circulations and the overlying synoptic-scale flows.

  11. Key wintertime meteorological features of the Grand Canyon and the Colorado Plateaus Basin

    SciTech Connect

    Whiteman, C.D.; Allwine, K.J.

    1992-06-01

    In the winter of 1989--1990 a major meteorological and air pollution experiment was conducted in the Colorado Plateaus Basin (Richards et al., 1991). The focus of the experiment, conducted by Arizona's Soft River Project, was to investigate the influence of three 750-MW coal-fired power plant units at the Navajo Generating Station near Page, Arizona, on visibility at Grand Canyon National Park. As part of the meteorological experiment, surface and upper air data were collected from multiple sites within the basin. This data set is the most comprehensive meteorological data set ever collected within the region, and the purpose of this paper is to briefly summarize the key wintertime meteorological features of the Colorado Plateaus Basin and the Grand Canyon, through which the basin drains, using analyses of the Winter Visibility Study data. Our analyses focused primarily on thermally driven circulations within the basin and the Grand Canyon, but we also investigated the surface energy budget that drives these circulations and the interactions between the thermal circulations and the overlying synoptic-scale flows.

  12. Ground water in the Springfield-Salem plateaus of southern Missouri and northern Arkansas

    USGS Publications Warehouse

    Harvey, Edward Joseph

    1980-01-01

    Average ground-water conditions have not changed significantly in the Springfield-Salem plateaus section of southern Missouri and northern Arkansas in the past 25 years except in the vicinity of well fields. The amount of ground water pumped is approximately 200 cubic feet per second, which is about 5 percent of the total discharge at the 80 percent point on flow-duration curves for major streams. Ground-water recharge is variable and occurs through sinkholes by infiltration in upland areas of good permeability, and through streambeds that lose flow. Main waterbearing zones lie in the Potosi Dolomite and the lower dolomite and sandstone of the Gasconade Dolomite. Cavernous connections from ground surface to depths as great as 1,500 feet occur in the West Plains area, Mo., and result in deep circulation of water. Municipal well-water in the area often becomes turbid after rainstorms, despite well depths of 1 ,500 feet and 950 to 1,000 feet of pressure-grouted casing. Ground-water movement is generaly north and south from the crest of the Springfield-Salem plateaus, which extend across southern Missouri from the St. Francois Mountains to the southwest. Interbasin diversion of surface- and ground-water flow is common. (USGS)

  13. Colorado Plateaus Ecoregion: Chapter 21 in Status and trends of land change in the Western United States--1973 to 2000

    USGS Publications Warehouse

    Stier, Michael P.

    2012-01-01

    The Colorado Plateaus Ecoregion covers approximately 129,617 km2 (50,045 mi2) within southern and eastern Utah, western Colorado, and the extreme northern part of Arizona (fig. 1). The terrain of this ecoregion is characterized by broad plateaus, ancient volcanoes, and deeply dissected canyons (Booth and others, 1999; fig. 2). The ecoregion is bounded on the east by the Wyoming Basin and Southern Rockies Ecoregions in Colorado and on the northwest by the Wasatch and Uinta Mountains Ecoregion in northern and central Utah. To the south, the ecoregion borders the Arizona/New Mexico Plateau Ecoregion, which has a higher elevation and more grasslands than the Colorado Plateaus Ecoregion (Omernik, 1987; U.S. Environmental Protection Agency, 1997).

  14. Hip instability.

    PubMed

    Smith, Matthew V; Sekiya, Jon K

    2010-06-01

    Hip instability is becoming a more commonly recognized source of pain and disability in patients. Traumatic causes of hip instability are often clear. Appropriate treatment includes immediate reduction, early surgery for acetabular rim fractures greater than 25% or incarcerated fragments in the joint, and close follow-up to monitor for avascular necrosis. Late surgical intervention may be necessary for residual symptomatic hip instability. Atraumatic causes of hip instability include repetitive external rotation with axial loading, generalized ligamentous laxity, and collagen disorders like Ehlers-Danlos. Symptoms caused by atraumatic hip instability often have an insidious onset. Patients may have a wide array of hip symptoms while demonstrating only subtle findings suggestive of capsular laxity. Traction views of the affected hip can be helpful in diagnosing hip instability. Open and arthroscopic techniques can be used to treat capsular laxity. We describe an arthroscopic anterior hip capsular plication using a suture technique. PMID:20473129

  15. Collective instabilities

    SciTech Connect

    K.Y. Ng

    2003-08-25

    The lecture covers mainly Sections 2.VIII and 3.VII of the book ''Accelerator Physics'' by S.Y. Lee, plus mode-coupling instabilities and chromaticity-driven head-tail instability. Besides giving more detailed derivation of many equations, simple interpretations of many collective instabilities are included with the intention that the phenomena can be understood more easily without going into too much mathematics. The notations of Lee's book as well as the e{sup jwt} convention are followed.

  16. Crustal Structure of the Ethiopian Rift and Adjacent Plateaus: Results of new integrated interpretation

    NASA Astrophysics Data System (ADS)

    Tadesse, K.; Keller, G. R.

    2006-12-01

    The Ethiopian rift is the large part of the East African Rift system, which represents an incipient divergent plate boundary. This important structure provides excellent opportunities to study the transition from continental to oceanic. As a result, geophysical data are becoming increasingly available but some results are contradictory. We used a newly enhanced gravity database and seismic information to produce an integrated interpretation of the crustal structure beneath the Ethiopian rift and the adjacent plateaus. Wide regions have been covered to assess the regional structures including the Kenyan and Ethiopian rifts and the area covered by the Ethiopian flood basalt. Broad negative Bouguer gravity anomalies are delineated over the Ethiopian Plateaus and the Kenyan dome. Residual gravity anomalies, which parallel the major border faults clearly highlight the segregation between the plateaus and the rift valleys. Results of other filtering techniques have clearly revealed individual volcanic centers within the rift valleys. Positive gravity anomalies outside the rift valleys may be associated with older structures, shield volcanoes, or structures that are related to the initiation and propagation of rifting. A long axial profile from the central part of Kenya to the Afar triple junction has been modeled to investigate along-axis crustal variation of the East African rift system, with emphasis on the Ethiopian rift. This modeling has been constrained using seismic refraction data from the Ethiopian Afar Geoscientific Lithospheric Experiment (EAGLE) and Kenya Rift International Seismic Project (KRISP) results. We are able to see a thin crust (~26 km) in the Afar triangle with a gradual thickening (~40 km) southwards towards the Main Central Ethiopian rift (MER). The crust thickness decreases towards Turkana rift (~22 km), and increases again towards the central eastern rift section in Kenya. Our profile model across the MER has revealed that the eastern rift

  17. NUMERICAL SIMULATIONS OF THE MAGNETIC RAYLEIGH-TAYLOR INSTABILITY IN THE KIPPENHAHN-SCHLUeTER PROMINENCE MODEL. I. FORMATION OF UPFLOWS

    SciTech Connect

    Hillier, Andrew; Isobe, Hiroaki; Shibata, Kazunari; Berger, Thomas

    2012-02-20

    The launch of the Hinode satellite led to the discovery of rising plumes, dark in chromospheric lines, that propagate from large ({approx}10 Mm) bubbles that form at the base of quiescent prominences. The plumes move through a height of approximately 10 Mm while developing highly turbulent profiles. The magnetic Rayleigh-Taylor instability was hypothesized to be the mechanism that drives these flows. In this study, using three-dimensional (3D) MHD simulations, we investigate the nonlinear stability of the Kippenhahn-Schlueter prominence model for the interchange mode of the magnetic Rayleigh-Taylor instability. The model simulates the rise of a buoyant tube inside the quiescent prominence model, where the interchange of magnetic field lines becomes possible at the boundary between the buoyant tube and the prominence. Hillier et al. presented the initial results of this study, where upflows of constant velocity (maximum found 6 km s{sup -1}) and a maximum plume width Almost-Equal-To 1.5 Mm which propagate through a height of approximately 6 Mm were found. Nonlinear interaction between plumes was found to be important for determining the plume dynamics. In this paper, using the results of ideal MHD simulations, we determine how the initial parameters for the model and buoyant tube affect the evolution of instability. We find that the 3D mode of the magnetic Rayleigh-Taylor instability grows, creating upflows aligned with the magnetic field of constant velocity (maximum found 7.3 km s{sup -1}). The width of the upflows is dependent on the initial conditions, with a range of 0.5-4 Mm which propagate through heights of 3-6 Mm. These results are in general agreement with the observations of the rising plumes.

  18. Landau-level dispersion and the quantum Hall plateaus in bilayer graphene

    SciTech Connect

    Zarenia, M.; Peeters, F. M.; Vasilopoulos, P.; Pourtolami, N.

    2013-12-04

    We study the quantum Hall effect (QHE) in bilayer graphene using the Kubo-Greenwood formula. At zero temperature the Hall conductivity σ{sub yx} is given by σ{sub yx} = 4(N+1)e{sup 2}/h with N the index of the highest occupied Landau level (LL). Including the dispersion of the LLs and their width, due to e.g. scattering by impurities, produces the plateau of the n = 0 LL in agreement with experimental results on doped samples and similar theoretical results on single-layer graphene plateaus widen with impurity concentration. Further, the evaluated resistivity ρ{sub xx} exhibits a strong, oscillatory dependence on the electron concentration. Explicit results are obtained for δ-function impurities.

  19. Do Children With Fragile X Syndrome Show Declines or Plateaus in Adaptive Behavior?

    PubMed

    Hahn, Laura J; Brady, Nancy C; Warren, Steven F; Fleming, Kandace K

    2015-09-01

    This study explores if children with fragile X syndrome (FXS) show advances, declines, or plateaus in adaptive behavior over time and the relationship of nonverbal cognitive abilities and autistic behavior on these trajectories. Parents of 55 children with FXS completed the Vineland Adaptive Behavior Scales ( Sparrow, Balla, & Cicchetti, 1984 ; Sparrow, Cicchetti, & Balla, 2005 ) between 3 and 6 times from 2 to 10 years of age. Using raw scores, results indicate that about half of the sample showed advances in adaptive behavior, whereas the other half showed declines, indicating a regression in skills. Children who were more cognitively advanced and had less autistic behaviors had higher trajectories. Understanding the developmental course of adaptive behavior in FXS has implications for educational planning and intervention, especially for those children showing declines. PMID:26322389

  20. Giant quantum Hall plateaus generated by charge transfer in epitaxial graphene.

    PubMed

    Alexander-Webber, J A; Huang, J; Maude, D K; Janssen, T J B M; Tzalenchuk, A; Antonov, V; Yager, T; Lara-Avila, S; Kubatkin, S; Yakimova, R; Nicholas, R J

    2016-01-01

    Epitaxial graphene has proven itself to be the best candidate for quantum electrical resistance standards due to its wide quantum Hall plateaus with exceptionally high breakdown currents. However one key underlying mechanism, a magnetic field dependent charge transfer process, is yet to be fully understood. Here we report measurements of the quantum Hall effect in epitaxial graphene showing the widest quantum Hall plateau observed to date extending over 50 T, attributed to an almost linear increase in carrier density with magnetic field. This behaviour is strong evidence for field dependent charge transfer from charge reservoirs with exceptionally high densities of states in close proximity to the graphene. Using a realistic framework of broadened Landau levels we model the densities of donor states and predict the field dependence of charge transfer in excellent agreement with experimental results, thus providing a guide towards engineering epitaxial graphene for applications such as quantum metrology. PMID:27456765

  1. 2D seismic data processing for straight lines in the loess plateaus in Fuxian of Shanbei

    NASA Astrophysics Data System (ADS)

    Li, Minjie; Chen, Yequan; Zhang, Hai; Pang, Shangming; Deng, Guozhen

    2005-01-01

    The crooked seismic lines along valleys were irregular previously in Fuxian of Shanbei, showing an irregular branch in plane, and hard to complete close grids. Therefore, it’s difficult to conduct reservoir inversion of 2D seismic data. In 2001, Zhongyuan Oilfield Company carried out the study on field acquisition methods and seismic processing technology in Fuxian. Straight lines were passing through plateaus and formed seismic grids by using flexible geometry with variable linear bins. Data processing involved model-inversion based refraction static correction, surface consistent amplitude compensation, deconvolution, and pre-stack noise attenuation. As the result, seismic data with a high fidelity was provided for the subsequent reservoir predictions, small-amplitude structure interpretation and integrative geologic study. Because all lines were jointed to form grids, comprehensive interpretation of reservoir inversion could be finally implemented by using the pseudo logging method to control lines without wells.

  2. Giant quantum Hall plateaus generated by charge transfer in epitaxial graphene

    NASA Astrophysics Data System (ADS)

    Alexander-Webber, J. A.; Huang, J.; Maude, D. K.; Janssen, T. J. B. M.; Tzalenchuk, A.; Antonov, V.; Yager, T.; Lara-Avila, S.; Kubatkin, S.; Yakimova, R.; Nicholas, R. J.

    2016-07-01

    Epitaxial graphene has proven itself to be the best candidate for quantum electrical resistance standards due to its wide quantum Hall plateaus with exceptionally high breakdown currents. However one key underlying mechanism, a magnetic field dependent charge transfer process, is yet to be fully understood. Here we report measurements of the quantum Hall effect in epitaxial graphene showing the widest quantum Hall plateau observed to date extending over 50 T, attributed to an almost linear increase in carrier density with magnetic field. This behaviour is strong evidence for field dependent charge transfer from charge reservoirs with exceptionally high densities of states in close proximity to the graphene. Using a realistic framework of broadened Landau levels we model the densities of donor states and predict the field dependence of charge transfer in excellent agreement with experimental results, thus providing a guide towards engineering epitaxial graphene for applications such as quantum metrology.

  3. Solitons versus parametric instabilities during ionospheric heating

    NASA Technical Reports Server (NTRS)

    Nicholson, D. R.; Payne, G. L.; Downie, R. M.; Sheerin, J. P.

    1984-01-01

    Various effects associated with ionospheric heating are investigated by numerically solving the modified Zakharov (1972) equations. It is shown that, for typical ionospheric parameters, the modulational instability is more important than the parametric decay instability in the spatial region of strongest heater electric field. It is concluded that the modulational instability leads to the formation of solitons, as originally predicted by Petviashvili (1976).

  4. A Fundamental Plane for Long Gamma-Ray Bursts with X-Ray Plateaus

    NASA Astrophysics Data System (ADS)

    Dainotti, M. G.; Postnikov, S.; Hernandez, X.; Ostrowski, M.

    2016-07-01

    A class of long gamma-ray bursts (GRBs) presenting light curves with an extended plateau phase in their X-ray afterglows obeys a correlation between the rest-frame end-time of the plateau, T a , and its corresponding X-ray luminosity, L a , (Dainotti et al). In this work we perform an analysis of a total sample of 176 Swift GRBs with known redshifts, exhibiting afterglow plateaus. By adding a third parameter that is the peak luminosity in the prompt emission, L peak, we discover the existence of a new three-parameter correlation. The scatter of data about this plane becomes smaller when a class-specific GRB sample is defined. This sample of 122 GRBs is selected from the total sample by excluding GRBs with associated supernovae (SNe), X-ray flashes and short GRBs with extended emission. With this sample the three-parameter correlation identifies a GRB “fundamental plane.” Moreover, we further limit our analysis to GRBs with light curves with good data coverage and almost flat plateaus, 40 GRBs forming our “gold sample.” The intrinsic scatter, {σ }{int}=0.27+/- 0.04, for the three-parameter correlation for this last sub-class is more than two times smaller than the value for the {L}a-{T}a one, making this the tightest three-parameter correlation that involves the afterglow plateau phase. Finally, we also show that a slightly less tight correlation is present between L peak and a proxy for the total energy emitted during the plateau phase, {L}a{T}a, confirming the existence of an energy scaling between the prompt and afterglow phases.

  5. Dislocation motion and instability

    NASA Astrophysics Data System (ADS)

    Zhu, Yichao; Chapman, Stephen Jonathan; Acharya, Amit

    2013-08-01

    The Peach-Koehler expression for the stress generated by a single (non-planar) curvilinear dislocation is evaluated to calculate the dislocation self stress. This is combined with a law of motion to give the self-induced motion of a general dislocation curve. A stability analysis of a rectilinear, uniformly translating dislocation is then performed. The dislocation is found to be susceptible to a helical instability, with the maximum growth rate occurring when the dislocation is almost, but not exactly, pure screw. The non-linear evolution of the instability is determined numerically, and implications for slip band formation and non-Schmid behavior in yielding are discussed.

  6. Jurassic plume-origin ophiolites in Japan: accreted fragments of oceanic plateaus

    NASA Astrophysics Data System (ADS)

    Ichiyama, Yuji; Ishiwatari, Akira; Kimura, Jun-Ichi; Senda, Ryoko; Miyamoto, Tsuyoshi

    2014-07-01

    later and in lower temperatures than the E- and D-type rocks. The Mikabu and Sorachi Plateaus were formed in a low-latitude region of the Late Jurassic Pacific Ocean possibly near a subduction zone, partially experienced high P/ T metamorphism during subduction, and then uplifted in association with (or without, in case of Mikabu) the supra-subduction zone ophiolite. The Mikabu and Sorachi Plateaus may be the Late Jurassic oceanic LIPs that could have been formed in brotherhood with the Shatsky Rise.

  7. 2014/2015 Investigations of the Ontong Java and Kerguelen Plateaus

    NASA Astrophysics Data System (ADS)

    Coffin, M. F.; Whittaker, J. M.

    2013-12-01

    The two largest oceanic plateaus, Ontong Java in the western Pacific, and Kerguelen in the southern Indian Ocean, will be the focus of scheduled multidisciplinary/interdisciplinary shipboard expeditions in 2014 and 2015. In mid-2014, scientists aboard the Schmidt Ocean Institute's RV Falkor will investigate the origin and evolution of two large atolls, Ontong Java and Nukumanu, surmounting the ca 122 Ma Ontong Java Plateau, as well how Kroenke Canyon, which deeply incises the plateau, formed and evolved. First-ever multibeam bathymetry and sub-bottom profiling data from the atolls and canyon will reveal their submarine and shallow sub-seafloor morphology, and, if combined with geochemical and geochronological analyses of potential igneous basement samples, will yield important information on their origin and evolution. The primary goals of this atoll and canyon project are: to test potential genetic relationships between a) the atolls and the OJP, and b) the atolls and Kroenke Canyon; to understand and model how atolls and canyons form and evolve on oceanic plateaus, isolated from terrestrial influences and subject to sea level fluctuations; and to contribute to understanding tsunami risk on low-lying atolls. In late 2014 and early 2015, researchers aboard Australia's new Marine National Facility, RV Investigator, will investigate active submarine hotspot volcanism on the Kerguelen Plateau and its consequences. The project's overall aim is to test the hypothesis that hydrothermal activity driven by active submarine magmatism fertilizes surface waters with iron that enhances primary biological productivity. Surmounting the Cretaceous plateau, Heard and McDonald Islands are among the world's most active hotspot volcanoes, and new multibeam bathymetry and sub-bottom profiling data will enable identification of candidate active submarine volcanoes, which we will sample. In the overlying water column, we will collect samples to test for the presence or absence of

  8. Magnetization plateaus of an easy-axis kagome antiferromagnet with extended interactions

    NASA Astrophysics Data System (ADS)

    Plat, X.; Alet, F.; Capponi, S.; Totsuka, K.

    2015-11-01

    We investigate the properties in finite magnetic field of an extended anisotropic X X Z spin-1/2 model on the kagome lattice, originally introduced by Balents, Fisher, and Girvin [Phys. Rev. B 65, 224412 (2002), 10.1103/PhysRevB.65.224412]. The magnetization curve displays plateaus at magnetization m =1 /6 and 1 /3 when the anisotropy is large. Using low-energy effective constrained models (quantum loop and quantum dimer models), we discuss the nature of the plateau phases, found to be crystals that break discrete rotation and/or translation symmetries. Large-scale quantum Monte Carlo simulations were carried out in particular for the m =1 /6 plateau. We first map out the phase diagram of the effective quantum loop model with an additional loop-loop interaction to find stripe order around the point relevant for the original model as well as a topological Z2 spin liquid. The existence of a stripe crystalline phase is further evidenced by measuring both standard structure factor and entanglement entropy of the original microscopic model.

  9. Regional Curves for Bankfull Channel Characteristics in the Appalachian Plateaus, West Virginia

    USGS Publications Warehouse

    Messinger, Terence

    2009-01-01

    Streams in the Appalachian Plateaus Physiographic Province in West Virginia were classified as a single region on the basis of bankfull characteristics. Regression lines for annual peak flow and drainage area measured at streamgages in the study area at recurrence intervals between 1.2 and 1.7 years fell within the 99-percent confidence interval of the regression line for bankfull flow. Channel characteristics were intermediate among those from surrounding states and regions where comparable studies have been done. The stream reaches that were surveyed were selected for apparent stability, and to represent gradients of drainage area, elevation, and mean annual precipitation. Profiles of high-water marks left by bankfull and near-bankfull peaks were surveyed, either as part of slope-area flow measurements at ungaged reaches, or to transfer known flow information to cross sections for gaged reaches. The slope-area measurements made it possible to include ungaged sites in the study, but still relate bankfull dimensions to peak flow and frequency.

  10. The Occurrence of Anomalous Conductance Plateaus and Spin Textures in Quantum Point Contacts

    NASA Astrophysics Data System (ADS)

    Wan, J.; Cahay, M.; Debray, P.; Newrock, R.

    2010-03-01

    Recently, we used a NEGF formalism [1] to provide a theoretical explanation for the experimentally observed 0.5G0 (G0=2e^2/h) plateau in the conductance of side-gated quantum point contacts (QPCs) in the presence of lateral spin-orbit coupling (LSOC) [2]. We showed that the 0.5G0 plateau appears in the QPCs without any external magnetic field as a result of three ingredients: an asymmetric lateral confinement, a LSOC, and a strong electron-electron (e-e) interaction. In this report, we present the results of simulations for a wide range of QPC dimensions and biasing parameters showing that the same physics predicts the appearance of other anomalous plateaus at non-integer values of G0, including the well-known 0.7G0 anomaly. These features are related to a plethora of spin textures in the QPC that depend sensitively on material, device, biasing parameters, temperature, and the strength of the e-e interaction. [1] J. Wan, M. Cahay, P. Debray, and R.S. Newrock, Phys. Rev. B 80, 155440 (2009). [2] P. Debray, S.M. Rahman, J. Wan, R.S. Newrock, M. Cahay, A.T. Ngo, S.E. Ulloa, S.T. Herbert, M. Muhammad, and M. Johnson, Nature Nanotech. 4, 759 (2009).

  11. A fundamental plane for Gamma Ray Bursts with X-Ray plateaus

    NASA Astrophysics Data System (ADS)

    Dainotti, Maria; Petrosian, Vahe'; Postnikov, Sergey; Hernandez, Xavier; Ostrowski, Michal

    2016-06-01

    A class of long Gamma-ray bursts (GRBs) presenting light curves with an extended plateau phase in their X-ray afterglows obey a two parameter correlation (Dainotti et al. 2008) betweenthe rest frame end time of the plateau, T_a, and its corresponding X-ray luminosity, L_{a}. This correlation is a useful tool towards cosmological parameter inferences and yields important constraints on physical GRB models.We here show through an analysis of 122 Swift long GRBs with known redshifts and presentingafterglow plateaus, that including also the peak luminosity in the prompt emission, L_{peak}, as a third parameter, a fundamental plane appears. A careful statistical analysis excludes the possibility of this correlation being due to detection thresholds or observational selection biases. The correlation is tightest when a class-specific GRB sample is defined excluding GRB-SNe, X-ray flashes and short GRBs with extended emission, and required also to have a good data coverage and an almost flat plateau. The intrinsic scatter, σ_{int}=0.27 ± 0.04, for the three-parameter correlation for this specific subclass (40 GRBs) is 53.5% smaller than the (L_{a}, T_a) one,making this the tighest three parameter correlation involving the plateau phase. This result will hence enhance cosmological parameter determinations and GRB model constraints. Finally, we also show that a slightly noisier correlation is also present between L_{peak} and a proxy of the total energy emitted during the plateau phase, the product L_{a}T_{a}, hinting at the presence of an energy scaling between prompt and afterglow phases.

  12. Carbon Abundance Plateaus among Carbon-Enhanced Metal-Poor Stars

    NASA Astrophysics Data System (ADS)

    Yoon, Jinmi; He, Siyu; Placco, Vinicius; Carollo, Daniela; Beers, Timothy C.

    2016-01-01

    A substantial fraction of low-metallicity stars in the Milky Way, the Carbon-Enhanced Metal-Poor (CEMP) stars, exhibit enhancements of their carbon-to-iron relative to the solar value ([C/Fe] > +0.7). They can be divided into several sub-classes, depending on the nature and degree of the observed enhancements of their neutron-capture elements, providing information on their likely progenitors. CEMP-s stars (which exhibit enhanced s-process elements) are thought to be enhanced by mass transfer from an evolved AGB companion, while CEMP-no stars (which exhibit no over-abundances of neutron-capture elements) appear to be associated with explosions of the very first generations of stars. High-resolution spectroscopic analyses are generally required in order to make these sub-classifications.Several recent studies have suggested the existence of bimodality in the distribution of absolute carbon abundances among CEMP stars -- most CEMP-no stars belong to a low-C band ((A(C) ˜ 6.5), while most CEMP-s stars reside on a high-C band (A(C) ˜ 8.25). The number of CEMP stars considered by individual studies is, however, quite small, so we have compiled all available high-resolution spectroscopic data for CEMP stars, in order to further investigate the existence of the claimed carbon bi-modality, and to consider what can be learned about the progenitors of CEMP-s and CEMP-no stars based on the observed distribution of A(C) on the individual plateaus.We acknowledge partial support from the grant PHY 14-30152; Physics Frontier Center/JINA Center for the Evolution of the Elements (JINA-CEE), awarded by the US National Science Foundation.

  13. Hybrid rocket instability

    NASA Technical Reports Server (NTRS)

    Greiner, B.; Frederick, R. A., Jr.

    1993-01-01

    The paper provides a brief review of theoretical and experimental studies concerned with hybrid rocket instability. The instabilities discussed include atomization and mixing instabilities, chuffing instabilities, pressure coupled combustion instabilities, and vortex shedding. It is emphasized that the future use of hybrid motor systems as viable design alternatives will depend on a better understanding of hybrid instability.

  14. Hybrid rocket instability

    NASA Astrophysics Data System (ADS)

    Greiner, B.; Frederick, R. A., Jr.

    1993-06-01

    The paper provides a brief review of theoretical and experimental studies concerned with hybrid rocket instability. The instabilities discussed include atomization and mixing instabilities, chuffing instabilities, pressure coupled combustion instabilities, and vortex shedding. It is emphasized that the future use of hybrid motor systems as viable design alternatives will depend on a better understanding of hybrid instability.

  15. Age and geochemistry of volcanic rocks from the Hikurangi and Manihiki oceanic Plateaus

    NASA Astrophysics Data System (ADS)

    Hoernle, Kaj; Hauff, Folkmar; van den Bogaard, Paul; Werner, Reinhard; Mortimer, Nick; Geldmacher, Jörg; Garbe-Schönberg, Dieter; Davy, Bryan

    2010-12-01

    Here we present the first radiometric age data and a comprehensive geochemical data set (including major and trace element and Sr-Nd-Pb-Hf isotope ratios) for samples from the Hikurangi Plateau basement and seamounts on and adjacent to the plateau obtained during the R/V Sonne 168 cruise, in addition to age and geochemical data from DSDP Site 317 on the Manihiki Plateau. The 40Ar/ 39Ar age and geochemical data show that the Hikurangi basement lavas (118-96 Ma) have surprisingly similar major and trace element and isotopic characteristics to the Ontong Java Plateau lavas (ca. 120 and 90 Ma), primarily the Kwaimbaita-type composition, whereas the Manihiki DSDP Site 317 lavas (117 Ma) have similar compositions to the Singgalo lavas on the Ontong Java Plateau. Alkalic, incompatible-element-enriched seamount lavas (99-87 Ma and 67 Ma) on the Hikurangi Plateau and adjacent to it (Kiore Seamount), however, were derived from a distinct high time-integrated U/Pb (HIMU)-type mantle source. The seamount lavas are similar in composition to similar-aged alkalic volcanism on New Zealand, indicating a second wide-spread event from a distinct source beginning ca. 20 Ma after the plateau-forming event. Tholeiitic lavas from two Osbourn seamounts on the abyssal plain adjacent to the northeast Hikurangi Plateau margin have extremely depleted incompatible element compositions, but incompatible element characteristics similar to the Hikurangi and Ontong Java Plateau lavas and enriched isotopic compositions intermediate between normal mid-ocean-ridge basalt (N-MORB) and the plateau basement. These younger (˜52 Ma) seamounts may have formed through remelting of mafic cumulate rocks associated with the plateau formation. The similarity in age and geochemistry of the Hikurangi, Ontong Java and Manihiki Plateaus suggest derivation from a common mantle source. We propose that the Greater Ontong Java Event, during which ˜1% of the Earth's surface was covered with volcanism, resulted from a

  16. Pleistocene periglacial cryogenic mounds (lithalsas) on basalt plateaus in the western Pannonian Basin

    NASA Astrophysics Data System (ADS)

    Sebe, Krisztina; Csillag, Gábor

    2015-04-01

    On some basalt plateaus of the western Pannonian Basin, Hungary, fields of circular depressions occur. They are traditionally called "basalt karst' and their formation has been attributed either to collapse over karstifying rocks or to anthropogenic action (quarrying); however, both of these theories are questionable. The depressions are situated between elevations of 350-500 m a.s.l. and are characteristically surrounded by circular raised rims or ramparts. They measure a few m-s (up to ~10 m) in diameter, the ramparts emerge 0.5-1.5 m above the surrounding level ground and encircle a depression of 1-2 (-3) m deep in the middle. Depressions cluster in well delineated, high-density groups, with individual fields containing several dozens of these forms. Neighbouring ramparts are tightly packed, often interfere and depressions can thus coalesce creating composite forms. The ramparts are composed of coarse (dm-sized) basalt blocks, whose material is identical to that of the surrounding terrain and seems to originate from the depression. Many of the depressions host ephemeral ponds. Raised rims exclude formation of these landforms by any karstic processes. The anthropogenic theory is opposed by the lack of the remains of any facilities (e.g. roads), of tools and by the very illogical distribution and geometry of depressions from the point of human use. On the contrary, we interpret these ramparted depressions as being of periglacial origin, remnants of cryogenic mounds. The central depression and the emergent rampart can be well explained by the ice core raising the overlying rock and by the radial downsliding of this material on the ice core to the margins. Within cryogenic mounds, clustering and size of the forms fits the characteristics of perennial frost mound without peat cover, i.e. lithalsas or minerogenic palsas. Cryogenic mounds are important paleoclimatic indicators. Based on modern analogs, these lithalsa scars indicate the former presence of discontinuous

  17. Surface Uplift History of the Central Andes: Implications for the Growth of Orogenic Plateaus

    NASA Astrophysics Data System (ADS)

    Garzione, C. N.; Hoke, G. D.; Libarkin, J. C.; MacFadden, B. J.; Withers, S.

    2007-05-01

    have occurred in the Andes in late Miocene time, including removal of high density lower lithosphere and redistribution of crust by erosion/sedimentation and flow of low density middle-lower crust, are likely mechanisms for building broad, flat, high elevation plateaus in convergent tectonic settings.

  18. Progressive Landslides in Uplifted Volcanic Plateaus: Persistent Loci of Channel Perturbation

    NASA Astrophysics Data System (ADS)

    Safran, E. B.; Anderson, S. W.; Mills-Novoa, M.; Othus, S.; Ely, L.; House, P. K.; O'Connor, J. E.; Grant, G.; Fenton, C.; Beebee, R. A.

    2007-12-01

    The semi-arid uplifted volcanic plateaus of the southern interior Columbia River basin contain over 300 large landslides or landslide complexes, ranging in area from several tenths of a km2 to several tens of km2. The distribution of these landslides is dominated by the outcropping of key stratigraphic contacts between coherent, volcanic cap rock atop weak sedimentary or volcaniclastic units in areas of >100 m local relief. The morphologies of many of these landslide complexes suggest a progressive mode of mass movement, with rubble-capped failure slices arrayed downslope at intervals of 10s to 100s of meters and deep tension cracks separating incipient failure blocks from the intact headscarp. Field evidence from the Owyhee River in southeastern Oregon indicates that individual landslide complexes can persist for millions of years. In one reach, for example, remnants of a 1.9 million year old intracanyon lava flow are inset against ancient landslide blocks. In the same location, cosmogenic isotope dating of boulders on a likely dam-burst flood deposit reveal a channel-blocking mass movement that may be as young as Holocene in age. The persistence of these landslide complexes has important implications for channel evolution, as it suggests that, in some environments, sediment supply may be chronically elevated at point sources. On the basis of GIS-based mapping of regional landslides and on field study of individual landslide complexes, we hypothesize that: 1) this sediment supply becomes increasingly dominated by fine-grained material as channels progressively incise into the weak units underlying coherent lava caps; and 2) the mass movements that impinge on the channels become correspondingly more earthflow-like. Loci of persistent landsliding are also subject to episodic variations in channel width due to physical constrictions caused by impinging failure masses. The discrete localization of large landslide complexes by particular stratigraphic and topographic

  19. Instability mechanisms in swirling flows

    NASA Astrophysics Data System (ADS)

    Gallaire, F.; Chomaz, J.-M.

    2003-09-01

    We investigate the stability of the screened Rankine vortex with added plug flow where the azimuthal velocity decreases abruptly outside the core of the vortex. The jump in circulation is known to induce centrifugal and azimuthal Kelvin-Helmholtz instabilities. Their effect on the stability of the different azimuthal wave number m is discussed using physical considerations associated with asymptotic expansions and numerical computations of the dispersion relation. It is shown that the axial shear and centrifugal instability are active for all m, and that modes with |m|⩾2 are also destabilized by azimuthal shear. In contrast, the bending modes m=±1 are stabilized by a coupling with Kelvin waves in the core. Effects of rotation on the absolute/convective transition are also discussed. The absolute instability of positive helical modes is seen to be promoted by centrifugal instability and azimuthal shear.

  20. Extensional instability in electro-osmotic microflows of polymer solutions

    NASA Astrophysics Data System (ADS)

    Bryce, R. M.; Freeman, M. R.

    2010-03-01

    Fluid transport in microfluidic systems typically is laminar due to the low Reynolds number characteristic of the flow. The inclusion of suspended polymers imparts elasticity to fluids, allowing instabilities to be excited when substantial polymer stretching occurs. For high molecular weight polymer chains we find that flow velocities achievable by standard electro-osmotic pumping are sufficient to excite extensional instabilities in dilute polymer solutions. We observe a dependence in measured fluctuations on polymer concentration which plateaus at a threshold corresponding to the onset of significant molecular crowding in macromolecular solutions; plateauing occurs well below the overlap concentration. Our results show that electro-osmotic flows of complex fluids are disturbed from the steady regime, suggesting potential for enhanced mixing and requiring care in modeling the flow of complex liquids such as biopolymer suspensions.

  1. Superresonant instability of a compressible hydrodynamic vortex

    NASA Astrophysics Data System (ADS)

    Oliveira, Leandro A.; Cardoso, Vitor; Crispino, Luís C. B.

    2016-06-01

    We show that a purely circulating and compressible system, in an adiabatic regime of acoustic propagation, presents superresonant instabilities. To show the existence these instabilities, we compute the quasinormal mode frequencies of this system numerically using two different frequency domain methods.

  2. X-ray flares, plateaus and chromatic breaks of GRB afterglows from up-scattered forward-shock emission

    NASA Astrophysics Data System (ADS)

    Panaitescu, A.

    2008-01-01

    Scattering of the forward-shock synchrotron emission by a relativistic outflow located behind the leading blast wave may produce an X-ray emission brighter than that coming directly from the forward shock and may explain four features displayed by Swift X-ray afterglows: flares, plateaus (slow decays), chromatic light-curve breaks and fast post-plateau decays. For a cold scattering outflow, the reflected flux overshines the primary one if the scattering outflow is nearly baryon-free and highly relativistic. These two requirements can be relaxed if the scattering outflow is energized by weak internal shocks, so that the incident forward-shock photons are also inverse-Compton scattered, in addition to bulk scattering. Sweeping-up of the photons left behind by the forward shock naturally yields short X-ray flares. Owing to the boost in photon energy produced by bulk scattering, the reflected emission is more likely to overshine that coming directly from the forward shock at higher photon energies, yielding light-curve plateaus and breaks that appear only in the X-ray. The brightness, shape and decay of the X-ray light-curve plateau depend on the radial distribution of the scatterer's Lorentz factor and mass flux. Chromatic X-ray light-curve breaks and sharp post-plateau decays cannot be accommodated by the direct forward-shock emission and argue in favour of the scattering-outflow model proposed here. On the other hand, the X-ray afterglows without plateaus, those with achromatic breaks and those with very long lived power-law decays are more naturally accommodated by the standard forward-shock model. Thus, the diversity of X-ray light curves arises from the interplay of the scattered and direct forward-shock emissions.

  3. New insights on lithospheric foundering from thermo-mechanically coupled numerical modelling

    NASA Astrophysics Data System (ADS)

    Pastor-Galán, Daniel; Thieulot, Cedric

    2015-04-01

    Earth's lithosphere is recycled into the mantle as required by global mass considerations. At least during the latest 1 G.y. the main mechanism of lithospheric foundering into the mantle has been subduction. Yet other mechanisms of mantle removal such as Rayleigh-Taylor-type instability or delamination have significant influence at present as revealed by mantle anomalies, and are considered to be likely candidates for the main recycling mechanisms of lithospheric during the Archean. Although lithospheric mantle removal has been geophysically imaged, e.g. Carpathians, Colorado Plateau, at many other locations geophysical and geological observations also seem to indicate that mantle lithosphere is anomalously thin or absent. Potential places where lithospheric mantle foundering processes took place are The Urals, the Variscides, underneath the Ibero Armorican Orocline in western Europe, and the Tibetan, Puna and Anatolian Plateaus. Lithospheric foundering has been blamed for, among others, cratonization processes, rapid surface uplift, generation of voluminous magmatism, changes in crustal stress from compression to extension and a long etc. However, its triggering mechanisms are not well studied, and a variety of possible explanations have been given for lithospheric foundering processes, including convective instability following orogenic thickening or some other perturbation of thermal boundary layers, development of eclogitic roots, erosion of the lithosphere by a flat-subducting slab or partial melting of the asthenosphere, and partial intruding pyroxenites into the base of lithosphere. To understand the mechanisms, causes and consequences of lithospheric foundering, we explored lithospheric foundering in an assortment of scenarios using the numerical code, ELEFANT, an user-friendly multipurpose geodynamics code. Preliminary results indicate that changes in geometry, thermal state and composition of the lithosphere, associated with mantle flow, can have a first

  4. Magnetohydrodynamic instability

    NASA Technical Reports Server (NTRS)

    Priest, E. R.; Cargill, P.; Forbes, T. G.; Hood, A. W.; Steinolfson, R. S.

    1986-01-01

    There have been major advances in the theory of magnetic reconnection and of magnetic instability, with important implications for the observations, as follows: (1) Fast and slow magnetic shock waves are produced by the magnetohydrodynamics of reconnection and are potential particle accelerators. (2) The impulsive bursty regime of reconnection gives a rapid release of magnetic energy in a series of bursts. (3) The radiative tearing mode creates cool filamentary structures in the reconnection process. (4) The stability analyses imply that an arcade can become unstable when either its height or twist of plasma pressure become too great.

  5. Quantification of subaerial and episodic subglacial erosion rates on high latitude upland plateaus: Cumberland Peninsula, Baffin Island, Arctic Canada

    NASA Astrophysics Data System (ADS)

    Margreth, Annina; Gosse, John C.; Dyke, Arthur S.

    2016-02-01

    Long-term rates of subaerial and episodic subglacial erosion by predominately cold-based ice cover are determined for tors on weathered plateaus on Cumberland Peninsula. By measuring terrestrial cosmogenic nuclide concentrations in differentially weathered surfaces on a given tor, we reconstruct the complex exposure and erosion history involving recurring cold-based ice cover. The style and rate of subaerial and subglacial erosion at multiple tor sites on Cumberland Peninsula are assessed with a Monte Carlo approach that computes plausible exposure histories based on a proxy record of global ice volume. Constant subaerial erosion rates by weathering are low (<2 mm ka-1), corroborated by nuclide concentrations measured on two tors located on coastal ridge crests that have likely never been glaciated. Summit plateaus intermittently covered by cold-based ice throughout the Quaternary have experienced episodic subglacial erosion by plucking of fractured bedrock with rates between 1 and 16 mm ka-1. Variation of episodic erosion rates is associated with topographic position of the sampled tors and bedrock fracture density. Most of the tors were last glacially plucked in pre-ultimate glaciations, not during the Wisconsinan glaciation. Furthermore, the new approach provides evidence for the extent of ice coverage during the late Wisconsinan, which is significant if no erratics are available for exposure dating. Despite late Pleistocene intervals of ice cover and glacial plucking, tor-studded landscapes of Cumberland Peninsula are of considerable antiquity.

  6. Subwavelength modulational instability and plasmon oscillons in nanoparticle arrays.

    PubMed

    Noskov, Roman E; Belov, Pavel A; Kivshar, Yuri S

    2012-03-01

    We study modulational instability in nonlinear arrays of subwavelength metallic nanoparticles and analyze numerically nonlinear scenarios of the instability development. We demonstrate that modulational instability can lead to the formation of regular periodic or quasiperiodic modulations of the polarization. We reveal that such nonlinear nanoparticle arrays can support long-lived standing and moving oscillating nonlinear localized modes--plasmon oscillons. PMID:22463637

  7. Late Cretaceous - Paleogene forearc sedimentation and accretion of oceanic plateaus and seamounts along the Middle American convergent margin (Costa Rica)

    NASA Astrophysics Data System (ADS)

    Baumgartner, Peter O.; Baumgartner-Mora, Claudia; Andjic, Goran

    2016-04-01

    The Late Cretaceous-Paleogene sedimentation pattern in space and time along the Middle American convergent margin was controlled by the accretion of Pacific plateaus and seamounts. The accretion of more voluminous plateaus must have caused the temporary extinction of the arc and tectonic uplift, resulting in short lived episodes of both pelagic and neritic biogenic sedimentation. By the Late Eocene, shallow carbonate environments became widespread on a supposed mature arc edifice, that is so far only documented in arc-derived sediments. In northern Costa Rica forearc sedimentation started during the Coniacian-Santonian on the Aptian-Turonian basement of the Manzanillo Terrane. The arrival and collision of the Nicoya Terrane (a CLIP-like, 139-83 Ma Pacific plateau) and the Santa Elena Terrane caused the extinction of the arc during late Campanian- Early Maastrichtian times, indicated by the change to pelagic limestone sedimentation (Piedras Blancas Formation) in deeper areas and shallow-water rudistid - Larger Benthic Foraminfera limestone on tectonically uplifted areas of all terranes. Arc-derived turbidite sedimentation resumed in the Late Maastrichtian and was again interrupted during the Late Paleocene - Early Eocene, perhaps due to the underplating of a yet unknown large seamount. The extinction of the arc resulted in the deposition of the siliceous pelagic Buenavista Formation, as well as the principally Thanetian Barra Honda carbonate platform on a deeply eroded structural high in the Tempisque area. In southern Costa Rica the basement is thought to be the western edge of the CLIP. It is Santonian-Campanian in age and is only exposed in the southwestern corner of Herradura. Cretaceous arc-forearc sequences are unknown, except for the Maastrichtian-Paleocene Golfito Terrane in southeastern Costa Rica. The distribution and age of shallow/pelagic carbonates vs. arc-derived detrital sediments is controlled by the history of accretion of Galápagos hot spot

  8. Instability of EDS maglev systems

    SciTech Connect

    Cai, Y.; Chen, S.S.

    1993-09-01

    Instabilities of an EDS maglev suspension system with 3 D.O.F. and 5 D.O.F. vehicles traveling on a double L-shaped set of guideway conductors have been investigated with various experimentally measured magnetical force data incorporated into the theoretical models. Divergence and flutter are obtained from both analytical and numerical solutions for coupled vibration of the 3 D.O.F. maglev vehicle model. Instabilities of five direction motions (heave, slip, rill, pitch and yaw) are observed for the 4 D.O.F. vehicle model. It demonstrates that system parameters, such as, system damping, vehicle geometry and coupling effects among five different motions play very important roles in the occurrence of dynamic instabilities of maglev vehicles.

  9. Tidal instability in exoplanetary systems evolution

    NASA Astrophysics Data System (ADS)

    Cébron, D.; Moutou, C.; Le Bars, M.; Le Gal, P.; Farès, R.

    2011-02-01

    A new element is proposed to play a role in the evolution of extrasolar planetary systems: the tidal (or elliptical) instability. It comes from a parametric resonance and takes place in any rotating fluid whose streamlines are (even slightly) elliptically deformed. Based on theoretical, experimental and numerical works, we estimate the growth rate of the instability for hot-jupiter systems, when the rotation period of the star is known. We present the physical process, its application to stars, and preliminary results obtained on a few dozen systems, summarized in the form of a stability diagram. Most of the systems are trapped in the so-called "forbidden zone", where the instability cannot grow. In some systems, the tidal instability is able to grow, at short timescales compared to the system evolution. Implications are discussed in the framework of misaligned transiting systems, as the rotational axis of the star would be unstable in systems where this elliptical instability grows.

  10. Combustion Instabilities Modeled

    NASA Technical Reports Server (NTRS)

    Paxson, Daniel E.

    1999-01-01

    NASA Lewis Research Center's Advanced Controls and Dynamics Technology Branch is investigating active control strategies to mitigate or eliminate the combustion instabilities prevalent in lean-burning, low-emission combustors. These instabilities result from coupling between the heat-release mechanisms of the burning process and the acoustic flow field of the combustor. Control design and implementation require a simulation capability that is both fast and accurate. It must capture the essential physics of the system, yet be as simple as possible. A quasi-one-dimensional, computational fluid dynamics (CFD) based simulation has been developed which may meet these requirements. The Euler equations of mass, momentum, and energy have been used, along with a single reactive species transport equation to simulate coupled thermoacoustic oscillations. A very simple numerical integration scheme was chosen to reduce computing time. Robust boundary condition procedures were incorporated to simulate various flow conditions (e.g., valves, open ends, and choked inflow) as well as to accommodate flow reversals that may arise during large flow-field oscillations. The accompanying figure shows a sample simulation result. A combustor with an open inlet, a choked outlet, and a large constriction approximately two thirds of the way down the length is shown. The middle plot shows normalized, time-averaged distributions of the relevant flow quantities, and the bottom plot illustrates the acoustic mode shape of the resulting thermoacoustic oscillation. For this simulation, the limit cycle peak-to-peak pressure fluctuations were 13 percent of the mean. The simulation used 100 numerical cells. The total normalized simulation time was 50 units (approximately 15 oscillations), which took 26 sec on a Sun Ultra2.

  11. Electrostatic ion cyclotron velocity shear instability

    NASA Technical Reports Server (NTRS)

    Lemons, D. S.; Winske, D.; Gary, S. P.

    1992-01-01

    A local electrostatic dispersion equation is derived for a shear flow perpendicular to an ambient magnetic field, which includes all kinetic effects and involves only one important parameter. The dispersion equation is cast in the form of Gordeyev integrals and is solved numerically. Numerical solutions indicate that an ion cyclotron instability is excited. The instability occurs roughly at multiples of the ion cyclotron frequency (modified by the shear), with the growth rate or the individual harmonics overlapping in the wavenumber. At large values of the shear parameter, the instability is confined to long wavelengths, but at smaller shear, a second distinct branch at shorter wavelengths also appears. The properties of the instability obtained are compared with those obtained in the nonlocal limit by Ganguli et al. (1985, 1988).

  12. Water-quality assessment of the Ozark Plateaus study unit, Arkansas, Kansas, Missouri, and Oklahoma- fish communities in streams of the Ozark Plateaus and their relations to selected environmental factors

    USGS Publications Warehouse

    Petersen, James C.

    1998-01-01

    Fish communities from 22 reaches at 18 stations in the Ozark Plateaus were sampled in 1993, 1994, and 1995. The 18 stations were chosen to represent selected combinations of major environmental factors (geology/physiographic area, land use, and basin size). Additional physical, chemical, and biological factors also were measured for each of the 22 reaches and the influence of these factors upon the fish communities was investigated. Fish community samples collected at the 22 reaches identified differences in these communities that can be attributed to differences in land use and related water-quality and habitat characteristics. Communities from agriculture reaches tended to have more species, increased relative abundance of stonerollers and members of the sucker family, and decreased relative abundance of members of the sunfish and darter families. Several groups of environmental factors (concentrations of nutrients, organic carbon, suspended sediment, and dissolved oxygen; measures related to ionic strength; measures related to riparian vegetation; measures related to substrate; and measures related to stream size) appear to be related to land-use differences and fish community differences. Three multivariate analysis techniques (two ordination techniques and a classification technique) yielded similar results when applied to the fish community data. Fish communities from reaches with more similar land use in their basins and with similar drainage areas generally were grouped closer together in the analysis. Water quality, substrate, stream morphology, and riparian measures appear to be affecting fish communities at these reaches. The relations between land use, stream size, and fish communities have implications for waterquality assessments of Ozark streams. Compared to other parts of the United States, many fish species live in the Ozark Plateaus. At least 19 of these species are endemic to the Ozarks area. Many of these species are intolerant of habitat or

  13. Computational methods for probability of instability calculations

    NASA Technical Reports Server (NTRS)

    Wu, Y.-T.; Burnside, O. H.

    1990-01-01

    This paper summarizes the development of the methods and a computer program to compute the probability of instability of a dynamic system than can be represented by a system of second-order ordinary linear differential equations. Two instability criteria based upon the roots of the characteristics equation or Routh-Hurwitz test functions are investigated. Computational methods based on system reliability analysis methods and importance sampling concepts are proposed to perform efficient probabilistic analysis. Numerical examples are provided to demonstrate the methods.

  14. On the convective-absolute nature of river bedform instabilities

    NASA Astrophysics Data System (ADS)

    Vesipa, Riccardo; Camporeale, Carlo; Ridolfi, Luca; Chomaz, Jean Marc

    2014-12-01

    River dunes and antidunes are induced by the morphological instability of stream-sediment boundary. Such bedforms raise a number of subtle theoretical questions and are crucial for many engineering and environmental problems. Despite their importance, the absolute/convective nature of the instability has never been addressed. The present work fills this gap as we demonstrate, by the cusp map method, that dune instability is convective for all values of the physical control parameters, while the antidune instability exhibits both behaviors. These theoretical predictions explain some previous experimental and numerical observations and are important to correctly plan flume experiments, numerical simulations, paleo-hydraulic reconstructions, and river works.

  15. Combustion instability modeling and analysis

    SciTech Connect

    Santoro, R.J.; Yang, V.; Santavicca, D.A.; Sheppard, E.J.

    1995-12-31

    It is well known that the two key elements for achieving low emissions and high performance in a gas turbine combustor are to simultaneously establish (1) a lean combustion zone for maintaining low NO{sub x} emissions and (2) rapid mixing for good ignition and flame stability. However, these requirements, when coupled with the short combustor lengths used to limit the residence time for NO formation typical of advanced gas turbine combustors, can lead to problems regarding unburned hydrocarbons (UHC) and carbon monoxide (CO) emissions, as well as the occurrence of combustion instabilities. The concurrent development of suitable analytical and numerical models that are validated with experimental studies is important for achieving this objective. A major benefit of the present research will be to provide for the first time an experimentally verified model of emissions and performance of gas turbine combustors. The present study represents a coordinated effort between industry, government and academia to investigate gas turbine combustion dynamics. Specific study areas include development of advanced diagnostics, definition of controlling phenomena, advancement of analytical and numerical modeling capabilities, and assessment of the current status of our ability to apply these tools to practical gas turbine combustors. The present work involves four tasks which address, respectively, (1) the development of a fiber-optic probe for fuel-air ratio measurements, (2) the study of combustion instability using laser-based diagnostics in a high pressure, high temperature flow reactor, (3) the development of analytical and numerical modeling capabilities for describing combustion instability which will be validated against experimental data, and (4) the preparation of a literature survey and establishment of a data base on practical experience with combustion instability.

  16. Turbine instabilities: Case histories

    NASA Technical Reports Server (NTRS)

    Laws, C. W.

    1985-01-01

    Several possible causes of turbine rotor instability are discussed and the related design features of a wide range of turbomachinery types and sizes are considered. The instrumentation options available for detecting rotor instability and assessing its severity are also discussed.

  17. TRANSVERSE INSTABILITIES IN RHIC.

    SciTech Connect

    Blaskiewicz, M; Cameron, P; Catalan-Lasheras, N; Dawson, C; Degen, C; Drees, K; Fischer, W; Koropsak, E; Michnoff, R; Montag, C; Roser, T

    2003-05-12

    The beam quality in RHIC can be significantly impacted by a transverse instability which can occur just after transition [1]. Data characterizing the instability are presented and analyzed. Techniques for ameliorating the situation are considered.

  18. Dust deposition during the Early Holocene on the loess plateaus of the Vojvodina region in Northern Serbia

    NASA Astrophysics Data System (ADS)

    Markovic, Slobodan; Timar-Gabor, Alida; Stevens, Thomas; Guo, Zhengtang; Hao, Qingzhen; Song, Yang; Hambach, Ulrich; Lehmkuhl, Frank; Peric, Zoran; Obreht, Igor; Zeeden, Christian; Veres, Daniel; Gavrilov, Milivoj

    2015-04-01

    The Northern Serbian province of Vojvodina is a lowland area encompassing the confluence of the Danube, Sava, Tisa (Tisza), Drava, Morava and Tamiš (Temes, Timiş) rivers, which separate several remnant loess plateaus. Loess sediments in the Vojvodina region are among the oldest and most complete loess-paleosol formations in Europe. These thick sequences contain a detailed paleoclimatic record since the Early Pleistocene. The better preservation of Serbian loess-paleosol sequences compared to other European loess records is most likely related to the persistence of much drier conditions in the region, coupled with "plateau-like" dust accumulation style. Recently and through detailed luminescence-based chronological investigations of accumulation derived from several loess sections we aimed at addressing the timing of the onset of Holocene soil (S0) formation in the wider region. So far, the chronological results demonstrate a lack of intensive pedogenesis coeval with the postulated Holocene onset (ie., 11.7 ka BP), and continuation of Aeolian dust deposition during the Early Holocene in some of the investigated sections. Lake sediment and speleothem records from the wider area also suggest that, at least regionally, the hydroclimatic characteristics of the Early Holocene differed markedly. This evidence leads to an important question about the validity of previously generalized direct stratigraphic correlations between regional terrestrial environmental archives and global marine and ice core records (direct synchronization of records vs. acknowledging leads/lags), that employ the Late Pleistocene/Holocene boundary at 11.7 as an absolute tie point.

  19. A Compilation of Provisional Karst Geospatial Data for the Interior Low Plateaus Physiographic Region, Central United States

    USGS Publications Warehouse

    Taylor, Charles J.; Nelson, Hugh L.

    2008-01-01

    Geospatial data needed to visualize and evaluate the hydrogeologic framework and distribution of karst features in the Interior Low Plateaus physiographic region of the central United States were compiled during 2004-2007 as part of the Ground-Water Resources Program Karst Hydrology Initiative (KHI) project. Because of the potential usefulness to environmental and water-resources regulators, private consultants, academic researchers, and others, the geospatial data files created during the KHI project are being made available to the public as a provisional regional karst dataset. To enhance accessibility and visualization, the geospatial data files have been compiled as ESRI ArcReader data folders and user interactive Published Map Files (.pmf files), all of which are catalogued by the boundaries of surface watersheds using U.S. Geological Survey (USGS) eight-digit hydrologic unit codes (HUC-8s). Specific karst features included in the dataset include mapped sinkhole locations, sinking (or disappearing) streams, internally drained catchments, karst springs inventoried in the USGS National Water Information System (NWIS) database, relic stream valleys, and karst flow paths obtained from results of previously reported water-tracer tests.

  20. Rogue Waves and Modulational Instability

    NASA Astrophysics Data System (ADS)

    Zakharov, V. E.; Dyachenko, A.

    2015-12-01

    The most plausible cause of rogue wave formation in a deep ocean is development of modulational instability of quasimonochromatic wave trains. An adequate model for study of this phenomenon is the Euler equation for potential flow of incompressible fluid with free surface in 2-D geometry. Numerical integration of these equations confirms completely the conjecture of rogue wave formation from modulational instability but the procedure is time consuming for determination of rogue wave appearance probability for a given shape of wave energy spectrum. This program can be realized in framework of simpler model using replacement of the exact interaction Hamiltonian by more compact Hamiltonian. There is a family of such models. The popular one is the Nonlinear Schrodinger equation (NLSE). This model is completely integrable and suitable for numerical simulation but we consider that it is oversimplified. It misses such important phenomenon as wave breaking. Recently, we elaborated much more reliable model that describes wave breaking but is as suitable as NLSE from the point of numerical modeling. This model allows to perform massive numerical experiments and study statistics of rogue wave formation in details.

  1. The plasmoid instability during asymmetric inflow magnetic reconnection

    NASA Astrophysics Data System (ADS)

    Murphy, Nicholas A.; Young, Aleida K.; Shen, Chengcai; Lin, Jun; Ni, Lei

    2013-06-01

    Theoretical studies of the plasmoid instability generally assume that the reconnecting magnetic fields are symmetric. We relax this assumption by performing two-dimensional resistive magnetohydrodynamic simulations of the plasmoid instability during asymmetric inflow magnetic reconnection. Magnetic asymmetry modifies the onset, scaling, and dynamics of this instability. Magnetic islands develop preferentially into the weak magnetic field upstream region. Outflow jets from individual X-points impact plasmoids obliquely rather than directly as in the symmetric case. Consequently, deposition of momentum by the outflow jets into the plasmoids is less efficient, the plasmoids develop net vorticity, and shear flow slows down secondary merging between islands. Secondary merging events have asymmetry along both the inflow and outflow directions. Downstream plasma is more turbulent in cases with magnetic asymmetry because islands are able to roll around each other after exiting the current sheet. As in the symmetric case, plasmoid formation facilitates faster reconnection for at least small and moderate magnetic asymmetries. However, when the upstream magnetic field strengths differ by a factor of 4, the reconnection rate plateaus at a lower value than expected from scaling the symmetric results. We perform a parameter study to investigate the onset of the plasmoid instability as a function of magnetic asymmetry and domain size. There exist domain sizes for which symmetric simulations are stable but asymmetric simulations are unstable, suggesting that moderate magnetic asymmetry is somewhat destabilizing. We discuss the implications for plasmoid and flux rope formation in solar eruptions, laboratory reconnection experiments, and space plasmas. The differences between symmetric and asymmetric simulations provide some hints regarding the nature of the three-dimensional plasmoid instability.

  2. The plasmoid instability during asymmetric inflow magnetic reconnection

    SciTech Connect

    Murphy, Nicholas A.; Young, Aleida K.; Shen, Chengcai; Lin, Jun; Ni, Lei

    2013-06-15

    Theoretical studies of the plasmoid instability generally assume that the reconnecting magnetic fields are symmetric. We relax this assumption by performing two-dimensional resistive magnetohydrodynamic simulations of the plasmoid instability during asymmetric inflow magnetic reconnection. Magnetic asymmetry modifies the onset, scaling, and dynamics of this instability. Magnetic islands develop preferentially into the weak magnetic field upstream region. Outflow jets from individual X-points impact plasmoids obliquely rather than directly as in the symmetric case. Consequently, deposition of momentum by the outflow jets into the plasmoids is less efficient, the plasmoids develop net vorticity, and shear flow slows down secondary merging between islands. Secondary merging events have asymmetry along both the inflow and outflow directions. Downstream plasma is more turbulent in cases with magnetic asymmetry because islands are able to roll around each other after exiting the current sheet. As in the symmetric case, plasmoid formation facilitates faster reconnection for at least small and moderate magnetic asymmetries. However, when the upstream magnetic field strengths differ by a factor of 4, the reconnection rate plateaus at a lower value than expected from scaling the symmetric results. We perform a parameter study to investigate the onset of the plasmoid instability as a function of magnetic asymmetry and domain size. There exist domain sizes for which symmetric simulations are stable but asymmetric simulations are unstable, suggesting that moderate magnetic asymmetry is somewhat destabilizing. We discuss the implications for plasmoid and flux rope formation in solar eruptions, laboratory reconnection experiments, and space plasmas. The differences between symmetric and asymmetric simulations provide some hints regarding the nature of the three-dimensional plasmoid instability.

  3. Granular Rayleigh-Taylor instability

    SciTech Connect

    Vinningland, Jan Ludvig; Johnsen, Oistein; Flekkoey, Eirik G.; Maaloey, Knut Joergen; Toussaint, Renaud

    2009-06-18

    A granular instability driven by gravity is studied experimentally and numerically. The instability arises as grains fall in a closed Hele-Shaw cell where a layer of dense granular material is positioned above a layer of air. The initially flat front defined by the grains subsequently develops into a pattern of falling granular fingers separated by rising bubbles of air. A transient coarsening of the front is observed right from the start by a finger merging process. The coarsening is later stabilized by new fingers growing from the center of the rising bubbles. The structures are quantified by means of Fourier analysis and quantitative agreement between experiment and computation is shown. This analysis also reveals scale invariance of the flow structures under overall change of spatial scale.

  4. Nonlinear evolution of Buneman instability

    NASA Astrophysics Data System (ADS)

    Ishihara, O.; Hirose, A.; Langdon, A. B.

    1981-03-01

    The nonlinear evolution of one-dimensional electron-ion two-stream instability in a field-free plasma is investigated analytically and by computer simulation. The instability is dominated by the fastest growing mode and its harmonics, provided that the initial fluctuation level is sufficiently small. A nonlinear dispersion relation is obtained and solved numerically, with allowance for the frequency and growth rate modulation, the electric field up to a specified order, and the renormalized particle distribution functions. It is shown that the model can explain computer simulation results, including the presence of an algebraic growth stage following the breakdown of the exponential linear growth, the appearance of harmonics, and the final saturation level.

  5. Ion-beam-driven electrostatic ion cyclotron instabilities

    SciTech Connect

    Miura, A.; Okuda, H.; Ashour-Abdalla, M.

    1982-10-01

    We present results of numerical simulations on the electrostatic ion cyclotron instabilities driven by the ion beam parallel to the magnetic field. For the beam speed exceeding the thermal speed of background ions and the beam temperature much lower than the background ion temperature, it is found that the instability results in strong perpendicular heating and slowing down of parallel drift of the beam ions, leading to the saturation of the instability. Applications to plasma heating and space plasma physics are discussed.

  6. Beam instabilities in the scale-free regime.

    PubMed

    Folli, V; DelRe, E; Conti, C

    2012-01-20

    The instabilities arising in a one-dimensional beam sustained by the diffusive photorefractive nonlinearity in out-of-equilibrium ferroelectrics are theoretically and numerically investigated. In the "scale-free model," in striking contrast with the well-known spatial modulational instability, two different beam instabilities dominate: a defocusing and a fragmenting process. Both are independent of the beam power and are not associated with any specific periodic pattern. PMID:22400741

  7. Reduced modeling of the magnetorotational instability

    NASA Astrophysics Data System (ADS)

    Jamroz, Ben F.

    2009-06-01

    well understood. Many recent numerical investigations of this problem are performed in a local domain, where the global cylindrical background state is projected onto a local Cartesian domain. The resulting system is then numerically modeled within a "shearing box" framework to obtain estimates of angular momentum transport and therefore accretion. However, the simplified geometry of the local domain, and the projection of global quantities leads to a model where the instability is able to grow unboundedly. Utilizing disparate characteristic scales, this thesis presents a reduced asymptotic model for the magnetorotational instability that allows a large scale feedback of local stresses (Reynolds, Maxwell and mixed) onto the projected background state. This system is investigated numerically to determine the impact of allowing this feedback on the saturated level of angular momentum transport.

  8. Baroclinic instability in stellar radiation zones

    SciTech Connect

    Kitchatinov, L. L.

    2014-03-20

    Surfaces of constant pressure and constant density do not coincide in differentially rotating stars. Stellar radiation zones with baroclinic stratification can be unstable. Instabilities in radiation zones are of crucial importance for angular momentum transport, mixing of chemical species, and, possibly, for magnetic field generation. This paper performs linear analysis of baroclinic instability in differentially rotating stars. Linear stability equations are formulated for differential rotation of arbitrary shape and then solved numerically for rotation nonuniform in radius. As the differential rotation increases, r- and g-modes of initially stable global oscillations transform smoothly into growing modes of baroclinic instability. The instability can therefore be interpreted as stability loss to r- and g-modes excitation. Regions of stellar parameters where r- or g-modes are preferentially excited are defined. Baroclinic instability onsets at a very small differential rotation of below 1%. The characteristic time of instability growth is about 1000 rotation periods. Growing disturbances possess kinetic helicity. Magnetic field generation by the turbulence resulting from baroclinic instability in differentially rotating radiation zones is therefore possible.

  9. Perturbation solutions of combustion instability problems

    NASA Technical Reports Server (NTRS)

    Googerdy, A.; Peddieson, J., Jr.; Ventrice, M.

    1979-01-01

    A method involving approximate modal analysis using the Galerkin method followed by an approximate solution of the resulting modal-amplitude equations by the two-variable perturbation method (method of multiple scales) is applied to two problems of pressure-sensitive nonlinear combustion instability in liquid-fuel rocket motors. One problem exhibits self-coupled instability while the other exhibits mode-coupled instability. In both cases it is possible to carry out the entire linear stability analysis and significant portions of the nonlinear stability analysis in closed form. In the problem of self-coupled instability the nonlinear stability boundary and approximate forms of the limit-cycle amplitudes and growth and decay rates are determined in closed form while the exact limit-cycle amplitudes and growth and decay rates are found numerically. In the problem of mode-coupled instability the limit-cycle amplitudes are found in closed form while the growth and decay rates are found numerically. The behavior of the solutions found by the perturbation method are in agreement with solutions obtained using complex numerical methods.

  10. Extracting dynamic topography from river profiles and cosmogenic nuclide geochronology in the Middle Atlas and the High Plateaus of Morocco

    NASA Astrophysics Data System (ADS)

    Pastor, Alvar; Babault, Julien; Owen, Lewis A.; Teixell, Antonio; Arboleya, María-Luisa

    2015-11-01

    The Moulouya river system has intensely eroded the Arhbalou, Missour, and Guercif Neogene foreland basins in northeastern Morocco, having changed from net aggradation during the Miocene-early Pliocene to net incision punctuated by alluvial fan deposition at late Pliocene or early Quaternary time. This region as a whole has experienced mantle-driven, surface uplift (dynamic topography) since the late Cenozoic, being locally affected by uplift due to crustal shortening and thickening of the Middle Atlas too. Knickpoints located along the major streams of the Moulouya fluvial network, appear on both the undeformed margins of the Missour and Guercif foreland basins (High Plateaus), as well as along the thrust mountain front of the southern Middle Atlas, where they reach heights of 800-1000 m. 500-550 m of the knickpoint vertical incision might be explained by long-wavelength mantle-driven dynamic surface uplift, whereas the remaining 450-500 m in the southern Middle Atlas front and 200-300 m in the northeastern Middle Atlas front seem to be thrust-related uplift of the Jebel Bou Naceur. Be-10 terrestrial cosmogenic nuclides have been used to date two Quaternary river terraces in the Chegg Ard valley at 62 ± 14 ka and 411 ± 55 ka. The dated terraces allow the incision rates associated with the frontal structures of the Middle Atlas to be estimated at ~ 0.3 mm yr- 1. Furthermore, these ages have served to evaluate mantle-driven regional surface uplift since the middle Pleistocene in the central Missour basin, yielding values of ~ 0.1-0.2 mm yr- 1.