Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational studies”

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,675 records · Page 93

Photoenzymatic Csp 3 –Csp 3 bond formation via enzyme-templated radical–radical coupling

Cross-couplings are essential reactions in modern chemical synthesis, enabling the rapid construction of complex molecules from simple precursors. Transition metal catalysts are prized for these transformations because their reactivity and selectivity can be tuned via judicious selection of the metal and ligand. Although enzymes offer analogous opportunities for tuning via protein engineering, their application to cross-coupling remains limited, as nature relies on alternative paradigms for building molecular complexity. Here, we report the cross-coupling of alkyl halides and benzylic carboxylic acids using an engineered flavin-dependent lactate monooxygenase—a photoenzyme. The enzyme achieves this feat by exploiting the redox versatility of the flavin cofactor. Stoichiometric experiments, ultrafast spectroscopy, and computational studies support a mechanism in which photoexcited flavin quinone initiates the reaction via oxidative decarboxylation to generate a benzylic radical. The resulting flavin semiquinone can reduce the alkyl halide to form a second organic radical within the protein active site, which rapidly engages in C(sp 3 )–C(sp 3 ) bond formation. A variant was engineered to control the stereochemical outcome of this radical–radical coupling event, highlighting the ability of the protein to alter the energetic barrier for a mechanistic step that is traditionally understood to be near barrierless. This work demonstrates that the scope for nonnative reaction mechanisms in biocatalysis far exceeds previously established bounds and has potential to solve a variety of reactivity challenges in cross-coupling chemistry.

biocatalysis↗

Synthesis and investigation into explosive sensitivity for a series of new picramide explosives

Tailoring the molecular properties that govern energetic material sensitivity is essential to improve safety and help develop new energetic materials. Despite this need, understanding the complex chemistry and physics of explosive initiation and propagation is still a challenge. Recent work by our group has reinforced the view that explosive sensitivity under sub-shock conditions is connected to the strength of the weakest covalent bond in the molecule, that is, its “trigger linkage.” These correlations have been observed with different classes of energetic molecules and indicate that “trigger linkage” bond breaking, and heat of explosion are good indicators for the sensitivity trends. Herein we report the synthesis of aliphatic energetic materials with ethane, propane and neopentane backbones. Experimental and computational studies show that the trigger linkage model, based on results from quantum molecular dynamics simulations, correctly predicts trends observed in the impact sensitivity of the molecules. However, while the model predicts the impact sensitivities of the ethane series, the neopentane series has higher impact sensitivities than predicted, which is presumably influenced by crystal packing effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CHARACTERIZING THE INFRARED SPECTRA OF SMALL, NEUTRAL, FULLY DEHYDROGENATED POLYCYCLIC AROMATIC HYDROCARBONS

We present the results of a computational study to investigate the infrared spectroscopic properties of a large number of polycyclic aromatic hydrocarbon (PAH) molecules and their fully dehydrogenated counterparts. We constructed a database of fully optimized geometries for PAHs that is complete for eight or fewer fused benzene rings, thus containing 1550 PAHs and 805 fully dehydrogenated aromatics. A large fraction of the species in our database have clearly non-planar or curved geometries. For each species, we determined the frequencies and intensities of their normal modes using density functional theory calculations. Whereas most PAH spectra are fairly similar, the spectra of fully dehydrogenated aromatics are much more diverse. Nevertheless, these fully dehydrogenated species show characteristic emission features at 5.2 μm, 5.5 μm, and 10.6 μm; at longer wavelengths, there is a forest of emission features in the 16-30 μm range that appears as a structured continuum, but with a clear peak centered around 19 μm. We searched for these features in Spitzer-IRS spectra of various positions in the reflection nebula NGC 7023. We find a weak emission feature at 10.68 μm in all positions except that closest to the central star. We also find evidence for a weak 19 μm feature at all positions that is not likely due to C{sub 60}. We interpret these features as tentative evidence for the presence of a small population of fully dehydrogenated PAHs, and discuss our results in the framework of PAH photolysis and the formation of fullerenes.

79 ASTRONOMY AND ASTROPHYSICS↗

Compression mechanisms of ferroelectric PbTiO 3 via high pressure neutron scattering

Switchable atomic displacements generate electric dipole moments in ferroelectric materials utilized in many contemporary devices. Lead titanate, a perovskite oxide with formula PbTiO 3 , has been referred to as a textbook example of a prototype displacive ferroelectric and is a testing platform of widely used models of piezoelectric response of complex solid-solutions. PbTiO 3 has been addressed by experimental and computational studies, often with apparently conflicting conclusions. To date, hydrostatic pressure experiments have been interpreted in terms of a model in which the dipole moments gradually diminish with increasing pressure until a transition to a cubic phase, characterized by a zero average dipole moment, occurs. Additionally, the model unrealistically assumes an even compression of the crystal. Here we show by high-pressure neutron powder diffraction measurements that a fast and slow shrinkage of 12-oxygen cages around Pb and octahedra around Ti, respectively, takes place. A phase diagram consolidating earlier and present results is given.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dynamic focusing of laser driven positron jets by self-generated fields

