Geometric transitions and D-term SUSY breaking (open access)

Geometric transitions and D-term SUSY breaking

We propose a new way of using geometric transitions to study metastable vacua in string theory and certain confining gauge theories. The gauge theories in question are N=2 supersymmetric theories deformed to N=1 by superpotential terms. We first geometrically engineer supersymmetry-breaking vacua by wrapping D5 branes on rigid 2-cycles in noncompact Calabi-Yau geometries, such that the central charges of the branes are misaligned. In a limit of slightly misaligned charges, this has a gauge theory description, where supersymmetry is broken by Fayet-Iliopoulos D-terms. Geometric transitions relate these configurations to dual Calabi-Yaus with fluxes, where H_RR, H_NS and dJ are all nonvanishing. We argue that the dual geometry can be effectively used to study the resulting non-supersymmetric, confining vacua
Date: November 5, 2007
Creator: Aganagic, Mina; Aganagic, Mina & Beem, Christopher
System: The UNT Digital Library
Search for Invisible Decays of a Light Scalar in Radiative Transitions Y(3S)->gamma A0 (open access)

Search for Invisible Decays of a Light Scalar in Radiative Transitions Y(3S)->gamma A0

We search for a light scalar particle produced in single-photon decays of the {Upsilon}(3S) resonance through the process {Upsilon}(3S) {yields} {gamma} + A{sup 0}, A{sup 0} {yields} invisible. Such an object appears in Next-to-Minimal Supersymmetric extensions of the Standard Model, where a light CP-odd Higgs boson naturally couples strongly to b-quarks. If, in addition, there exists a light, stable neutralino, decays of A{sup 0} could be preferentially to an invisible final state. We search for events with a single high-energy photon and a large missing mass, consistent with a 2-body decay of {Upsilon}(3S). We find no evidence for such processes in a sample of 122 x 10{sup 6} {Upsilon}(3S) decays collected by the BABAR collaboration at the PEP-II B-factory, and set 90% C.L. upper limits on the branching fraction {Beta}({Upsilon}(3S) {yields} {gamma}A{sup 0}) x {Beta}(A{sup 0} {yields} invisible) at (0.7-31) x 10{sup -6} in the mass range m{sub A{sup 0}} {le} 7.8 GeV. The results are preliminary.
Date: November 5, 2008
Creator: Aubert, B.
System: The UNT Digital Library
Evidence for B Semileptonic Decays into the Lambda_c Charm Baryon (open access)

Evidence for B Semileptonic Decays into the Lambda_c Charm Baryon

We present the first evidence for B semileptonic decays into the charmed baryon {Lambda}{sub c}{sup +} based on 420 fb{sup -1} of data collected at the {Upsilon}(4S) resonance with the BABAR detector at the PEP-II e{sup +}e{sup -} storage rings. Events are tagged by fully reconstructing one of the B mesons in a hadronic decay mode. We measure the relative branching fraction {Beta}({bar B} {yields} {Lambda}{sub c}{sup +} X{ell}{sup -}{bar {nu}}{sub {ell}})/{Beta}({bar B} {yields} {Lambda}{sub c}{sup +}/{bar {Lambda}}{sub c}{sup -}X) = (3.2 {+-} 0.9{sub stat.} {+-} 0.9{sub syst.})%. The significance of the signal including the systematic uncertainty is 4.9 standard deviations.
Date: November 5, 2008
Creator: Aubert, Bernard; Bona, M.; Karyotakis, Y.; Lees, J. P.; Poireau, V.; Prencipe, E. et al.
System: The UNT Digital Library
Measurement of the Branching Fractions of Bbar -> D** l^- nubar_l Decays in Events Tagged by a Fully Reconstructed B Meson (open access)

Measurement of the Branching Fractions of Bbar -> D** l^- nubar_l Decays in Events Tagged by a Fully Reconstructed B Meson

We report a measurement of the branching fractions of {bar B} {yields} D** {ell}{sup -}{bar {nu}}{sub {ell}} decays based on 417 fb{sup -1} of data collected at the {Upsilon}(4S) resonance with the BABAR detector at the PEP-II e{sup +}e{sup -} storage rings. Events are selected by fully reconstructing one of the B mesons in a hadronic decay mode. A fit to the invariant mass differences m(D{sup (*)})-m(D{sup (*)}) is performed to extract the signal yields of the different D** states. We observe the {bar B} {yields} D**{ell}{sup -}{bar {nu}}{sub {ell}} decay modes corresponding to the four D** states predicted by Heavy Quark Symmetry with a significance greater than six standard deviations including systematic uncertainties.
Date: November 5, 2008
Creator: Aubert, Bernard; Bona, M.; Karyotakis, Y.; Lees, J. P.; Poireau, V.; Prencipe, E. et al.
System: The UNT Digital Library
Asymptotic Diffusion-Limit Accuracy of Sn Angular Differencing Schemes (open access)

