Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “observation geometry”

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 379 records · Page 21

Investigation of magnetic fluctuations in L-H and H-L transition dynamics on DIII-D

The dynamics of the L-H transition is not fully understood, with many parameters changing the threshold power to enter H-mode and the self-regulation between zonal flows and turbulence in the plasma edge. This paper is primarily a presentation of experimental results for DIII-D L-H and H-L transitions and speculation on the observations made. Power threshold analysis and measurements of pedestal temperatures for these transitions are presented. A comparison is made between an L-H transition and H-L transition of comparable Psep exhibiting oscillatory behaviour, showing symmetry between forward and backward transition dynamics. This paper shows the first observations of magnetic fluctuations during L-H and H-L transitions on DIII-D, and shows that L-H and H-L transitions have similar magnetic fluctuation dynamics. Information geometry analysis has been performed on measurements of plasma density fluctuations, perpendicular plasma velocity fluctuations, and magnetic field fluctuations to investigate the self-regulation and evolution of these variables during the transitions. Perpendicular flow evolution is shown to dominate the transition dynamics in both directions, but self-regulation behaviour is observed between all three variables. A strong correlation between magnetic fluctuation information rate and density fluctuation information rate for these two shots shows the strong influence of magnetic behaviour on both the L-H and H-L transition, and that these transition dynamics necessarily include electromagnetic effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Anatomy of an Unusual Edge-on Protoplanetary Disk. I. Dust Settling in a Cold Disk

As the earliest stage of planet formation, massive, optically thick, and gas-rich protoplanetary disks provide key insights into the physics of star and planet formation. When viewed edge-on, high-resolution images offer a unique opportunity to study both the radial and vertical structures of these disks and relate this to vertical settling, radial drift, grain growth, and changes in the midplane temperatures. In this work, we present multi-epoch Hubble Space Telescope and Keck scattered light images, and an Atacama Large Millimeter/submillimeter Array 1.3 mm continuum map for the remarkably flat edge-on protoplanetary disk SSTC2DJ163131.2–242627, a young solar-type star in ρ Ophiuchus. We model the 0.8 μm and 1.3 mm images in separate Markov Chain Monte Carlo (MCMC) runs to investigate the geometry and dust properties of the disk using the MCFOST radiative transfer code. In scattered light, we are sensitive to the smaller dust grains in the surface layers of the disk, while the submillimeter dust continuum observations probe larger grains closer to the disk midplane. An MCMC run combining both data sets using a covariance-based log-likelihood estimation was marginally successful, implying insufficient complexity in our disk model. The disk is well characterized by a flared disk model with an exponentially tapered outer edge viewed nearly edge-on, though some degree of dust settling is required to reproduce the vertically thin profile and lack of apparent flaring. A colder than expected disk midplane, evidence for dust settling, and residual radial substructures all point to a more complex radial density profile to be probed with future, higher-resolution observations.

47 OTHER INSTRUMENTATION↗

Spacecraft Geometry Effects on Kinetic Impactor Missions

Abstract The DART (Double Asteroid Redirection Test) mission will impact a spacecraft on the secondary (Dimorphos) of the binary asteroid system Didymos in 2022 September, with the goal of altering the orbital period of Dimorphos about Didymos sufficiently to be observed from ground-based observations. Numerical impact modeling is a crucial component in understanding the outcome of the DART experiment, and while many have investigated the effects of target properties, such as material strength and porosity (which remain unknown), an often overlooked factor is the importance of accurately representing the spacecraft itself in such models. Most impact modeling to date has considered simple impactor geometries such as a solid uniform sphere, but in reality the spacecraft is a complex shape full of different components, open spaces, and thin walled structures. At a minimum, a simple solid representation underestimates the surface area of the impact: for a small body such as Dimorphos (approximately 160 m in diameter), the difference between a spacecraft spanning 20 m (including solar arrays) impacting and a sub-1 m idealized shape may be important. In this paper, we compare models impacting high-fidelity models of the spacecraft based on the CAD geometry with various simplified impactors, in order to assess the potential importance of this effect. We find that the difference between the simplest impactor geometries (such as a uniform sphere) and the real spacecraft is measurable, and has an interesting dependence on the material properties of the asteroid itself.

