Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “COMPUTER SIMULATION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

Phonon-magnon interactions in body centered cubic iron: A combined molecular and spin dynamics study

Combining an atomistic many-body potential with a classical spin Hamiltonian parameterized by first principles calculations, molecular-spin dynamics computer simulations were performed to investigate phonon-magnon interactions in body centered cubic iron. Results obtained for spin-spin and density-density dynamic structure factors show that noticeable softening and damping of magnon modes occur due to the presence of lattice vibrations. Furthermore, as a result of the phonon-magnon coupling, additional longitudinal spin wave excitations are observed, with the same frequencies as the longitudinal phonon modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Rare Earth and Transition Metal Containing Glasses

Transition metal (TM) and rare earth (RE) ions have been incorporated into many glass systems such as silicate, phosphate, and borosilicate-based oxide glasses, as well as in halide and chalcogenide glasses, that find applications ranging from optical, photonic, and magnetic devices, solid-state battery, to nuclear waste disposal. Understanding the structural role of RE and TM in these glasses can help to develop glass compositions for targeted applications with either high-optical emission efficiency, electrical conductivity, or chemical durability. In this chapter, we first provide a general introduction of the applications and structural features of RE and TM in glasses, then the critical aspects of molecular dynamics (MD) simulations of these glasses such as interatomic potentials, structural analysis tools to study RE and TM ions in glasses and their clustering behaviors, Quantitative Structure–Property Analysis (QSPR), diffusion and dynamic property calculations, and electronic structure calculations to understand electronic defects such as charge trapping and radiation effects are introduced. Three representative case studies are presented: the first one is on MD simulations of erbium- and europium-doped silica and silicate glasses, as well as cerium doped aluminophosphate glasses, that revealed the effect of glass composition on RE ion local structure and clustering behavior. Electronic structure calculations of cerium-doped glass show how the existence of multioxidation states help to mediate radiation-induced damages caused by excited electron–hole pairs was also discussed. The second one focuses on alkali vanadophosphate glasses where the existence of two vanadium oxidation states help to provide electronic conduction in the glasses while alkali ions provide ionic conduction. MD simulations were used to understand vanadium environments and other structural aspects in the phosphate glasses, as well as the ionic transport behaviors of alkali ions. The third case study is on zirconium-containing borosilicate and aluminosilicate glasses that find wide applications in nuclear waste disposal. MD simulations help to provide structural details of zirconium ions that are validated by diffraction and EXAFS spectra. The structural information was used to interpret changes of mechanical properties and chemical durability by using QSPR and other analyses-based MD-generated structure models.

Du, Jincheng↗

Deep subsurface drip irrigation using coal-bed sodic water: Part I. Water and solute movement

Water co-produced with coal-bed methane (CBM) in the semi-arid Powder River Basin of Wyoming and Montana commonly has relatively low salinity and high sodium adsorption ratios that can degrade soil permeability where used for irrigation. Nevertheless, a desire to derive beneficial use from the water and a need to dispose of large volumes of it have motivated the design of a deep subsurface drip irrigation (SDI) system capable of utilizing that water. Drip tubing is buried 92 cm deep and irrigates at a relatively constant rate year-round, while evapotranspiration by the alfalfa and grass crops grown is seasonal. We use field data from two sites and computer simulations of unsaturated flow to understand water and solute movements in the SDI fields. Combined irrigation and precipitation exceed potential evapotranspiration by 300–480 mm annually. Initially, excess water contributes to increased storage in the unsaturated zone, and then drainage causes cyclical rises in the water table beneath the fields. Native chloride and nitrate below 200 cm depth are leached by the drainage. Some CBM water moves upward from the drip tubing, drawn by drier conditions above. Chloride from CBM water accumulates there as root uptake removes the water. Year over year accumulations indicated by computer simulations illustrate that infiltration of precipitation water from the surface only partially leaches such accumulations away. Field data show that 7% and 27% of added chloride has accumulated above the drip tubing in an alfalfa and grass field, respectively, following 6 years of irrigation. Maximum chloride concentrations in the alfalfa field are around 45 cm depth but reach the surface in parts of the grass field, illustrating differences driven by crop physiology. Deep SDI offers a means of utilizing marginal quality irrigation waters and managing the accumulation of their associated solutes in the crop rooting zone.

