Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence rate”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Extension during continental convergence, with application to the Tibetan Plateau

An explanation in terms of the thermal evolution of thickened continental lithosphere is offered to explain the transition, in the late Tertiary to Quaternary, from north-south compression to east-west extension in the strain rate field of the Tibetan Plateau. The lower part of the lithosphere consists of a thermal boundary layer which, when thickened by horizontal shortening, is colder and denser than its surroundings. Convective instability of the thickened thermal boundary layer and its replacement by hot asthenosphere would rapidly raise the surface elevation and gravitational potential energy of the overlying part of the lithosphere. The convective instability would happen in a time brief compared with the collision time scale but would only occur after there had already been substantial thickening of the lithosphere. The increase in surface height and of potential energy are sufficient for east-west extension to replace north-south compression as the dominant feature of the stress field.

England, Philip↗
DOE OSTI2023

Unraveling Anisotropic and Pulsating Etching of ZnO Nanorods in Hydrochloric Acid via Correlative Electron Microscopy

Despite much technical progress achieved so far, the exact surface and shape evolution during wet chemical etching is less unraveled, especially in ionically bonded ceramics. Herein, by using in situ liquid cell transmission electron microscopy, a repeated two-stage anisotropic and pulsating periodic etching dynamic is discovered during the pencil shape evolution of a single crystal ZnO nanorod in aqueous hydrochloric acid. Specifically, the nanopencil tip shrinks at a slower rate along [0001̅] than that along the $\langle$101̅0$\rangle$ directions, resulting in a sharper ZnO pencil tip. Afterward, rapid tip dissolution happens due to accelerated etching rates along various crystal directions. Concurrently, the vicinal base region of the original nanopencil tip emerges as a new tip followed by the repeated sequence of tip shrinking and removal. The high-index surfaces, such as {101̅m} (m = 0, 1, 2, or 3) and {21̅1̅n} (n = 0, 1, 2, or 3), are found to preferentially expose in different ratios. Our 3D electron tomography, convergent beam electron diffraction, middle-angle bright-field STEM, and XPS results indicate the dissociative Cl – species were bound to the Zn-terminated tip surfaces. Furthermore, DFT calculation suggests the preferential Cl – passivation over the {101̅1} and (0001) surfaces of lower energy than others, leading to preferential surface exposures and the oscillatory variation of different facet etching rates. The boosted reactivity due to high-index nanoscale surface exposures is confirmed by comparatively enhanced chemical sensing and CO 2 hydrogenation activity. In conclusion, these findings provide an in-depth understanding of anisotropic wet chemical etching of ionic nanocrystals and offer a design strategy for advanced functional materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Pleistocene Colorado River terraces in the canyonlands region (Utah, USA) record unsteady, transient incision and growth of the Cataract Canyon knickzone by salt tectonics

Colorado River terraces in the canyonlands region of southeastern Utah can be dated and analyzed to address the controls of incision and nature of the Colorado Plateau's largest channel-steepness anomaly, Cataract Canyon. Field correlations supported by luminescence and cosmogenic-nuclide ages on strath terraces along Meander Canyon, upstream of Cataract Canyon, reveal a complex record of unsteady incision over the past ~340 k.y. at an average rate of ~0.4 mm/yr. Both an upstream progression of rapid incision and a unique sigmoidal long-profile pattern of terraces indicate incision in response to episodic baselevel fall. Also, terraces converge downstream with the anomalously low-gradient modern channel above Cataract Canyon. Here we interpret these results as indicating that growth of the Cataract Canyon knickzone is due to an erosion-salt tectonics feedback since at least the Mid Pleistocene, which has imparted unsteady, and currently elevated, local baselevel. More broadly, the canyonlands region is marked by rapid and unsteady incision that is complicated by local geologic controls even while being absent of any regional or mantle-driven uplift.

58 GEOSCIENCES↗
NASA NTRS1987

Accuracy of derivatives of control performance using a reduced structural model

