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,639 records · Page 91

Three-dimensional potential flows from functions of a 3D complex variable

Potential, or ideal, flow velocities can be found from the gradient of an harmonic function. An ordinary complex valued analytic function can be written as the sum of two real valued functions, both of which are harmonic. Thus, 2D complex valued functions serve as a source of functions that describe two-dimensional potential flows. However, this use of complex variables has been limited to two-dimensions. Recently, a new system of three-dimensional complex variables has been developed at the NASA Ames Research Center. As a step toward application of this theory to the analysis of 3D potential flow, several functions of a three-dimensional complex variable have been investigated. The results for two such functions, the 3D exponential and 3D logarithm, are presented in this paper. Potential flows found from these functions are investigated. Important characteristics of these flows fields are noted.

Kelly, Patrick↗

A study of the topology of dissipating motions in direct numerical simulations of time-developing compressible and incompressible mixing layers

A preliminary investigation of the geometry of flow patterns in numerically simulated compressible and incompressible mixing layers was carried out using 3-D critical point methodology. Motions characterized by high rates of kinetic energy dissipation and/or high enstrophy were of particular interest. In the approach the partial derivatives of the velocity field are determined at every point in the flow. These are used to construct the invariants of the velocity gradient tensor and the rate-of-strain tensor (P, Q, R, and P(sub s), Q(sub s), R(sub s) respectively). For incompressible flow the first invariant is zero. For the conditions of the compressible simulation, the first invariant is found to be everywhere small, relative to the second and third invariants, and so in both cases the local topology at a point is mainly determined by the second and third invariants. The data at every grid point is used to construct scatter plots of Q versus R and Q(sub s) versus R(sub s). Most points map to a cluster near the origin in Q-R space. However, fine scale motions, that is motions which are characterized by velocity derivatives which scale with the square root of R(sub delta), tend to map to regions which lie far from the origin. Definite trends are observed for motions characterized by high enstrophy and/or high dissipation. The observed trends suggest that, for these motions, the second and third invariants of the velocity gradient and rate-of-strain tensors are strongly correlated. The second and third invariants of the rate-of-strain tensor are related by K(-Q(sub s))(exp 3/2), which is consistent with the above scaling of velocity derivatives. The quantity K appears to depend on Reynolds number with an upper limit K = 2(the square root of 3)/9 corresponding to locally axisymmetric flow. For both the compressible and incompressible mixing layer, regions corresponding to high rates of dissipation are found to be characterized by comparable magnitudes of R(sub ij)R(sub ij) and S(sub ij)S(sub ij). For the incompressible mixing layer, regions characterized by the highest values of enstrophy are found to have relatively low strain rates.

Chen, J. H.↗

Microbial Fossils from Terrestrial Subsurface Hydrothermal Environments: Examples and Implications for Mars

The recognition of biological signatures in ancient epithermal deposits has special relevance for studies of early blaspheme evolution and in exploring for past life on Mars. Recently, proposals for the existence of an extensive subsurface blaspheme on Earth, dominated by chemoautotrophic microbial life, has gained prominence. However, reports of fossilized microbial remains, or biosedimentary structures (e.g. stromatolites) from the deposits of ancient subsurface systems, are rare. Microbial preservation is favoured where high population densities co-exist with rapid mineral precipitation. Near-surface epithetical systems with strong gradients in temperature and redox are good candidates for the abundant growth and fossilization of microorganisms, and are also favorable environments for the precipitation of ore minerals. Therefore, we might expect microbial remain, to be particularly well preserved in various kinds of hydrothermal and diagenetic mineral precipitates that formed below the upper temperature limit for life (approx. 120 C).

Hofmann, Beda A.↗

Using Magnetic Fields to Control Convection during Protein Crystallization: Analysis and Validation Studies