Asymptotic Diffusion-Limit Accuracy of Sn Angular Differencing Schemes

In a previous paper, Morel and Montry used a Galerkin-based diffusion analysis to define a particular weighted diamond angular discretization for S{sub n}n calculations in curvilinear geometries. The weighting factors were chosen to ensure that the Galerkin diffusion approximation was preserved, which eliminated the discrete-ordinates flux dip. It was also shown that the step and diamond angular differencing schemes, which both suffer from the flux dip, do not preserve the diffusion approximation in the Galerkin sense. In this paper we re-derive the Morel and Montry weighted diamond scheme using a formal asymptotic diffusion-limit analysis. The asymptotic analysis yields more information than the Galerkin analysis and demonstrates that the step and diamond schemes do in fact formally preserve the diffusion limit to leading order, while the Morel and Montry weighted diamond scheme preserves it to first order, which is required for full consistency in this limit. Nonetheless, the fact that the step and diamond differencing schemes preserve the diffusion limit to leading order suggests that the flux dip should disappear as the diffusion limit is approached for these schemes. Computational results are presented that confirm this conjecture. We further conjecture that preserving the Galerkin diffusion approximation is equivalent to preserving the …
Date: November 5, 2009
Creator: Bailey, T S; Morel, J E & Chang, J H
System: The UNT Digital Library
Determination of Dark Matter Properties at High-Energy Colliders (open access)

Determination of Dark Matter Properties at High-Energy Colliders

If the cosmic dark matter consists of weakly-interacting massive particles, these particles should be produced in reactions at the nextgeneration of high-energy accelerators. Measurements at these accelerators can then be used to determine the microscopic properties of the dark matter. From this, we can predict the cosmic density, the annihilation cross sections, and the cross sections relevant to direct detection. In this paper, we present studies in supersymmetry models with neutralino dark matter that give quantitative estimates of the accuracy that can be expected. We show that these are well matched to the requirements of anticipated astrophysical observations of dark matter. The capabilities of the proposed International Linear Collider (ILC) are expected to play a particularly important role in this study.
Date: November 5, 2006
Creator: Baltz, Edward A.; Battaglia, Marco; Peskin, Michael E. & Wizansky, Tommer
System: The UNT Digital Library
Holographic Systematics of D-brane Inflation (open access)

Holographic Systematics of D-brane Inflation

We provide a systematic treatment of possible corrections to the inflaton potential for D-brane inflation in the warped deformed conifold. We consider the D3-brane potential in the presence of the most general possible corrections to the throat geometry sourced by coupling to the bulk of a compact Calabi-Yau space. This corresponds to the potential on the Coulomb branch of the dual gauge theory, in the presence of arbitrary perturbations of the Lagrangian. The leading contributions arise from perturbations by the most relevant operators that do not destroy the throat geometry. We find a generic contribution from a non-chiral operator of dimension {Delta} = 2 associated with a global symmetry current, resulting in a negative contribution to the inflaton mass-squared. If the Calabi-Yau preserves certain discrete symmetries, this is the dominant correction to the inflaton potential, and fine-tuning of the inflaton mass is possible. In the absence of such discrete symmetries, the dominant contribution comes from a chiral operator with {Delta} = 3/2, corresponding to a {phi}{sup 3/2} term in the inflaton potential. The resulting inflationary models are phenomenologically identical to the inflection point scenarios arising from specific D7-brane embeddings, but occur under far more general circumstances. Our strategy extends immediately …
Date: November 5, 2008
Creator: Baumann, Daniel; /Harvard U., Phys. Dept. /Princeton U.; Dymarsky, Anatoly; /Stanford U., Phys. Dept.; Kachru, Shamit; /Stanford U., Phys. Dept. /SLAC et al.
System: The UNT Digital Library
The Emma Commissioning Procedure (open access)

The Emma Commissioning Procedure

