Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient Limiters”

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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Experimental demonstration of externally driven millimeter-wave particle accelerator structure

We report the experimental demonstration of a mm-wave electron accelerating structure powered by a high-power rf source. We demonstrate reliable coupling of an unprecedented rf power—up to 575kW into the mm-wave accelerator structure using a quasi-optical setup. This standing wave accelerating structure consists of a single-cell copper cavity and a Gaussian to TM 01 mode converter. The accelerator structure is powered by 110 GHz, 10-ns long rf pulses. These pulses are chopped from 3 microsecond pulses from a gyrotron oscillator using a laser-driven silicon switch. We show an unprecedented high gradient up to 230MV/m that corresponds to a peak surface electric field of more than 520 MV/m. We have achieved these results after conditioning the cavity with more than 100,000 pulses. We also report preliminary measurements of rf breakdown rates, which are important for understanding rf breakdown physics in the millimeter-wave regime. Furthermore, these results open up many frontiers for applications not only limited to the next generation particle accelerators but also x-ray generation, probing material dynamics, and nonlinear light-matter interactions at mm-wave frequency.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Anti-Ultralocality and Plateau Models of Inflation

Anti-ultralocality refers to the growth of spatial gradient terms relative to velocity terms in the coupled Einstein--scalar field equations. It is a characteristic feature of decelerated expansion before the onset of inflation. Previous numerical relativity studies have shown that anti-ultralocality prevents the onset of inflation in models with power-law inflaton potentials. In this paper, we show that models with plateau-shaped inflaton potentials, which are considered to be the simplest way to generate a tensor-to-scalar ratio below current observational upper limits, are especially vulnerable to anti-ultralocality effects. The reasons are the flatness of the plateau and the energy density gap of $\sim 10$ orders of magnitude between the Planck density and the plateau potential energy. To study the problem, we develop a protocol for assessing the viability of inflationary models in general, and we apply it to a plateau potential using a previously validated numerical relativity code. We find that, starting from generic initial conditions, the growth of gradient terms in the Einstein equations relative to non-gradient terms either prevents inflation from lasting for enough $e$-folds or triggers a phase of quantum runaway. We show that the fine-tuning of initial conditions necessary to avoid these issues becomes more severe as the energy scale of inflation is made smaller, disfavoring common approaches for reducing the tensor-to-scalar ratio.

FOS: Physical sciences↗
DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Predicting Kyasanur forest disease in resource-limited settings using event-based surveillance and transfer learning

In recent years, the reports of Kyasanur forest disease (KFD) breaking endemic barriers by spreading to new regions and crossing state boundaries is alarming. Effective disease surveillance and reporting systems are lacking for this emerging zoonosis, hence hindering control and prevention efforts. We compared time-series models using weather data with and without Event-Based Surveillance (EBS) information, i.e., news media reports and internet search trends, to predict monthly KFD cases in humans. We fitted Extreme Gradient Boosting (XGB) and Long Short-Term Memory models at the national and regional levels. We utilized the rich epidemiological data from endemic regions by applying Transfer Learning (TL) techniques to predict KFD cases in new outbreak regions where disease surveillance information was scarce. Overall, the inclusion of EBS data, in addition to the weather data, substantially increased the prediction performance across all models. The XGB method produced the best predictions at the national and regional levels. The TL techniques outperformed baseline models in predicting KFD in new outbreak regions. Novel sources of data and advanced machine-learning approaches, e.g., EBS and TL, show great potential towards increasing disease prediction capabilities in data-scarce scenarios and/or resource-limited settings, for better-informed decisions in the face of emerging zoonotic threats.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Developing ductile and isotropic Ti alloy with tailored composition for laser powder bed fusion

Laser powder bed fusion (LPBF) is an advanced technology to create metallic components with complex geometry. Ti-6Al-4V (Ti64) is one of the most frequently used materials for LPBF. The intrinsic high cooling rates and ultra-high directional thermal gradient of melt, however, lead to a hierarchical structure composed of fine martensitic alpha' within columnar prior-beta grains in LPBF Ti64. This microstructure results in poor ductility and strong mechanical anisotropy in as-built Ti64 components, which largely limit their applications. Here, we report that the above shortcomings can be overcome by tailoring the alloy composition. Specifically, by reducing the aluminum content from 6 wt% to 4 wt%, the resultant LPBF Ti-4Al-4V (Ti44) alloy exhibits substantially improved ductility and mechanical isotropy. These improvements are attributed to grain refinement during solidification and the activation of multiple slip modes and twinning in the as-built material during deformation. This work offers a simple strategy to design new titanium alloys for LPBF with significantly improved ductility and isotropy.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Dual color optogenetic control of neural populations using low-noise, multishank optoelectrodes

