Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic profiling”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Influence of the Dielectric Constant on the Ionic Current Rectification of Bipolar Nanopores

In this paper, we investigate how the dielectric constant, ϵ, of an electrolyte solvent influences the current rectification characteristics of bipolar nanopores. It is well recognized that bipolar nanopores with two oppositely charged regions rectify current when exposed to an alternating electric potential difference. Here, we consider dilute electrolytes with NaCl only and with a mixture of NaCl and charged nanoparticles. These systems are studied using two levels of description, all-atom explicit water molecular dynamics (MD) simulations and coarse-grained implicit solvent MD simulations. The charge density and electric potential profiles and current-voltage relationship predicted by the implicit solvent simulations with ϵ = 11.3 show good agreement with the predictions from the explicit water simulations. Under nonequilibrium conditions, the predictions of the implicit solvent simulations with a dielectric constant closer to the one of bulk water are significantly different from the predictions obtained with the explicit water model. Further, these findings are closely aligned with experimental data on the dielectric constant of water when confined to nanometric spaces, which suggests that ϵ decreases significantly compared to its value in the bulk. Moreover, the largest electric current rectification is observed in systems containing nanoparticles when ϵ = 78.8. Using enhanced sampling, we have shown that this larger rectification arises from the presence of a significantly deeper minimum in the free energy of the system with a larger ϵ, and when a negative voltage bias is applied. Since implicit solvent models and mean-field continuum theories are often used to design Janus membranes based on bipolar nanopores, this work highlights the importance of properly accounting for the effects of confinement on the dielectric constant of the electrolyte solvent. The results presented here indicate that the dielectric constant in implicit solvent simulations may be used as an adjustable parameter to approximately account for the effects of nanometric confinement on aqueous electrolyte solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

DETECTION OF EQUATORWARD MERIDIONAL FLOW AND EVIDENCE OF DOUBLE-CELL MERIDIONAL CIRCULATION INSIDE THE SUN

Meridional flow in the solar interior plays an important role in redistributing angular momentum and transporting magnetic flux inside the Sun. Although it has long been recognized that the meridional flow is predominantly poleward at the Sun's surface and in its shallow interior, the location of the equatorward return flow and the meridional flow profile in the deeper interior remain unclear. Using the first 2 yr of continuous helioseismology observations from the Solar Dynamics Observatory/Helioseismic Magnetic Imager, we analyze travel times of acoustic waves that propagate through different depths of the solar interior carrying information about the solar interior dynamics. After removing a systematic center-to-limb effect in the helioseismic measurements and performing inversions for flow speed, we find that the poleward meridional flow of a speed of 15 m s{sup -1} extends in depth from the photosphere to about 0.91 R{sub Sun }. An equatorward flow of a speed of 10 m s{sup -1} is found between 0.82 and 0.91 R{sub Sun} in the middle of the convection zone. Our analysis also shows evidence of that the meridional flow turns poleward again below 0.82 R{sub Sun }, indicating an existence of a second meridional circulation cell below the shallower one. This double-cell meridional circulation profile with an equatorward flow shallower than previously thought suggests a rethinking of how magnetic field is generated and redistributed inside the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Electron trapping in evolving coronal structures during a large gradual hard X-ray/radio burst

Gradual hard X-ray/radio bursts are characterized by their long duration, smooth time profile, time delays between peaks at different hard X-ray energies and microwaves, and radiation from extended sources in the low and middle corona. Their characteristic properties have been ascribed to the dynamic evolution of the accelerated electrons in coronal magnetic traps or to the separate acceleration of high-energy electrons in a 'second step' process. The information available so far was drawn from quality considerations of time profiles or even only from the common occurrence of emissions in different spectral ranges. This paper presents model computations of the temporal evolution of hard X-ray and microwave spectra, together with a qualitative discussion of radio lightcurves over a wide spectral range, and metric imaging observations. The basis hypothesis investigated is that the peculiar 'gradual' features can be related to the dynamical evolution of electrons injected over an extended time interval in a coronal trap, with electrons up to relativistic energies being injected simultaneously. The analyzed event (26 April. 1981) is particularly challenging to this hypothesis because of the long time delays between peaks at different X-ray energies and microwave frequencies. The observations are shown to be consistent with the hypothesis, provided that the electrons lose their energy by Coulomb collisions and possibly betatron deceleration. The access of the electrons to different coronal structures varies in the course of the event. The evolution and likely destabilization of part of the coronal plasma-magnetic field configuration is of crucial influence in determining the access to these structures and possibly the dynamical evolution of the trapped electrons through betatron deceleration in the late phase of the event.

