SEARCH · Engineering Papers
Results for “dynamic models”
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.
Recent Advances to Estimation of Fixed-Interface Modal Models Using Dynamic Substructuring
In 2010, Allen & Mayes proposed to estimate the fixed-interface modes of a structure by measuring the modes of the structure bolted to a fixture and then applying constraints to the fixture using the transmission simulator method. While the method proved useful, and has indeed been used in studies since that point, a few peculiarities were noted. First, in some cases the estimated fixed-base natural frequencies were observed to converge very slowly to the true values (in simulated experiments) as the number of constraints was increased. To formulate these constraints, prior studies used only the free-interface modes of the fixture or the measured modes of the assembly. This work extends that to consider other sets of constraints, showing improved results. Furthermore, in some prior studies it has been observed that there were errors of more than 10% in the natural frequencies even when the fixture motion was hundreds of times smaller than the motion of the structure of interest (and so it had presumably been removed). This work explores this phenomenon, seeking to use the strain energy in the fixture, to the extent that it can be estimated using a test-analysis model for the fixture, as a metric to predict frequency error. The proposed methods are explored by applying them to simulated measurements from a beam and from the NASA Space Launch System coupled to the Mobile Launcher.
High Reynolds number calculations using the dynamic subgrid-scale stress model
The dynamic subgrid-scale eddy viscosity model was used in the large-eddy simulation of the turbulent flow in a plane channel for Reynolds numbers based on friction velocity and channel half-width ranging between 200 and 2000, a range including values significantly higher than in previous simulations. The computed wall stress, mean velocity, and Reynolds stress profiles compare very well with experimental and direct simulation data. Comparison of higher moments is also satisfactory. Although the grid in the near-wall region is fairly coarse, the results are quite accurate: the turbulent kinetic energy peaks at y(+) of about 12, and the near-wall behavior of the resolved stresses is captured accurately. The model coefficient is o(0.001) in the buffer layer and beyond, where the cutoff wave numbers are in the decaying region of the spectra; in the near-wall region, the cutoff wave numbers are nearer the energy-containing range, and the resolved turbulent stresses become a constant fraction of the resolved stresses. This feature is responsible for the correct near-wall behavior of the model coefficient.
Oxide driven strength evolution of silicon surfaces
Previous experiments have shown a link between oxidation and strength changes in single crystal silicon nanostructures but provided no clues as to the mechanisms leading to this relationship. Using atomic force microscope-based fracture strength experiments, molecular dynamics modeling, and measurement of oxide development with angle resolved x-ray spectroscopy we study the evolution of strength of silicon (111) surfaces as they oxidize and with fully developed oxide layers. We find that strength drops with partial oxidation but recovers when a fully developed oxide is formed and that surfaces intentionally oxidized from the start maintain their high initial strengths. MD simulations show that strength decreases with the height of atomic layer steps on the surface. These results are corroborated by a completely separate line of testing using micro-scale, polysilicon devices, and the slack chain method in which strength recovers over a long period of exposure to the atmosphere. Combining our results with insights from prior experiments we conclude that previously described strength decrease is a result of oxidation induced roughening of an initially flat silicon (1 1 1) surface and that this effect is transient, a result consistent with the observation that surfaces flatten upon full oxidation.
Searching for dark matter annihilation in the Smith high-velocity cloud
Recent observations suggest that some high-velocity clouds may be confined by massive dark matter halos. In particular, the proximity and proposed dark matter content of the Smith Cloud make it a tempting target for the indirect detection of dark matter annihilation. We argue that the Smith Cloud may be a better target than some Milky Way dwarf spheroidal satellite galaxies and use γ-ray observations from the Fermi Large Area Telescope to search for a dark matter annihilation signal. No significant γ-ray excess is found coincident with the Smith Cloud, and we set strong limits on the dark matter annihilation cross section assuming a spatially extended dark matter profile consistent with dynamical modeling of the Smith Cloud. Notably, these limits exclude the canonical thermal relic cross section (∼ 3 × 10{sup –26} cm{sup 3} s{sup –1}) for dark matter masses ≲ 30 GeV annihilating via the b b-bar or τ{sup +}τ{sup –} channels for certain assumptions of the dark matter density profile; however, uncertainties in the dark matter content of the Smith Cloud may significantly weaken these constraints.
Strengthening of the Coordination Shell by Counter Ions in Aqueous Th 4+ Solutions
The presence of counter ions in solutions containing highly charged metal cations can trigger processes such as ion-pair formation, hydrogen bond breakages and subsequent reformation, and ligand exchanges. In this work, it is shown how halide (Cl-, Br-) and perchlorate (ClO4-) anions affect the strength of the primary solvent coordination shells around Th4+ using explicit solvent and finite temperature ab initio molecular dynamics modeling methods. The 9-fold solvent geometry was found to be the most stable hydration structure in each aqueous solution. Relative to the dilute aqueous solution, the presence of the counter ions did not significantly alter the geometry of the primary hydration shell. However, the free energy analyses indicated that the 10-fold hydrated states were thermodynamically accessible in dilute and bromide aqueous solutions within 1 kcal/mol. Analysis of the results showed that the hydrogen bond lifetimes were longer and solvent exchange energy barriers were larger in solutions with counter ions in comparison with the solution with no counter ions. This implies that the presence of the counter ions induces a strengthening of the Th4+ hydration shell.
Conducting thermal energy to the membrane/water interface for the enhanced desalination of hypersaline brines using membrane distillation
Membrane distillation (MD) is a membrane-based thermal desalination process capable of treating hypersaline brines. Standard MD systems rely on preheating the feed to drive the desalination process. However, relying on the feed to carry thermal energy is limited by a decline of the thermal driving force as the water moves across the membrane, and temperature polarization. In contrast, supplying heat directly into the feed channel, either through the membrane or other channel surfaces, has the potential of minimizing temperature polarization, increasing single-pass water recoveries, and decreasing the number of heat exchangers in the system. When solar thermal energy can be utilized, particularly if the solar heat is optimally delivered to enhance water evaporation and process performance, MD processes can potentially be improved in terms of energy efficiency, environmental sustainability, or operating costs. Here we describe an MD process using layered composite membranes that include a high-thermal-conductivity layer for supplying heat directly to the membrane-water interface and the flow channel. Here, the MD system showed stable performance with water flux up to 9 L/m 2 /hr, and salt rejection >99.9% over hours of desalinating hypersaline feed (100 g/L NaCl). In addition to bench-scale system, we developed a computational fluid dynamics model that successfully described the transport phenomena in the system.
VLTI-GRAVITY Measurements of Cool Evolved Stars I. Variable Photosphere and Extended Atmosphere of the Mira Star R Peg
Context. Dynamic model atmospheres of Mira stars predict variabilities in the photospheric radius and in atmospheric molecular layers which are not yet strongly constrained by observations. Aim. Here we measure the variability of the oxygen-rich Mira star R Peg in near-continuum and molecular bands. Methods. We used near-infrared K-band spectro-interferometry with a spectral resolution of about 4000 obtained at four epochs between post-maximum and minimum visual phases employing the newly available GRAVITY beam combiner at the Very Large Telescope Interferometer (VLTI). Results. Our observations show a continuum radius that is anti-correlated with the visual lightcurve. Uniform disc (UD) angular diameters at a near-continuum wavelength of 2.24 m are steadily increasing with values of 8.7 0.1 mas, 9.4 0.1 mas, 9.8 0.1 mas, and 9.9 0.1 mas at visual phases of 0.15, 0.36, 0.45, 0.53, respectively. UD diameters at a bandpass around 2.05 m, dominated by water vapour, follow the near-continuum variability at larger UD diameters between 10.7 mas and 11.7 mas. UD diameters at the CO 2-0 bandhead, instead, are correlated with the visual lightcurve and anti-correlated with the near-continuum UD diameters, with values between 12.3 mas and 11.7 mas. Conclusions. The observed anti-correlation between continuum radius and visual lightcurve is consistent with an earlier study of the oxygen-rich Mira S Lac, and with recent 1D CODEX dynamic model atmosphere predictions. The amplitude of the variation is comparable to the earlier observations of S Lac, and smaller than predicted by CODEX models. The wavelength-dependent visibility variations at our epochs can be reproduced by a set of CODEX models at model phases between 0.3 and 0.6. The anti-correlation of water vapour and CO contributions at our epochs suggests that these molecules undergo different processes in the extended atmosphere along the stellar cycle. The newly available GRAVITY instrument is suited to conducting longer time series observations, which are needed to provide strong constraints on the model-predicted intra- and inter-cycle variability.
An Overview of the Waste-to-Energy System Simulation (WESyS) Model
The leveraging of waste streams for energy and chemical production could add revenue to waste disposal operations, and it presents opportunities for addressing a variety of economic and environmental objectives at the local, state, and national levels. The Waste-to-Energy System Simulation (WESyS) model is a system dynamics model that was created to simulate the development of the U.S. waste-to-energy industry over time. For each of the three primary waste resources modeled (landfills, concentrated animal feeding operations, and publicly owned treatment works), WESyS simulates technically feasible scenarios for use of the waste, including direct conversion to fuels, and anaerobic digestion followed by flaring, electricity generation, combined heat and power, cleanup and compression to compressed natural gas, and cleanup and injection into an existing pipeline. The model allows users to explore numerous plausible future scenarios for the development of the U.S. waste-to-energy industry. This report provides an overview of the WESyS model and documents the key assumptions, equations, and data sources used to create the model.
An elastography method based on the scanning contact resonance of a piezoelectric cantilever
Purpose: Most tissues may become significantly stiffer than their normal states when there are lesions inside. The tissue's modulus can then act as an identification parameter for clinic diagnosis of tumors or fibrosis, which leads to elastography. This study introduces a novel elastography method that can be used for modulus imaging of superficial organs. Methods: This method is based on the scanning contact-resonance of a unimorph piezoelectric cantilever. The cantilever vibrates in its bending mode with the tip pressed tightly on the sample. The contact resonance frequency of the cantilever-sample system is tracked at each scanning point, from which the sample's modulus can be derived based on a beam dynamic model and a contact mechanics model. Scanning is performed by a three-dimensional motorized stage and the whole system is controlled by a homemade software program based on LabVIEW. Results: Testing onin vitro beef tissues indicates that the fat and the muscle can be easily distinguished using this system, and the accuracy of the modulus measurement can be comparable with that of nanoindentation. Imaging on homemade gelatin phantoms shows that the depth information of the abnormalities can be qualitatively obtained by varying the pressing force. The detection limit of this elastography method is specially examined both experimentally and numerically. Results show that it can detect the typical lesions in superficial organs with the depth of several centimeters. The lateral resolution of this elastography method/system is better than 0.5 mm, and could be further enhanced by using more scanning points. Conclusions: The proposed elastography system can be regarded as a sensitive palpation robot, which may be very promising in early diagnosis of tumors in superficial organs such as breast and thyroid.
Dynamic neutron scattering from conformational dynamics. I. Theory and Markov models
The dynamics of complex molecules can be directly probed by inelastic neutron scattering experiments. However, many of the underlying dynamical processes may exist on similar timescales, which makes it difficult to assign processes seen experimentally to specific structural rearrangements. Here, we show how Markov models can be used to connect structural changes observed in molecular dynamics simulation directly to the relaxation processes probed by scattering experiments. For this, a conformational dynamics theory of dynamical neutron and X-ray scattering is developed, following our previous approach for computing dynamical fingerprints of time-correlation functions [F. No , S. Doose, I. Daidone, M. L llmann, J. Chodera, M. Sauer, and J. Smith, Proc. Natl. Acad. Sci. U.S.A.108, 4822 (2011)]. Markov modeling is used to approximate the relaxation processes and timescales of the molecule via the eigenvectors and eigenvalues of a transition matrix between conformational substates. Furthermore, this procedure allows the establishment of a complete set of exponential decay functions and a full decomposition into the individual contributions, i.e., the contribution of every atom and dynamical process to each experimental relaxation process.
Trade Studies of Space Launch Architectures using Modular Probabilistic Risk Analysis
A top-down risk assessment in the early phases of space exploration architecture development can provide understanding and intuition of the potential risks associated with new designs and technologies. In this approach, risk analysts draw from their past experience and the heritage of similar existing systems as a source for reliability data. This top-down approach captures the complex interactions of the risk driving parts of the integrated system without requiring detailed knowledge of the parts themselves, which is often unavailable in the early design stages. Traditional probabilistic risk analysis (PRA) technologies, however, suffer several drawbacks that limit their timely application to complex technology development programs. The most restrictive of these is a dependence on static planning scenarios, expressed through fault and event trees. Fault trees incorporating comprehensive mission scenarios are routinely constructed for complex space systems, and several commercial software products are available for evaluating fault statistics. These static representations cannot capture the dynamic behavior of system failures without substantial modification of the initial tree. Consequently, the development of dynamic models using fault tree analysis has been an active area of research in recent years. This paper discusses the implementation and demonstration of dynamic, modular scenario modeling for integration of subsystem fault evaluation modules using the Space Architecture Failure Evaluation (SAFE) tool. SAFE is a C++ code that was originally developed to support NASA s Space Launch Initiative. It provides a flexible framework for system architecture definition and trade studies. SAFE supports extensible modeling of dynamic, time-dependent risk drivers of the system and functions at the level of fidelity for which design and failure data exists. The approach is scalable, allowing inclusion of additional information as detailed data becomes available. The tool performs a Monte Carlo analysis to provide statistical estimates. Example results of an architecture system reliability study are summarized for an exploration system concept using heritage data from liquid-fueled expendable Saturn V/Apollo launch vehicles.
Coupling Thermal-Hydraulics with Reactor Pressure Vessel Fracture Models
Because of the importance of maintaining the integrity of the reactor pressure vessel (RPV) in light-water reactors, significant effort has been devoted over the last several decades to understand the potential failure mechanisms, characterize how material properties evolve due to environmental conditions, and develop tools to assess the probability of fracture in aged RPVs during potential transient events. One important aspect of this problem that has not yet been fully addressed is the effects of nonuniform cooling on the RPV inner wall due to the way low-temperature water injected into the primary coolant system is distributed within the RPV during potential loss-of-coolant accident transients. Colder temperatures are expected in the regions near the inlets, causing a “cold-plume” effect. Rapid cooling of the RPV wall increases the likelihood of fracture initiation, and localized cooling due to cold-plume effects could potentially exacerbate this problem. The Grizzly code can perform a three-dimensional simulation of the RPV thermomechanical response as well as the probabilistic assessment of fracture for a population of pre-existing flaws, which is necessary to account for spatial variations in coolant temperature. This report documents a first effort to couple Grizzly RPV fracture models with computational fluid dynamics models performed using the Cardinal code, which allows for a realistic representation of the temperature distribution. This report demonstrates this coupling approach on two test cases: a simplified proof-of-concept scenario and a realistic small-break loss-of-coolant accident scenario. In both cases, the effects of the spatially varying temperature are evident; although, because the small-break loss-of-coolant accident scenario was far from challenging the RPV’s margins, it is still inconclusive how important cold-plume effects are in realistic accident scenarios. However, this effort represents an important step toward performing realistic analyses to determine to what extent spatially nonuniform cooling might affect RPV integrity under more aggressive accident scenarios.
Bond Length Alternation and Internal Dynamics in Model Aromatic Substituents of Lignin
In this report broadband microwave spectra were recorded over the 2-18 GHz frequency range for a series of four model aromatic components of lignin; namely, guaiacol (ortho-methoxy phenol, G), syringol (2,6-dimethoxy phenol, S), 4-methyl guaiacol (MG), and 4-vinyl guaiacol (VG), under jet-cooled conditions in the gas phase. Using a combination of 13 C isotopic data and electronic structure calculations, distortions of the phenyl ring by the substituents on the ring are identified. In all four molecules, the r C(1)-C(6) bond between the two substituted C-atoms lengthens, leading to clear bond alternation that reflects an increase in the phenyl ring resonance structure with double bonds at r C(1)-C(2) , r C(3)-C(4) and r C(5)-C(6) . Syringol, with its symmetric methoxy substituents, possesses a microwave spectrum with tunneling doublets in the a-type transitions associated with H-atom tunneling. These splittings were fit to determine a barrier to hindered rotation of the OH group of 1975 cm -1 , a value nearly 50% greater than that in phenol, due to the presence of the intramolecular OH…OCH 3 H-bonds at the two equivalent planar geometries. In 4-methyl guaiacol, methyl rotor splittings are observed and used to confirm and refine an earlier measurement of the three-fold barrier V 3 = 67 cm -1 . Finally, 4-vinyl guaiacol shows transitions due to two conformers differing in the relative orientations of the vinyl and OH groups.
Assessing dynamic vegetation model parameter uncertainty across Alaskan arctic tundra plant communities
Abstract As the Arctic region moves into uncharted territory under a warming climate, it is important to refine the terrestrial biosphere models (TBMs) that help us understand and predict change. One fundamental uncertainty in TBMs relates to model parameters, configuration variables internal to the model whose value can be estimated from data. We incorporate a version of the Terrestrial Ecosystem Model (TEM) developed for arctic ecosystems into the Predictive Ecosystem Analyzer (PEcAn) framework. PEcAn treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We examined how variation in 21 parameters in the equation for gross primary production influenced model sensitivity and uncertainty in terms of two carbon fluxes (net primary productivity and heterotrophic respiration) and two carbon (C) pools (vegetation C and soil C). We set up different parameterizations of TEM across a range of tundra types (tussock tundra, heath tundra, wet sedge tundra, and shrub tundra) in northern Alaska, along a latitudinal transect extending from the coastal plain near Utqiaġvik to the southern foothills of the Brooks Range, to the Seward Peninsula. TEM was most sensitive to parameters related to the temperature regulation of photosynthesis. Model uncertainty was mostly due to parameters related to leaf area, temperature regulation of photosynthesis, and the stomatal responses to ambient light conditions. Our analysis also showed that sensitivity and uncertainty to a given parameter varied spatially. At some sites, model sensitivity and uncertainty tended to be connected to a wider range of parameters, underlining the importance of assessing tundra community processes across environmental gradients or geographic locations. Generally, across sites, the flux of net primary productivity (NPP) and pool of vegetation C had about equal uncertainty, while heterotrophic respiration had higher uncertainty than the pool of soil C. Our study illustrates the complexity inherent in evaluating parameter uncertainty across highly heterogeneous arctic tundra plant communities. It also provides a framework for iteratively testing how newly collected field data related to key parameters may result in more effective forecasting of Arctic change.
Smoky Mountain Data Challenge 2020: An Open Call to Solve Data Problems in the Areas of Neutron Science, Material Science, Urban Modeling and Dynamics, Geophysics, and Biomedical Informatics
The 2020 Smoky Mountains Computational Sciences and Engineering Conference enlists research scientists from across Oak Ridge National Laboratory (ORNL) to be data sponsors and help create data analytics challenges for eminent data sets at the laboratory. This work describes the significance of each of the seven data sets and their as- sociated challenge questions. The challenge questions for each data set were required to cover multiple difficulty levels. An international call for participation was sent to students, and researchers asking them to form teams of up to four people to apply novel data analytics techniques to these data sets.
Study on the Fluid Flow Characteristics of Coherent Jets with CO{sub 2} and O{sub 2} Mixed Injection in Electric Arc Furnace Steelmaking Processes
As an efficient oxygen supplying technology, coherent jets are widely applied in electric arc furnace (EAF) steelmaking processes to strengthen chemical energy input, speed up smelting rhythm, and promote the uniformity of molten bath temperature and compositions. Recently, the coherent jet with CO{sub 2} and O{sub 2} mixed injection (COMI) was proposed and demonstrated great application potentiality in reducing the dust production in EAF steelmaking. In the present study, based on the eddy dissipation concept model, a computational fluid dynamics model of coherent jets with COMI was built with the overall and detailed chemical kinetic mechanisms (GRI-Mech 3.0). Compared with one-step combustion reaction, GRI-Mech 3.0 consists of 325 elementary reactions with 53 components and can predict more accurate results. The numerical simulation results were validated by the combustion experiment data. The jet behavior and the fluid flow characteristics of coherent jets with COMI under 298 K and 1700 K (25 °C and 1427 °C) were studied and the results showed that for coherent jets with COMI, the chemical effect of CO{sub 2} significantly weakened the shrouding combustion reactions of CH{sub 4} and the relative importance of the chemical effect of CO{sub 2} increases with CO{sub 2} concentration increasing. The potential core length of coherent jet decreases with the volume fraction of CO{sub 2} increasing. Moreover, it also can be found that the potential core length of coherent jets was prolonged with higher ambient temperature.