Optogenetics allows for optical manipulation of neuronal activity and has been increasingly combined with intracellular and extracellular electrophysiological recordings. Genetically-identified classes of neurons are optically manipulated, though the versatility of optogenetics would be increased if independent control of distinct neural populations could be achieved on a sufficient spatial and temporal resolution. We report a scalable multisite optoelectrode design that allows simultaneous optogenetic control of two spatially intermingled neuronal populations in vivo. We describe the design, fabrication, and assembly of low-noise, multisite/multicolor optoelectrodes. Each shank of the four-shank assembly is monolithically integrated with 8 recording sites and a dual-color waveguide mixer with a 7 × 30 μm cross-section, coupled to 405 nm and 635 nm injection laser diodes (ILDs) via gradient-index (GRIN) lenses to meet optical and thermal design requirements. To better understand noise on the recording channels generated during diode-based activation, we developed a lumped-circuit modeling approach for EMI coupling mechanisms and used it to limit artifacts to amplitudes under 100 μV upto an optical output power of 450 μW. We implanted the packaged devices into the CA1 pyramidal layer of awake mice, expressing Channelrhodopsin-2 in pyramidal cells and ChrimsonR in paravalbumin-expressing interneurons, and achieved optical excitation of each cell type using sub-mW illumination. We highlight the potential use of this technology for functional dissection of neural circuits.

42 ENGINEERING↗
DOE OSTI2017

Comparative first-principles studies of prototypical ferroelectric materials by LDA, GGA, and SCAN meta-GGA

Originating from a broken spatial inversion symmetry, ferroelectricity is a functionality of materials with an electric dipole that can be switched by external electric fields. Spontaneous polarization is a crucial ferroelectric property, and its amplitude is determined by the strength of polar structural distortions. Density functional theory (DFT) is one of the most widely used theoretical methods to study ferroelectric properties, yet it is limited by the levels of approximations in electron exchange-correlation. On the one hand, the local density approximation (LDA) is considered to be more accurate for the conventional perovskite ferroelectrics such as BaTiO 3 and PbTiO 3 than the generalized gradient approximation (GGA),which suffers from the so-called super-tetragonality error. On the other hand, GGA is more suitable for hydrogen-bonded ferroelectrics than LDA, which largely overestimates the strength of hydrogen bonding in general. We demonstrate here that the recently developed general-purpose strongly constrained and appropriately normed (SCAN) meta-GGA functional significantly improves over the traditional LDA/GGA for structural, electric, and energetic properties of diversely bonded ferroelectric materials with a comparable computational effort and thus enhances largely the predictive power of DFT in studies of ferroelectric materials. We also address the observed system-dependent performances of LDA and GGA for ferroelectrics from a chemical bonding point of view.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Energy efficient spark plasma sintering: Breaking the threshold of large dimension tooling energy consumption

An energy efficient spark plasma sintering method enabling the densification of large size samples assisted by very low electric current levels is developed. In this method, the electric current is concentrated in the graphite foils around the sample. High-energy dissipation is then achieved in this area enabling the heating and full densification of large (alumina) parts (Ø 40 mm) at relatively low currents (800 A). The electrothermal mechanical simulation reveals that the electric current needed to heat the large samples is 70% lower in the energy efficient configuration compared to the traditional configuration. The presence of thermal and densification gradients is also revealed for the larger size samples. Potential solutions for this problem are discussed. The experiments confirm the possibility of full densification (96%-99%) of large alumina samples. This approach allows using small (and low cost) SPS devices (generally limited to 10-15 mm samples) for large-size samples (40-50 mm). Finally, the developed technique enables also an optimized energy consumption by large-scale SPS systems.

36 MATERIALS SCIENCE↗
DOE OSTI2016