36 MATERIALS SCIENCE↗

Drift kinetic electrostatic simulations of the edge localized mode heat pulse

In the present work, electrostatic drift kinetic simulations of parallel plasma transport within the tokamak scrape-off layer (SOL) are conducted using the COGENT code. The SOL configuration is represented in one-dimensional slab geometry, incorporating a heat source localized in the midplane. The heat source parameters correspond to those characterizing edge-localized modes observed in the Joint European Torus (JET) tokamak. The numerical model includes kinetic treatment of both ions and electrons, a simplified model for the gyrokinetic Poisson equation that allows one to step over short time scales associated with fast electrostatic shear Alfvèn waves, and the logical sheath boundary condition (LSBC) that enforces global system quasineutrality. A third-order accurate LSBC is derived to be consistent with the third-order accurate upwind advection scheme utilized in the code, and it was shown to noticeably impact the simulation results, especially parallel heat flux at the target plate. The findings of this study are in agreement with results from preceding fluid and kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards Multilayered Coatings of 304L Stainless Steels Using Friction Surfacing

Abstract The objective of this work is to study friction surfacing process variability when depositing multilayered coatings. This is motivated by the need to maintain deposition quality when depositing multiple friction surfacing layers, whether for repair, remanufacturing, or new part creation using this solid-state metal additive manufacturing process. In this study, 10-mm-diameter 304L stainless steel rods were used to create up to five layers of 40-mm-long coatings on 304L substrates using a constant set of processing parameters. In-process measurement of forces (X, Y, Z), flash temperature, flash geometry, layer temperature, and post-process measurement of layer geometry, microhardness, and microstructure are used to characterize changes in the friction surfacing process as more layers are deposited. It was observed that with increasing layers: layer thickness and deposition efficiency decrease; offsetting of the deposition towards the retreating side, and temperature in the deposited layer increase; and flash temperature does not change. Metallurgical analyses of friction-surfaced cross-sections revealed fine grain refinement and transformation of base austenite to strain-induced martensite. It is concluded that the process parameters need to be adjusted even after the second or third layer is deposited, corrections to the tool path are required after a couple of layers, and the measured process forces, as well as deposited layer temperature, may be useful to monitor and control the process and its instabilities.

Engineering↗

A Parallel Cut-Cell Algorithm for the Free-Boundary Grad--Shafranov Problem

A parallel cut-cell algorithm is described to solve the free-boundary problem of the Grad--Shafranov equation. The algorithm reformulates the free-boundary problem in an irregular bounded domain and its important aspects include a searching algorithm for the magnetic axis and separatrix, a surface integral along the irregular boundary to determine the boundary values, an approach to optimize the coil current based on a targeting plasma shape, Picard iterations with Aitken's acceleration for the resulting nonlinear problem, and a Cartesian grid embedded boundary method to handle the complex geometry. Here the algorithm is implemented in parallel using a standard domain-decomposition approach and a good parallel scaling is observed. Numerical results verify the accuracy and efficiency of the free-boundary Grad--Shafranov solver.

97 MATHEMATICS AND COMPUTING↗

Algorithm 1049: The Delaunay Density Diagnostic

Accurate approximation of a real-valued function depends on two aspects of the available data: the density of inputs within the domain of interest and the variation of the outputs over that domain. There are few methods for assessing whether the density of inputs is sufficient to identify the relevant variations in outputs—i.e., the “geometric scale” of the function—despite the fact that sampling density is closely tied to the success or failure of an approximation method. In this article, we introduce a general purpose, computational approach to detecting the geometric scale of real-valued functions over a fixed domain using a deterministic interpolation technique from computational geometry. The algorithm is intended to work on scalar data in moderate dimensions (2–10). Our algorithm is based on the observation that a sequence of piecewise linear interpolants will converge to a continuous function at a quadratic rate (in L 2 norm) if and only if the data are sampled densely enough to distinguish the feature from noise (assuming sufficiently regular sampling). We present numerical experiments demonstrating how our method can identify feature scale, estimate uncertainty in feature scale, and assess the sampling density for fixed (i.e., static) datasets of input–output pairs. Finally, we include analytical results in support of our numerical findings and have released lightweight code that can be adapted for use in a variety of data science settings.

