Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sharp features”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2014

Design and Expected Performance of GISMO-2, a Two Color Millimeter Camera for the IRAM 30 m Telescope

We present the main design features for the GISMO-2 bolometer camera, which we build for background-limited operation at the IRAM 30 m telescope on Pico Veleta, Spain. GISMO-2 will operate simultaneously in the 1 and 2 mm atmospherical windows. The 1 mm channel uses a 32 × 40 TES-based backshort under grid (BUG) bolometer array, the 2 mm channel operates with a 16 × 16 BUG array. The camera utilizes almost the entire full field of view provided by the telescope. The optical design of GISMO-2 was strongly influenced by our experience with the GISMO 2mm bolometer camera, which is successfully operating at the 30 m telescope. GISMO is accessible to the astronomical community through the regularIRAMcall for proposals.

Two-Color mm camera↗
DOE OSTI2022

Prediction of DIII-D Pedestal Structure from Externally Controllable Parameters

The sharp increase of pressure at the edge of a high confinement mode (H-mode) plasma, the pedestal, strongly impacts overall plasma performance. Predicting the pedestal is a necessity to control and optimize tokamak operations. An experimental data-driven machine learning (ML) approach is presented that predicts the pedestal heights and widths of electron density (ne) and electron temperature (Te) profiles as well as the separatrix ne from externally controllable parameters such as the plasma shape, heating method and power, and gas puff rate and integrated gas puff. The OMFIT framework was used with DIII-D data to efficiently, robustly, and automatically build a database of pedestal parameters to train machine learning models. Database creation was enabled by the search engine tool for DIII-D data, TokSearch, which parallelizes data fetching, enabling fast searches through basic signals of thousands of DIII-D shots and selection of relevant time intervals. Principal Component Analysis (PCA) separated the database into three clusters that represent classes of plasma shapes that are regularly used in DIII-D. The most important parameters for setting the pedestal structure were plasma current (Ip), toroidal magnetic field (Bφ), neutral beam heating power (PNBI) and shaping quantities. The Deep Jointly Informed Neural Networks (DJINN) algorithm was applied to identify suitable neural network (NN) architectures that appropriately capture the features of the pedestal database. Separate NNs were implemented for each pedestal parameter, and ensembling methods were used to improve the prediction accuracy and allowed estimation of the prediction uncertainty. The pedestal predictions of the test dataset lie within the measurement uncertainties of the pedestal parameters. The NN outperformed simple Linear Regression (LR) analysis, indicating non-linear dependencies in the pedestal structure. The presented achievements illustrate a promising path for future research, using feature extraction to infer experimental trends and thereby improve pedestal models as well as deploying NN for a fast pedestal prediction in DIII-D scenario development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Proposed Mars Surveyor Landing Sites in Northern Meridiani Sinus, Southern Elysium Planitia, and Argyre Planitia