The author begins with a brief review of the goals of the EMMA experiment. He then describe two stages of EMMA commissioning. The first stage is simply to get the beam to circulate a full turn in the ring, and is done only once during the course of the experiment. The second stage will be repeated several times, at least once for each lattice configuration, and involves two parts: setting the required values for the machine parameters, and determining the tunes and time of flight as a function of energy.
Date: November 5, 2007
Creator: Berg, J. S.
System: The UNT Digital Library
A HAMILTONIAN FORMULATION FOR SPIRAL-SECTOR ACCELERATORS. (open access)

A HAMILTONIAN FORMULATION FOR SPIRAL-SECTOR ACCELERATORS.

I develop a formulation for Hamiltonian dynamics in an accelerator with magnets whose edges follow a spiral. I demonstrate using this Hamiltonian that a spiral FFAG can be made perfectly 'scaling'. I examine the effect of tilting an RF cavity with respect a radial line from the center of the machine, potentially with a different angle than the spiral of the magnets.
Date: November 5, 2007
Creator: Berg, J. S.
System: The UNT Digital Library
Effusive-flow characterization of arbitrary size and geometry target/vapor transport systems: radioactive ion beam applications (open access)

Effusive-flow characterization of arbitrary size and geometry target/vapor transport systems: radioactive ion beam applications

The principal factors that severely limit intensities of short-lived radioactive ion beams produced by the Isotope Separator On-Line (ISOL) technique are time delays due to diffusion of radioactive species from solid or liquid target materials and their effusive-flow transport to the ion source. Although diffusion times can be reduced by proper design of short diffusion length, highly refractory targets, effusive-flow times are more difficult to assess. After diffusion from the target material, the species must travel through the target material and vapor transport system to the ion source. The time required for effusive-flow transport to the ion source depends on the conduction path, chemical reactions between the species and target material and materials of construction as well as the physical size and geometry of the transport system. We have developed a fast-valve (1 ms closing time) for introducing gaseous or vapor-state species into the target/vapor transport/ion source/system th at permits measurement of effusive-flow times for any gaseous or vaporous species (chemically active or chemically inactive) through any vapor transport system, independent of size and geometry. Characteristic times are determined from the exponential decay of the momentum analyzed ion beam intensity for the species.during effusive-flow through the vapor transport system under …
Date: November 5, 2001
Creator: Bilheux, J.-C.
System: The UNT Digital Library
Progress in Infrared Pyrometry Measurements of Shocked Solids (open access)

Progress in Infrared Pyrometry Measurements of Shocked Solids

Temperature measurement is one of the grand challenges still facing experimental shock physics. A shock experiment fundamentally measures E({sigma}{sub x}, {var_epsilon}{sub 11}) which is an incomplete equation of state since temperature (or entropy) remains unspecified. Ideally, one would like to experimentally determine a free energy F(T, {var_epsilon}{sub ij}) from which all other thermo-mechanical properties might be derived. In practice, temperature measurement would allow direct comparison with theory/simulation since T and {var_epsilon}{sub 11} are in most theories the underlying variables. Temperature is a sensitive measure of energy partitioning, knowledge of which would increase our understanding phase boundaries and thermally activated processes (such as chemical reactivity (including dissociation and ionization)). Temperature measurement would also allow a thermodynamically consistent coupling of hydrodynamic equations of state to the material's constitutive (deformation) behavior. The measurement of the temperature of a material that has undergone severe strains at small time-scales is extremely difficult, and we are developing a method using infrared reflectance and pyrometry. The emitted power from a warm surface is measured over a range of wavelengths using a multi-channel IR detector with a response time of {approx}0.1 {micro}s. Each channel of the detector passes the radiation from a selected wavelength interval into a detector. …
Date: November 5, 2001
Creator: Cazamias, J U; Hare, D E & Poulsen, P
System: The UNT Digital Library
Time-Reversal Analysis for Scatterer Characterization (open access)

Time-Reversal Analysis for Scatterer Characterization

Only the simplest monopole scattering behavior has usually been treated in previous time-reversal analyses. A new application of time-reversal processing of wave scattering data permits characterization of scatterers by analyzing the number and nature of the singular functions (or eigenfunctions) associated with individual scatterers when they have multiple contributions from monopole, dipole and/or quadrupole scattering terms. We discuss acoustic, elastic, and electromagnetic scattering problems for low frequencies (ka < 1, k being the wavenumber and a the radius of the scatterer). Specific examples for electromagnetic scattering from one of a number of small conducting spheres show that each sphere can have up to six distinct time-reversal eigenfunctions associated with it.
Date: November 5, 2003
Creator: Chambers, D H & Berryman, J G
System: The UNT Digital Library
A PIV Methodology for High-Resolution Measurement of Flow Statistics (open access)

