Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak scaling”

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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Alternating Direction Decomposition with Strong Bounding and Convexification (ADDSBC) for Solving Security Constrained AC Unit Commitment Problems

This project aims to develop efficient and robust computational methods for solving the security-constrained unit commitment and alternating current optimal power flow problem (SC-UC-ACOPF). The SC-UC-ACOPF problem is at the center of the short-term operation of the U.S. Power Grid. It is solved every week, every day, and every 10 minutes to plan for the optimal action of electricity generation and consumption by minimizing the generation cost and maintaining power system reliability against potential disruptions of equipment failures. In mathematical terms, SC-UC-ACOPF is a challenging large-scale mixed-integer nonlinear optimization model. This means that the decisions involve both discrete variables, e.g. the turning on and off of generators and switching of transmission lines and transformers, and continuous decisions, e.g. the amount of energy generated by each generator and the power flows in the power grid. The physics of the power flow is described by nonlinear equations involving real and reactive power and bus voltages. Another key feature is the large number of contingencies, i.e. the system needs to stay reliable in face of failure of any one equipment, such as transmission lines and generators. The U.S. power grids are extremely complicated and large scale with more than 5,000 generators, 50,000 buses, and 100,000 high-voltage transmission lines, making the SC-UC-ACOPF a very large-scale computation challenge. The research developed in this project aims to solve the SC-UC-ACOPF problems in the three timescales, i.e. weekly, daily, and every 10-min. The proposed computational methods are built on a principled algorithmic approach of decomposition and penalization. More specifically, the algorithm develops spatial and temporal decomposition by exploiting the strong temporal coupling and weak spatial coupling of the UC problem and the complementary feature, i.e. weak temporal coupling and strong spatial coupling of the ACOPF problem. The algorithm also leverages recent progresses in strong convex relaxation of ACOPF. A unique feature of the proposed approach is that it generates a valid, global upper bound on the optimal maximum profit. In this way, a global optimality gap is available to measure the quality of the solution. To further speed up computation, the research team has developed a plethora of effective heuristics to strengthen the iterative penalty-based decomposition framework. For instance, a heuristic is developed to construct inner approximations of the time coupling constraints within the time decoupled problems. Contingencies are pre-screened and low-rank matrix computation is exploited to find the almost unique solution to each contingency. A novel heuristic for line switching is proposed and tested with positive impacts on instances where line switching is beneficial. Taking a systematic approach and carefully handling every detail of the problem pays off. The TIM-GO’s performance throughout the trials and the final event was stellar. TIM-GO garnered the second highest total prize money and is ranked in the top three positions across all categories of comparison.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Weak Form Scientific Machine Learning: Test Function Construction for System Identification

Weak form Scientific Machine Learning (WSciML) is a recently developed framework for data-driven modeling and scientific discovery. It leverages the weak form of equation error residuals to provide enhanced noise robustness in system identification via convolving model equations with test functions, reformulating the problem to avoid direct differentiation of data. The performance, however, relies on wisely choosing a set of compactly supported test functions. In this work, we mathematically motivate a novel data-driven method for constructing Single-scale-Local reference functions for creating the set of test functions. Our approach numerically approximates the integration error introduced by the quadrature and identifies the support size for which the error is minimal, without requiring access to the model parameter values. Through numerical experiments across various models, noise levels, and temporal resolutions, we demonstrate that the selected supports consistently align with regions of minimal parameter estimation error. We also compare the proposed method against the strategy for constructing Multi-scale-Global (and orthogonal) test functions introduced in our prior work, demonstrating the improved computational efficiency.

FOS: Computer and information sciences↗
NASA NTRS1982

Mesospheric water vapor

Water vapor in the earth's mesosphere has been observed at the frequency of 22.235 GHz as an absorption against the sun by utilizing a ground-based radio telescope. The H2O mixing ratio of 4.3 + or - 1 ppm is obtained in the height 40-70 km for a constant distribution model. The data are equally well represented by a photochemical model of Crutzen (1974), which was scaled by a factor of 0.93 yielding a peak mixing ratio of 6 + or - 1 ppm at 55 km. These results are 3-day averages during November 1979, where the averaging periods cover the about + or - 4 hours of the sun's transit of our meridian. The measurements are insensitive to H2O above 70 km, where the absorption is very weak. It is not possible to estimate the mixing ratios below 50 km, since the pressure broadening in the line exceeds the bandwidth of our spectrometer at the lower altitudes.