01 COAL, LIGNITE, AND PEAT↗

Simulation of computed radiography with imaging plate detectors

Computed radiography (CR) using phosphor imaging plate detectors is taking an increasing place in Radiography Testing. CR uses similar equipment as conventional radiography except that the classical X-ray film is replaced by a numerical detector, called image plate (IP), which is made of a photostimulable layer and which is read by a scanning device through photostimulated luminescence. Such digital radiography has already demonstrated important benefits in terms of exposure time, decrease of source energies and thus reduction of radioprotection area besides being a solution without effluents. This paper presents a model for the simulation of radiography with image plate detectors in CIVA together with examples of validation of the model. The study consists in a cross comparison between experimental and simulation results obtained on a step wedge with a classical X-ray tube. Results are proposed in particular with wire Image quality Indicator (IQI) and duplex IQI.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

How Does Vertical Wind Shear Influence Updraft Characteristics and Hydrometeor Distributions in Supercell Thunderstorms?

Abstract Vertical wind shear is known to affect supercell thunderstorms by displacing updraft hydrometeor mass downshear, thereby facilitating the storms’ longevity. Shear also impacts the size of supercell updrafts, with stronger shear leading to wider, less dilute, and stronger updrafts with likely greater hydrometeor production. To more clearly define the role of shear across different vertical layers on hydrometeor concentrations and displacements relative to supercell updrafts, a suite of idealized numerical model simulations of supercells was conducted. Shear magnitudes were systematically varied across the 0–1, 1–6, and 6–12 km AGL layers, while the thermodynamic environment was held fixed. Simulations show that as shear magnitude increases, especially from 1 to 6 km, updrafts become wider and less dilute with an increase in hydrometeor loading, along with an increase in the low-level precipitation area/rate and total precipitation accumulation. Even with greater updraft hydrometeor loading amid stronger shear, updrafts are more intense in stronger shear simulations due to larger thermal buoyancy owing to wider, less dilute updraft cores. Furthermore, downshear hydrometeor displacements are larger in environments with stronger 1–6-km shear. In contrast, there is relatively less sensitivity of hydrometeor concentrations and displacements to variations in either 0–1- or 6–12-km shear. Results are consistent across free tropospheric relative humidity sensitivity simulations, which show an increase in updraft size and hydrometeor mass with increasing free tropospheric relative humidity owing to a reduction in entrainment-driven dilution for wider updrafts in moister environments. Significance Statement Rotating thunderstorms, known as supercells, are able to persist for multiple hours. One common explanation is that large changes in wind speed and/or direction with height, or shear, transport rain/hail away from supercell updrafts, supporting their maintenance. The strong shear within supercell environments, however, may also lead to greater rail/hail amounts, thereby leading to weaker storms due to this extra mass of water/ice within updrafts. Furthermore, the impact of shear across different height layers on supercell rain/hail characteristics has not been thoroughly investigated. In this study, computer simulations of supercells were conducted to determine that shear occurring between 1 and 6 km above ground level has a large impact on rain/hail distribution in supercells and that stronger shear in this layer leads to wider/stronger supercells with greater rain/hail accumulations at the surface. Additionally, some of the extra mass of water/ice is transported farther away from updrafts due to the stronger environmental storm-relative winds.

Meteorology & Atmospheric Sciences↗

Collaborative Research: Tomographic imaging of laser-plasma structures