A PIV Methodology for High-Resolution Measurement of Flow Statistics

Particle-image velocimetry (PIV) is a flow-diagnostic technique that provides velocity fields from a comparison of images of particulate-laden flow. We have developed a PIV processing methodology that extracts measurements of the particle-displacement histogram from a flow video or ensemble of flow-image pairs. Single-pixel measurement of mean velocity can be obtained from an ensemble of {Omicron}10{sup 3} images. Measurements of higher-order moments of the velocity histogram require spatial averaging (i.e., lower spatial resolution), larger ensembles of images, or a combination of the two. We present single-pixel-resolution PIV measurements of a steady microflow and high-resolution measurements of the velocity histogram of a stationary turbulent flow. This methodology has applications in quantifying velocity statistics in other stochastic flows, e.g., bulk and near-wall boiling.
Date: November 5, 2000
Creator: Cummings, Eric B.; Schefer, Robert W. & Chung, Jacob N.
System: The UNT Digital Library
Initial test results of an ionization chamber shower detector for a LHC luminosity monitor (open access)

Initial test results of an ionization chamber shower detector for a LHC luminosity monitor

A novel, segmented, multi-gap, pressurized gas ionization chamber is being developed for optimization of the luminosity of the LHC. The ionization chambers are to be installed in the front quadrupole and zero degree neutral particle absorbers in the high luminosity IRs and sample the energy deposited near the maxima of the hadronic/electromagnetic showers in these absorbers. The ionization chambers are instrumented with low noise, fast, pulse shaping electronics to be capable of resolving individual bunch crossings at 40 MHz. In this paper we report the initial results of our second test of this instrumentation in an SPS external proton beam. Single 300 GeV protons are used to simulate the hadronic/electromagnetic shower produced by the forward collision products from the interaction regions of the LHC. The capability of instrumentations to measure the luminosity of individual bunches in a 40 MHz bunch train is demonstrated.
Date: November 5, 2002
Creator: Datte, P.; Beche, J.-F.; Haguenauer, M.; Manfredi, P.F.; Manghisoni, M.; Millaud, J. et al.
System: The UNT Digital Library
Optimization of signal extraction and front-end design in a fast,multigap ionization chamber (open access)

Optimization of signal extraction and front-end design in a fast,multigap ionization chamber

This paper discusses the criteria that have been adopted tooptimize the signal processing in a shower detector to be employed as LHCbeam luminosity monitor. The original aspect ofthis instrument is itsablility to operate on a bunch-by-bunch basis. This means that it mustperform accurate charge measurements at a repetition rate of 40 MHz. Thedetector must withstand an integrated dose of 100 Grad, that is, two tothree orders of magnitude beyond those expected in the experiments. Tomeet the above requirements, an ionization chamber consisting of severalgaps of thickness 0.5 mm, filled with a gas that is expected to beradiation resistant, has been designed. Crucial in the development of thesystem is the signal processing, as the electronic noise may set thedominant limitation to the accuracy of the measurement. This is relatedto two aspects. One is the short time available for the chargemeasurement. The second one is the presence of a few meter cable betweenthe detector and the preamplifier, as this must be located out of theregion of highest radiation field. Therefore the optimization of thesignal-to-noise ratio requires that the best configuration of the chambergaps be determined under the constraint of the presence of a cable ofnon-negligible length between detector and preamplifier. The remoteplacement …
Date: November 5, 2001
Creator: Datte, P.S.; Manfredi, P.F.; Millaud, J.E.; Placidi, M.; Ratti,L.; Speziali, V. et al.
System: The UNT Digital Library
ACCUMULATION OF RADIOCESIUM BY MUSHROOMS IN THE ENVIRONMENT: A LITERATURE REVIEW AND IMAGE GALLERY (open access)

ACCUMULATION OF RADIOCESIUM BY MUSHROOMS IN THE ENVIRONMENT: A LITERATURE REVIEW AND IMAGE GALLERY