Deguchi, S.↗
NASA NTRS1987

Magnetic braking in weakly ionized media

The combined magnetic braking-ambipolar diffusion problem in weakly ionized, rigidly rotating disks is studied. An analytical solution is presented for a disk whose angular velocity and magnetic yield vectors are aligned with the symmetry axis, illustrating the effects of the relative azimuthal drift of neutrals and ions. The effects of radial drift are added, commenting on the ratio of the characteristic ambipolar diffusion and magnetic braking time scales in high-mass and low-mass disks. A numerical calculation is used to show the combined action of these two processes.

Konigl, Arieh↗
NASA NTRS1993

Reactive nitrogen and its correlation with ozone in the lower stratosphere and upper troposphere

High resolution in situ measurements of reactive nitrogen (NO(y)) and O3 were made in the upper troposphere and lower stratosphere at a variety of latitudes and seasons. In the lower stratosphere, NO(y) and O3 are very highly positively correlated at all times and spatial scales sampled. The ratio NO(y)/NO3 is much less variable than either species measured separately. The ratio has a much weaker gradient with altitude than the mixing ratios of O3 or NO(y). The ratio is smaller and decreases more rapidly with altitude in the tropics than at midlatitudes. In the upper troposphere NO(y) and O3 are only weakly correlated. Their ratio in the tropical upper troposphere is about 0.005-0.025 and the ratio in the midlatitude upper troposphere is about 0.004-0.010. The NO(y) in the upper troposhere is probably partly due to lightning.

Murphy, D. M.↗
NASA NTRS1996

Turbulent Compressible Convection with Rotation: Flow Structure and Evolution - Part 1

The effects of Coriolis forces on compressible convection are studied using three-dimensional numerical simulations carried out within a local modified f-plane model. The physics is simplified by considering a perfect gas occupying a rectilinear domain placed tangentially to a rotating sphere at various latitudes, through which a destabilizing heat flux is driven. The resulting convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers, evaluating conditions where the influence of rotation is both weak and strong. Given the computational demands of these high-resolution simulations, the parameter space is explored sparsely to ascertain the differences between laminar and turbulent rotating convection. The first paper in this series examines the effects of rotation on the flow structure within the convection, its evolution, and some consequences for mixing. Subsequent papers consider the large-scale mean shear flows that are generated by the convection, and the effects of rotation on the convective energetics and transport properties. It is found here that the structure of rotating turbulent convection is similar to earlier nonrotating studies, with a laminar, cellular surface network disguising a fully turbulent interior punctuated by vertically coherent structures. However, the temporal signature of the surface flows is modified by inertial motions to yield new cellular evolution patterns and an overall increase in the mobility of the network. The turbulent convection contains vortex tubes of many scales, including large-scale coherent structures spanning the full vertical extent of the domain involving multiple density scale heights. Remarkably, such structures align with the rotation vector via the influence of Coriolis forces on turbulent motions, in contrast with the zonal tilting of streamlines found in laminar flows. Such novel turbulent mechanisms alter the correlations which drive mean shearing flows and affect the convective transport properties. In contrast to this large-scale anisotropy, small-scale vortex tubes at greater depths are randomly orientated by the rotational mixing of momentum, leading to an increased degree of isotropy on the medium to small scales of motion there. Rotation also influences the thermodynamic mixing properties of the convection. In particular, interaction of the larger coherent vortices causes a loss of correlation between the vertical velocity and the temperature leaving a mean stratification which is not isentropic.

Brummell, Nicholas H.↗
NASA NTRS1998

Computer Modeling of the Thermal Conductivity of Cometary Ice

The thermal conductivity was found to be only weakly dependent on the microstructure of the amorphous ice. In general, the amorphous ices were found to have thermal conductivities of the same order of magnitude as liquid water. This is in contradiction to recent experimental estimates of the thermal conductivity of amorphous ice, and it is suggested that the extremely low value obtained experimentally is due to larger-scale defects in the ice, such as cracks, but is not an intrinsic property of the bulk amorphous ice.

Bunch, Theodore E.↗
NASA NTRS1998

Computer Modeling of the Thermal Conductivity of Cometary Ice