Bruggmann, G.↗
DOE OSTI2022

Mechanism of tethered agonist-mediated signaling by polycystin-1

Mutations of polycystin-1 (PC1) are the major cause (85% of cases) of autosomal dominant polycystic kidney disease (ADPKD), which is the fourth leading cause of kidney failure. PC1 is thought to function as an atypical G protein-coupled receptor, yet the mechanism by which PC1 regulates G-protein signaling remains poorly understood. A significant portion of ADPKD mutations of PC1 encode a protein with defects in maturation or reduced function that may be amenable to functional rescue. In this work, we have combined complementary biochemical and cellular assay experiments and accelerated molecular simulations, which revealed an allosteric transduction pathway in activation of the PC1 C-terminal fragment. Our findings will facilitate future rational drug design efforts targeting the PC1 signaling function. Polycystin-1 (PC1) is an important unusual G protein-coupled receptor (GPCR) with 11 transmembrane domains, and its mutations account for 85% of cases of autosomal dominant polycystic kidney disease (ADPKD). PC1 shares multiple characteristics with Adhesion GPCRs. These include a GPCR proteolysis site that autocatalytically divides these proteins into extracellular, N-terminal, and membrane-embedded, C-terminal fragments (CTF), and a tethered agonist (TA) within the N-terminal stalk of the CTF that is suggested to activate signaling. However, the mechanism by which a TA can activate PC1 is not known. Here, we have combined functional cellular signaling experiments of PC1 CTF expression constructs encoding wild type, stalkless, and three different ADPKD stalk variants with all-atom Gaussian accelerated molecular dynamics (GaMD) simulations to investigate TA-mediated signaling activation. Correlations of residue motions and free-energy profiles calculated from the GaMD simulations correlated with the differential signaling abilities of wild type and stalk variants of PC1 CTF. They suggested an allosteric mechanism involving residue interactions connecting the stalk, Tetragonal Opening for Polycystins (TOP) domain, and putative pore loop in TA-mediated activation of PC1 CTF. Key interacting residues such as N3074–S3585 and R3848–E4078 predicted from the GaMD simulations were validated by mutagenesis experiments. Together, these complementary analyses have provided insights into a TA-mediated activation mechanism of PC1 CTF signaling, which will be important for future rational drug design targeting PC1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

Constructing Synoptic Maps of Stratospheric Column Ozone from HALOE, SAGE and Balloonsonde Data Using Potential Vorticity Isentropic Coordinate Transformations

In this study we utilize potential vorticity - isentropic (PVI) coordinate transformations as a means of combining ozone data from different sources to construct daily, synthetic three-dimensional ozone fields. This methodology has been used successfully to reconstruct ozone maps in particular regions from aircraft data over the period of the aircraft campaign. We expand this method to create high-resolution daily global maps of profile ozone data, particularly in the lower stratosphere, where high-resolution ozone data are sparse. Ozone climatologies in PVI-space are constructed from satellite-based SAGE II and UARS/HALOE data, both of which-use solar occultation techniques to make high vertical resolution ozone profile measurements, but with low spatial resolution. A climatology from ground-based balloonsonde data is also created. The climatologies are used to establish the relationship between ozone and dynamical variability, which is defined by the potential vorticity (in the form of equivalent latitude) and potential temperature fields. Once a PVI climatology has been created from data taken by one or more instruments, high-resolution daily profile ozone field estimates are constructed based solely on the PVI fields, which are available on a daily basis from NCEP analysis. These profile ozone maps could be used for a variety of applications, including use in conjunction with total ozone maps to create a daily tropospheric ozone product, as input to forecast models, or as a tool for validating independent ozone measurements when correlative data are not available. This technique is limited to regions where the ozone is a long-term tracer and the flow is adiabatic. We evaluate the internal consistency of the technique by transforming the ozone back to physical space and comparing to the original profiles. Biases in the long-term average of the differences are used to identify regions where the technique is consistently introducing errors. Initial results show the technique is useful in the lower stratosphere at most latitudes throughout the year,and in the winter hemisphere in the middle stratosphere. The results are problematic in the summer hemisphere middle stratosphere due to increased ozone photochemistry and weak PV gradients. Alternate techniques in these regions will be discussed. An additional limitation is the quality and resolution of the meteorological data.