The interaction of intense short laser pulses with ionized gases, or plasmas, underlies many applications such as acceleration of elementary particles, production of energy by laser fusion, generation of x-ray and far-infrared “terahertz” pulses for medical and materials probing, remote sensing of explosives and pollutants, and generation of guide stars. Such laser-plasma interactions create tiny electron density structures (analogous to the wake behind a boat) inside the plasma in the shape of waves, bubbles and filaments that move at the speed of light, and evolve as they propagate. Prior to recent work by the PI of this proposal, detailed knowledge of such structures came exclusively from intensive computer simulations. Now “snapshots” of these elusive, light-velocity structures can be taken in the laboratory using dynamic variant of holography, the technique used to produce ID cards and DVDs, and dynamic variant of tomography, the technique used in medicine to image internal bodily organs. These fast visualization techniques are important for understanding, improving and scaling the above-mentioned applications of laser-plasma interactions. In this project, we accomplished three things: 1) We took holographic pictures of a laser-driven plasma-wave in the act of accelerating electrons to high energy, and used computer simulations to understand the pictures. 2) Using results from this experiment to optimize the performance of the accelerator, and the brightness of x-rays that it emits. These x-rays will be useful for medical and materials science applications. 3) We made technical improvements to the holographic technique that enables us to see finer details in the recorded pictures. Four refereed journal papers were published, and two students earned PhDs and moved on to scientific careers in US National Laboratories based on their work under this project.

43 PARTICLE ACCELERATORS↗

Study of a macrodefect in a silicon carbid single crystal by means of X-ray phase contrast

The morphology of a macrodefect in a single-crystal silicon carbide wafer has been investigated by the computer simulation of an experimental X-ray phase-contrast image. A micropipe, i.e., a long cavity with a small (elliptical in the general case) cross section, in a single crystal has been considered as a macrodefect. A far-field image of micropipe has been measured with the aid of synchrotron radiation without a monochromator. The parameters of micropipe elliptical cross section are determined based on one projection in two directions: parallel and perpendicular to the X-ray beam propagation direction, when scanning along the pipe axis. The results demonstrate the efficiency of the phase contrast method supplemented with computer simulation for studying such macrodefects when the defect position in the sample volume is unknown beforehand.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

An Experimental and Computational Study of a Single Diesel Droplet Impinging on a Dry Surface

This work presents experiments and computational simulations of a single diesel droplet impinging on a dry surface to provide a better understanding of the dynamics of the impingement process. The experimental work was carried out at several substrate surface temperatures where a single diesel droplet was dropped from a precision syringe pump to provide a specific droplet diameter and impact velocity. The effect of the surface temperature on the temporal evolution of the droplet spreading factor, dynamic contact angle, and contact line velocity was studied. The experimental study on the temporal evolution of the droplet-wall interaction was used to develop a dynamic contact angle model for computational simulations. Using an in-house multiphase flow solver, simulations of the single diesel droplet impingement were performed. The simulation results of the spreading factor were compared to the experimental results and the effects of liquid temperature and contact angle model parameters were investigated. The development of a dynamic contact angle model for diesel droplets enhances the ability to study diesel droplet impingement, an important factor in diesel fuel injection. With an accurate dynamic contact angle model, future studies can be performed at scales and velocities representative of the conditions present in internal combustion engines.

Markt, D.↗

Speckle contrast diffuse correlation tomography of complex turbid medium flow

Purpose: Developed herein is a three-dimensional (3D) flow contrast imaging system leveraging advancements in the extension of laser speckle contrast imaging theories to deep tissues along with our recently developed finite-element diffuse correlation tomography (DCT) reconstruction scheme. This technique, termed speckle contrast diffuse correlation tomography (scDCT), enables incorporation of complex optical property heterogeneities and sample boundaries. When combined with a reflectance-based design, this system facilitates a rapid segue into flow contrast imaging of larger, in vivo applications such as humans. Methods: A highly sensitive CCD camera was integrated into a reflectance-based optical system. Four long-coherence laser source positions were coupled to an optical switch for sequencing of tomographic data acquisition providing multiple projections through the sample. This system was investigated through incorporation of liquid and solid tissue-like phantoms exhibiting optical properties and flow characteristics typical of human tissues. Computer simulations were also performed for comparisons. A uniquely encountered smear correction algorithm was employed to correct point-source illumination contributions during image capture with the frame-transfer CCD and reflectance setup. Results: Measurements with scDCT on a homogeneous liquid phantom showed that speckle contrast-based deep flow indices were within 12% of those from standard DCT. Inclusion of a solid phantom submerged below the liquid phantom surface allowed for heterogeneity detection and validation. The heterogeneity was identified successfully by reconstructed 3D flow contrast tomography with scDCT. The heterogeneity center and dimensions and averaged relative flow (within 3%) and localization were in agreement with actuality and computer simulations, respectively. Conclusions: A custom cost-effective CCD-based reflectance 3D flow imaging system demonstrated rapid acquisition of dense boundary data and, with further studies, a high potential for translatability to real tissues with arbitrary boundaries. A requisite correction was also found for measurements in the fashion of scDCT to recover accurate speckle contrast of deep tissues.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Roll-up of validation results to a target application.