During the last 50 years, a large amount of information on radionuclide accumulators or 'sentinel-type' organisms in the environment has been published. Much of this work focused on the risks of food-chain transfer of radionuclides to higher organisms such as reindeer and man. However, until the 1980's and 1990's, there has been little published data on the radiocesium ({sup 134}Cs and {sup 137}Cs) accumulation by mushrooms. This presentation will consist of a review of the published data for {sup 134,137}Cs accumulation by mushrooms in nature. The review will consider the time of sampling, sample location characteristics, the radiocesium source term and other aspects that promote {sup 134,137}Cs uptake by mushrooms. This review will focus on published data for mushrooms that demonstrate a large propensity for use in the environmental biomonitoring of radiocesium contamination. It will also provide photographs and descriptions of habitats for many of these mushrooms to facilitate their collection for biomonitoring.
Date: November 5, 2006
Creator: Duff, M & Mary Ramsey, M
System: The UNT Digital Library
Long Term Corrosion Potential and Corrosion Rate of Creviced Alloy 22 in Chloride Plus Nitrate Brines (open access)

Long Term Corrosion Potential and Corrosion Rate of Creviced Alloy 22 in Chloride Plus Nitrate Brines

Alloy 22 is a nickel base alloy highly resistant to all forms of corrosion. In conditions where tight crevices exist in hot chloride containing solutions and at anodic potentials, Alloy 22 may suffer crevice corrosion, a form of localized attack. The occurrence (or not) of crevice corrosion in a given environment (e.g. salt concentration and temperature), is governed by the values of the critical potential (E{sub crit}) for crevice corrosion and the corrosion potential (E{sub corr}) that the alloy may establish in the studied environment. If E{sub corr} is equal or higher than E{sub crit}, crevice corrosion may be expected. In addition, it is generally accepted that as Alloy 22 becomes passive in a certain environment, its E{sub corr} increases and its corrosion rate (CR) decreases. This paper discusses the evolution of E{sub corr} and corrosion rate (CR) of creviced Alloy 22 specimens in six different mixtures of sodium chloride (NaCl) and potassium nitrate (KNO{sub 3}) at 100 C. The effect of immersion time on the value of E{sub crit} was also determined. Two types of specimens were used, polished as-welded (ASW) and as-welded plus solution heat-treated (ASW+SHT). The latter contained the black annealing oxide film on the surface. Results …
Date: November 5, 2005
Creator: Evans, K. J.; Stuart, M. L.; Etien, R. A.; Hust, G. A.; Estill, J. C. & Rebak, R. B.
System: The UNT Digital Library
Constrained Walks and Self-Avoiding Walks: Implications for Protein Structure Determination (open access)

Constrained Walks and Self-Avoiding Walks: Implications for Protein Structure Determination

None
Date: November 5, 2000
Creator: FAULON,JEAN-LOUP; RINTOUL,MARK DANIEL & YOUNG,MALIN M.
System: The UNT Digital Library
AIR PERMIT COMPLIANCE FOR WASTE RETRIEVAL OEPRATIONS INVOLVING MULTI-UNIT OPERATIONS (open access)

AIR PERMIT COMPLIANCE FOR WASTE RETRIEVAL OEPRATIONS INVOLVING MULTI-UNIT OPERATIONS

Since 1970, approximately 38,000 suspect-transuranic and transuranic waste containers have been placed in retrievable storage on the Hanford Site in the 200 Areas burial grounds. Hanford's Waste Retrieval Project is retrieving these buried containers and processing them for safe storage and disposition. Container retrieval activities require an air emissions permit to account for potential emissions of radionuclides. The air permit covers the excavation activities as well as activities associated with assaying containers and installing filters in the retrieved transuranic containers lacking proper venting devices. Fluor Hanford, Inc. is required to track radioactive emissions resulting from the retrieval activities. Air, soil, and debris media contribute to the emissions and enabling assumptions allow for calculation of emissions. Each of these activities is limited to an allowed annual emission (per calendar year) and .contributes to the overall total emissions allowed for waste retrieval operations. Tracking these emissions is required to ensure a permit exceedance does not occur. A tracking tool was developed to calculate potential emissions in real time sense. Logic evaluations are established within the tracking system to compare real time data against license limits to ensure values are not exceeded for either an individual activity or the total limit. Data input …
Date: November 5, 2007
Creator: FM, SIMMONS
System: The UNT Digital Library
Far-from-equilibrium measurements of thermodynamic length (open access)

Far-from-equilibrium measurements of thermodynamic length

Thermodynamic length is a path function that generalizes the notion of length to the surface of thermodynamic states. Here, we show how to measure thermodynamic length in far-from-equilibrium experiments using the work fluctuation relations. For these microscopic systems, it proves necessary to define the thermodynamic length in terms of the Fisher information. Consequently, the thermodynamic length can be directly related to the magnitude of fluctuations about equilibrium. The work fluctuation relations link the work and the free energy change during an external perturbation on a system. We use this result to determine equilibrium averages at intermediate points of the protocol in which the system is out-of-equilibrium. This allows us to extend Bennett's method to determine the potential of mean force, as well as the thermodynamic length, in single molecule experiments.
Date: November 5, 2008
Creator: Feng, Edward H. & Crooks, Gavin E.
System: The UNT Digital Library
Magnetic Detachment and Plume Control in Escaping Magnetized Plasma (open access)