Hollandsworth, Stacey M.↗
DOE OSTI2020

Dynamic optimization and economic evaluation of flexible heat integration in a hybrid concentrated solar power plant

Hybridization of concentrated solar power (CSP) plants provides flexibility in operation that can drastically improve the solar-to-electric (STE) efficiency and levelized cost of electricity (LCOE) relative to standalone CSP plants. Flexible heat integration (FHI) is a novel concept where the collection and integration of CSP within a power plant is modified relative to the amount of solar energy available. FHI improves the thermal efficiency of a hybrid solar tower steam Rankine cycle power plant but leads to increased pumping needs due to continuously elevated molten salt flow rates through the collection system, which can negatively impact STE efficiency. The present work is carried out to maximize the STE efficiency of a hybrid CSP plant utilizing FHI by employing a dynamic optimization framework where a genetic algorithm optimizes the operation of the plant over a given solar irradiance profile. The study concerns a plant hypothetically located in Salt Lake City, Utah. Here, the optimization results confirm the accuracy of a predictive heuristic where the preferred operation of the plant can be estimated relative to local peaks in the incident power generated by the heliostat collection field. The optimized FHI operation demonstrates a yearly STE efficiency of 13.8%, whereas the equivalent base-level hybrid and solar-only plants exhibit solar efficiencies of 13.4% and 11.2%, respectively. Economic analysis shows that FHI reduces yearly natural gas costs, leading to a $\$0.5$/MWh reduction in LCOE relative to the base-level hybrid configuration. Overall, the results show that hybrid FHI schemes exhibit economic benefits along with observed thermodynamic improvements.

14 SOLAR ENERGY↗
NASA NTRS1993

Unsteady Flow Simulations in Support of the SSME HEX Turning Vane Cracking Investigation with the ATD HPOTP

Unsteady flow computations are being performed with the P&W (ATD) and the Rocketdyne baseline configurations of the SSME LO2 turbine turnaround duct (TAD) and heat exchanger (HEX). The work is in support of the HEX inner turning vane cracking investigation. Fatigue cracking has occurred during hot firings with the P&W configuration on the HEX inner vane, and it appears the fix will involve changes to the TAD splitter vane position and to the TAD inner wall curvature to reduce the dynamic loading on the inner vane. Unsteady flow computations on the P&W baseline and fix and on the Rocketdyne baseline reference follow steady-flow screening computations done by MSFC/ED32 on several trial configurations arriving at the fix. The P&W TAD inlet velocity profile has a strong radial velocity component that directs the flow toward the inner wall and raises the local velocity a factor of two and the dynamic pressure a factor, of four. The fix is intended to redistribute the flow more evenly across the HEX inner and outer vanes like the Rocketdyne baseline reference. Vane buffeting at frequencies around 4,000 Hz is the leading suspected cause of the problem. Our simulations (work in progress) are being done with the USA 2D axisymmetric code approximating the flow as axisymmetric u+v 2D (axial, u, and radial, v, components only). The HEX coils are included in the model to make sure the fix does not adversely affect the HEX environment. Turbulent kinetic energy, k, levels where k = 1/2 v' rms2 are locally as high as 10,000 ft2/sec2 for the P&W baseline at the engine interface (between the TAD and HEX) at the HEX inner vane location. However, k is less than 8,000 on the HEX outer vane and only about 4,500 on the HEX inner vane for the Rocketdyne baseline. Unsteady turbulence intensity, v'rms/v, and pressure, p', are being computed in the present computations to compare with steady-flow Reynolds-averaged computations where p'rms = const (pk) for overall rms random turbulence from 0.1 to 12,000 Hz frequency. Random overall static, p'rms fluctuations as large as 1.7 psi are estimated from k on the HEX inner vane for the P&W baseline configuration but only about 0.7 psi for the Rocketdyne configuration.

Dougherty, N. S.↗
NASA NTRS1972

Recent satellite measurements of the morphology and dynamics of the plasmasphere.