The effect of convection during the crystallization of proteins is not very well understood. In a gravitational field, convection is caused by crystal sedimentation and by solutal buoyancy induced flow and these can lead to crystal imperfections. While crystallization in microgravity can approach diffusion limited growth conditions (no convection), terrestrially strong magnetic fields can be used to control fluid flow and sedimentation effects. In this work, we develop the analysis for magnetic flow control and test the predictions using analog experiments. Specifically, experiments on solutal convection in a paramagnetic fluid were conducted in a strong magnetic field gradient using a dilute solution of Manganese Chloride. The observed flows indicate that the magnetic field can completely counter the settling effects of gravity locally and are consistent with the theoretical predictions presented. This phenomenon suggests that magnetic fields may be useful in mimicking the microgravity environment of space for some crystal growth ana biological applications where fluid convection is undesirable.

Ramachandran, N.↗

Robustness Against Interfraction Prostate Movement in Scanned Ion Beam Radiation Therapy

Purpose: To assess the robustness of scanned ion beam treatment plans against the interfraction internal target motion and evaluate the limits of validity of target-based isocenter realignment for prostate cancer radiation therapy. Methods and Materials: For 12 prostate patients, scanned beam carbon ion treatment plans were prepared using 2 lateral opposed beams and the raster scanning technique with different clinical target volume (CTV) to planning target volume (PTV) margins (2-10 mm). Internal target motion of 2-15 mm in anteroposterior (AP), superoinferior (SI), and left-to-right (LR) directions was simulated by displacing the CTV contours with respect to the computed tomography data. The plans were recalculated with and without target-based isocenter realignment and the CTV coverage was assessed. Results: For CTV shifts within the applied planning margin, the CI-98% is greater than 98.0%, both with and without isocenter realignment. Without realignment, because of the sharp lateral gradients, the CI-98% shows rapid fall as soon as the target shift exceeds the applied planning margin for all displacement directions. With isocenter realignment, the coverage improves notably for shifts in AP and SI directions and the CI-98% is restored to >95.0% for plans optimized with a 2-mm margin, >97% with a 3-mm margin, and >98% with larger margins. For large corrections, predominately in AP direction, cold spots in the CTV may occur. Their magnitude is dependent on the patients' individual anatomies. Conclusions: Within the physiological limits of internal prostate movement, target-based isocenter realignment results in improved CTV coverage for shifts in AP and SI directions exceeding the applied planning margin. Assuming optimal patient setup reproducibility (eg, immobilization, setup error correction, patient preparation protocols), hence negligible interfraction bone and soft-tissue variations, changes in traversed densities resulting from target-based realignment applied with a reduced planning margin do not induce significant dose deterioration in most of the cases.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

The effects of advanced spectral line shapes on atmospheric carbon dioxide retrievals

Line shape effects such as speed-dependence and collisional (Dicke) narrowing are manifest in high signal-to-noise ratio laboratory spectra for a wide range of molecular species and sample conditions. However, these effects are usually neglected in spectroscopic models used in remote sensing of atmospheric greenhouse gases. Here, in this study, we examine the effects of line profile choice on laboratory spectra and XCO 2 values retrieved from Total Carbon Column Observing Network (TCCON) spectra. To this end, we utilize limiting cases of the Hartmann-Tran profile (HTP) with shape parameters recently measured by cavity ring-down spectroscopy. Significant differences highlighting the importance of including speed-dependence effects in forward models were observed both in scaling retrievals as well as profile retrievals. However, because of uncertainties associated with the atmospheric temperature gradient, it was not possible to extract accurate carbon column profiles even with inclusion of these new line shape parameters and advanced line profiles.

54 ENVIRONMENTAL SCIENCES↗

Geology and elevation shape bacterial assembly in Antarctic endolithic communities