Abstract Focusing effect of laser-driven positron jets by self-generated target sheath fields has been observed for the first time experimentally and the results are supported by the computational studies. In the experiment, OMEGA EP short-pulse (0.7 ps, 500 J) irradiates mm-size gold targets with a concave back surface and reference flat-surface targets. Both targets exhibited positrons with quasi-monoenergetic energy peaks while targets with concave curvature also showed increased number of positrons at the detector. The data is consistent with hybrid-PIC simulations confirming that the time-varying electric fields driven by electrons escaping from the target significantly change the trajectories of positrons. These simulations show a small radius of curvature on the rear side increases the relative focusing effect and the positrons to electrons ratio in the escaping plasma. For the smallest radius of curvature, positron jets that are up to 10 times denser can be achieved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dependence of alpha-particle-driven Alfvén eigenmode linear stability on device magnetic field strength and consequences for next-generation tokamaks

Recently-proposed tokamak concepts use magnetic fields up to 12 T, far higher than in conventional devices, to reduce size and cost. Theoretical and computational study of trends in plasma behavior with increasing field strength is critical to such proposed devices. This paper considers trends in Alfvén eigenmode (AE) stability. Energetic particles, including alphas from deuterium–tritium fusion, can destabilize AEs, possibly causing loss of alpha heat and damage to the device. AEs are sensitive to device magnetic field via the field dependence of resonances, alpha particle beta, and alpha orbit width. We describe the origin and effect of these dependences analytically and by using recently-developed numerical techniques (Rodrigues et al 2015 Nucl. Fusion 55 083003). The work suggests high-field machines where fusion-born alphas are sub-Alfvénic or nearly sub-Alfvénic may partially cut off AE resonances, reducing growth rates of AEs and the energy of alphas interacting with them. High-field burning plasma regimes have non-negligible alpha particle beta and AE drive, but faster slowing down time, provided by high electron density, and higher field strength reduces this drive relative to low-field machines with similar power densities. The toroidal mode number of the most unstable modes will tend to be higher in high magnetic field devices. Finally, the work suggests that high magnetic field devices have unique, and potentially advantageous, AE instability properties at both low and high densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A novel measurement of marginal Alfvén eigenmode stability during high power auxiliary heating in JET

The interaction of Alfvén eigenmodes (AEs) and energetic particles is one of many important factors determining the success of future tokamaks. In JET, eight in-vessel antennas were installed to actively probe stable AEs with frequencies ranging 25–250 kHz and toroidal mode numbers |n| < 20. During the 2019–2020 deuterium campaign, almost 7500 resonances and their frequencies f 0 , net damping rates γ < 0, and toroidal mode numbers were measured in almost 800 plasma discharges. From a statistical analysis of this database, continuum and radiative damping are inferred to increase with edge safety factor, edge magnetic shear, and when including non-ideal effects. Both stable AE observations and their associated damping rates are found to decrease with |n|. Active antenna excitation is also found to be ineffective in H-mode as opposed to L-mode; this is likely due to the increased edge density gradient’s effect on accessibility and ELM-related noise’s impact on mode identification. A novel measurement is reported of a marginally stable, edge-localized ellipticity-induced AE probed by the antennas during high-power auxiliary heating (ion cyclotron resonance heating and neutral beam injection) up to 25 MW. NOVA-K kinetic-MHD simulations show good agreement with experimental measurements of f 0 , γ, and n, indicating the dominance of continuum and electron Landau damping in this case. Similar experimental and computational studies are planned for the recent hydrogen and ongoing tritium campaigns, in preparation for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the generation of horizontal shear waves by underground explosions in jointed rocks

This paper describes computational studies of non-spherical ground motions generated by spherical explosions in a heavily jointed granite formation. Various factors affecting the shear wave generation are considered, including joint spacing, orientation, persistence and properties. Simulations are performed both in 2D for a single joint set to elucidate the basic response mechanisms, and in 3D for multiple joint sets to realistically represent in situ conditions in a realistic geologic setting. The joints are modeled explicitly using both contact elements and weakness planes in the material. Simulations are performed both deterministically and stochastically to quantify the effects of geologic uncertainties on near field ground motions. The mechanical properties of the rock and the joints as well as the joint spacing and orientation are taken from experimental test data and geophysical logs corresponding to the Climax Stock granitic outcrop, which is the geologic setting of the Source Physics Experiment (SPE). Agreement between simulation results and near field wave motion data from SPE enables newfound understanding of the origin and extent of non-spherical motions associated with underground explosions in fractured geologic media.