Suites of experiments are preformed over a validation hierarchy to test computational simulation models for complex applications. Experiments within the hierarchy can be performed at different conditions and configurations than those for an intended application, with each experiment testing only part of the physics relevant for the application. The purpose of the present work is to develop methodology to roll-up validation results to an application, and to assess the impact the validation hierarchy design has on the roll-up results. The roll-up is accomplished through the development of a meta-model that relates validation measurements throughout a hierarchy to the desired response quantities for the target application. The meta-model is developed using the computation simulation models for the experiments and the application. The meta-model approach is applied to a series of example transport problems that represent complete and incomplete coverage of the physics of the target application by the validation experiments.

Hills, Richard Guy↗

Decomposition of toluene in a steady-state atmospheric-pressure glow discharge

Results are presented from experimental studies of decomposition of toluene (C{sub 6}H{sub 5}CH{sub 3}) in a polluted air flow by means of a steady-state atmospheric pressure glow discharge at different water vapor contents in the working gas. The experimental results on the degree of C{sub 6}H{sub 5}CH{sub 3} removal are compared with the results of computer simulations conducted in the framework of the developed kinetic model of plasma chemical decomposition of toluene in the N{sub 2}: O{sub 2}: H{sub 2}O gas mixture. A substantial influence of the gas flow humidity on toluene decomposition in the atmospheric pressure glow discharge is demonstrated. The main mechanisms of the influence of humidity on C{sub 6}H{sub 5}CH{sub 3} decomposition are determined. The existence of two stages in the process of toluene removal, which differ in their duration and the intensity of plasma chemical decomposition of C{sub 6}H{sub 5}CH{sub 3} is established. Based on the results of computer simulations, the composition of the products of plasma chemical reactions at the output of the reactor is analyzed as a function of the specific energy deposition and gas flow humidity. The existence of a catalytic cycle in which hydroxyl radical OH acts a catalyst and which substantially accelerates the recombination of oxygen atoms and suppression of ozone generation when the plasma-forming gas contains water vapor is established.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling of polymer photodegradation for solar cell modules

It was shown that many of the experimental observations in the photooxidation of hydrocarbon polymers can be accounted for with a computer simulation using an elementary mechanistic model with corresponding rate constants for each reaction. For outdoor applications, however, such as in photovoltaics, the variation of temperature must have important effects on the useful lifetimes of such materials. The data bank necessary to replace the isothermal rate constant values with Arrhenius activation parameters: A (the pre-exponential factor) and E (the activation energy) was searched. The best collection of data assembled to data is summarized. Note, however, that the problem is now considerably enlarged since from a theoretical point of view, with 51 of the input variables replaced with 102 parameters. The sensitivity of the overall scheme is such that even after many computer simulations, a successful photooxidation simulation with the expanded variable set was not completed. Many of the species in the complex process undergo a number of competitive pathways, the relative importance of each being often sensitive to small changes in the calculated rate constant values.

Somersall, A. C.↗

Typical motions in multiple systems

In very old times, people counted - one, two, many. The author wants to show that they were right. Consider the motions of isolated bodies: (1) N = 1 - simple motion; (2) N = 2 - Keplerian orbits; and (3) N = 3 - this is the difficult problem. In general, this problem can be studied only by computer simulations. The author studied this problem over many years (see, e.g., Agekian and Anosova, 1967; Anosova, 1986, 1989 a,b). The principal result is that two basic types of dynamics take place in triple systems. The first special type is the stable hierarchical systems with two almost Keplerian orbits. The second general type is the unstable triple systems with complicated motions of the bodies. By random choice of the initial conditions, by the Monte-Carlo method, the stable systems comprised about approx. 10% of the examined cases; the unstable systems comprised the other approx. 90% of cases under consideration. In N greater than 3, the studies of dynamics of such systems by computer simulations show that we have in general also the motions roughly as at the cases 1 - 3 with the relative negative or positive energies of the bodies. In the author's picture, the typical trajectories of the bodies in unstable triple systems of the general type of dynamics are seen. Such systems are disrupted always after close triple approaches of the bodies. These approaches play a role like the gravitational slingshot. Often, the velocities of escapers are very large. On the other hand, the movie also shows the dynamical processes of a formation, dynamical evolution and disruption of the temporary wide binaries in triples and a formation of final hard massive binaries in the final evolution of triples.