Our objective is to propose two landing sites that the Mars Surveyor 2001 Lander and Athena Rover could go to on Mars that should meet the safety requirements of the spacecraft landing system and optimize surface operations (chiefly driven by power and communications requirements). An additional site within Argyre Planitia, initially proposed by Parker to the Mars Surveyor Landing Site program, is also proposed for potential consideration for post-2001 missions to Mars, as it is well outside the current latitude limits for the Athena Rover. All three sites are designed to be situated as close to a diversity of geologic units within a few kilometers of the landing site so that diversity can be placed in a geologic context. This objective is very different from the Mars Pathfinder requirement to land at a site with a maximum chance for containing a diversity of rocks within a few tens of meters of the lander. That requirement was driven by the Sojourner mobility limit of a few tens of meters. It can be argued that the Athena project, with its much larger mobility capability, might actually want to avoid such a site, because placing collected samples in geologic context would be difficult. While it has been argued, both before and after the Mars Pathfinder landing, that the provenance for local blocks may be determined by orbiter spectra, primarily from the MGS TES instrument, our ability to do so has yet to be demonstrated. Indeed, several months after conclusion of the Pathfinder mission, we have yet to reach a consensus on the composition of local materials. Our primary data set for selecting a landing site within the latitude and elevation constraints of the 2001 mission is the Viking Orbiter image archive. The site must be selected to place the landing ellipse so as to avoid obvious hazards, such as steep slopes, large or numerous craters, or abundant large knobs. For this purpose, we chose a resolution limit of better than 50 m/pixel. This necessarily excludes from the present study images from current and future orbiter spacecraft, until such data does become readily available. Within each proposed region, it may be possible to identify additional sites once these data become available. Second, the fine-component thermal inertia data, should be greater than about 5 or 6 cgs Units (10(exp -3) cal/sq cm s(exp -0.5)/K). Low thermal inertias imply dusty environments, which could pose a mobility hazard. Similarly, the albedo of the site should not be particularly high, which would also suggest dusty surfaces. Low albedos are preferred, as they often coincide with low Viking red:violet ratios and indicate less dusty surfaces. Next, the Modeled Block Abundance should also not be too high or too low. Based on the Viking Lander and Mars Pathfinder experiences, percentages of blocks should be on the order of 5-25%. Too many blocks could pose a hazard to the landing and mobility. Too few blocks could also indicate a dusty surface. Primary Landing Site: Northern Meridiani Sinus (Proposed by T. J. Parker and K., S. Edgett) Vital Statistics: (1) Latitude, Longitude: 0-3 N, 350-2 W. *Elevation (Viking): about0.5-1.5 Ian. (2) Viking Orbiter Image coverage: Excellent coverage by 15 - 25 m/pixel images (orbits 709A and 410B). Possible stereo coverage in region where two orbits overlap (probably small parallax angle, as these orbits are not listed in NASA Contractor Report 3501) (3) Albedo: about .18 -.26 (4) Block Abundance: 5-26% (5)Fine-Component Thermal Inertia: 5-9 cgs units This region consists of bright deposits similar to those described by Edgett et al, that also lie within a prominent dark albedo region. These deposits are flat-lying, to such a degree that they ramp against topography rather than draping over it. This led Edgett and Parker to suggest that they may be subaqueous sediments, possibly lacustrine or marine evaporites, laid down sometime from the late Noachian to middle Hesperian (age determination pending crater counts). A contact between this material and elevated, dissected highlands to the south was identified , and is described by Edgett et al. Our desire in proposing this landing site is to sample the edge of this deposit where it has been exposed through etching, presumably eolian deflation (the deposit, though in the highlands, is itself only lightly to moderately cratered). This should enable access to in situ stratigraphy. The actual landing site will be selected where slopes are not expected to be steep, such that the rover itself should be able to traverse them and sample layered materials on the way, either up or down the slope. Perhaps due to uncertainties at this time as to the friability or meter-scale roughness of the deposit, it might make sense to place the landing ellipse on the exhumed highland surface adjacent to the deflated margin of the deposit and plan on driving to the deposit rather than landing on it and driving downslope. This should also enable imaging the margin for evidence of layering should it prove too difficult to climb. A target ellipse on the highland surface should also allow Athena access to ancient Noachian highland materials, particularly if placed near crater ejecta or an inlier of knobby material. Secondary Landing Site: Southern Elysium Planitia (Proposed by T. J. Parker) Vital Statistics: (1) Latitude, Longitude: 1.5-3.5 S, 195-198 W. (2) Elevation (Viking): -1.0 km. (3) Viking Orbiter Image coverage: Excellent coverage by 15 - 25 m/pixel images (orbit 725). Possible stereo coverage between images from beginning and end of orbit that overlap (probably small parallax angle) (4) Albedo: about .27-.28 (5) Block Abundance: 4-7% (6) Fine-Component Thermal Inertia: about 3 cgs units This region consists of eroded knobby material, probably of Noachian age, though much of the crater population has been destroyed, that is onlapped at a sharp contact by an extensive plains unit in southern Elysium Planitia that is Amazonian in age. The plains materials have been attributed to unusually low-viscosity flood lavas from fissures south of the Elysium volcanic rise, or to lacustrine materials associated with a large, Amazonian lake at the source of Marte Vallis. Parker and Schenk presented evidence in support of the latter interpretation, though they attributed the putative shore morphology to an embayment of a northern plains ocean into the southern Elysium region. Detailed examination of the margin of the deposit, showing erosion, not simply burial, of small crater rims and fluidized ejecta blankets, also points to lacustrine or marine sedimentation rather than volcanic plains burial. The plains surface exhibits a "crusty" appearance that many researchers have attributed to pressure ridges in lava flows. In a lacustrine context, they also resemble pressure ridges in desiccated evaporite deposits and salt-rimmed pools (now dry) similar in scale and morphology to spectacular, hundred meter-scale pool rims in alkaline Lake Natron, East African Rift. The eroded highland margin surface adjacent to these plains appears to be fairly smooth, even at 15 m/pixel. Isolated knob inliers are scattered from a few kilometers to several tens of "kilometers apart. Heights of the knobs have not been measured yet but, based on experience with similar features in the Pathfinder landing ellipse, are probably typically on the order of several tens of meters high and smaller, though some of the largest knobs in the region are probably up to a few hundred meters high. Two craters larger than a kilometer in diameter, with fluidized deposits, lie nearby the proposed landing site. Very high-resolution images from MOC should help to determine whether a landing site navigable by the Athena rover could be placed in this region. The space between knobs and craters is large enough to enable placement of a target landing ellipse between them but still provide access to one or more of them and to the margin of the Elysium plains material. Post-2001 Mars Surveyor Landing Site: Argyre Planitia (Proposed by T. J. Parker) Vital Statistics: (1) Latitude, Longitude: 55-56 S, 41-43 W. (2) Elevation (Viking): 1.0 km. (3) Viking Orbiter Image coverage: Excellent coverage by 40 m/pixel images (orbits 567B, 568B, and 569B). Excellent stereo coverage with large parallax angles over the entire landing site region, and much of central and southern Argyre. (4) Albedo: about .23-.24 (5) Block Abundance: No data (6) Fine-Component Thermal Inertia: No data The floors of both the Argyre and Hellas basins contain etched layered materials that are probably thick accumulations of channel or lacustrine sediments. The deposits in Hellas are much more eroded than those in Argyre, and Hellas lacks a channel outlet. Argyre is unique in that Uzboi Vallis flowed out of the basin, requiring overflow of a standing body of water within Argyre. This makes it the largest impact basin on Mars with channels both draining into it and flowing out from it. Hellas' channels may be catastrophic flood channels, whereas Argyre was fed by modest-scale valley networks, though the outlet at Uzboi Vallis was a catastrophic flood Highland craters and basins of this kind should be high-priority landing targets for missions intended to focus on the search for either prebiotic organic materials or even simple fossil microorganisms. Basins with internally-draining valley networks should be preferred over flood channels, as they could have provided the long-term influx of water favorable to the origin of life. (Catastrophic floods are not conducive to fossil preservation, due to their very short durations and high transportation energies). They also afford an opportunity to study the evolution of the planet's climate and volatiles during the period of time between the late Noachian and early Hesperian, when a drastic change from a proposed early warm, wet climate to one more closely resembling the modern environment is thought to have occurred. Large basin