Ice free areas of continental Antarctica are among the coldest and driest environments on Earth, and yet, they support surprisingly diverse and highly adapted microbial communities. Endolithic growth is one of the key adaptations to such extreme environments and often represents the dominant life-form. Despite growing scientific interest, little is known of the mechanisms that influence the assembly of endolithic microbiomes across these harsh environments. Here, we used metagenomics to examine the diversity and assembly of endolithic bacterial communities across Antarctica within different rock types and over a large elevation range. While granite supported richer and more heterogeneous communities than sandstone, elevation had no apparent effect on taxonomic richness, regardless of rock type. Conversely, elevation was clearly associated with turnover in community composition, with the deterministic process of variable selection driving microbial assembly along the elevation gradient. The turnover associated with elevation was modulated by geology, whereby for a given elevation difference, turnover was consistently larger between communities inhabiting different rock types. Overall, selection imposed by elevation and geology appeared stronger than turnover related to other spatially-structured environmental drivers. Our findings indicate that at the cold-arid limit of life on Earth, geology and elevation are key determinants of endolithic bacterial heterogeneity. This also suggests that warming temperatures may threaten the persistence of such extreme-adapted organisms.

54 ENVIRONMENTAL SCIENCES↗

Making a Superbolt: Reconciling Observations of the Optically Brightest Lightning on Earth From Different Satellites

We previously documented geographic distributions of the optically brightest lightning on Earth—known as “superbolts”—using two space-based instruments: the photodiode detector (PDD) on the Fast On-orbit Recording of Transient Events (FORTE) satellite and the Geostationary Lightning Mapper (GLM) on NOAA's Geostationary Operational Environmental Satellites. In this study, we further examine the superbolts identified by the PDD and GLM to reconcile the differences between their geographic distributions. We find that both the physical extent of the parent flash and the development speed of its leaders are important for making a superbolt. The oceanic PDD superbolts tend to occur early in flashes that rapidly expand laterally into long horizontal “megaflashes.” The top GLM superbolts occur over land at later times in particularly large megaflashes. These land-based flashes grow more slowly until they extend over multiple hundreds of kilometers. The FORTE PDD missed these delayed superbolts due to limitations in its triggering. Coincident Tropical Rainfall Measuring Mission measurements show that the warm season megaflash superbolts detected by Lightning Imaging Sensor/GLM and wintertime oceanic superbolts observed by the PDD occur in otherwise similar thunderstorm environments. Both are marked by: low storm heights (<10 km), widespread precipitation near the surface, small infrared brightness temperature gradients, and low flash rates. We suggest that the vertically compact, stratiform nature of these clouds provides favorable conditions for superbolt production.

54 ENVIRONMENTAL SCIENCES↗

Identifying geologic characteristics and operational decisions to meet global carbon sequestration goals

Geologic carbon sequestration is the process of injecting and storing CO 2 in subsurface reservoirs and is an essential technology for global environmental security (e.g., climate change mitigation) and economic security (e.g., CO 2 tax credits). To meet energy, economic, and environmental goals, society will have to identify vast volumes of high-capacity, low-cost, and viable storage reservoirs for sequestering CO 2 . In turn, this requires understanding how major geologic characteristics (such as reservoir depth, thickness, permeability, porosity, and temperature) and design and operational decisions (such as injection well spacing) impact CO 2 injection rates, storage capacity, and economics. Although many numerical simulation tools exist, they cannot repeat the required thousands or millions of simulations to identify ideal reservoir properties and the sensitivity and interaction between geologic parameters and operational decisions. Here, we use SCO 2 T—a fast-running, reduced-order modeling framework—to explore the sensitivity of major geologic parameters and operational decisions to engineering (CO 2 injection rates, plume dimensions, and storage capacities and effectiveness) and costs. Our results show, for the first time, benefits and impacts such as allowing CO 2 plumes to overlap, how different well spacing patterns affect CO 2 sequestration, the effects on costs of including brine treatment and disposal, and the effect of restricting injection rates to 1 MtCO 2 per y based on well limitations. We reveal multiple novel and unintuitive findings including: (i) deeper reservoirs have reduced carbon sequestration costs until injection rates reach 1 MtCO 2 per y, at which point deeper reservoirs become more expensive, (ii) thicker formations allow for increased injection rates and storage capacity, but thickness barely impacts plume areas, (iii) higher geothermal gradients result in reduced sequestration costs, unless brine treatment/disposal costs are included, at which point reservoirs having lower geothermal gradients are more economical because they produce less brine for each unit of injected CO 2 , and (iv) allowing plumes to overlap has a significantly positive impact of increasing storage capacities but has only a small influence on reducing sequestration costs. Altogether, our results illustrate new scientific conclusions to help identify suitable sites to inject and store CO 2 , to help understand the complex interaction between geology and resulting costs, and to help support the pursuit of meeting global sequestration targets.