SU-F-J-97: A Joint Registration and Segmentation Approach for Large Bladder Deformations in Adaptive Radiotherapy

Purpose: To improve results of deformable image registration (DIR) in adaptive radiotherapy for large bladder deformations in CT/CBCT pelvis imaging. Methods: A variational multi-modal DIR algorithm is incorporated in a joint iterative scheme, alternating between segmentation based bladder matching and registration. Using an initial DIR to propagate the bladder contour to the CBCT, in a segmentation step the contour is improved by discrete image gradient sampling along all surface normals and adapting the delineation to match the location of each maximum (with a search range of +−5/2mm at the superior/inferior bladder side and step size of 0.5mm). An additional graph-cut based constraint limits the maximum difference between neighboring points. This improved contour is utilized in a subsequent DIR with a surface matching constraint. By calculating an euclidean distance map of the improved contour surface, the new constraint enforces the DIR to map each point of the original contour onto the improved contour. The resulting deformation is then used as a starting guess to compute a deformation update, which can again be used for the next segmentation step. The result is a dense deformation, able to capture much larger bladder deformations. The new method is evaluated on ten CT/CBCT male pelvis datasets, calculating Dice similarity coefficients (DSC) between the final propagated bladder contour and a manually delineated gold standard on the CBCT image. Results: Over all ten cases, an average DSC of 0.93±0.03 is achieved on the bladder. Compared with the initial DIR (0.88±0.05), the DSC is equal (2 cases) or improved (8 cases). Additionally, DSC accuracy of femoral bones (0.94±0.02) was not affected. Conclusion: The new approach shows that using the presented alternating segmentation/registration approach, the results of bladder DIR in the pelvis region can be greatly improved, especially for cases with large variations in bladder volume. Fraunhofer MEVIS received funding from a research grant by Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Boride-based Ceramic Super-high Temperature Thermocouples in Harsh Environments (Final Scientific/Technical Report)

An electromotive force (emf) can be generated along a temperature gradient between the cold end and hot end of a thermoelectric material, termed the Seebeck effect. Based on the Seebeck effect, metallic alloys have been extensively employed to detect temperatures for centuries, named thermocouples. However, commercially available thermocouple alloys suffer from limitations, such as oxidation, chemical degradation, and poor long-term stability under high-temperature harsh environments. This DOE-funded project aimed to develop high-temperature, chemically tolerant thermocouples suitable for operation in extreme environments relevant to semiconducting thermoelectric materials. The research focused on boride-based semiconducting thermoelectric compounds as candidates for next-generation thermocouples with enhanced oxidation resistance, chemical stability, and thermal robustness under conditions representative of charcoal-fired electricity facilities. During the funded years, boride materials were synthesized using an arc-plasma technique under ambient air and argon atmospheres, enabling scalable and cost-effective production compared with conventional boride fabrication methods. The synthesized borides were processed into nanostructured powders, followed by consolidation into dense bulk materials using a spark plasma sintering (SPS) bottom-up approach. Comprehensive characterization was performed, including microstructural analysis, electrical transport measurements, and optical and thermal property evaluation. Both p-type and n-type boride electric legs were fabricated and integrated into boride-based thermocouples. The thermal and irradiation stabilities of the boride nanomaterials and bulk thermoelectric materials were systematically evaluated to assess suitability for long-term operation in harsh environments. Additionally, 12 students were broadly hands-on trained spanning the full research workflow, including word processing and technical editing (e.g., LATEX for manuscript and poster preparation), data collection and analysis (using Python and related libraries and hardware interfaces), sample preparation (including arc-plasma synthesis and spark plasma sintering), and advanced characterization techniques (such as X-ray diffraction, UV–vis spectroscopy, electron microscopy, differential thermal analysis (DTA), and Seebeck coefficient measurements, etc). Overall, this project demonstrated the feasibility of boride-based thermoelectric materials as durable high-temperature thermocouples, providing a promising pathway toward robust temperature sensing technologies aligned with DOE energy infrastructure and extreme-environment monitoring needs.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Jipole: A Differentiable ipole-based Code for Radiative Transfer in Curved Spacetimes