97 MATHEMATICS AND COMPUTING↗

Field Results from New Tensor Borehole Optical Fiber Strainmeter Installations in Oklahoma and Utah

The time evolving strain field contains a wealth of information that can be used to interpret subsurface behavior. For example, injecting or removing fluids from reservoirs or aquifers causes deformation that can be used as a diagnostic signal in some cases, while it can interfere with geodetic interpretations in other cases. We've previously demonstrated the feasibility of measuring the strain tensor at a depth of 30m caused by injection into a reservoir at 530m. The observed strain signals were interpreted using four independent analytic and numerical methods that resulted in estimates of the poroelastic properties and geometry of the reservoir that was consistent with data from well logs. However, studies like these are only possible if these deformations can be reliably measured. Years of lab and field work has culminated in the development of a novel borehole strainmeter capable of resolving multiple components of strain using embedded optical fibers configured as Michelson interferometers. It features four horizontal gauges separated by 45° to resolve the horizontal strain tensor as well as a vertical strain gauge and a sixth null component for state-of-health monitoring. The downhole sensing package also includes an open pipe through its center for grout circulation during single-trip deployments and a fully welded stainless steel exterior for robustness and longevity. These instruments have a resolution of 2x10-13 strain that can easily measure the solid earth tides. Preliminary data are available from four strainmeters in shale at our Oklahoma site and four in compacted sand and gravel in Utah. These are deployed from 40-60m, except one of the strainmeters in Oklahoma is deployed at 500m. The data include strains from the initial grout curing, comparisons to predicted earth tide models and in-situ calibration results, barometric pressure admittances and spectral analyses as well as signals from underground injections and surface waves from teleseismic events. Preliminary analyses indicate behavior consistent with other strainmeter deployments, and comparison to data from a Gladwin strainmeter at the Oklahoma site validate the performance of the new design. Analyses from a suite of six well tests at the Oklahoma site show for the first time how the strain tensor field varies with location during well testing.

DeWolf, Scott↗

Study of Magnetic Field and Turbulence in the TeV Halo around the Monogem Pulsar

Magnetic fields are ubiquitous in the interstellar medium, including extended objects such as supernova remnants and diffuse halos around pulsars. Its turbulent characteristics govern the diffusion of cosmic rays and the multiwavelength emission from pulsar wind nebulae (PWNe). However, the geometry and turbulence nature of the magnetic fields in the ambient region of PWN is still unknown. Recent gamma-ray observations from HAWC and synchrotron observations suggest a highly suppressed diffusion coefficient compared to the mean interstellar value. In this study, we present the first direct observational evidence of the orientation of the mean magnetic field and turbulent characteristics by employing a recently developed statistical parameter "Y turb " in the extended halo around the Monogem pulsar. Our study points to two possible scenarios: nearly aligned toward the line of sight (LOS) with compressible modes dominance or high inclination angle toward the LOS and characterized by Alfvénic turbulence. The first scenario appears consistent with other observational signatures. Furthermore, we report that the magnetic field has an observed correlation length of approximately 3 ± 0.6 pc in the Monogem halo. Our study highlights the pivotal role of magnetic field and turbulence in unraveling the physical processes in TeV halos and cosmic-ray transport.

79 ASTRONOMY AND ASTROPHYSICS↗

A mathematical model of solid-state dewetting of barium thin films on W (112)