58 GEOSCIENCES↗

Finite system size effects on drift wave stability

Unstable electrostatic resistive modes, driven by density gradients, are identified in a bounded sheared slab. The boundary conditions play a crucial role and are shown to change the nature of the problem, which is related to so called “universal” mode. The dispersion relation and the structure of the eigenmodes of the instability are derived and are shown to depend on a limited set of dimensionless parameters. The occurrence and possible impact of these modes on numerical simulations and actual plasmas are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structural properties of small rhodium clusters

We report a systematic study of the structural properties of rhodium clusters at the atomistic level. A novel global-minimum search algorithm, known as parallel tempering multicanonical basin hopping plus genetic algorithm (PTMBHGA), is used to obtain the geometrical structures with lowest minima at the semi-empirical level where Gupta potential is used to describe the atomic interaction among the rhodium atoms. These structures are then re-optimized at the density functional theory (DFT) level with exchange-correlation energy approximated by Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA). The structures are optimized for different spin multiplicities. The ones with lowest energies will be taken as ground-state structures. In most cases, we observe only minor changes in the geometry and bond length of the clusters as a result of DFT-level re-optimization. Only in some limited cases, the initial geometries obtained from the PTMBHGA are modified by the re-optimization. The variation of structural properties, such as ground-state geometry, symmetry and binding energy, with respect to the cluster size is studied and agreed well with other results available in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Generalized kinetic equation for tokamak plasma equilibrium distribution function

A generalised kinetic equation for the equilibrium distribution function in a finite beta, arbitrary tokamak plasma is derived. The equation is correct to second order in ρ/L (ρ is the particle Larmor radius and L is the system size). Resolving finite Larmor radius length scales with no restriction on the ratio of poloidal to total equilibrium magnetic field, B $\vartheta$ /B, it generalises the drift kinetic theory of [Hazeltine Phys. Plasmas 15 (1973) 77] to the limit of B $\vartheta$ /B ~ 1 (e.g. to ensure validity for spherical tokamaks). Two cases are considered. The first provides the equilibrium distribution function, consistent with the generalised gyrokinetic formalism of [Dudkovskaia et al. Plasma Phys. Control. Fusion 65 (2023) 045010], derived specifically to capture neoclassical equilibrium currents in gyrokinetic stability analyses in strong gradient regions. The second assumes short length scales in the direction perpendicular to the magnetic field, which can occur as a result of small coherent magnetic structures in the plasma, such as neoclassical tearing mode magnetic islands close to threshold. This then extends the drift island equations of [Dudkovskaia et al. Nucl. Fusion 63 (2023) 016020] for the plasma response to magnetic islands to a spherical tokamak plasma configuration. In conclusion, resolving ρ ~ ρ $\vartheta$ (or B $\vartheta$ ~ B), where ρ $\vartheta$ is the particle poloidal Larmor radius, is also expected to influence calculations of the magnetic island propagation frequency and the associated contributions to the island onset conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Topological Guided Detection of Extreme Wind Phenomena: Implications for Wind Energy