Magnetic Detachment and Plume Control in Escaping Magnetized Plasma

The model of two-fluid, axisymmetric, ambipolar magnetized plasma detachment from thruster guide fields is extended to include plasmas with non-zero injection angular velocity profiles. Certain plasma injection angular velocity profiles are shown to narrow the plasma plume, thereby increasing exhaust efficiency. As an example, we consider a magnetic guide field arising from a simple current ring and demonstrate plasma injection schemes that more than double the fraction of useful exhaust aperture area, more than halve the exhaust plume angle, and enhance magnetized plasma detachment.
Date: November 5, 2008
Creator: Fisch, P. F. Schmit and N. J.
System: The UNT Digital Library
Electronic and magnetic properties of zincblende half-metal superlattices (open access)

Electronic and magnetic properties of zincblende half-metal superlattices

Zincblende half-metallic compounds such as CrAs, with large magnetic moments and high Curie temperatures, are promising materials for spintronic applications. They explore layered materials, consisting of alternating layers of zincblende half-metals, by first principles calculations, and find that superlattices of (CrAs){sub 1}(MnAs){sub 1} and (CrAs){sub 2}(MnAs){sub 2} are half-metallic with magnetic moments of 7.0{mu}{sub B} and 14.0{mu}{sub B} per unit cell, respectively. They discuss the nature of the bonding and half-metallicity in these materials and, based on the understanding acquired, develop a simple expression for the magnetic moment in such materials. They explore the range of lattice constants over which half-metallicity is manifested, and suggest corresponding substrates for growth in thin film form.
Date: November 5, 2003
Creator: Fong, C Y; Qian, M C; Pask, J; Yang, L H & Dag, S
System: The UNT Digital Library
Raman spectroscopy under extreme conditions (open access)

Raman spectroscopy under extreme conditions

We report the results of Raman measurements of various materials under simultaneous conditions of high temperature and high pressure in the diamond anvil cell (DAC). High temperatures are generated by laser heating or internal resistive (ohmic) heating or a combination of both. We present Raman spectra of cubic boron nitride (cBN) to 40 GPa and up to 2300 K that show a continuous pressure and temperature shift of the frequency of the transverse optical mode. We have also obtained high-pressure Raman spectra from a new noble metal nitride, which we synthesized at approximately 50 GPa and 2000 K. We have obtained high-temperature spectra from pure nitrogen to 39 GPa and up to 2000 K, which show the presence of a hot band that has previously been observed in CARS measurements. These measurements have also allowed us to constrain the melting curve and to examine changes in the intramolecular potential with pressure.
Date: November 5, 2004
Creator: Goncharov, A F & Crowhurst, J C
System: The UNT Digital Library
A Quantitative Comparison of Numerical Methods for the Compressible Euler Equations: Fifth-order WENO and Piecewise-Linear Godunov (open access)

A Quantitative Comparison of Numerical Methods for the Compressible Euler Equations: Fifth-order WENO and Piecewise-Linear Godunov

A numerical study is undertaken comparing a fifth-order version of the weighted essentially non-oscillatory numerical (WENO5) method to a modern piecewise-linear, second-order, version of Godunov's (PLMDE) method for the compressible Euler Equations. A series of one-dimensional test problems are examined beginning with classical linear problems and ending with complex shock interactions. The problems considered are: (1) linear advection of a Gaussian pulse in density, (2) Sod's shock tube problem, (3) the ''peak'' shock tube problem, (4) a version of the Shu and Osher shock entropy wave interaction and (5) the Woodward and Colella interacting shock wave problem. For each problem and method, run times, density error norms and convergence rates are reported for each method as produced from a common code test-bed. The linear problem exhibits the advertised convergence rate for both methods as well as the expected large disparity in overall error levels; WENO5 has the smaller errors and an enormous advantage in overall efficiency (in accuracy per unit CPU time). For the nonlinear problems with discontinuities, however, we generally see both first-order self-convergence of error as compared to an exact solution, or when an analytic solution is not available, a converged solution generated on an extremely fine grid. …
Date: November 5, 2003
Creator: Greenough, J A & Rider, W J
System: The UNT Digital Library