The characteristic morphology and dynamics of the plasmasphere vary with local time and with geomagnetic conditions. On the nightside the plasmapause position changes predictably with changing magnetic activity. Once established at a specific L-shell value, the steep density gradient on the nightside corotates into the dayside, where filling from the ionosphere takes place. In the duskside bulge region the characteristic density profile inside the plasmapause displays a smooth decrease proportional to 1/R to the fourth power where R is radial distance. Plasmasphere morphology and dynamics can be understood in terms of a time-varying convection electric-field model of the magnetosphere that includes the bulge region as part of the main circulation pattern of the plasmasphere.

Chappell, C. R.↗
NASA NTRS2013

Conductance of Ion Channels - Theory vs. Experiment

Transmembrane ion channels mediate a number of essential physiological processes in a cell ranging from regulating osmotic pressure to transmission of neural signals. Kinetics and selectivity of ion transport is of critical importance to a cell and, not surprisingly, it is a subject of numerous experimental and theoretical studies. In this presentation we will analyze in detail computer simulations of two simple channels from fungi - antiamoebin and trichotoxin. Each of these channels is made of an alpha-helical bundle of small, nongenomically synthesized peptides containing a number of rare amino acids and exhibits strong antimicrobial activity. We will focus on calculating ionic conductance defined as the ratio of ionic current through the channel to applied voltage. From molecular dynamics simulations, conductance can be calculated in at least two ways, each involving different approximations. Specifically, the current, given as the number of charges transferred through the channel per unit of time, can be obtained from the number of events in which ions cross the channel during the simulation. This method works well for large currents (high conductance values and/or applied voltages). If the number of crossing events is small, reliable estimates of current are difficult to achieve. Alternatively, conductance can be estimated assuming that ion transport can be well approximated as diffusion in the external potential given by the free energy profile. Then, the current can be calculated by solving the one-dimensional diffusion equation in this external potential and applied voltage (the generalized Nernst-Planck equation). To do so three ingredients are needed: the free energy profile, the position-dependent diffusion coefficient and the diffusive flux of ions into the channel. All these quantities can be obtained from molecular dynamics simulations. An important advantage of this method is that it can be used equally well to estimating large and small currents. In addition, once the free energy profile becomes available the full current-voltage dependence can be readily obtained. For both channels we carried out calculations using both approaches. We also tested the main assumptions underlying the diffusive model, such as uncorrelated nature of individual crossing events and Fickian diffusion. The accuracy and consistency of different methods will be discussed. Finally we will discuss how comparisons between calculated and measured ionic conductance and selectivity of transport can be used for determining structural models of the channels.

Pohorille, Andrew↗
NASA NTRS1994

(abstract) Monitoring the Freeze/Thaw Transitions in Taiga Forests Using ERS-1 SAR

Automated recording stations have been installed at the Bonanza Creek Experimental Forest, a Long Term Ecological Research (LTER) site located near Fairbanks, Alaska, in a forest stand of the Tanana River floodplain underlain by discontinuous permafrost. These stations provide a continuous record of dielectric constant and temperature of tree trunks, and soil moisture and temperature profiles down to the root zone. Along with the weather stations deployed at the same location, these measurements provide a continuous record of the environmental and phenologic conditions of the forest during a complete seasonal cycle. At the same time, ERS-1 SAR imaged the study site repeatedly from space to provide radar backscatter measurements of the forest approximately three times a month. Here, we examine the temporal dynamic of ERS-1 SAR measurements in relation with the changing environmental and phenologic state of the forest canopy and of the forest ground layers during the winter/spring and fall/winter transitions of 1992 and 1993. During these transitions, we examine whether changes in radar backscatter observed by ERS-1 may be related to freezing or thawing of the soil and vegetation in order to determine the start and end of the growing season for the forest. The results of this analysis are used in turn to determine whether similar changes are observed over larger regions. Mosaics of SAR data generated along three different North-South Alaskan ERS-1 transects that intercept with our study site are used in combination with hourly air temperature and daily precipitation rates gathered at airport weather stations by the National Weather Service. Results obtained using ERS-1 data collected from January 1992 to mid-1993 will be discussed.

permafrost temporal dynamics Alaska ecology radar ↗
DOE OSTI2016

Linking density functional and mode coupling models for supercooled liquids