Extreme wind phenomena play a crucial role in the efficient operation of wind farms for renewable energy generation. However, existing detection methods are computationally expensive, limited to specific coordinate. In real-world scenarios, understanding the occurrence of these phenomena over a large area is essential. Therefore, there is a significant demand for a fast and accurate approach to forecast such events. In this paper, we propose a novel method for detecting wind phenomena using topological analysis, leveraging the gradient of wind speed or critical points in a topological framework. By extracting topological features from the wind speed profile within a defined region, we employ topological distance to identify extreme wind phenomena. Our results demonstrate the effectiveness of utilizing topological features derived from regional wind speed profiles. We validate our approach using high-resolution simulations with the Weather Research and Forecasting model (WRF) over a month in the US East Coast.

MATHEMATICS AND COMPUTING,WIND ENERGY↗

Surface Water Quality Data from Beaver-Impacted Streams; Trail Creek and East River, Colorado 2025

This data package contains surface water chemistry measurements collected in 2025 to evaluate how beaver damming and low-tech process-based stream restoration influence water quality and metal mobility in mountainous headwater systems of the Upper Colorado River Basin. Sampling was conducted at Trail Creek (Taylor Park watershed, Colorado), a tributary undergoing restoration through installation of low-tech process-based structures (i.e., beaver dam analogs), and at off-channel beaver ponds within the East River floodplain (East River watershed, Colorado). Samples were collected along longitudinal transects spanning upstream control reaches, beaver-influenced ponded reaches, and downstream segments. Additional samples were collected from near-surface pore waters within a beaver dam seepage face. The dataset includes concentrations of major and trace elements measured by inductively coupled plasma–mass spectrometry (ICP-MS) and inductively coupled plasma–optical emission spectrometry (ICP-OES), major anions measured by ion chromatography (IC), and dissolved organic carbon (DOC; reported as non-purgeable organic carbon, NPOC). Samples were size-fractionated at 0.45 micrometers (µm), 0.22 µm, and 0.02 µm to distinguish particulate (>0.45 µm), colloidal (0.22–0.02 µm), and dissolved (<0.02 µm) fractions. The data package consists of comma-separated value (.csv) files containing tabulated chemical concentration data, sample metadata (site identifiers, geographic coordinates, sampling dates, fraction type), and quality control flags. All files are provided in open, non-proprietary formats that can be accessed using standard data analysis software such as Microsoft Excel, R, Python, MATLAB, or other programs capable of reading .csv files. Units, detection limits, and analytical methods are documented in accompanying metadata files. The dataset is designed to support analyses of (1) how beaver impoundment and restoration structures alter elemental partitioning and transport, (2) the role of iron and organic carbon in mediating trace metal mobility, and (3) reach-scale changes in water quality across restoration gradients. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. Part of this work was performed at SLAC Accelerator Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-76SF00515.

Anions↗

Pressure Barrier Applicability to Polymer Flood Design

Effective oil displacement from a reservoir requires adequate and properly directed pressure gradients in areas of high oil saturation. If the polymer bank is too large or too viscous during a polymer flood, the pressure drop from the injection well to the polymer front may act as a pressure barrier by usurping most of the downstream driving force for oil displacement. Polymer injection pressures must be limited. The maximum allowable injection pressure is commonly constrained by caprock integrity, injection equipment, and/or regulations, even though fractures can be beneficial to polymer injectivity (and even sweep efficiency in some cases). This paper examines when the pressure-barrier concept limits the size and viscosity of the polymer bank during a polymer flood. Both analytical and numerical methods are used to address this issue. We examine the relevance of the pressure barrier concept for a wide variety of circumstances, including oil viscosities ranging from 10-cp (like at Daqing, China) to 1650-cp (like at Pelican Lake, Alberta), vertical wells (like at Tambaredjo, Suriname) versus horizontal wells (like at Milne Point, Alaska), single versus multiple layered reservoirs, permeability contrast, and with versus with crossflow between layers. We also examine the relation between the pressure-barrier concept and fractures and fracture extension during polymer injection. We demonstrate that in reservoirs with single layers, the pressure-barrier concept only limits the optimum viscosity of the injected polymer if the mobility of the polymer bank is less than the mobility of the displaced oil bank. The same is true for multi-zoned reservoirs with no crossflow between layers. Thus, for these cases, the optimum polymer viscosity is likely to be dictated by the mobility of the oil bank, unless other factors (like fracture extension) intervene. For multi-zoned reservoirs with free crossflow between layers, the situation is different. A compromise must be reached between injected polymer viscosity and the efficiency of oil recovery. The relevance of our findings is applied to operations for several existing polymer floods. This work is particularly relevant to viscous-oil reservoirs (like Pelican Lake and others) where the injected polymer viscosities are substantially lower than the oil viscosity