The sensitivity of control system performance to structural changes is calculated for a multi-span beam with direct-rate feedback vibration control. Reduced models based on the natural modes of the structure are formed and derivatives of the damping ratios of the closed-loop eigenvalues are calculated. The convergence of the derivatives of the damping ratios with increasing number of modes is shown to be slower than the convergence of the damping ratios themselves. In particular, in some cases the convergence of finite-element approximations to the derivatives is much faster than the convergence of the modal approximations. The results indicate that the use of reduced models based on natural vibration modes may be ill-advised for calculating the sensitivity of control system performance to changes in the controlled structure.

Sandridge, Chris A.↗
NASA NTRS2016

Numerical Evaluation of an Ejector-Enhanced Resonant Pulse Combustor with a Poppet Inlet Valve and a Converging Exhaust Nozzle

A computational investigation of a pressure-gain combustor system for gas turbine applications is presented. The system consists of a valved pulse combustor and an ejector, housed within a shroud. The study focuses on two enhancements to previous models, related to the valve and ejector components. First, a new poppet inlet valve system is investigated, replacing the previously used reed valve configuration. Secondly, a new computational approach to approximating the effects of choked turbine inlet guide vanes present immediately downstream of the Ejector-Enhanced Resonant Pulse Combustor (EERPC) is investigated. Instead of specifying a back pressure at the EERPC exit boundary (as was done in previous studies) the new model adds a converging-diverging (CD) nozzle at the exit of the EERPC. The throat area of the CD nozzle can be adjusted to obtain the desired back pressure level and total mass flow rate. The results presented indicate that the new poppet valve configuration performs nearly as well as the original reed valve system, and that the addition of the CD nozzle is an effective method to approximate the exit boundary effects of a turbine present downstream of the EERPC. Furthermore, it is shown that the more acoustically reflective boundary imposed by a nozzle as compared to a constant pressure surface does not significantly affect operation or performance.

