7,147 Matching Results

Search Results

Advanced search parameters have been applied.

Simulating relativistic beam and plasma systems using an optimal boosted frame

Description: It was shown recently that it may be computationally advantageous to perform computer simulations in a Lorentz boosted frame for a certain class of systems. However, even if the computer model relies on a covariant set of equations, it was pointed out that algorithmic difficulties related to discretization errors may have to be overcome in order to take full advantage of the potential speedup. In this paper, we summarize the findings, the difficulties and their solutions, and review the applications of the technique that have been performed to date.
Date: May 1, 2009
Creator: Vay, J. L.; Bruhwiler, D. L.; Geddes, C. G. R.; Fawley, W. M.; Martins, S. F.; Cary, J. R. et al.
Partner: UNT Libraries Government Documents Department

Manifold methods for methane combustion

Description: Objective is to develop a new method for studying realistic chemistry in turbulent methane combustion with NO{sub x} mechanism. The realistic chemistry is a simplification to a more detailed chemistry based on the manifold method; accuracy is determined by interaction between the transport process and the chemical reaction. In this new (tree) method, probability density function or partially stirred reactor calculations are performed. Compared with the reduced mechanism, manifold, and tabulation methods, the new method overcomes drawbacks of the reduced mechanism method and preserves the advantages of the manifold method. Accuracy is achieved by specifying the size of the cell.
Date: December 31, 1995
Creator: Yang, B. & Pope, S.B.
Partner: UNT Libraries Government Documents Department

Development of an advanced system identification technique for comparing ADAMS analytical results with modal test data for a MICON 65/13 wind turbine

Description: This work uses the theory developed in NREL/TP--442-7110 to analyze simulated data from an ADAMS (Automated Dynamic Analysis of Mechanical Systems) model of the MICON 65/13 wind turbine. The Observer/Kalman Filter identification approach is expanded to use input-output time histories from ADAMS simulations or structural test data. A step by step outline is offered on how the tools developed in this research, can be used for validation of the ADAMS model.
Date: July 1, 1995
Creator: Bialasiewicz, J.T.
Partner: UNT Libraries Government Documents Department

Using EnergyPlus for California Title-24 compliancecalculations

Description: For the past decade, the non-residential portion of California's Title-24 building energy standard has relied on DOE-2.1E as the reference computer simulation program for development as well as compliance. However, starting in 2004, the California Energy Commission has been evaluating the possible use of Energy Plus as the reference program in future revisions of Title-24. As part of this evaluation, the authors converted the Alternate Compliance Method (ACM) certification test suite of 150 DOE-2 files to Energy Plus, and made parallel DOE-2 and Energy Plus runs for this extensive set of test cases. A customized version of DOE-2.1E named doe2ep was developed to automate the conversion process. This paper describes this conversion process, including the difficulties in establishing an apples-to-apples comparison between the two programs, and summarizes how the DOE-2 and Energy Plus results compare for the ACM test cases.
Date: August 26, 2006
Creator: Huang, Joe; Bourassa, Norman; Buhl, Fred; Erdem, Ender & Hitchcock, Rob
Partner: UNT Libraries Government Documents Department

Computer-Assisted Discovery and Proof

Description: With the advent of powerful, widely-available mathematical software, combined with ever-faster computer hardware, we are approaching a day when both the discovery and proof of mathematical facts can be done in a computer-assisted manner. his article presents several specific examples of this new paradigm in action.
Date: December 10, 2007
Creator: Bailey, David H. & Borwein, Jonathan M.
Partner: UNT Libraries Government Documents Department

Playing With Sandpiles.

Description: The Bak-Tang-Wiesenfeld sandpile model provides a simple and elegant system with which to demonstrate self-organized criticality. This model has rather remarkable mathematical properties first elucidated by Dhar. I demonstrate some of these properties graphically with a simple computer simulation.
Date: January 21, 2003
Creator: Creutz, M.
Partner: UNT Libraries Government Documents Department

ALE advantage in hypervelocity impact calculations

Description: The ALE3D code is used to model experiments relevant to hypervelocity impact lethality, carried out in the 4-5 km/s velocity range. The code is run in the Eulerian and ALE modes. Zoning in the calculations is refined beyond the level found in most lethality calculations, but still short of convergence. The level of zoning refinement that produces equivalent results in uniformly zoned Eulerian calculations and ALE ones utilizing specialized zoning, weighting and relaxation techniques is established. It takes 11 times fewer zones and about 60% as many cycles when ALE capabilities are used. Calculations are compared to experimental results.
Date: October 1, 1998
Creator: Gerassimenko, M. & Rathkopf, J.
Partner: UNT Libraries Government Documents Department