58 GEOSCIENCES↗

Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂⁡(10 2 ) ⁢GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂⁡(10 0 ) ⁢GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Monte Carlo methods↗

Modeling phase transformations in Mn-rich disordered rocksalt cathodes with machine-learning interatomic potentials

Mn-rich disordered rocksalt (DRX) cathode materials exhibit a phase transformation from a disordered to a partially disordered spinel-like structure (δ-phase) during electrochemical cycling. Here, in this computational study, we use charge-informed molecular dynamics with a fine-tuned CHGNet foundation potential to investigate the phase transformation in LixMn 0.8 Ti 0.1 O 1.9 F 0.1 . Our results indicate that transition metal migration occurs and reorders to form the spinel-like ordering in an FCC anion framework. The transformed structure contains a higher concentration of nontransition metal (0-TM) face-sharing channels, which are known to improve Li transport kinetics. Analysis of the Mn valence distribution suggests that the appearance of tetrahedral Mn 2+ is a consequence of spinel-like ordering, rather than the trigger for cation migration as previously suggested. Calculated equilibrium intercalation voltage profiles demonstrate that the δ-phase, unlike the ordered spinel, exhibits solid-solution signatures at low voltage. A higher Li capacity is obtained than in the DRX phase. This study provides atomic insights into solid-state phase transformation and its relation to experimental electrochemistry, highlighting the potential of machine-learning interatomic potentials for understanding complex oxide materials.