The main objective of this research was to estimate the thermal conductivity of cometry ices from computer simulations of model amorphous ices. This was divided into four specific tasks: (1) Generating samples of amorphous water ices at different microporosities; (2) Comparing the resulting molecular structures of the ices with experimental results, for those densities where data was available; (3) Calculating the thermal conductivities of liquid water and bulk amorphous ices and comparing these results with experimentally determined thermal conductivities; and (4) Investigating how the thermal conductivity of amorphous ice depends upon the microscopic porosity of the samples. The thermal conductivity was found to be only weakly dependent on the microstructure of the amorphous ice. In general, the amorphous ices were found to have thermal conductivities of the same order of magnitude as liquid water. This is in contradiction to recent experimental estimates of the thermal conductivity of amorphous ice, and it is suggested that the extremely low value obtained experimentally is due to larger-scale defects in the ice, such as cracks, but it is not an intrinsic property of the bulk amorphous ice.

Bunch, Theodore E.↗
NASA NTRS2001

Structure of the Jovian Magnetodisk Current Sheet: Initial Galileo Observations

The ten-degree tilt of the Jovian magnetic dipole causes the magnetic equator to move back and forth across Jupiter's rotational equator and tile Galileo orbit that lies therein. Beyond about 24 Jovian radii, the equatorial current sheet thins and tile magnetic structure changes from quasi-dipolar into magnetodisk-like with two regions of nearly radial but antiparallel magnetic field separated by a strong current layer. The magnetic field at the center of the current sheet is very weak in this region. Herein we examine tile current sheet at radial distances from 24 55 Jovian radii. We find that the magnetic structure very much resembles tile structure seen at planetary magnetopause and tail current sheet crossings. Tile magnetic field variation is mainly linear with little rotation of the field direction, At times there is almost no small-scale structure present and the normal component of the magnetic field is almost constant through the current sheet. At other times there are strong small-scale structures present in both the southward and northward directions. This small-scale structure appears to grow with radial distance and may provide the seeds for tile explosive reconnection observed at even greater radial distances oil tile nightside. Beyond about 40 Jovian radii, the thin current sheet also appears to be almost constantly in oscillatory motion with periods of about 10 min. The amplitude of these oscillations also appears to grow with radial distance. The source of these fluctuations may be dynamical events in tile more distant magnetodisk.

Russell, C. T.↗
DOE OSTI2020

CCCP and MENeaCS: (updated) weak-lensing masses for 100 galaxy clusters

ABSTRACT Large area surveys continue to increase the samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at 0.05 < z < 0.15, observed as part of the Multi Epoch Nearby Cluster Survey, and re-evaluate the mass estimates for 52 clusters from the Canadian Cluster Comparison Project. Updated high-fidelity photometric redshift catalogues of reference deep fields are used in combination with advances in shape measurements and state-of-the-art cluster simulations, yielding an average systematic uncertainty in the lensing signal below 5 per cent, similar to the statistical uncertainty for our cluster sample. We derive a scaling relation with Planck measurements for the full sample and find a bias in the Planck masses of 1 − b = 0.84 ± 0.04 (stat) ±0.05 (syst). We find no statistically significant trend of the mass bias with redshift or cluster mass, but find that different selections could change the bias by up to 0.07. We find a gas fraction of 0.139 ± 0.014 (stat) for eight relaxed clusters in our sample, which can also be used to infer cosmological parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Monsoon Depression Amplification by Horizontal Shear and Humidity Gradients: A Shallow Water Perspective

Transient, synoptic-scale vortices produce a large fraction of total rainfall in most monsoon regions and are often associated with extreme precipitation. However, the mechanism of their amplification remains a topic of active research. For monsoon depressions, which are the most prominent synoptic-scale vortex in the Asian–Australian monsoon, recent work has suggested that meridional gradients in zonal wind in the vortex environment may produce growth through barotropic instability, while meridional gradients in environmental humidity have also been proposed to cause amplification through coupling with precipitating convection. Here, a two-dimensional shallow water model on a sphere with parameterized precipitation is used to examine the relative role played by these two environmental gradients. By systematically varying the meridional moisture gradient and meridional wind shear for both weak, quasi-linear waves and finite-amplitude isolated vortices, we show that rotational winds in the initial vortex are amplified most strongly by meridional shear of the environmental zonal wind, while vortex precipitation rates are most sensitive to environmental moisture gradients. The growth rate in the presence of both gradients is less than the sum of growth rates in the presence of isolated gradients, as the phase relation between moisture and vorticity anomalies becomes distorted with increasing shear. In conclusion, these results suggest that background meridional gradients in both zonal wind and environmental humidity can contribute to the amplification of vortices to monsoon depression strength, but with some degree of decoupling of the dry rotational flow and the moist convection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Linear simulation of kinetic electromagnetic instabilities in a tokamak plasma with weak magnetic shear