Simulations of implosions with a 3D, parallel, unstructured-grid, radiation-hydrodynamics code

Description: An unstructured-grid, radiation-hydrodynamics code is used to simulate implosions. Although most of the problems are spherically symmetric, they are run on 3D, unstructured grids in order to test the code�s ability to maintain spherical symmetry of the converging waves. Three problems, of increasing complexity, are presented. In the first, a cold, spherical, ideal gas bubble is imploded by an enclosing high pressure source. For the second, we add non-linear heat conduction and drive the implosion with twelve laser beams centered on the vertices of an icosahedron. In the third problem, a NIF capsule is driven with a Planckian radiation source.
Date: December 28, 1998
Creator: Kaiser, T. B.; Milovich, J. L.; Prasad, M. K.; Rathkopf, J. & Shestakov, A. I.
Partner: UNT Libraries Government Documents Department

New computational method for non-LTE, the linear response matrix

Description: We investigate non-local thermodynamic equilibrium atomic kinetics using nonequilibrium thermodynamics and linear response theory. This approach gives a rigorous general framework for exploiting results from non-LTE kinetic calculations and offers a practical data-tabulation scheme suitable for use in plasma simulation codes. We describe how this method has been implemented to supply a fast and accurate non-LTE option in Lasnex.
Date: October 1, 1998
Creator: Fournier, K. B.; Graziani, F. R.; Harte, J. A.; Libby, S. B.; More, R. M.; Rathkopf, J. et al.
Partner: UNT Libraries Government Documents Department

Dynamic simulation of the Hanford tank waste remediation system