Wang, Dongmei↗

Utilization of Synthetic Near-Infrared Spectra via Generative Adversarial Network to Improve Wood Stiffness Prediction

Near-infrared (NIR) spectroscopy is widely used as a nondestructive evaluation (NDE) tool for predicting wood properties. When deploying NIR models, one faces challenges in ensuring representative training data, which large datasets can mitigate but often at a significant cost. Machine learning and deep learning NIR models are at an even greater disadvantage because they typically require higher sample sizes for training. In this study, NIR spectra were collected to predict the modulus of elasticity (MOE) of southern pine lumber (training set = 573 samples, testing set = 145 samples). To account for the limited size of the training data, this study employed a generative adversarial network (GAN) to generate synthetic NIR spectra. The training dataset was fed into a GAN to generate 313, 573, and 1000 synthetic spectra. The original and enhanced datasets were used to train artificial neural networks (ANNs), convolutional neural networks (CNNs), and light gradient boosting machines (LGBMs) for MOE prediction. Overall, results showed that data augmentation using GAN improved the coefficient of determination (R 2 ) by up to 7.02% and reduced the error of predictions by up to 4.29%. ANNs and CNNs benefited more from synthetic spectra than LGBMs, which only yielded slight improvement. All models showed optimal performance when 313 synthetic spectra were added to the original training data; further additions did not improve model performance because the quality of the datapoints generated by GAN beyond a certain threshold is poor, and one of the main reasons for this can be the size of the initial training data fed into the GAN. LGBMs showed superior performances than ANNs and CNNs on both the original and enhanced training datasets, which highlights the significance of selecting an appropriate machine learning or deep learning model for NIR spectral-data analysis. The results highlighted the positive impact of GAN on the predictive performance of models utilizing NIR spectroscopy as an NDE technique and monitoring tool for wood mechanical-property evaluation. Further studies should investigate the impact of the initial size of training data, the optimal number of generated synthetic spectra, and machine learning or deep learning models that could benefit more from data augmentation using GANs.

59 BASIC BIOLOGICAL SCIENCES↗

Diffusion flame in homologous turbulent shear flows.

The simplifed statistical theory developed previously is employed to analyze the equilibrium and near-equilibrium combustion of initially unmixed reactants. It is found that the flame zone in the limit of large Damkohler numbers is very thick and is of the order of the local integral scale of turbulence. This is in contrast to the existing phenomenological theories which predict the infinitesimally thin flame sheet, in the same limit, as it is with the laminar diffusion flame. Qualitative agreements with the available experimental results are shown. It is found that singularities exist at the edges of the flame which are removed as Damkohler number is reduced. Also, it is found that the heat transfer may take place against the local mean temperature gradient in certain regions within the flame.

Chung, P. M.↗

Effects of vane-induced rotation on diffuser flow distortion in an axisymmetric mixed-compression inlet

An investigation of vane-induced flow rotation to modify distorted steady-state total-pressure patterns in the subsonic diffuser of a supersonic mixed-compression inlet was conducted. Radial static-pressure gradients generated by the rotation was the mechanism used to modify the total-pressure distributions. Significant redistribution of circumferential distortion patterns into more compatible radial patterns was realized, but flow problems near the duct walls reduced the general effectiveness of the technique. Total-pressure losses associated with the swirl vanes were slight. Limited turbulence data indicated that vane istallation resulted in reduced turbulence levels.

Meleason, E. T.↗