Anosova, Joanna P.↗

Test-Free Fracture Toughness

Computational simulation results can give the prediction of damage growth and progression and fracture toughness of composite structures. The experimental data from literature provide environmental effects on the fracture behavior of metallic or fiber composite structures. However, the traditional experimental methods to analyze the influence of the imposed conditions are expensive and time consuming. This research used the CODSTRAN code to model the temperature effects, scaling effects and the loading effects of fiberbraided composite specimens with and without fiber-optic sensors on the damage initiation and energy release rates. The load-displacement relationship and fracture toughness assessment approach is compared with the test results from literature and it is verified that the computational simulation, with the use of established material modeling and finite element modules, adequately tracks the changes of fracture toughness and subsequent fracture propagation for any fiberbraided composite structure due to the change of fiber orientations, presence of large diameter optical fibers, and any loading conditions.

Minnetyan, Levon↗

Test-Free Fracture Toughness

Computational simulation results can give the prediction of damage growth and progression and fracture toughness of composite structures. The experimental data from literature provide environmental effects on the fracture behavior of metallic or fiber composite structures. However, the traditional experimental methods to analyze the influence of the imposed conditions are expensive and time consuming. This research used the CODSTRAN code to model the temperature effects, scaling effects and the loading effects of fiber/braided composite specimens with and without fiber-optic sensors on the damage initiation and energy release rates. The load-displacement relationship and fracture toughness assessment approach is compared with the test results from literature and it is verified that the computational simulation, with the use of established material modeling and finite element modules, adequately tracks the changes of fracture toughness and subsequent fracture propagation for any fiber/braided composite structure due to the change of fiber orientations, presence of large diameter optical fibers, and any loading conditions.

Minnetyan, Levon↗

Controlling Nanoparticle Orientations in the Self-Assembly of ‘Patchy’ Quantum Dot-Gold Heterostructural Nanocrystals

Self-assembly of nanocrystals is a promising route for creating macroscale materials that derive function from the properties of their nanoscale building blocks. While much progress has been made assembling nanocrystals into different superlattices, controlling the relative orientations of nanocrystals in those lattices remains a challenge. Here, we combine experiments with computer simulations to study the self-assembly of patchy heterostructural nanocrystals (HNCs), consisting of near-spherical quantum dots decorated with regular arrangements of small gold satellites, into close-packed superlattices with pronounced orientational alignment of HNCs. Our simulations indicate that the orientational alignment is caused by van der Waals interactions between gold patches and is sensitive to the interparticle distance in the superlattice. We demonstrate experimentally that the degree and type of orientational alignment can be controlled by changing ligand populations on HNCs. This study provides guidance for the design and fabrication of nanocrystal superlattices with enhanced structural control.

ZHu, Hua↗

The sheath structure around a negatively charged rocket payload

The sheath structure around a rocket payload charged up to 460 V negative relative to the ambient ionospheric plasma is investigated experimentally and by computer simulations. The experimental results come from the Charge 2 sounding rocket experiment in which the payload was split into two separate sections (mother and daughter) connected with a conducting, insulated tether. In one of the experimental modes, the voltage between the payloads was increased linearly from 0 to 460 V in 2.5 s. A floating probe array was mounted on the mother with probes located 25, 50, 75, and 100 cm from the rocket surface. The internal impedance of the array was smaller than the probe/plasma impedance, which influenced the potential measurements. The measurements contain signatures, resulting from the outward expansion of the ion sheath with increasing negative mother potential. This conclusion is substantiated by computer simulations of space charge limited flow.

Neubert, T.↗