Recent imaging of supermassive black holes by the Event Horizon Telescope has relied on exhaustive parameter-space searches, matching observations to large, precomputed libraries of theoretical models. As observational data become increasingly precise, the limitations of this computationally expensive approach grow more acute, creating a pressing need for more efficient methods. In this work, we present Jipole, an automatically differentiable (AD), ipole-based code for radiative transfer in curved spacetimes, designed to compute image gradients with respect to underlying model parameters. These gradients quantify how parameter changes—such as the black hole’s spin or the observer’s inclination—affect the image, enabling more efficient parameter estimation and reducing the number of required images. We validate Jipole against ipole in two analytical tests and then compare pixelwise intensity derivatives from AD with those from finite-difference methods. We then demonstrate the utility of these gradients by performing parameter recovery for an analytical model in three increasingly complex cases for the injected image: ideal, blurred, and blurred with added noise. In most cases, high-accuracy fits are obtained in only a few optimization steps, failing only in cases with extremely low signal-to-noise ratios. These results highlight the potential of AD-based methods to accelerate robust, high-fidelity model-data comparisons in current and future black hole imaging efforts.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Molecules and chains in a strong magnetic field - Statistical treatment

A Thomas-Fermi-Dirac-Weizsaecker statistical model is developed and employed to investigate diatomic molecules and infinite molecular chains in strong magnetic fields. The standard magnetic Thomas-Fermi-Dirac kinetic, potential, and exchange energy functionals are supplemented by a gradient correction to the kinetic energy. The numerical method used for solving this system in two spatial dimensions is detailed. Numerical solutions for a wide range of magnetic strengths and elements are presented to demonstrate the robustness, as well as the limitations, of the statistical approach. These calculations qualitatively reproduce many of the results of detailed quantum mechanical treatments. For example, the fractional binding energy is greatest for low atomic numbers and for strong magnetic fields.

Abrahams, Andrew M.↗
NASA NTRS2013

Recalibration of the Shear Stress Transport Model to Improve Calculation of Shock Separated Flows

The Menter Shear Stress Transport (SST) k . turbulence model is one of the most widely used two-equation Reynolds-averaged Navier-Stokes turbulence models for aerodynamic analyses. The model extends Menter s baseline (BSL) model to include a limiter that prevents the calculated turbulent shear stress from exceeding a prescribed fraction of the turbulent kinetic energy via a proportionality constant, a1, set to 0.31. Compared to other turbulence models, the SST model yields superior predictions of mild adverse pressure gradient flows including those with small separations. In shock - boundary layer interaction regions, the SST model produces separations that are too large while the BSL model is on the other extreme, predicting separations that are too small. In this paper, changing a1 to a value near 0.355 is shown to significantly improve predictions of shock separated flows. Several cases are examined computationally and experimental data is also considered to justify raising the value of a1 used for shock separated flows.

Georgiadis, Nicholas J.↗
NASA NTRS2016

Least Squares Shadowing Sensitivity Analysis of Chaotic Flow Around a Two-Dimensional Airfoil

Gradient-based sensitivity analysis has proven to be an enabling technology for many applications, including design of aerospace vehicles. However, conventional sensitivity analysis methods break down when applied to long-time averages of chaotic systems. This breakdown is a serious limitation because many aerospace applications involve physical phenomena that exhibit chaotic dynamics, most notably high-resolution large-eddy and direct numerical simulations of turbulent aerodynamic flows. A recently proposed methodology, Least Squares Shadowing (LSS), avoids this breakdown and advances the state of the art in sensitivity analysis for chaotic flows. The first application of LSS to a chaotic flow simulated with a large-scale computational fluid dynamics solver is presented. The LSS sensitivity computed for this chaotic flow is verified and shown to be accurate, but the computational cost of the current LSS implementation is high.

Blonigan, Patrick J.↗
DOE OSTI2015

Magnetic design constraints of helical solenoids

Helical solenoids have been proposed as an option for a Helical Cooling Channel for muons in a proposed Muon Collider. Helical solenoids can provide the required three main field components: solenoidal, helical dipole, and a helical gradient. In general terms, the last two are a function of many geometric parameters: coil aperture, coil radial and longitudinal dimensions, helix period and orbit radius. In this paper, we present design studies of a Helical Solenoid, addressing the geometric tunability limits and auxiliary correction system.

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