Parker, T. J.↗
NASA NTRS1998

The Coronae of Moderate-Mass Giants in the Hertzsprung Gap and the Clump

We have used the Roentgensatellit (ROSAT), the Extreme Ultraviolet Explorer (EUVE), and the Hubble Space Telescope (HST) to measure X-ray and ultraviolet emissions of moderate-mass (APPROX. 2-3 solar mass) giants in the Hertzsprung gap (spectral types early F to mid-G) and the post-helium flash "clump" (approx. G8-K0). Our motivation was to document the evolution of hot coronae (T greater than 10(exp 6)K) along the post-main-sequence trajectories traveled by such stars in order to gain insight concerning the "X-ray deficiency" of the F-GO giants and the strong braking of stellar rotation at the red edge of the Hertzsprung gap. With few exceptions, Hertzsprung gap and clump giants observed by ROSAT PSPC show hot (T approx. 10(exp 7)K) coronal energy distributions, regardless of any X-ray deficiency, EUVE spectra of gap star 31 Com (G0 111) indicate a broad coronal emission measure hump at approx. 10(exp 7.2)K, while the active clump giant beta Ceti (K0 III) displays a sharp peak at approx. 10(exp 6.8)K, as seen previously in the mixed clump/gap binary Capella (alpha Aur: G8 III + G0 III). The gap giants upsilon Peg (F8 III) and 24 UMa (G4 III) have EUV emissions of intermediate temperature (approx. 10(exp 7.0)K). The stars 31 Com, psi(sup 3) Psc (G0 III), and beta Cet exhibit redshifted transition zone (TZ: approx. 10(exp 5)K) lines in HST GHRS spectra, as reported earlier in Procyon (alpha CMi: F5 IV-V) and Capella G0. Such redshifts on the Sun are thought to signify flows in magnetic loops. beta Cas (F2 III)-a rare soft coronal source among the gap stars-displays blueshifts of C iv and 0 iv], although emissions at cooler and hotter temperatures are near the photospheric velocity. The remarkably broad line profiles of the fastest rotating gap giants suggest that the 10(exp 5) K "subcoronal " emission zones extend to h approx. R(sub *) above the photosphere, about 50 scale heights. In contrast to the TZ line redshifts, the upper chromospheric emissions (e.g., Mg II and Si III) of 31 Com and upsilon(sup 3) Psc have blueshifted cores. Blue-asymmetric peaks in the solar Mg a lines are thought to indicate dynamical heating in the chromosphere. Observations of the H(sub I) Ly(alpha) feature of 31 Com taken 9 months apart reveal striking profile changes, reminiscent of those noted previously in the Ly(alpha) blue peak of the Capella G0 star. We used the far-ultraviolet diagnostics, in combination with ROSAT X-ray photometry and EUVE high-excitation line strengths, to constrain physical models of the stellar outer atmospheres. Quasi-static magnetic loops can simulate the empirical coronal emission measures of the giant stars, but the inferred pressures for sensible loop lengths conflict with direct measurements of subcoronal densities. Furthermore, the high rate of emission at approx. 10(exp 5) K cannot be explained by thermal conduction down the legs of hot quasi-static loops. On the other hand, the possible existence of elongated (l - R(sub *) emission structures on the gap giants leads to a speculative scenario to explain the X-ray deficiency. It is based on the increased importance of the dynamical filling phase ("explosive evaporation") of the loop life cycle; conductive cooling, yielding TZ emissions at the footpoints, when the heating is interrupted; and the possibility for transitions between " hot " and " cool " energy balance solutions owing to dynamical suspension and centrifugal trapping of the cooling gas. The long loops might represent a vestigial global " magnetosphere " inherited from the main-sequence phase, which ultimately is disrupted near approx. G0 by the deepening convective envelope and growth of a more solar-like dynamo. Coronal emissions might be boosted temporarily as the X-ray deficiency is removed but soon would be quenched by wind braking previously inhibited by the magnetospheric "dead zone."

Ayres, Thomas R.↗
DOE OSTI2021

Prediction of DIII-D Pedestal Structure From Externally Controllable Parameters

The sharp increase of pressure at the edge of a high confinement mode (H-mode) plasma, the pedestal, strongly impacts overall plasma performance. Predicting the pedestal is a necessity to control and optimize tokamak operations. Here, an experimental data-driven machine learning (ML) approach is presented that predicts the pedestal heights and widths of electron density (n e ) and electron temperature (T e ) profiles as well as the separatrix ne from externally controllable parameters such as the plasma shape, heating method and power, and gas puff rate and integrated gas puff. The OMFIT framework was used with DIII-D data to efficiently, robustly, and automatically build a database of pedestal parameters to train machine learning models. Database creation was enabled by the search engine tool for DIII-D data, TokSearch, which parallelizes data fetching, enabling fast searches through basic signals of thousands of DIII-D shots and selection of relevant time intervals. Principal Component Analysis (PCA) separated the database into three clusters that represent classes of plasma shapes that are regularly used in DIII-D. The most important parameters for setting the pedestal structure were plasma current (I p ), toroidal magnetic field (B Φ ), neutral beam heating power (P NBI ) and shaping quantities. The Deep Jointly Informed Neural Networks (DJINN) algorithm was applied to identify suitable neural network (NN) architectures that appropriately capture the features of the pedestal database. Separate NNs were implemented for each pedestal parameter, and ensembling methods were used to improve the prediction accuracy and allowed estimation of the prediction uncertainty. The pedestal predictions of the test dataset lie within the measurement uncertainties of the pedestal parameters. The NN outperformed simple Linear Regression (LR) analysis, indicating non-linear dependencies in the pedestal structure. The presented achievements illustrate a promising path for future research, using feature extraction to infer experimental trends and thereby improve pedestal models as well as deploying NN for a fast pedestal prediction in DIII-D scenario development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2012

Global Patterns of Tectonism on Titan from Mountain Chains and Virgae

This research is based on the exploration of tectonic patterns on Titan from a global perspective. Several moons in the outer solar system display patterns of surface tectonic features that imply global stress fields driven or modified by global forces. Patterns such as these are seen in Europa's tidally induced fracture patterns, Enceladus's tiger stripes, and Ganymede's global expansion induced normal fault bands. Given its proximity to Saturn, as well as its eccentric orbit, tectonic features and global stresses may be present on Titan as well. Titan displays possible tectonic structures, such as mountain chains along its equator (Radebaugh et al. 2007), as well as the unexplored dark linear streaks termed virgae by the IAU. Imaged by Cassini with the RADAR instrument, mountain chains near the equator are observed with a predominante east-west orientation (Liu et al. 2012, Mitri et al. 2010). Orientations such as these can be explained by modifications in the global tidal stress field induced by global contraction followed by rotational spin-up. Also, due to Titan's eccentric orbit, its current rotation rate may be in an equilibrium between tidal spin-up near periapsis and spin-down near apoapsis (Barnes and Fortney 2003). Additional stress from rotational spin-up provides an asymmetry to the stress field. This, combined with an isotropic stress from radial contraction, favors the formation of equatorial mountain chains in an east-west direction. The virgae, which have been imaged by Cassini with both the Visual and Infrared Mapping Spectrometer (VIMS) and Imaging Science Subsystem (ISS) instruments, are located predominately near 30 degrees latitude in either hemisphere. Oriented with a pronounced elongation in the east-west direction, all observed virgae display similar characteristics: similar relative albedos as the surrounding terrain however darkened with an apparent neutral absorber, broken-linear or rounded sharp edges, and connected, angular elements with distinct, linear edges. Virgae imaged during northern latitude passes are oriented with their long dimensions toward Titan's antiSaturn point. If the virgae are of tectonic origin, for instance if the turn out to be i.e. grabens, they could serve as markers to Titan's global stress field. Using them in this way allows for a mapping of global tectonic patterns. These patterns will be tested for consistency against the various sources of global stress and orientations of mountain chains. By determining what drives Titan's tectonics globally, we will be able to place Titan within the context of the other outer planet icy satellites.

Cook, C.↗
DOE OSTI2012

Synthesis and morphology of Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals with tungsten bronze structure in RE{sub 2}O{sub 3}-BaO-Nb{sub 2}O{sub 5}-B{sub 2}O{sub 3} glasses (RE: Sm, Eu, Gd, Dy, Er)

Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals with a tetragonal tungsten bronze (TTB) structure are synthesized using a conventional glass crystallization technique in 2.3RE{sub 2}O{sub 3}-27.4BaO-34.3Nb{sub 2}O{sub 5}-36B{sub 2}O{sub 3} (mol%) (RE=Sm, Eu, Gd, Dy, and Er) glasses. One sharp crystallization peak is observed at {approx}670 Degree-Sign C in both powdered and bulk glasses, and the formation of Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals with unit cell parameters of a{approx}1.24 nm and c{approx}0.39 nm was confirmed. It is found from high resolution transmission electron microscope observations that the morphology of Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals is ellipsoidal. Their average particle size is in the range of 15-60 nm and decreases with decreasing ionic radius of RE{sup 3+} being present in the precursor glasses. The optical transparent crystallized glass (bulk) shows the total photoluminescence (PL) quantum yield of 53% in the visible region of Eu{sup 3+} ions, suggesting a high potential of Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals as PL materials. - Graphical abstract: This figure shows a TEM photograph for the heat-treated (667 Degree-Sign C, 3 h) sample of 2.3Dy{sub 2}O{sub 3}-27.4BaO-34.3Nb{sub 2}O{sub 5}-36B{sub 2}O{sub 3}. An ellipsoidal-shaped Ba{sub 1-x}Dy{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystal with diameters of 17 and 28 nm is observed. The ellipsoidal morphology is a common feature in Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals synthesized by the crystallization of 2.3RE{sub 2}O{sub 3}-27.4BaO-34.3Nb{sub 2}O{sub 5}-36B{sub 2}O{sub 3} glasses. Highlights: Black-Right-Pointing-Pointer Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals with a tetragonal tungsten bronze structure are synthesized. Black-Right-Pointing-Pointer A glass crystallization technique was applied. Black-Right-Pointing-Pointer The morphology of Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals is ellipsoidal. Black-Right-Pointing-Pointer The average particle size with 15-60 nm decreases with decreasing ionic radius of RE{sup 3+}. Black-Right-Pointing-Pointer Ba{sub 1-x}RE{sub 2x/3}Nb{sub 2}O{sub 6} nanocrystals have a high potential as photoluminescence materials.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Coherent Vibrational Wavepacket Dynamics in Platinum(II) Dimers and Their Implications

Vibrational coherence in the metal-metal-to-ligand-charge transfer (MMLCT) excited state of cyclometalated platinum dimers with a pseudo C 2 symmetry was investigated where two nearly degenerate transitions from HOMO (metal-metal σ* orbital) to higher energy ligand π* orbitals could be simultaneously induced. We observed oscillatory features in femtosecond degenerate transient absorption (TA) signals from complexes [(ppy)Pt(μ- t Bu 2 pz)] 2 (1) and anti-[(ppy)Pt(μ-pyt)] 2 (2), which are attributed to coherent nuclear motions that modulate the HOMO (antibonding σ*) energy level, and hence the energy for the MMLCT transition. The characteristics of such coherent nuclear motions, such as the oscillatory frequency and the dephasing time, differ between 1 and 2 and are explained by mainly two structural factors that could influence the vibrational coherence, the Pt-Pt distance (2.97 Å for 1 vs. 2.85 Å for 2) and molecular shape (1 in an “A” frame vs. 2 in an “H” frame). Because the electronic coupling between the two Pt atoms determines the energy splitting of the bonding σ and antibonding σ* orbital, the Pt-Pt stretching mode coupled to the MMLCT transition changes the inter Pt distance from that of the ground state. Interestingly, while 1 shows a single PtPt stretching frequency of 120 cm-1 in the MMLCT state, 2 exhibits multiple downshifted frequencies (80 and 105 cm-1) in the MMLCT state along with a shorter vibrational dephasing time than 1. Based on the ground state optimized structures and Raman calculations, the changes evident in the vibrational wavepaket dynamics in 2 are closely correlated with the “H” framed geometry in 2 compared to the “A” frame in 1, leading to the sharp increase in - interaction between ppy ligands. Although the TA experiments do not directly reveal the ultrafast intersystem crossing (ISC) because of a strong coherent spike at early timescales, the dependence of the vibrational wavepacket dynamics on molecular geometry can be understood based on previously proposed potential energy surfaces (PES) as function of PtPt distance, suggesting that the interaction between the cyclometalating ligands can be a key factor in determining the PtPt shortening and the related energy relaxation dynamics in the Pt(II) dimers. Further experiments using femtosecond broadband TA and X-ray scattering spectroscopy are planned to investigate directly the ISC and PtPt contraction, to support the relationship between ground state molecular geometry and photoinduced structural changes in the Pt(II) dimers.

Kim, Pyosang↗
DOE OSTI2016

SU-C-BRA-04: Automated Segmentation of Head-And-Neck CT Images for Radiotherapy Treatment Planning Via Multi-Atlas Machine Learning (MAML)

Purpose: Accurate image segmentation is a crucial step during image guided radiation therapy. This work proposes multi-atlas machine learning (MAML) algorithm for automated segmentation of head-and-neck CT images. Methods: As the first step, the algorithm utilizes normalized mutual information as similarity metric, affine registration combined with multiresolution B-Spline registration, and then fuses together using the label fusion strategy via Plastimatch. As the second step, the following feature selection strategy is proposed to extract five feature components from reference or atlas images: intensity (I), distance map (D), box (B), center of gravity (C) and stable point (S). The box feature B is novel. It describes a relative position from each point to minimum inscribed rectangle of ROI. The center-of-gravity feature C is the 3D Euclidean distance from a sample point to the ROI center of gravity, and then S is the distance of the sample point to the landmarks. Then, we adopt random forest (RF) in Scikit-learn, a Python module integrating a wide range of state-of-the-art machine learning algorithms as classifier. Different feature and atlas strategies are used for different ROIs for improved performance, such as multi-atlas strategy with reference box for brainstem, and single-atlas strategy with reference landmark for optic chiasm. Results: The algorithm was validated on a set of 33 CT images with manual contours using a leave-one-out cross-validation strategy. Dice similarity coefficients between manual contours and automated contours were calculated: the proposed MAML method had an improvement from 0.79 to 0.83 for brainstem and 0.11 to 0.52 for optic chiasm with respect to multi-atlas segmentation method (MA). Conclusion: A MAML method has been proposed for automated segmentation of head-and-neck CT images with improved performance. It provides the comparable result in brainstem and the improved result in optic chiasm compared with MA. Xuhua Ren and Hao Gao were partially supported by the NSFC (#11405105), the 973 Program (#2015CB856000), and the Shanghai Pujiang Talent Program (#14PJ1404500).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Onset of normal and inverse homoclinic bifurcation in a double plasma system near a plasma fireball

Plasma fireballs are generated due to a localized discharge and appear as a luminous glow with a sharp boundary, which suggests the presence of a localized electric field such as electrical sheath or double layer structure. The present work reports the observation of normal and inverse homoclinic bifurcation phenomena in plasma oscillations that are excited in the presence of fireball in a double plasma device. The controlling parameters for these observations are the ratio of target to source chamber (n{sub T}/n{sub S}) densities and applied electrode voltage. Homoclinic bifurcation is noticed in the plasma potential fluctuations as the system evolves from narrow to long time period oscillations and vice versa with the change of control parameter. The dynamical transition in plasma fireball is demonstrated by spectral analysis, recurrence quantification analysis (RQA), and statistical measures, viz., skewness and kurtosis. The increasing trend of normalized variance reflects that enhancing n{sub T}/n{sub S} induces irregularity in plasma dynamics. The exponential growth of the time period is strongly indicative of homoclinic bifurcation in the system. The gradual decrease of skewness and increase of kurtosis with the increase of n{sub T}/n{sub S} also reflect growing complexity in the system. The visual change of recurrence plot and gradual enhancement of RQA variables DET, L{sub max}, and ENT reflects the bifurcation behavior in the dynamics. The combination of RQA and spectral analysis is a clear evidence that homoclinic bifurcation occurs due to the presence of plasma fireball with different density ratios. However, inverse bifurcation takes place due to the change of fireball voltage. Some of the features observed in the experiment are consistent with a model that describes the dynamics of ionization instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2017

Shed Vortex Structure and Phase-Averaged Velocity Statistics in Symmetric/Asymmetric Turbulent Flat Plate Wakes

The near wake of a flat plate is investigated via direct numerical simulations (DNS). Many earlier experimental investigations have used thin plates with sharp trailing edges and turbulent boundary layers to create the wake. This results in large theta divided by D (sub TE) values (theta is the boundary layer momentum thickness towards the end of the plate and D (sub TE) is the trailing edge thickness). In the present study the emphasis is on relatively thick plates with circular trailing edges (CTE) resulting in theta divided by D values less than one (D is the plate thickness and the diameter of the CTE), and vigorous vortex shedding. The Reynolds numbers based on the plate length and D are 1.255 x 10 (sup 6) and 10,000, respectively. Two cases are computed; one with turbulent boundary layers on both the upper and lower surfaces of the plate (statistically the same, symmetric wake, Case TT) and, a second with turbulent and laminar boundary layers on the upper and lower surfaces, respectively (asymmetric case, Case TL). The data and understanding obtained is of considerable engineering interest, particularly in turbomachinery where the pressure side of an airfoil can remain laminar or transitional because of a favorable pressure gradient and the suction side is turbulent. Shed-vortex structure and phase-averaged velocity statistics obtained in the two cases are compared here. The upper negative shed vortices in Case TL (turbulent separating boundary layer) are weaker than the lower positive ones (laminar separating boundary layer) at inception (a factor 1.27 weaker in terms of peak phase-averaged spanwise vorticity at first appearance of a peak). The upper vortices weaken rapidly as they travel downstream. A second feature of interest in Case TL is a considerable increase in the peak phase-averaged, streamwise normal intensity (random component) with increasing streamwise distance (x divided by D) that occurs nears the positive vortex cores. This behavior is observed for a few diameters in the near wake. This is counter to Case TT where the peak value essentially decreases with increasing x divided by D. Both these effects are examined in detail and the important contributors are identified.

direct numerical simulations↗
DOE OSTI2019

Signatures of tunable superconductivity in a trilayer graphene moiré superlattice

Understanding the mechanism of high-transition-temperature (high-T c ) superconductivity is a central problem in condensed matter physics. It is often speculated that high-Tc superconductivity arises in a doped Mott insulator1 as described by the Hubbard model 2,3,4 . An exact solution of the Hubbard model, however, is extremely challenging owing to the strong electron–electron correlation in Mott insulators. Therefore, it is highly desirable to study a tunable Hubbard system, in which systematic investigations of the unconventional superconductivity and its evolution with the Hubbard parameters can deepen our understanding of the Hubbard model. Here we report signatures of tunable superconductivity in an ABC-trilayer graphene (TLG) and hexagonal boron nitride (hBN) moiré superlattice. Unlike in ‘magic angle’ twisted bilayer graphene, theoretical calculations show that under a vertical displacement field, the ABC-TLG/hBN heterostructure features an isolated flat valence miniband associated with a Hubbard model on a triangular superlattice 5,6 where the bandwidth can be tuned continuously with the vertical displacement field. Upon applying such a displacement field we find experimentally that the ABC-TLG/hBN superlattice displays Mott insulating states below 20 kelvin at one-quarter and one-half fillings of the states, corresponding to one and two holes per unit cell, respectively. Upon further cooling, signatures of superconductivity (‘domes’) emerge below 1 kelvin for the electron- and hole-doped sides of the one-quarter-filling Mott state. The electronic behaviour in the ABC-TLG/hBN superlattice is expected to depend sensitively on the interplay between the electron–electron interaction and the miniband bandwidth. By varying the vertical displacement field, we demonstrate transitions from the candidate superconductor to Mott insulator and metallic phases. Our study shows that ABC-TLG/hBN heterostructures offer attractive model systems in which to explore rich correlated behaviour emerging in the tunable triangular Hubbard model.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Quantifying Interannual Variability of the UTLS Ozone Using Assimilation of Satellite Data

An accurate representation of spatial and temporal variability of the Upper Troposphere Lower Stratosphere (UTLS) ozone is essential for understanding both the tropospheric ozone budget and ozone s contribution to radiative forcing. The complex, dynamically driven structure of trace gas fields in the UTLS presents a challenge to data-based and modelling studies. Small features are not fully resolved in data from limb-sounding instruments such as the Microwave Limb Sounder on EOS-Aura (the EOS-MLS), but are captured in assimilation of those data as vertical structure is added from the assimilated meteorology. This will be demonstrated using a multi-year assimilation of EOS-MLS observations in the Goddard Earth Observing System, Version 5 (GEOS-5) data assimilation system. The results demonstrate the realism of the seasonal and year to year variability of laminar structures in the mid-latitudinal ozone field between years 2005-2007, for which independent validation data are available from the HIRDLS instrument. The analysis is done in the context of the underlying large scale dynamics. The lifetimes of most research instruments are too short for them to be used throughout the duration of long-term (at least 3 decades) reanalyses. For example, the EOS-MLS instrument has operated since mid-2004 until present. By contrast, Solar Backscatter Ultra Violet (SBUV) measurements provide continuous data since late 1978, but their vertical resolution is insufficient to represent the profile shape in the UTLS. Assimilation of these SBUV/2 observations in the GEOS-5 data assimilation system has hitherto not captured a realistic ozone structure in the UTLS, even though transport studies using GEOS-5 wind fields do show such structures. We show that careful construction of the background error covariance structure in GEOS-5 can lead to more realistic UTLS ozone fields when assimilating SBUV/2 observations. The reasoning behind this will be discussed, emphasizing the need to retain the sharp gradient of ozone concentrations across the tropopause. We analyze the UTLS ozone distributions in multi-year SBUV/2 assimilation experiments, comparing the results against the independent HIRDLSdataset and, for a longer period, with the MLS assimilation and discuss the consequences for tropospheric ozone and radiative forcing.

Wargan, K.↗
DOE OSTI2025

Sharp detection of low-dimensional structure in probability measures via dimensional logarithmic Sobolev inequalities

Identifying low-dimensional structure in high-dimensional probability measures is an essential pre-processing step for efficient sampling. To identify this structure, we approximate the target measure as a perturbation of an arbitrary reference measure along a few directions in $\mathbb{R}^{d}$. These directions are determined by minimizing an upper bound on the Kullback–Leibler (KL) divergence between the target and its approximation. Our contribution improves upon previous works by leveraging dimensional logarithmic Sobolev inequalities to refine the bound on the KL divergence. These inequalities lead to a uniformly tighter bound on the KL divergence, thereby enhancing the identification of the most significant perturbation directions. In particular, when the target and reference are both Gaussian, minimizing the resulting bound is equivalent to minimizing the KL divergence. We further demonstrate the applicability of this analysis to the squared Hellinger distance, where analogous reasoning shows that the dimensional Poincaré inequality offers improved bounds.

Bayesian inference↗
DOE OSTI2024

A Thermo‐Flow‐Mechanics‐Fracture Model Coupling a Phase‐Field Interface Approach and Thermo‐Fluid‐Structure Interaction

This work proposes a novel approach for coupling non-isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid-filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface-capturing phase-field fracture method and interface-tracking thermo-fluid-structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase-field approach to represent fractures and reconstruct the geometry to frame a thermo-fluid-structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase-field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid-filled fracture with the phase-field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase-field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.

fracture↗
DOE OSTI2020

Tunable correlated Chern insulator and ferromagnetism in a moiré superlattice

Studies of two-dimensional electron systems in a strong magnetic field revealed the quantum Hall effect, a topological state of matter featuring a finite Chern number C and chiral edge states. Haldane later theorized that Chern insulators with integer quantum Hall effects could appear in lattice models with complex hopping parameters even at zero magnetic field. The ABC-trilayer graphene/hexagonal boron nitride (ABC-TLG/hBN) moiré superlattice provides an attractive platform with which to explore Chern insulators because it features nearly flat moiré minibands with a valley-dependent, electrically tunable Chern number. In this work, we report the experimental observation of a correlated Chern insulator in an ABC-TLG/hBN moiré superlattice. We show that reversing the direction of the applied vertical electric field switches the moiré minibands of ABC-TLG/hBN between zero and finite Chern numbers, as revealed by large changes in magneto-transport behaviour. For topological hole minibands tuned to have a finite Chern number, we focus on quarter filling, corresponding to one hole per moiré unit cell. The Hall resistance is well quantized at h/2e 2 (where h is Planck’s constant and e is the charge on the electron), which implies C = 2, for a magnetic field exceeding 0.4 tesla. The correlated Chern insulator is ferromagnetic, exhibiting substantial magnetic hysteresis and a large anomalous Hall signal at zero magnetic field. Finally, our discovery of a C = 2 Chern insulator at zero magnetic field should open up opportunities for discovering correlated topological states, possibly with topological excitations, in nearly flat and topologically nontrivial moiré minibands.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Models for predicting objective function weights in prostate cancer IMRT

Purpose: To develop and evaluate the clinical applicability of advanced machine learning models that simultaneously predict multiple optimization objective function weights from patient geometry for intensity-modulated radiation therapy of prostate cancer. Methods: A previously developed inverse optimization method was applied retrospectively to determine optimal objective function weights for 315 treated patients. The authors used an overlap volume ratio (OV) of bladder and rectum for different PTV expansions and overlap volume histogram slopes (OVSR and OVSB for the rectum and bladder, respectively) as explanatory variables that quantify patient geometry. Using the optimal weights as ground truth, the authors trained and applied three prediction models: logistic regression (LR), multinomial logistic regression (MLR), and weighted K-nearest neighbor (KNN). The population average of the optimal objective function weights was also calculated. Results: The OV at 0.4 cm and OVSR at 0.1 cm features were found to be the most predictive of the weights. The authors observed comparable performance (i.e., no statistically significant difference) between LR, MLR, and KNN methodologies, with LR appearing to perform the best. All three machine learning models outperformed the population average by a statistically significant amount over a range of clinical metrics including bladder/rectum V53Gy, bladder/rectum V70Gy, and dose to the bladder, rectum, CTV, and PTV. When comparing the weights directly, the LR model predicted bladder and rectum weights that had, on average, a 73% and 74% relative improvement over the population average weights, respectively. The treatment plans resulting from the LR weights had, on average, a rectum V70Gy that was 35% closer to the clinical plan and a bladder V70Gy that was 29% closer, compared to the population average weights. Similar results were observed for all other clinical metrics. Conclusions: The authors demonstrated that the KNN and MLR weight prediction methodologies perform comparably to the LR model and can produce clinical quality treatment plans by simultaneously predicting multiple weights that capture trade-offs associated with sparing multiple OARs.

60 APPLIED LIFE SCIENCES↗