Computational Fluid Dynamics (CFD↗
DOE OSTI2021

Predicting nonresonant pressure-driven MHD modes in equilibria with low magnetic shear

Nonresonant internal modes can be difficult to anticipate as there is no resonant surface in the plasma. However, equilibria that are unstable to multiple nonresonant magnetohydrodynamic (MHD) modes may be more prone to global loss of confinement since these instabilities generate spatially extended linear displacements, potentially enhancing magnetic field line chaos via nonlinear interactions. Here, we successfully predict the unstable nonresonant pressure-driven modes for equilibria with zero shear in the plasma core, irrational q on axis, and a central pressure gradient, which is consistent with pre-crash profiles in sawtoothing tokamak plasmas in the large-aspect-ratio limit. A criterion for identifying nonresonant modes most likely to be unstable is developed from the convergents of the continued fraction representation of q 0 . A higher-order analysis of the standard Energy Principle reveals the conditions under which these modes are expected to dominate. Linear growth rate spectra, as a function of toroidal mode number (up to n = 30), calculated using the initial-value extended-MHD code, M3D-C1, recover the characteristic dependence observed for ideal infernal modes. Nonresonant modes have also been invoked in some ideal sawtooth crash models. Finally, this work provides a mechanism to predict the mode numbers of infernal modes and, potentially, the width of some post-sawtooth-crash profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Blockchain for Fault-Tolerant Grid Operations

Radial topology and vast geographic coverage make distribution systems prone to widespread power outages upon the failure of a single (or multiple) upstream component. Fault-handling algorithms depend heavily on correct estimations of the system’s state to effectively isolate the affected area and reduce the number of affected customers while maintaining operational safety. The work described here leverages the core features of distributed, consensus-based decision-making processes and the immutability of blockchain, and demonstrates their value in improving fault-tolerant grid operations. In this work, blockchain was used to create a trusted data-sharing platform that enables independent actors to reconstruct the system state; this enables distributed resources to make intelligent decisions with limited knowledge. Although the process requires data sharing, its algorithms have been designed to limit the amount of private information that is exchanged, which helps preserve business-sensitive data and maintain customer privacy. In addition, by reducing the information that must be shared, the communication requirements are also reduced; (however, an in-depth analysis of the communication requirements is beyond the scope of this project). The proposed use cases are intended to represent a foundational basis for third parties to develop functional solutions that can eventually be deployed in the field. To further provide guidance, the envisioned use cases have incorporated design requirements that consider the blockchain characteristics and a need to limit information from surrounding resources, which preserve the assumption and the possibility that such resources could belong to different entities. This report presents a detailed design of the three use cases with the tools needed to enable the analysis being tested. The implemented gross error detection method can detect mismatches when the error exceeds 3.8 times the sensor’s rated accuracy. Detection of the circuit breaker state successfully identified the correct states across all simulation tests. A distribution-system power-flow solution in the simulator OpenDSS generally possesses a convergency tolerance of 0.01% on the voltage magnitude. The evaluation of possible reconnection using voltage magnitude—preserving the data ownership—has a voltage magnitude difference smaller than 0.001% from the OpenDSS result. The results preserving data ownership have a difference within the expected power flow tolerance with full knowledge of the system, which surpasses expectations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1989

Closure of the Africa-Eurasia-North America plate motion circuit and tectonics of the Gloria fault

The current motions of the African, Eurasian, and North American plates are examined. The problems addressed include whether there is resolvable motion of a Spitsbergen microplate, the direction of motion between the African and North American plates, whether the Gloria fault is an active transform fault, and the implications of plate circuit closures for rates of intraplate deformation. Marine geophysical data and magnetic profiles are used to construct a model which predicts about 4 mm/yr slip across the Azores-Gibraltar Ridge, and west-northwest convergence near Gibraltar. The analyzed data are consistent with a rigid plate model with the Gloria fault being a transform fault.

Argus, Donald F.↗
DOE OSTI2021

Excited State Intramolecular Proton Transfer with Nuclear-Electronic Orbital Ehrenfest Dynamics

The recent development of the Ehrenfest dynamics approach in the nuclear-electronic orbital (NEO) framework provides a promising way to simulate coupled nuclear-electronic dynamics. Our previous study showed that the NEO-Ehrenfest approach with a semi- classical traveling proton basis method yields accurate predictions of molecular vibrational frequencies. In this work, we provide a more thorough analysis of the semiclassical traveling proton basis method to elucidate its validity and convergence behavior. We also conduct NEO-Ehrenfest dynamics simulations to study an excited state intramolecular proton transfer process. These simulations reveal that nuclear quantum effects influence the predictions of proton transfer reaction rates and kinetic isotope effects due to the intrinsic delocalized nature of the quantum nuclear wave function. This work illustrates the importance of nuclear quantum effects in coupled nuclear-electronic dynamical processes and shows that the NEO-Ehrenfest approach can be a powerful tool to provide insights and predictions for these processes.

Zhao, Luning↗
NASA NTRS1979

Structural analysis of hollow blades: Torsional stress analysis of hollow fan blades for aircraft jet engines

A torsional stress analysis of hollow fans blades by the finite element method is presented. The fans are considered to be double circular arc blades, hollowed 30 percent, and twisted by a component of the centrifugal force by the rated revolution. The effects of blade hollowing on strength and rigidity are discussed. The effects of reinforcing webs, placed in the hollowed section in varying numbers and locations, on torsional rigidity and the convergence of stresses, are reported. A forecast of the 30 percent hollowing against torsional loadings is discussed.

Ogawa, A.↗
NASA NTRS1982

Aerothermal modeling

The objectives, approach, and status of a program to develop the computational fluid dynamics tools needed to improve combustor design and analysis are outlined. The calculation procedure selected consists of a finite difference solution of the time averaged, steady state, primitive variable, elliptic form of the Reynolds equations. Standard TEACH type numerics are used to solve the resulting equations. These include hybrid differencing, SIMPLE algorithm for the pressure field, line by line iterative solution using the ADI method and the tridiagonal matrix algorithm (TDMA). Convergence is facilitated by using under relaxation. The physical processes are modeled by a two equation eddy viscosity model for turbulence; combustion is represented by a simple, irreversible, one step chemical reaction whose rate is influenced only by the time scale of the turbulence. The model evaluation procedure is also described.

Sturgess, G. J.↗
DOE OSTI2025

Inertial collapse of a gas bubble in a shear flow near a rigid wall

Despite the extensive research on bubble collapse near rigid walls, the bubble collapse dynamics in the presence of shear flow near a rigid wall is poorly understood. We conduct direct simulations of the Navier–Stokes equations to explore the bubble dynamics and pressures during bubble collapse near a rigid, flat wall under linear shear flow conditions. We examine the dependence of the bubble collapse morphology and wall pressures on the initial bubble location and shear rate. We find that shear distorts the bubble, generating two re-entrant jets – one developing from the side opposite to the mean flow and the other from the far end toward the wall. Upon impact of the jet on the opposite side of the bubble, water-hammer shocks are produced, which propagate outward and interact with the convoluted bubble shape. The shock stretches the bubble towards the wall, resulting in a closer impact location for the jet originating from the far end compared with the case with no shear flow. The water-hammer pressure location can be approximated as the theoretical distance travelled by a particle initialised at the bubble centre with the corresponding constant shear flow velocity. The maximum wall pressures can thus be predicted by considering the distance between the far jet impingement location and the wall along the wall-normal direction. As the shear rate is increased, the maximum wall pressure increases, although only marginally. We determine the critical initial stand-off distance from the wall at which the bubble morphology is shear dominated, i.e. characterised by converging re-entrant jets.

Bhola, Sahil (ORCID:0000000236766302)↗
DOE OSTI2026

4th TDAMM Workshop White Paper

Time-Domain and Multi-Messenger Astrophysics (TDAMM) is entering a fundamentally new phase characterized by an unprecedented increase in the rate and diversity of astrophysical transient detections. The community is transitioning from a discovery-limited to a follow-up-limited era, driven by major investments across electromagnetic, gravitational-wave, and neutrino observatories. Upcoming facilities such as the Vera C. Rubin Observatory, the Nancy Grace Roman Space Telescope, and wide-field survey instruments will produce a deluge of time-domain alerts, reaching millions of events per night. Simultaneously, upgrades to the gravitational-wave network (LVK O5 and beyond) and neutrino observatories (IceCube Gen2) will significantly increase the detection rates of non-electromagnetic messengers. New high-energy missions and expansions of the InterPlanetary Network (IPN) will further enhance discovery capabilities across the gamma-ray and X-ray regimes. This convergence of capabilities represents a transformative opportunity: for the first time, the community will routinely detect rare and high-impact events across multiple messengers. However, the scientific return from these discoveries will depend critically on the ability to rapidly identify, prioritize, and coordinate follow-up observations across a heterogeneous and globally distributed set of facilities.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

Characteristic-based algorithms for flows in thermo-chemical nonequilibrium

A generalized finite-rate chemistry algorithm with Steger-Warming, Van Leer, and Roe characteristic-based flux splittings is presented in three-dimensional generalized coordinates for the Navier-Stokes equations. Attention is placed on convergence to steady-state solutions with fully coupled chemistry. Time integration schemes including explicit m-stage Runge-Kutta, implicit approximate-factorization, relaxation and LU decomposition are investigated and compared in terms of residual reduction per unit of CPU time. Practical issues such as code vectorization and memory usage on modern supercomputers are discussed.

Walters, Robert W.↗
NASA NTRS1990

Ignition transient analysis of solid rocket motor

To predict pressure-time and thrust-time behavior of solid rocket motors, a one-dimensional numerical model is developed. The ignition phase of solid rocket motors (time less than 0.4 sec) depends critically on complex interactions among many elements, such as rocket geometry, heat and mass transfer, flow development, and chemical reactions. The present model solves the mass, momentum, and energy equations governing the transfer processes in the rocket chamber as well as the attached converging-diverging nozzle. A qualitative agreement with the SRM test data in terms of head-end pressure gradient and the total thrust build-up is obtained. Numerical results show that the burning rate in the star-segmented head-end section and the erosive burning are two important parameters in the ignition transient of the solid rocket motor (SRM).

Han, Samuel S.↗
NASA NTRS2007

Computational Study of an Axisymmetric Dual Throat Fluidic Thrust Vectoring Nozzle for a Supersonic Aircraft Application

A computational investigation of an axisymmetric Dual Throat Nozzle concept has been conducted. This fluidic thrust-vectoring nozzle was designed with a recessed cavity to enhance the throat shifting technique for improved thrust vectoring. The structured-grid, unsteady Reynolds- Averaged Navier-Stokes flow solver PAB3D was used to guide the nozzle design and analyze performance. Nozzle design variables included extent of circumferential injection, cavity divergence angle, cavity length, and cavity convergence angle. Internal nozzle performance (wind-off conditions) and thrust vector angles were computed for several configurations over a range of nozzle pressure ratios from 1.89 to 10, with the fluidic injection flow rate equal to zero and up to 4 percent of the primary flow rate. The effect of a variable expansion ratio on nozzle performance over a range of freestream Mach numbers up to 2 was investigated. Results indicated that a 60 circumferential injection was a good compromise between large thrust vector angles and efficient internal nozzle performance. A cavity divergence angle greater than 10 was detrimental to thrust vector angle. Shortening the cavity length improved internal nozzle performance with a small penalty to thrust vector angle. Contrary to expectations, a variable expansion ratio did not improve thrust efficiency at the flight conditions investigated.

Deere, Karen A.↗
DOE OSTI2026

Characterization and automated optimization of laser-driven proton beams from converging liquid sheet jet targets

Compact, stable, and versatile laser-driven ion sources hold great promise for applications ranging from medicine to materials science and fundamental physics. While single-shot sources have demonstrated favorable beam properties, including the peak fluxes necessary for several applications, high-repetition-rate operation will be necessary to generate and sustain the high average flux needed for many of the most exciting applications of laser-driven ion sources. Further, to navigate through the high-dimensional space of laser and target parameters toward experimental optima, it is essential to develop ion acceleration platforms compatible with machine learning techniques and capable of autonomous real-time optimization. Here, we present a multi-Hz ion acceleration platform employing a liquid sheet jet target. We characterize the laser-plasma interaction and the laser-driven proton beam across a variety of key parameters governing the interaction using an extensive suite of online diagnostics. We also demonstrate real-time, closed-loop optimization of the ion beam maximum energy by tuning the laser wave front using a Bayesian optimization scheme. This approach increased the maximum proton energy by 11% compared to a manually optimized wave front by enhancing the energy concentration within the laser focal spot, demonstrating the potential for closed-loop optimization schemes to tune future ion accelerators for robust high-repetition-rate operation.

Glenn, G. D. [SLAC National Accelerator Laboratory↗
DOE OSTI2021

Growth of Mesoscale Convective Systems in Observations and a Seasonal Convection-Permitting Simulation over Argentina

Abstract A 6.5-month, convection-permitting simulation is conducted over Argentina covering the Remote Sensing of Electrification, Lightning, And Mesoscale/Microscale Processes with Adaptive Ground Observations and Clouds, Aerosols, and Complex Terrain Interactions (RELAMPAGO-CACTI) field campaign and is compared with observations to evaluate mesoscale convective system (MCS) growth prediction. Observed and simulated MCSs are consistently identified, tracked, and separated into growth, mature, and decay stages using top-of-the-atmosphere infrared brightness temperature and surface rainfall. Simulated MCS number, lifetime, seasonal and diurnal cycles, and various cloud-shield characteristics including growth rate are similar to those observed. However, the simulation produces smaller rainfall areas, greater proportions of heavy rainfall, and faster system propagations. Rainfall area is significantly underestimated for long-lived MCSs but not for shorter-lived MCSs, and rain rates are always overestimated. These differences result from a combination of model and satellite retrieval biases, in which simulated MCS rain rates are shifted from light to heavy, while satellite-retrieved rainfall is too frequent relative to rain gauge estimates. However, the simulation reproduces satellite-retrieved MCS cloud-shield evolution well, supporting its usage to examine environmental controls on MCS growth. MCS initiation locations are associated with removal of convective inhibition more than maximized low-level moisture convergence or instability. Rapid growth is associated with a stronger upper-level jet (ULJ) and a deeper northwestern Argentinean low that causes a stronger northerly low-level jet (LLJ), increasing heat and moisture fluxes, low-level vertical wind shear, baroclinicity, and instability. Sustained growth corresponds to similar LLJ, baroclinicity, and instability conditions but is less sensitive to the ULJ, large-scale vertical motion, or low-level shear. Growth sustenance controls MCS maximum extent more than growth rate.

54 ENVIRONMENTAL SCIENCES↗