Description: Cleaning up and disposing of approximately 50 years of nuclear waste is the main mission at the U.S. Department of Energy`s Hanford Nuclear Reservation, located in the southeastern part of the state of Washington. A major element of the total cleanup effort involves retrieving, processing, and disposing of radioactive and hazardous waste stored in 177 underground storage tanks. This effort, referred to as the Tank Waste Remediation System (TWRS), is expected to cost billions of dollars and take approximately 25 years to complete. Several computer simulations of this project are being created, focusing on both programmatic and detailed engineering issues. This paper describes one such simulation activity, using the ithink(TM)computer simulation software. The ithink(TM) simulation includes a representation of the complete TWRS cleanup system, from retrieval of waste through intermediate processing and final vitrification of waste for disposal. Major issues addressed to date by the simulation effort include the need for new underground storage tanks to support TWRS activities, and the estimated design capacities for various processing facilities that are required to support legally mandated program commitment dates. This paper discusses how the simulation was used to investigate these questions.
Date: May 3, 1996
Creator: Harmsen, R. W.
Partner: UNT Libraries Government Documents Department

Bundle Binding in Polyelectrolyte Solutions

Description: Stiff polyelectrolytes are found to spontaneously form oriented bundles. Conditions under which bundling occurs are found. Molecular dynamics simulations show that divalent counterions are necessary, and the chains must be sufficiently long and stiff. No aggregation occurs for monovalent counterions. For flexible or short chains aggregation occurs, but bundle formation does not. Due to dynamical constraints the systems tend to order into a network of connected bundles, not a single bundle.
Date: January 21, 1999
Creator: Stevens, M.J.
Partner: UNT Libraries Government Documents Department


Description: We present results obtained with the Odyssey simulation code. Odyssey is a 1, 2, and 3 dimensional AMR code using cartesian, cylindrical, and spherical coordinates. The results provide an interesting snapshot of Odyssey at this point in its development. Results include parallel performance and scaling, Eulerian hydrodynamics algorithm comparisons, ADI based diffusion solvers on hierarchical meshes, ECB treatment of material interfaces in diffusion solves.
Date: October 1, 1998
Creator: Braddy, D.; Brown, S.; Cook, G.; Kueny, C.; Lambert, M.; Peters, D. et al.
Partner: UNT Libraries Government Documents Department

Multi-Domain Surety Modeling and Analysis for High Assurance Systems

Description: Engineering systems are becoming increasingly complex as state of the art technologies am incorporated into designs. Surety modeling and analysis is an emerging science that permits an engineer to qualitatively and quantitatively predict and assess the completeness and predictability of a design. Surety is a term often used in the Department of Defense (DoD) and Department of Energy (DOE) communities, which refers to the integration of safety, security, reliability and performance aspects of design. Current risk assessment technologies for analyzing complex systems fail to adequately describe the problem, thus making assessment fragmented and non-integrated. To address this problem, we have developed a methodology and extensible software tool set to address model integration and complexity for high consequence systems. The MultiGraph Architecture (MGA) facilitates multi-domain, model-integrated modeling and analyses of complex, high-assurance systems. The MGA modeling environment allows the engineer to customize the modeling environment to match a design paradigm representative of the actual design. Previous modeling tools have a point-defined model space that forms the modeler to work in less than optimal environments. Current approaches for the problem to be bounded and constrained by requirements of the modeling tool and not the actual design problem. In some small cases, this is only maximally adequate MM facilitates the implementation of a surety methodology, which is used to represent high assurance systems with respect to safety and reliability. Formal mathematical models am used to correctly describe design safety and reliability functionality and behavioral fictional and behavioral representations of the design w then analyzed using commercial-off-the-shelf tools.
Date: January 27, 1999
Creator: Davis, J.; Martinez, M.; Scott, J. & Sztipanovits, J.
Partner: UNT Libraries Government Documents Department

Axions from string decay

Description: We have studied numerically the evolution and decay of axion strings. These global defects decay mainly by axion emission and thus contribute to the cosmological axion energy density. The relative importance of this source relative to misalignment production of axions depends on the spectrum. Radiation spectra for various string loop configurations are presented. They support the contention that the string decay contribution is of the same order of magnitude as the contribution from misalignment.
Date: July 9, 1998
Creator: Hagmann, C., LLNL
Partner: UNT Libraries Government Documents Department

New Development in WARP: Progress Toward End-to-End Simulation

Description: The development of a high current, heavy-ion beam driver for inertial confinement fusion requires a detailed understanding of the behavior of the beam, including effects of the strong self-fields. The necessity of including the self-fields of the beam makes particle-in-cell (PIC) simulation techniques ideal, and for this reason, the multi-dimensional PIC/accelerator code WARP has been developed. WARP has been used extensively to study the creation and propagation of ion beams both in experiments and for the understanding of basic beam physics. An overview of the structure of the code will be presented along with a discussion of features that make the code an effective tool in the understanding of space-charge dominated beam behavior. Much development has been done on WARP increasing its flexibility and generality. Major additions include a generalized field description, an efficient steady-state modelling technique, a transverse slice model with a bending algorithm, further improvement of the parallel processing version, and capabilities for linking to chamber transport codes. With these additions, the capability of modeling a large scale accelerator from end-to-end comes closer to reality.
Date: November 7, 1997
Creator: Grote, D.P.; Friedman, A.; Haber, I.; Fawley, W. & Luc Vay, J.
Partner: UNT Libraries Government Documents Department

Mixing of Isotactic and Syndiotactic Polypropylenes in the Melt

Description: The miscibility of polypropylene (PP) melts in which the chains differ only in stereochemical composition has been investigated by two different procedures. One approach used detailed local information from a Monte Carlo simulation of a single chain, and the other approach takes this information from a rotational isomeric state model devised decades ago, for another purpose. The first approach uses PRISM theory to deduce the intermolecular packing in the polymer blend, while the second approach uses a Monte Carlo simulation of a coarse-grained representation of independent chains, expressed on a high-coordination lattice. Both approaches find a positive energy change upon mixing isotactic PP (iPP) and syndiotactic polypropylene (sPP) chains in the melt. This conclusion is qualitatively consistent with observations published recently by Muelhaupt and coworkers. The size of the energy chain on mixing is smaller in the MC/PRISM approach than in the RIS/MC simulation, with the smaller energy change being in better agreement with the experiment. The RIS/MC simulation finds no demixing for iPP and atactic polypropylene (aPP) in the melt, consistent with several experimental observations in the literature. The demixing of the iPP/sPP blend may arise from attractive interactions in the sPP melt that are disrupted when the sPP chains are diluted with aPP or iPP chains.
Date: July 14, 2000
Partner: UNT Libraries Government Documents Department


Description: Significant progress continued to be made during the past reporting quarter on both major technical tasks. During the reporting period at OSU, computational investigations were conducted of addition vs. abstraction reactions of H, O(3 P), and OH with monocyclic aromatic hydrocarbons. The potential energy surface for more than 80 unique reactions of H, O ( 3 P), and OH with aromatic hydrocarbons were determined at the B3LYP/6-31G(d) level of theory. The calculated transition state barriers and reaction free energies indicate that the addition channel is preferred at 298K, but that the abstraction channel becomes dominant at high temperatures. The thermodynamic preference for reactivity with aromatic hydrocarbons increases in the order O(3 P) < H < OH. Abstraction from six-membered aromatic rings is more facile than abstraction from five-membered aromatic rings. However, addition to five-membered rings is thermodynamically more favorable than addition to six-membered rings. The free energies for the abstraction and addition reactions of H, O, and OH with aromatic hydrocarbons and the characteristics of the respective transition states can be used to calculate the reaction rate constants for these important combustion reactions. Experimental work at Brown University on the effect of reaction on the structural evolution of different chars (i.e., phenolic resin char and chars produced from three different coals) have been investigated in a TGA/TPD-MS system. It has been found that samples of different age of these chars appeared to lose their "memory" concerning their initial structures at high burn-offs. During the reporting period, thermal desorption experiments of selected samples were conducted. These spectra show that the population of low temperature oxygen surface complexes, which are primarily responsible for reactivity, are more similar for the high burn-off than for the low burn-off samples of different ages; i.e., the population of active sites are more similar for the �younger� ...
Date: January 13, 1999
Partner: UNT Libraries Government Documents Department

Numerical simulations of magnetic reversal in layered spring magnets.

Description: This report summarizes the results of numerical investigations of magnetic reversal in layered spring magnets. A one-dimensional model is used of a film consisting of several atomic layers of soft material on top of several atomic layers of hard material. Each atomic layer is taken to be uniformly magnetized, and spatial inhomogeneities within an atomic layer are neglected. The state of such a system is described by a chain of magnetic spin vectors. Each spin vector behaves like a spinning top driven locally by the effective magnetic field and subject to damping (Landau-Lifshitz-Gilbert equation). A numerical integration scheme for the LLG equation is presented that is unconditionally stable and preserves the magnitude of the magnetization vector at all times. The results of numerical investigations for a bilayer in a rotating in-plane magnetic field show hysteresis with a basic period of 2{pi} at moderate fields and hysteresis with a basic period of {pi} (or any multiple thereof) at strong fields.
Date: January 24, 2001
Creator: Jiang, J.S.; Kaper, H.G. & Leaf, G.K.
Partner: UNT Libraries Government Documents Department

Massively Parallel Direct Simulation of Multiphase Flow

Description: The authors understanding of multiphase physics and the associated predictive capability for multi-phase systems are severely limited by current continuum modeling methods and experimental approaches. This research will deliver an unprecedented modeling capability to directly simulate three-dimensional multi-phase systems at the particle-scale. The model solves the fully coupled equations of motion governing the fluid phase and the individual particles comprising the solid phase using a newly discovered, highly efficient coupled numerical method based on the discrete-element method and the Lattice-Boltzmann method. A massively parallel implementation will enable the solution of large, physically realistic systems.
Date: August 10, 2000
Partner: UNT Libraries Government Documents Department

Study of the transverse beam motion in the DARHT Phase II accelerator

Description: The accelerator for the second-axis of the Dual Axis Radiographic Hydrodynamic Test (DARHT) facility will accelerate a 4-kA, 3-MeV, 2--µs long electron current pulse to 20 MeV. The energy variation of the beam within the flat-top portion of the current pulse is (plus or equal to) 0.5%. The performance of the DARHT Phase II radiographic machine requires the transverse beam motion to be much less than the beam spot size which is about 1.5 mm diameter on the x-ray converter. In general, the leading causes of the transverse beam motion in an accelerator are the beam breakup instability (BBU) and the corkscrew motion. We have modeled the transverse beam motion in the DARHT Phase II accelerator with various magnetic tunes and accelerator cell configurations by using the BREAKUP code. The predicted sensitivity of corkscrew motion and BBU growth to different tuning algorithms will be presented.
Date: August 20, 1998
Creator: Chen, Yu-Jiuan; Fawley, W M & Houck, T L
Partner: UNT Libraries Government Documents Department

Advanced Techniques for Reservoir Simulation and Modeling of Non-Conventional Wells

Description: Research results for the second year of this project on the development of improved modeling techniques for non-conventional (e.g., horizontal, deviated or multilateral) wells were presented. The overall program entails the development of enhanced well modeling and general simulation capabilities. A general formulation for black-oil and compositional reservoir simulation was presented.
Date: August 23, 2001
Creator: Durlofsky, Louis J. & Aziz, Khalid
Partner: UNT Libraries Government Documents Department

Radio Frequency Plasma Applications for Space Propulsion

Description: Recent developments in solid-state radio frequency (RF) power technologies allow for the practical consideration of RF heated plasmas for space propulsion. These technologies permit the use of any electrical power source, de-couple the power and propellant sources, and allow for the effcient use of both the propellant mass and power. Effcient use of the propellant is obtained by expelling the rocket exhaust at the highest possible velocity, which can be orders of magnitude higher than those achieved in chemical rockets. Handling the hot plasma exhaust requires the use of magnetic nozzles, and the basic physics of ion detachment from the magnetic #12;eld is discussed. The plasma can be generated by RF using helicon waves to heat electrons. Further direct heating of the ions helps to reduce the line radiation losses, and the magnetic geometry is tailored to allow ion cyclotron resonance heating. RF #12;eld and ion trajectory calculations are presented to give a reasonably self-consistent picture of the ion acceleration process.
Date: September 13, 1999
Creator: Baity, F. W., Jr.; Barber, G. C.; Carter, M. D.; Chang-Diaz, F. R.; Goulding, R. H.; Ilin, A. V. et al.
Partner: UNT Libraries Government Documents Department

Meteorological Simulations of Ozone Episode Case Days during the 1996 Paso del Norte Ozone Study

Description: Meteorological simulations centered around the border cities of El Paso and Ciudad Juarez have been performed during an ozone episode that occurred on Aug. 13,1996 during the 1996 Paso del Norte Ozone Study field campaign. Simulations were petiormed using the HOTMAC mesoscale meteorological model using a 1,2,4, and 8 km horizontal grid size nested mesh system. Investigation of the vertical structure and evolution of the atmospheric boundary layer for the Aug. 11-13 time period is emphasized in this paper. Comparison of model-produced wind speed profiles to rawirisonde and radar profiler measurements shows reasonable agreement. A persistent upper-level jet was captured in the model simulations through data assimilation. In the evening hours, the model was not able to produce the strong wind direction shear seen in the radar wind profiles. Based on virtual potential temperature profile comparisons, the model appears to correctly simulate the daytime growth of the convective mixed layer. However, the model underestimates the cooling of the surface layer at night. We found that the upper-level jet significantly impacted the turbulence structure of the boundary layer, leading to relatively high turbulent kinetic energy (tke) values aloft at night. The model indicates that these high tke values aloft enhance the mid-morning growth of the boundary layer. No upper-level turbulence measurements were available to verify this finding, however. Radar profiler-derived mixing heights do indicate relatively rapid morning growth of the mixed layer.
Date: February 1, 1999
Creator: Brown, M. J.; Costigan, K.; Muller, C. & Wang, G.
Partner: UNT Libraries Government Documents Department

Longitudinal impedance measurement of an RK-TBA induction accelerating gap

Description: Induction accelerating gap designs are being studied for Relativistic Klystron Two-Beam Accelerator (RK-TBA) applications. The accelerating gap has to satisfy the following major requirements: hold-off of the applied accelerating voltage pulse, low transverse impedance to limit beam breakup, low longitudinal impedance at the beam-modulation frequency to minimize power loss. Various gap geometries, materials and novel insulating techniques were explored to optimize the gap design. We report on the experimental effort to evaluate the rf properties of the accelerating gaps in a simple pillbox cavity structure. The experimental cavity setup was designed using the AMOS, MAFIA and URMEL numerical codes. Longitudinal impedance measurements above beam-tube cut-off frequency using a single-wire measuring system are presented.
Date: May 1, 1997
Creator: Eylon, S.; Henestroza, E.; Kim, J.-S.; Houck, T.L.; Westenskow, G.A. & Yu, S.S.
Partner: UNT Libraries Government Documents Department