Zhong, Peichen [University of California, Berkeley↗

Evidence of topological charge polarization at active-passive interfaces in acoustically powered active liquid crystals

Theoretical and computational studies predict topological charge separation at active-passive interfaces in active nematics, but reliable experimental validation has been lacking. We utilize a synthetic acoustically energized liquid crystal to experimentally investigate the dynamics of topological defects across an active-passive interface. The interference pattern induced by the acoustic wave inside the experimental cell creates a spatial distribution of activity, resulting in the formation of active-passive interfaces with the liquid-crystalline director field aligned parallel to those interfaces. The activity gradient drives the reorientation of positively charged topological defects along its direction, causing them to migrate into the passive zone and form a layer with net topological charge. That layer “discharges” upon cessation of activity through annihilation of positive and negative topological charges, resembling the discharge of a capacitor.

Active defects↗

Atomistic three-dimensional coherent x-ray imaging of nonbiological systems

We computationally study the resolution limits for three-dimensional coherent x-ray diffractive imaging of heavy, nonbiological systems using Ar clusters as a prototype. We treat electronic and nuclear dynamics on an equal footing and remove the frozen-lattice approximation often used in electronic damage studies. We explore the achievable resolution as a function of pulse parameters (fluence level, pulse duration, and photon energy) and particle size. The contribution of combined lattice and electron dynamics is not negligible even for 2 fs pulses, and the Compton scattering is less deleterious than in biological systems for atomic-scale imaging. Although free-electron scattering represents a significant background, we find that recovery of the original structure is in principle possible with 3 angstrom resolution for particles of 11 nm diameter.

47 OTHER INSTRUMENTATION↗

Topological surface states in strained Dirac semimetal thin films

We computationally study the Fermi arc states in a Dirac semimetal, both in a semi-infinite slab and in the thin-film limit. We use Cd 3 As 2 as a model system, and include perturbations that break the C-4 symmetry and inversion symmetry while mirror symmetries remain. The surface states are protected by the mirror symmetries present in the bulk states and thus survive these perturbations. Additionally, the Fermi arc states persist down to very thin films, thinner than presently measured experimentally, but are affected by breaking the symmetry of the Hamiltonian. Recent experimental observations of transport in Cd 3 As 2 films strongly suggest transport by surface states. Our findings are consistent with those observations in that we show that Fermi arc states can survive to rather thin films, and our results also suggest that symmetry-breaking terms that preserve the Fermi arc states nevertheless can have a profound effect in the thin film limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evolution of surface structure in laser-preheated, perturbed materials

Here, we report an experimental and computational study investigating the effects of laser preheat on the hydrodynamic behavior of a material layer. In particular, we find that perturbation of the surface of the layer results in a complex interaction, in which the bulk of the layer develops density, pressure, and temperature structure and in which the surface experiences instability-like behavior, including mode coupling. A uniform one-temperature preheat model is used to reproduce the experimentally observed behavior, and we find that this model can be used to capture the evolution of the layer, while also providing evidence of complexities in the preheat behavior. Lastly, this result has important consequences for inertially confined fusion plasmas, which can be difficult to diagnose in detail, as well as for laser hydrodynamics experiments, which generally depend on assumptions about initial conditions in order to interpret their results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three-dimensional simulations of reshocked inclined Richtmyer-Meshkov instability: Effects of initial perturbations

The effect of initial perturbations on the evolution of the inclined Richtmyer-Meshkov turbulent mixing layer before and after reshock initiated by a shock wave with Mach number 1.55 is investigated through three-dimensional (3D) simulations using the flash code. The 3D simulations aim to reproduce both predominantly single-mode and multimode interfaces between light and heavy gases (N 2 -CO 2 , Atwood number, A≈0.22; amplitude to wavelength ratio of 0.088) which were created in an inclined shock tube facility to analyze the effects of initial conditions on mixing development in the entire flow field. The two-dimensional center slices of 3D simulations are compared with the experimental results to validate the computational code. Mixing width, mixed mass, mixed-mass thickness, and circulation in addition to concentration fields are shown to be in good agreement with the experimental data. The three-dimensional density and vorticity fields are first presented to qualitatively describe the flow behavior before and after reshock. Several measured density/velocity-related quantities indicate that the growth of the mixing material is strongly dependent on initial conditions. Before reshock and at early times after reshock, flow is clearly maintaining the memory of initial perturbations. However, at late time after reshock, although the large wavelength feature still dominates the flow motion, and the morphology of the two different interfaces indicates several differences, by breakdown of large-scale coherent structures to much finer scales, the memory of small scales of the multimode initial perturbation is not as clear as pre-reshock. Regarding three-dimensionality of the flow, before reshock in the multimode case, the baroclinic vorticity production, circulation, turbulent kinetic energy, and turbulent mass flux suggest that the small-scale roll-up features along the large inclined wavelength quickly evolves in all three dimensions. The coherent vortex tubes break down to smaller wormlike vortex structures, and turbulent fluctuations in the out-of-plane dimension are comparable to the spanwise direction. After reshock, this three-dimensionality of mixing growth was observed in the flow for both initial conditions. The results of this work represent a significant extension of previous computational studies performed on this specific topic. A different code with a different numerical method is validated through comparison with the experimental data. The initial perturbations are directly measured from the experimental results. Moreover, the entire three-dimensional experimental shock tube domain is simulated, and more quantities are investigated to understand the mixing mechanism and instability evolution in all three dimensions.

42 ENGINEERING↗

Controlling Rayleigh-Taylor Instabilities in Magnetically Driven Solid Metal Shells by Means of a Dynamic Screw Pinch

Magnetically driven implosions of solid metal shells are an effective vehicle to compress materials to extreme pressures and densities. Rayleigh-Taylor instabilities (RTI) are ubiquitous, yet typically undesired features in all such experiments where solid materials are rapidly accelerated to high velocities. In cylindrical shells (“liners”), the magnetic field driving the implosion can exacerbate the RTI. Here, we suggest an approach to implode solid metal liners enabling a remarkable reduction in the growth of magnetized RTI (MRTI) by employing a magnetic drive with a tilted, dynamic polarization, forming a dynamic screw pinch. Our calculations, based on a self-consistent analytic framework, demonstrate that the cumulative growth of the most deleterious MRTI modes may be reduced by as much as 1 to 2 orders of magnitude. One key application of this technique is to generate increasingly stable, higher-performance implosions of solid metal liners to achieve fusion [M. R. Gomez et al., Phys. Rev. Lett. 113, 155003 (2014)]. Finally, we weigh the potentially dramatic benefits of the solid liner dynamic screw pinch against the experimental tradeoffs required to achieve the desired drive field history and identify promising designs for future experimental and computational studies.

36 MATERIALS SCIENCE↗

Effect of Chain Length Dispersity on the Mobility of Entangled Polymers

While nearly all theoretical and computational studies of entangled polymer melts have focused on uniform samples, polymer synthesis routes always result in some dispersity, albeit narrow, of distribution of molecular weights ( Ð M = M w / M n ~ 1.02 – 1.04 ). Here, the effects of dispersity on chain mobility are studied for entangled, disperse melts using a coarse-grained model for polyethylene. Polymer melts with chain lengths set to follow a Schulz-Zimm distribution for the same average M w = 36 kg / mol with Ð M = 1.0 to 1.16, were studied for times of 600 – 800 μ s using molecular dynamics simulations. This time frame is longer than the time required to reach the diffusive regime. We find that dispersity in this range does not affect the entanglement time or tube diameter. However, while there is negligible difference in the average mobility of chains for the uniform distribution Ð M = 1.0 and Ð M = 1.02 , the shortest chains move significantly faster than the longest ones offering a constraint release pathway for the melts for larger Ð M .

36 MATERIALS SCIENCE↗