Restoration of the Potosi Dynamic Model 2010

In topical Report DOE/FE0002068-1 [2] technical performance evaluations on the Cambrian Potosi Formation were performed through reservoir modeling. The data included formation tops from mud logs, well logs from the VW1 and the CCS1 wells, structural and stratigraphic formation from three dimensional (3D) seismic data, and field data from several waste water injection wells for Potosi Formation. Intention was for two million tons per annum (MTPA) of CO2 to be injected for 20 years. In this Task the 2010 Potosi heterogeneous model (referred to as the "Potosi Dynamic Model 2010" in this report) was re-run using a new injection scenario; 3.2 MTPA for 30 years. The extent of the Potosi Dynamic Model 2010, however, appeared too small for the new injection target. It was not sufficiently large enough to accommodate the evolution of the plume. Also, it might have overestimated the injection capacity by enhancing too much the pressure relief due to the relatively close proximity between the injector and the infinite acting boundaries. The new model, Potosi Dynamic Model 2013a, was built by extending the Potosi Dynamic Model 2010 grid to 30 miles x 30 miles (48 km by 48 km), while preserving all property modeling workflows and layering. This model was retained as the base case. Potosi Dynamic Model 2013.a gives an average CO2 injection rate of 1.4 MTPA and cumulative injection of 43 Mt in 30 years, which corresponds to 45% of the injection target. This implies that according to this preliminary model, a minimum of three (3) wells could be required to achieve the injection target. The injectivity evaluation of the Potosi formation will be revisited in topical Report 15 during which more data will be integrated in the modeling exercise. A vertical flow performance evaluation could be considered for the succeeding task to determine the appropriate tubing size, the required injection tubing head pressure (THP) and to investigate whether the corresponding well injection rate falls within the tubing erosional velocity limit. After 30 years, the plume extends 15 miles (24 km) in E-W and 14 miles (22 km) in N-S directions. After injection is completed, the plume continues to migrate laterally, mainly driven by the remaining pressure gradient. After 100 years post-injection, the plume extends 17 miles (27 km) in E-W and 15 miles (24 km) in N-S directions. The increase of reservoir pressure at the end of injection is approximately 370 psia around the injector and gradually decreases away from the well. The reservoir pressure increase is less than 30 psia beyond 14 miles (22 km) away from injector. The initial reservoir pressure is restored after approximately 20 years post-injection. This result, however, is associated with uncertainties on the boundary conditions, and a sensitivity analysis could be considered for the succeeding tasks. It is important to remember that the respective plume extent and areal pressure increase corresponds to an injection of 43 Mt CO2. Should the targeted cumulative injection of 96 Mt be achieved; a much larger plume extent and areal pressure increase could be expected. Re-evaluating the permeability modeling, vugs and heterogeneity distributions, and relative permeability input could be considered for the succeeding Potosi formation evaluations. A simulation using several injectors could also be considered to determine the required number of wells to achieve the injection target while taking into account the pressure interference.

Adushita, Yasmin↗
DOE OSTI2015

Direct identification of predator-prey dynamics in gyrokinetic simulations

The interaction between spontaneously formed zonal flows and small-scale turbulence in nonlinear gyrokinetic simulations is explored in a shearless closed field line geometry. It is found that when clear limit cycle oscillations prevail, the observed turbulent dynamics can be quantitatively captured by a simple Lotka-Volterra type predator-prey model. Fitting the time traces of full gyrokinetic simulations by such a reduced model allows extraction of the model coefficients. Scanning physical plasma parameters, such as collisionality and density gradient, it was observed that the effective growth rates of turbulence (i.e., the prey) remain roughly constant, in spite of the higher and varying level of primary mode linear growth rates. The effective growth rate that was extracted corresponds roughly to the zonal-flow-modified primary mode growth rate. It was also observed that the effective damping of zonal flows (i.e., the predator) in the parameter range, where clear predator-prey dynamics is observed, (i.e., near marginal stability) agrees with the collisional damping expected in these simulations. This implies that the Kelvin-Helmholtz-like instability may be negligible in this range. The results imply that when the tertiary instability plays a role, the dynamics becomes more complex than a simple Lotka-Volterra predator prey.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