We compare predictions from two familiar models of the metastable supercooled liquid, respectively, constructed with thermodynamic and dynamic approaches. In the so called density functional theory the free energy F[ρ] of the liquid is a functional of the inhomogeneous density ρ(r). The metastable state is identified as a local minimum of F[ρ]. The sharp density profile characterizing ρ(r) is identified as a single particle oscillator, whose frequency is obtained from the parameters of the optimum density function. On the other hand, a dynamic approach to supercooled liquids is taken in the mode coupling theory (MCT) which predict a sharp ergodicity-non-ergodicity transition at a critical density. The single particle dynamics in the non-ergodic state, treated approximately, represents a propagating mode whose characteristic frequency is computed from the corresponding memory function of the MCT. The mass localization parameters in the above two models (treated in their simplest forms) are obtained, respectively, in terms of the corresponding natural frequencies depicted and are shown to have comparable magnitudes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Tracing the outflow kinematics in Type 2 active galactic nuclei

Here, we used a sample of 577 spectra of active galactic nuclei Type 1.8-2 (z < 0.25) taken from the Sloan Digital Sky Survey to trace the influence of the outflow kinematics on the profiles of different emission lines (Hβ, [O III], Hα, [N II], [S II]). All considered lines were fitted with two Gaussian components: one that fits the core of the line, and another that fits the wings. We provide a procedure for decomposition of the Hα+[N II] wavelength band for spectra where these lines overlap. The influence of the gravitational and non-gravitational kinematics on the line components is investigated by comparing the dispersions of the line components with stellar velocity dispersion. We find that wing components of all the considered emission lines have pure non-gravitational kinematics. The core components are consistent with gravitational kinematics for the Hα, [N II], and [S II] lines, while in the [O III] there is evidence for contribution from non-gravitational kinematics. We adopted the wing components as a proxy for the outflow contribution and investigated the outflow kinematics by analysing the correlations between the widths and shifts of the wing components of different lines. For this purpose, we used the subsets in which wing components are detected in both compared lines, and can be fitted independently. We find strong correlations between wing component shifts, as well as between wing component widths of all considered lines, with the exception of the Hβ wing component width. These correlations indicate that outflow dynamics systemically affect all emission lines in the spectrum. However, it reflects with a different strength in their profiles, which is observed as different widths of the wing components. This is investigated by comparison of the mean widths of the wing components in subsets where wing components are present in all lines. The strongest outflow signature is observed in the [O III] lines, which have the broadest wing components; weaker outflow signatures are found in Hα and [N II], and the weakest is found for [S II]. These results imply that the considered lines arise in different parts of an outflowing region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Assessing the Dynamics and Complexity of Disease Pathogenicity Using 4-Dimensional Immunological Data

Investigating disease pathogenesis and personalized prognostics are major biomedical needs. Because patients sharing the same diagnosis can experience different outcomes, such as survival or death, physicians need new personalized tools, including those that rapidly differentiate several inflammatory phases. To address these topics, a pattern recognition-based method (PRM) that follows an inverse problem approach was designed to assess, in <10 min, eight concepts: synergy , pleiotropy , complexity , dynamics , ambiguity , circularity , personalized outcomes , and explanatory prognostics (pathogenesis). By creating thousands of secondary combinations derived from blood leukocyte data, the PRM measures synergic, pleiotropic, complex and dynamic data interactions, which provide personalized prognostics while some undesirable features—such as false results and the ambiguity associated with data circularity-are prevented. Here, this method is compared to Principal Component Analysis (PCA) and evaluated with data collected from hantavirus-infected humans and birds that appeared to be healthy. When human data were examined, the PRM predicted 96.9 % of all surviving patients while PCA did not distinguish outcomes. Demonstrating applications in personalized prognosis, eight PRM data structures sufficed to identify all but one of the survivors. Dynamic data patterns also distinguished survivors from non-survivors, as well as one subset of non-survivors, which exhibited chronic inflammation. When the PRM explored avian data, it differentiated immune profiles consistent with no, early, or late inflammation. Yet, PCA did not recognize patterns in avian data. Findings support the notion that immune responses, while variable, are rather deterministic: a low number of complex and dynamic data combinations may be enough to, rapidly, unmask conditions that are neither directly observable nor reliably forecasted.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Large-scale atomistic simulations of low-energy helium implantation into tungsten single crystals