Gyrokinetic simulation and eigenvalue calculation of electromagnetic instabilities are carried out for an experimentally observed low-n mode in weak magnetic shear discharge. With different magnetic shear values, the Gyrokinetic Toroidal Code simulation of the ion temperature gradient mode is consistent with the eigenvalue code (HD7) calculation. Due to the sensitivity of the kinetic ballooning mode (KBM) to global equilibrium, the simulation of the KBM deviates from the eigenvalue results, for the ballooning representation used in HD7 satisfies its spatial scale separations. Under a flat safety factor profile, the KBM is more unstable and its mode structure tends to move with the peak of the ion temperature drive. Further simulation of the KBM in an HL-2A-like equilibrium shows that the β excitation threshold of the mode is lower than 0.2% and the dominant toroidal mode number is n = 4, which is consistent with the measured experimental spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗
DOE OSTI2021

Smart Inverter Stability Enhancement in Weak Grids Using Adaptive Virtual-Inductance

An adaptive virtual-inductance feedforward scheme is presented in this article to enhance the stable operation of smart inverters in weak grids. In the developed scheme, the virtual inductance is updated in real time according to an adaptation law to ensure that the inverter follows the dynamic response of a robust reference model, thereby ensuring stability without the need for any grid impedance estimation. In this article, the developed scheme is validated with experimental results obtained from testing a small-scale two-level, 208-V, 2-kW inverter.

42 ENGINEERING↗
NASA NTRS2012

Magnetic Field Strength Fluctuations in the Heliosheath: Voyager 1 Observations During 2009

We analyze the microscale fluctuations of the magnetic field strength B on a scale of several hours observed by Voyager1 (V1) in the heliosheath during 2009. The microscale fluctuations of B range from coherent to stochastic structures. The amplitude of microscale fluctuations of B during 1 day is measured by the standard deviation (SD) of 48 s averages of B. The distribution of the daily values of SD is lognormal. SD(t) from day of year (DOY) 1 to 331, 2009, is very intermittent. SD(t) has a 1/f or "pink noise" spectrum on scales from 1 to 100 days, and it has a broad multifractal spectrum f(alpha) with 0.57 less than or equal to alpha less than or equal to 1.39. The time series of increments SD(t + tau) -- SD(t) has a pink noise spectrum with alpha(1) = 0.88 plus or minus 0.14 on scales from 1 to 100 days. The increments have a Tsallis (q-Gaussian) distribution on scales from 1 to 165 days, with an average q = 1.75 plus or minus 0.12. The skewness S and kurtosis K have Gaussian and lognormal distributions, respectively. The largest spikes in K(t) and S(t) are often associated with a change in B across a data gap and with identifiable physical structures. The "turbulence" observed by V1 during 2009 was weakly compressible on average but still very intermittent, highly variable, and highly compressible at times. The turbulence observed just behind the termination shock by Voyager 2 was twice as strong. These observations place strong constraints on any model of turbulence in the heliosheath.

Brulaga, L. F.↗
NASA NTRS2016

Warming of the Global Ocean: Spatial Structure and Water-Mass Trends

This study investigates the multidecadal warming and interannual-to-decadal heat content changes in the upper ocean (0-700 m), focusing on vertical and horizontal patterns of variability. These results support a nearly monotonic warming over much of the World Ocean, with a shift toward Southern Hemisphere warming during the well-observed past decade. This is based on objectively analyzed gridded observational datasets and on a modeled state estimate. Besides the surface warming, a warming climate also has a subsurface effect manifesting as a strong deepening of the midthermocline isopycnals, which can be diagnosed directly from hydrographic data. This deepening appears to be a result of heat entering via subduction and spreading laterally from the high-latitude ventilation regions of subtropical mode waters. The basin-average multidecadal warming mainly expands the subtropical mode water volume, with weak changes in the temperature-salinity (u-S) relationship (known as ''spice'' variability). However, the spice contribution to the heat content can be locally large, for example in Southern Hemisphere. Multidecadal isopycnal sinking has been strongest over the southern basins and weaker elsewhere with the exception of the Gulf Stream/North Atlantic Current/subtropical recirculation gyre. At interannual to decadal time scales, wind-driven sinking and shoaling of density surfaces still dominate ocean heat content changes, while the contribution from temperature changes along density surfaces tends to decrease as time scales shorten.