Solid state dewetting occurs when a thin solid film is heated. The temperature of dewetting depends on the thickness of the film; dewetting can be observed in the range of 1/3 to 1/2 of the bulk melting temperature. While remaining solid, the film behaves in a manner similar to liquids dewetting and agglomerating to forming islands or droplets. One of the possible mechanisms is the conversion of a metastable thin film geometry into a more stable form. Heating the metastable film gives the film atoms higher mobility, and the films retract, dewetting the surface. This atomic motion can be restricted due to surface anisotropy. Herein we present in situ emission microscopy observations of barium thin films deposited onto W (112) by thermal evaporation. From the modeling viewpoint, the evolution of the film in this system could be divided in four stages: (i) the nucleation and growth of the thin film as a simply connected region; (ii) formation of droplets/islands/hillocks; (iii) nucleation of holes; (iv) evolution of the components of the disconnected film to their equilibrium state. We present a continuum model that is qualitatively consistent with the evolution of the film observed at the initial stage of the experiment and discuss the later stages of the evolution of surface structures.

36 MATERIALS SCIENCE↗

Excessive charge, beam loading, and impedance collapse thresholds for a velvet emitter

Excess emission has been observed from velvet cathodes with total emission times >150 ns. In this diode experiment, we produce a sub-relativistic electron beam with the ability to consistently change γ from 1.2 to 1.5 and β = 0.5–0.75. Electron emission in this particular diode geometry requires electric fields >40 kV/cm. The current increases at steady rates >0.05 A/ns after the head of the pulse, indicating an expansion of the emission surface and reduction in the effective AK gap. Small transients of excess emission (or arcs) are consistently observed for current pulses exceeding 150 ns. The excess emission results in beam loaded levels ≥10 kV on the diode voltage. The effects described here are compounded as diode voltage is increased. The principal objective of these experiments is to quantify electric field emission thresholds, current ramps, excess emission delays, effective beam loading, and impedance collapse.

43 PARTICLE ACCELERATORS↗

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗

E × B flow driven electron temperature bifurcation in a closed slot divertor with ion B × ∇ B away from the X-point in the DIII-D tokamak

An electron temperature bifurcation is observed in the small angle slot divertor, which has been developed to enhance neutral cooling across the divertor target by coupling a closed slot structure with appropriate target shaping. Experiments in the DIII-D tokamak and associated SOLPS-ITER modeling with full drifts find a strong interplay between drifts and divertor geometry on divertor dissipation. The coupling of divertor geometry and drift flows can strongly affect the path towards divertor detachment onset as the plasma density is raised. With the strike point on the inner slanted surface and ion B × ∇B away from the magnetic X-point, bifurcative transitions were observed with sharp decrease of T e towards detachment onset both experimentally and computationally. This differs from the situation for the open divertor where the T e cliff was only observed for ion B × ∇B towards the X-point. SOLPS-ITER modeling with full drifts demonstrates that the magnitude of the E × B drift flow is comparable with the main plasma flow. The reversal of both the poloidal and radial E × B flows near the strike point leads to rapid density accumulation right near the separatrix, which results in bifurcative step transition of divertor conditions with cold plasma across the entire divertor target plate. Furthermore, these results indicate that the interplay between geometry and drifts should be fully taken into account in future fusion reactor divertor designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Behavior of thermally quenched topological defects in quasicrystal artificial spin ices

We report that understanding the formation of topological defects in connected magnetic nanowires such as in artificial spin ices (ASIs) is critical to controlling their behavior. Moreover, emergent frustration behavior of these defects can be expected due to complex underlying lattice geometry of the ASIs. We have explored in situ thermalization and quenching of quasicrystalline ASIs, which exhibit additional frustration due to aperiodicity from near the Curie temperature, and we observe the formation of two sets of topological vortex defects: within the magnetic bars and within the vertices at which the bars meet. By varying the laser fluence, we show that the number of defects in the vertices varies as the quenching rate changes following a power law relation, which we interpret as being related to the ferromagnetic-to-paramagnetic phase transition. Conversely, the defects formed within the bars are metastable states related to the spontaneous magnetization reversal of the magnetic bars and are therefore strongly dependent on underlying lattice geometry. These results provide insight into the emergence and control of topological defects in confined frustrated magnetic systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantitative study of the effect of grain boundary parameters on the slip system level Hall-Petch slope for basal slip system in Mg-4Al