Large-scale molecular dynamics simulations of post-implantation helium behavior in plasma-facing tungsten single crystals reveal orientation-dependent depth profiles, surface evolution patterns, and other crystallographic and diffusion-related characteristics of helium behavior in tungsten during the first microsecond. Here, the flux of implanted helium atoms studied, Γ ≈ 4 × 10 25 m -2 s -1 , is about one order of magnitude larger than that expected in ITER, the experimental fusion reactor currently being constructed in France. With simulation times on the order of 1 μs, these results serve to discover the mechanisms involved in surface evolution as well as to serve as benchmarks for coarse-grained simulations such as kinetic Monte Carlo and continuum-scale drift–reaction–diffusion cluster dynamics simulations. The findings of our large-scale simulations are significant due to diminished finite-size effects and the longer times reached (corresponding to higher fluences). Specifically, our findings are drastically different from findings published previously in the literature for (001) surfaces under a helium flux of Γ ~ 10 28 m -2 s -1 , which is typical of smaller size and shorter time atomistic simulations. In particular, this study highlights the atomic-scale materials processes relevant to helium entrapment and transport in metals, which have implications not only for nuclear fusion–relevant processes, but also helium-induced embrittlement in irradiated materials such as hospital equipment and fission reactor materials.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Mean Dynamic Topography of the Arctic Ocean

ICESat and Envisat altimetry data provide measurements of the instantaneous sea surface height (SSH) across the Arctic Ocean, using lead and open water elevation within the sea ice pack. First, these data were used to derive two independent mean sea surface (MSS) models by stacking and averaging along-track SSH profiles gathered between 2003 and 2009. The ICESat and Envisat MSS data were combined to construct the high-resolution ICEn MSS. Second, we estimate the 5.5-year mean dynamic topography (MDT) of the Arctic Ocean by differencing the ICEn MSS with the new GOCO02S geoid model, derived from GRACE and GOCE gravity. Using these satellite-only data we map the major features of Arctic Ocean dynamical height that are consistent with in situ observations, including the topographical highs and lows of the Beaufort and Greenland Gyres, respectively. Smaller-scale MDT structures remain largely unresolved due to uncertainties in the geoid at short wavelengths.

Arctic Ocean↗
DOE OSTI2013

Report on the Implementation of Homogeneous Nucleation Scheme in MARMOT-based Phase Field Simulation

In this report, we summarized our effort in developing mesoscale phase field models for predicting precipitation kinetics in alloys during thermal aging and/or under irradiation in nuclear reactors. The first part focused on developing a method to predict the thermodynamic properties of critical nuclei such as the sizes and concentration profiles of critical nuclei, and nucleation barrier. These properties are crucial for quantitative simulations of precipitate evolution kinetics with phase field models. Fe-Cr alloy was chosen as a model alloy because it has valid thermodynamic and kinetic data as well as it is an important structural material in nuclear reactors. A constrained shrinking dimer dynamics (CSDD) method was developed to search for the energy minimum path during nucleation. With the method we are able to predict the concentration profiles of the critical nuclei of Cr-rich precipitates and nucleation energy barriers. Simulations showed that Cr concentration distribution in the critical nucleus strongly depends on the overall Cr concentration as well as temperature. The Cr concentration inside the critical nucleus is much smaller than the equilibrium concentration calculated by the equilibrium phase diagram. This implies that a non-classical nucleation theory should be used to deal with the nucleation of Cr precipitates in Fe-Cr alloys. The growth kinetics of both classical and non-classical nuclei was investigated by the phase field approach. A number of interesting phenomena were observed from the simulations: 1) a critical classical nucleus first shrinks toward its non-classical nucleus and then grows; 2) a non-classical nucleus has much slower growth kinetics at its earlier growth stage compared to the diffusion-controlled growth kinetics. 3) a critical classical nucleus grows faster at the earlier growth stage than the non-classical nucleus. All of these results demonstrated that it is critical to introduce the correct critical nuclei into phase field modeling in order to correctly capture the kinetics of precipitation. In most alloys the matrix phase and precipitate phase have different concentrations as well as different crystal structures. For example, Cu precipitates in FeCu alloys have fcc crystal structure while the matrix Fe-Cu solid solution has bcc structure at low temperature. The WBM model and KimS model, where both concentrations and order parameters are chosen to describe the microstructures, are commonly used to model precipitations in such alloys. The WBM and KimS models have not been implemented into Marmot yet. In the second part of this report, we focused on implementing the WBM and KimS models into Marmot. The Fe-Cu alloys, which are important structure materials in nuclear reactors, was taken as the model alloys to test the models.

Fe Alloys↗
DOE OSTI2017

Idealized Simulations of a Squall Line from the MC3E Field Campaign Applying Three Bin Microphysics Schemes: Dynamic and Thermodynamic Structure