Global Ocean↗
DOE OSTI2025

Hydroxy-substituted electron deficient Pd porphyrin cofactors illuminate ultrafast proton transfer reactions

Coordinated electron and proton movement drive bioenergetic functions. Relative to electron transfer reactions, tracking proton transport over fast-to-ultrafast time scales is challenging. Optical resolution of proton transfer dynamics can take advantage of chromophoric photoacids that not only trigger proton migration upon photoexcitation, but produce distinct spectroscopic changes associated with protonation/deprotonation. Here, in this work, we report the design of a hydroxy-substituted electron deficient Pd porphyrin, PPd(C 6 F 5 ) 3 OH; upon photoexcitation, the acidity constant of this weak acid (pKa = 6.49) dramatically drops (pKa* = 0.97). Electronic excitation of PPd(C 6 F 5 ) 3 OH triggers an ultrafast proton transfer reaction (PPd(C 6 F 5 ) 3 OH +:B + hυ → 1[PPd(C 6 F 5 ) 3 OH]* +:B → 1 [PPd(C 6 F 5 ) 3 O]−* + HB + ; τ PT = 342 fs) to a H-bonded base (B) in solution. Both PPd(C 6 F 5 ) 3 OH and its conjugate anion PPd(C 6 F 5 ) 3 O − exhibit distinct vis-NIR spectral features for their respective ground and excited states. Because the electroni-cally excited triplet lifetimes of these species exceed tens of microseconds, the PPd(C 6 F 5 ) 3 OH photoacid defines an ideal cofactor to probe light-triggered proton release and track long-range proton migration in protein environments.

Palladium↗
DOE OSTI2022

Energy-Storing Cryogenic Carbon Capture™ for Utility and Industrial-scale Processes (Final Report)

This project leverages work completed under funding from previous U.S. Department of Energy (DOE) projects and other sources of funding around the Cryogenic Carbon Capture Process™ (CCC). The purpose of this project is to further develop a perturbation of the CCC process that stores energy called the Cryogenic Carbon Capture™ Energy Storing (CCC-ES) process. The successful completion of the tasks in this project has prepared the CCC-ES process for scaling and developed initial plans and a feasibility analysis for an engineering-scale or small-commercial scale system that can store at least 10 MWh of energy. The specific areas of work in this project are: (1) Develop a conceptual, site-specific integrated process flow diagram: Sustainable Energy Solutions (SES) developed a conceptual, site-specific plan including a process flow diagram for integrating CCC-ES at a specific location. The proposed system was based on a typical engineering-scale CCC system designed to capture 30-50 metric tonnes of CO 2 per day (TPD), with modifications and additional changes and equipment put in place to add the energy storage element to the system. The 30-50 TPD scale corresponds to the required 10 MWh of stored energy. This task includes updating simulation software for transient analysis. SES also coordinated with the Jim Bridger Power Plant in Point of Rocks, Wyoming to get site-specific data. Using these data, SES completed site-specific preliminary process flow diagrams, estimated capital and operating costs for the CCC-ES system, and used the upgraded software package to estimate the primary figures of merit. (2) Complete a technoeconomic analysis (TEA): SES completed a full-scale TEA of the CCC-ES system to meet specific success criteria. This TEA used process information from the host site selected above. (3) Technology and data gap assessment for the CCC-ES process: As part of the technology gap assessment, SES performed an assessment of the best alternative energy storage solutions including figures of merit and operational strengths and weaknesses. SES then compared the proposed system with a discussion of how it addresses the weaknesses of alternative systems, showing for example a cost of stored energy at less than $50/MWh and with a roundtrip efficiency of greater than 95%. This assessment includes technical and other risk analyses and technology gaps with research and development needs to commercialize the process by 2030. Finally, work includes a commercialization roadmap/development pathway outlining the major milestones to move to a commercially operating CCC-ES system. (4) Phase II pre-FEED project plan: SES completed a project plan for the Phase II pre-FEED analysis using the same host site as in the Phase I conceptual study.

03 NATURAL GAS↗