Several theoretical studies have reported that the geometry and structure of grain boundaries in polycrystalline materials could impose a significant effect on the Hall-Petch slope. However, experimental observations are primarily limited by the ability of the techniques to accurately quantify the grain boundary strength and validate these theoretical models. Using high-resolution electron backscatter diffraction (HR-EBSD), the local stress tensor ahead of a slip band blocked by a grain boundary was quantified and coupled with a continuum dislocation pile-up model to assess the barrier strength of specific grain boundaries to specific slip systems, referred to as micro-Hall-Petch coefficient. For basal slip system in a deformed Mg-4Al alloy, the micro-Hall-Petch coefficient ($k^{basal}_{μ}$) varied significantly, from 0.054 to 0.184 MPa - m 1/2 for nine different grain boundaries. These results were correlated with geometric descriptors of the respective grain boundaries, with three-dimensional GB profile additionally measured via focused ion beam milling. It was found that the angle between the two slip plane traces on the grain boundary plane was the most sensitive parameter affecting $k^{basal}_{μ}$, followed by the angle between the slip directions. A functional form for calculation of $k^{basal}_{μ}$ depending on these two angles is proposed to augment crystal plasticity constitutive models with slip resistance dependent on some measure of the grain size. The method allows a new pathway to calibrate grain size strengthening parameters in crystal plasticity models, allowing further computational investigations of the interrelationship between texture, grain morphology, and the Hall Petch effect.

36 MATERIALS SCIENCE↗

Contribution of the catalytic dyad of SARS-CoV-2 main protease to binding covalent and noncovalent inhibitors

The effect of mutations of the catalytic dyad residues of SARS-CoV-2 main protease (MPro WT ) on the thermodynamics of binding of covalent inhibitors comprising nitrile [nirmatrelvir (NMV), NBH2], aldehyde (GC373), and ketone (BBH1) warheads to MPro is examined together with room temperature X-ray crystallography. When lacking the nucleophilic C145, NMV binding is ~400-fold weaker corresponding to 3.5 kcal/mol and 13.3 °C decrease in free energy (ΔG) and thermal stability (T m ), respectively, relative to MPro WT . The H41A mutation results in a 20-fold increase in the dissociation constant (K d ), and 1.7 kcal/mol and 1.4 °C decreases in ΔG and T m , respectively. Increasing the pH from 7.2 to 8.2 enhances NMV binding to MPro H41A , whereas no significant change is observed in binding to MPro WT . Structures of the four inhibitor complexes with MPro 1-304/C145A show that the active site geometries of the complexes are nearly identical to that of MPro WT with the nucleophilic sulfur of C145 positioned to react with the nitrile or the carbonyl carbon. These results support a two-step mechanism for the formation of the covalent complex involving an initial non-covalent binding followed by a nucleophilic attack by the thiolate anion of C145 on the warhead carbon. Noncovalent inhibitor ensitrelvir (ESV) exhibits a binding affinity to MPro WT that is similar to NMV but differs in its thermodynamic signature from NMV. The binding of ESV to MPro C145A also results in a significant, but smaller, increase in K d and decrease in ΔG and T m , relative to NMV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigating property-porosity relationships for micro-architected lattice structures

Micro-architected structures are increasingly valued for their light weight and tunable mechanical properties; this class of material includes sheet structures like the Schöen gyroid and beam lattices like the octet. For design purposes, it is critical to understand how to tune the relative density (RD) to obtain desired mechanical properties. This study investigates the mechanical response of Ti-6Al-4V gyroid structures, spanning a broad range of RDs (0.03–0.90), unit cell sizes (1–4 cm), and sheet thicknesses (0.2–7.8 mm). Results demonstrate that the classical Gibson-Ashby power law scaling between RD and modulus and yield stress does not adequately capture the response over a wide range of RDs, nor does it extrapolate correctly to the fully dense solid. Analytical models, in concert with experimental results and high-fidelity finite element calculations, show that the deviation from Gibson-Ashby reflects a transition from structure-dominated to material-dominated behavior. Distinctions are drawn between scaling relationships and property-porosity models, and an analytical model is proposed to better capture the evolution of mechanical properties with RD. These observations are mirrored in other architected topologies, like the octet, and highlight the importance of understanding the relationship between mechanical properties and geometry for design purposes.

Materials science↗