The squall line event on May 20, 2011, during the Midlatitude Continental Convective Clouds (MC3E) field campaign has been simulated by three bin (spectral) microphysics schemes coupled into the Weather Research and Forecasting (WRF) model. Semi-idealized three-dimensional simulations driven by temperature and moisture profiles acquired by a radiosonde released in the pre-convection environment at 1200 UTC in Morris, Oklahoma show that each scheme produced a squall line with features broadly consistent with the observed storm characteristics. However, substantial differences in the details of the simulated dynamic and thermodynamic structure are evident. These differences are attributed to different algorithms and numerical representations of microphysical processes, assumptions of the hydrometeor processes and properties, especially ice particle mass, density, and terminal velocity relationships with size, and the resulting interactions between the microphysics, cold pool, and dynamics. This study shows that different bin microphysics schemes, designed to be conceptually more realistic and thus arguably more accurate than bulk microphysics schemes, still simulate a wide spread of microphysical, thermodynamic, and dynamic characteristics of a squall line, qualitatively similar to the spread of squall line characteristics using various bulk schemes. Future work may focus on improving the representation of ice particle properties in bin schemes to reduce this uncertainty and using the similar assumptions for all schemes to isolate the impact of physics from numerics.

Xue, Lulin↗
DOE OSTI2023

Evaluation of 6-OxP-CD, an Oxime-based cyclodextrin as a viable medical countermeasure against nerve agent poisoning: Experimental and molecular dynamic simulation studies on its inclusion complexes with cyclosarin, soman and VX

The ability of the cyclodextrin-oxime construct 6-OxP-CD to bind and degrade the nerve agents Cyclosarin (GF), Soman (GD) and S -[2-[Di(propan-2-yl)amino]ethyl] O -ethyl methylphosphonothioate (VX) has been studied using 31 P-nuclear magnetic resonance (NMR) under physiological conditions. While 6-OxP-CD was found to degrade GF instantaneously under these conditions, it was found to form an inclusion complex with GD and significantly improve its degradation (t 1/2 ~ 2 hrs) relative over background (t 1/2 ~ 22 hrs). Consequently, effective formation of the 6-OxP-CD:GD inclusion complex results in the immediate neutralization of GD and thus preventing it from inhibiting its biological target. In contrast, NMR experiments did not find evidence for an inclusion complex between 6-OxP-CD and VX, and the agent’s degradation profile was identical to that of background degradation (t 1/2 ~ 24 hrs). As a complement to this experimental work, molecular dynamics (MD) simulations coupled with Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) calculations have been applied to the study of inclusion complexes between 6-OxP-CD and the three nerve agents. These studies provide data that informs the understanding of the different degradative interactions exhibited by 6-OxP-CD with each nerve agent as it is introduced in the CD cavity in two different orientations (up and down). For its complex with GF, it was found that the oxime in 6-OxP-CD lies in very close proximity (P GF …O Oxime ~ 4–5 Å) to the phosphorus center of GF in the ‘down GF ’ orientation for most of the simulation accurately describing the ability of 6-OxP-CD to degrade this nerve agent rapidly and efficiently. Further computational studies involving the center of masses (COMs) for both components (GF and 6-OxP-CD) also provided some insight on the nature of this inclusion complex. Distances between the COMs (ΔCOM) lie closer in space in the ‘down GF ’ orientation than in the ‘up GF ’ orientation; a correlation that seems to hold true not only for GF but also for its congener, GD. In the case of GD, calculations for the ‘down GD ’ orientation showed that the oxime functional group in 6-OxP-CD although lying in close proximity (P GD …O Oxime ~ 4–5 Å) to the phosphorus center of the nerve agent for most of the simulation, adopts another stable conformation that increase this distance to ~ 12–14 Å, thus explaining the ability of 6-OxP-CD to bind and degrade GD but with less efficiency as observed experimentally (t 1/2 ~ 4 hr. vs. immediate). Lastly, studies on the VX:6-OxP-CD system demonstrated that VX does not form a stable inclusion complex with the oxime-bearing cyclodextrin and as such does not interact in a way that is conducive to an accelerated degradation scenario. Collectively, these studies serve as a basic platform from which the development of new cyclodextrin scaffolds based on 6-OxP-CD can be designed in the development of medical countermeasures against these highly toxic chemical warfare agents.

60 APPLIED LIFE SCIENCES↗