SEARCH · Search NASA

Results for “time scale separation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

NASA NTRS2018

Magnetic Flux Cancelation as the Buildup and Trigger Mechanism for CME-Producing Eruptions in Two Small Active Regions

We follow two small, magnetically isolated CME-producing solar active regions (ARs) from the time of their emergence until several days later, when their core regions erupt to produce the CMEs. In both cases, magnetograms show: (a) following an initial period where the poles of the emerging regions separate from each other, the poles then reverse direction and start to retract inward; (b) during the retraction period, flux cancelation occurs along the main neutral line of the regions, (c) this cancelation builds the sheared core field/flux rope that eventually erupts to make the CME. In the two cases, respectively 30% and 50% of the maximum flux of the region cancels prior to the eruption. Recent studies indicate that solar coronal jets frequently result from small-scale filaments eruptions, with those "minifilament" eruptions also being built up and triggered by cancelation of magnetic flux. Together, the small-AR eruptions here and the coronal jet results suggest that isolated bipolar regions tend to erupt when some threshold fraction, perhaps in the range of 50%, of the region's maximum flux has canceled. Our observed erupting filaments/flux ropes form at sites of flux cancelation, in agreement with previous observations. Thus, the recent finding that minifilaments that erupt to form jets also form via flux cancelation is further evidence that minifilaments are small-scale versions of the long-studied full-sized filaments.

Sterling, Alphonse C.↗

DOE OSTI2026

Denudation, solute export, landscape evolution modeling, and geographic information system data for the East River watershed, Colorado, USA (2020-2024)

This data package contains geographic information system (GIS) layers and tabular datasets associated with the study of lithologic controls on denudation, solute export, carbon-scaling relationships, and transient landscape evolution in the East River watershed near Crested Butte, Colorado, USA. The package includes GIS layers used to produce the Figure 2 map, including drainage, hillshade, lithology, sample locations, and basin polygons, together with comma-separated value (CSV) tables and matching CSV data dictionaries. One group of tables reports sample-level and catchment-level information for river-sediment samples analyzed for in situ-produced cosmogenic beryllium-10 (10Be), including sample names, outlet elevations, geographic coordinates, upstream drainage area, rock-type classes, production-rate scaling scheme, analyzed nuclide, catchment-averaged denudation rates, and associated lower and upper analytical uncertainties. Sample and catchment attributes provide the basis for comparing denudation rates across intrusive, shale, sedimentary, and mixed-lithology settings. A second group of tables reports supporting information for landscape-evolution modeling and the mapped geologic framework of the study area. Included files list parameter values and definitions for the two-phase landscape-evolution simulations, summarize full-domain model erosion fluxes and topographic metrics for different simulation configurations, provide a fixed-area carbon-model scaling table, and summarize mapped geologic units within the East River study domain, including geologic code, formation name, lithologic description, mapped area, and lithologic class grouping. Model outputs and geologic summaries support interpretation of transient landscape behavior and its relation to the mapped distribution of shale, intrusive, sedimentary, and surficial units. A third group of tables reports hydrologic and hydrochemical information used to quantify dissolved export from the watershed. Included files provide site-level values for drainage area, mean annual solute export, standard error of annual export, area-normalized solute yield, and equivalent weathering rate for five East River monitoring sites, along with metadata describing the number, sampling cadence, and date range of discharge records and partial and full total dissolved solids observations used in the solute-yield analyses. The package also contains a supplementary daily ion-load time series with daily mean discharge, discharge observation counts, dissolved concentrations, and daily loads for calcium, magnesium, sodium, potassium, chloride, sulfate, nitrate, fluoride, dissolved silica, charge-balance bicarbonate, and total dissolved solids. The package contains GIS files, comma-separated value files (.csv), CSV data dictionaries, a file-level metadata table, a package-tree text file, and a readme text file.

10Be↗

DOE OSTI2026

Microcanonical Kinetics of Water-Mediated Proton Transfer in 4ABAH + ·(H 2 O) n = 4–6 Clusters (ABA = Aminobenzoic Acid): A Model System for Size-Dependent Relaxation to Ergodic Behavior

Here, we leverage the unique properties of the 4ABAH + · (H 2 O) n clusters (ABA = 4-aminobenzoic acid, n = 4−6) to quantitatively address how a finite, isolated system evolves into an ergodic condition starting from localized arrangements in configuration space. This system adopts two distinct structural isomers in which water molecules cluster around the cationic centers of its two protomers with widely separated positive charge centers. These isomers arise from excess proton attachment to either the acid (O) or amino (N) group on opposite sides of the benzene ring. Both forms are captured and kinetically trapped using cryogenic ion methods and then selectively vibrationally excited through their mutually exclusive IR bands involving NH and OH stretching fundamentals. Because the IR excitation lies below the water binding energy, the system can evolve to explore slow, rare events that lead to the interconversion between the two isomers. The rates of these intracluster reactions are determined by using a pump−probe scheme involving ∼5 ns IR pump and UV probe lasers. The rates occur on the microsecond time scale, leading to steady state populations of the isomers, thus revealing the cluster size-dependent fractionation between the two species at microcanonical equilibrium. The steady state distributions are correlated with the expected trend in the cluster size-dependent reaction energetics, which are in turn consistent with changes in the relative densities of states of the two species. These results thus provide an unusually clear example in which complex, protic-solvent-mediated chemical transformations are captured within a finite system at a precisely determined internal energy.

Rana, Abhijit [Yale Univ., New Haven, CT (United S↗

DOE OSTI2024

Membrane-based Carbon Capture Process Optimization using CFD Modeling

Carbon capture is a promising option to mitigate CO2 emissions from existing coal-fired power plants, cement and steel industries, and petrochemical complexes. Among the available technologies, membrane-based carbon capture presents the lowest energy consumption, operating costs, and carbon footprint. In addition, membrane processes have important operational flexibility and response times. On the other hand, the major challenges to widespread application of this technology are related to reducing capital costs and improving membrane stability and durability. To upscale the technology into stacked flat sheet configurations, high fidelity computational fluid dynamics (CFD) that describes the separation process accurately are required. High fidelity simulations have been shown to be effective in studying the complex transport phenomena in membrane systems. In addition, obtaining high CO2 recovery percentages and product purity re-quires a multi-stage membrane process, where the optimal network configuration of the mem-brane modules must be studied in a systematic way. In order to address the design problem at process scale, we formulate a superstructure for the membrane-based carbon capture, including up to three separation stages. In the formulation of the optimization problem, we include reduced models, based on rigorous CFD simulations of the membrane modules.

Pedrozo, Hector A.↗

DOE OSTI2026

Improving ProtoDUNE pion cross-section measurements with NuGraph Michel-electron tagging

Understanding hadron-argon interactions is essential for precise neutrino energy reconstruction and final-state interaction modeling in liquid-argon time projection chamber (LArTPC) experiments such as DUNE. In particular, pion absorption and charge-exchange processes constitute significant sources of systematic uncertainty in neutrino oscillation measurements. ProtoDUNE-SP, a large-scale LArTPC prototype operated at the CERN Neutrino Platform and exposed to charged-particle test beams in the few-GeV range, enables direct measurements of these processes. This work focuses on the measurement of differential cross sections for pion absorption and charge exchange using the 2 GeV/c pion beam data from the ProtoDUNE-SP run. A key component of this analysis is the identification of Michel electrons from $\pi \rightarrow \mu \rightarrow e$ decay chains, which helps separate different interaction topologies and improves background rejection. Michel electron identification will also assist in reliably calibrating the electromagnetic response in ProtoDUNE-SP data and for the future DUNE detectors. In this analysis, we apply NuGraph to identify Michel electrons. NuGraph is a graph neural network that models detector hits as nodes connected by spatial and temporal edges for particle and topology classification in LArTPC detectors. We first benchmark NuGraph’s Michel electron classification performance using ICEBERG data, a small-scale LArTPC prototype used for DUNE electronics and reconstruction development, and then transfer the approach to ProtoDUNE-SP. This poster presents the analysis strategy, NuGraph-based classification studies, and discusses how these developments are expected to improve the pion cross-section measurement.

Razafinime, Soamasina Herilala [Cincinnati U.] (OR↗

NASA NTRS1992

Full-scale transmission testing to evaluate advanced lubricants

Experimental tests were performed on the OH-58A helicopter main rotor transmission in the NASA Lewis 500 hp helicopter transmission test stand. The testing was part of a lubrication program. The objectives are to develop and show a separate lubricant for gearboxes with improved performance in life and load carrying capacity. The goal was to develop a testing procedure to fail certain transmission components using a MIL-L-23699 based reference oil and then to run identical tests with improved lubricants and show improved performance. The tests were directed at parts that failed due to marginal lubrication from Navy field experience. These failures included mast shaft bearing micropitting, sun gear and planet bearing fatigue, and spiral bevel gear scoring. A variety of tests were performed and over 900 hrs of total run time accumulated for these tests. Some success was achieved in developing a testing procedure to produce sun gear and planet bearing fatigue failures. Only marginal success was achieved in producing mast shaft bearing micropitting and spiral bevel gear scoring.

Lewicki, David G.↗

NASA NTRS1993

A magnetic suspension system with a large angular range

The paper describes a small-scale laboratory system, called the Large-Angle Magnetic Suspension Test Fixture (LAMSTF), constructed at NASA Langley Research Center in order to explore and develop technology required for the magnetic suspension of objects over large ranges of orientation. The LAMSTF hardware comprises five electromagnets in a circular arrangement, each driven from a separate bipolar power amplifier. The suspended element is a cylindrical axially magnetized permanent magnet core, within an aluminum tube. The element, which is 'levitated' by repulsive forces, is stabilized in five degrees-of-freedom, with rotation about the cylinder axis not controlled. The controller accommodates the changes in magnetic coupling between the electromagnets and the suspended element by real-time adaptation of a decoupling matrix. The paper presents performance measurements demonstrating that the major design objective of the 360 deg rotation was accomplished.

Britcher, Colin P.↗

NASA NTRS2011

Cluster Observations of Multiple Dipolarization Fronts

We present Cluster observations of a series of dipolarization fronts (DF 1 to 6) at the central current sheet in Earth's magnetotail. The velocities of fast earthward flow following behind each DF 1-3, are comparable to the Alfven velocity, indicating that the flow bursts might have been generated by bursty reconnection that occurred tailward of the spacecraft. Based on multi-spacecraft timing analysis, DF normals are found to propagate mainly earthward at $160-335$ km/s with a thickness of 900-1500 km, which corresponds to the ion inertial length or gyroradius scale. Each DF is followed by significant fluctuations in the $x$ and $y$ components of the magnetic field whose peaks are found 1-2 minutes after the DF passage. These $(B_{x},B_{y} )$-fluctuations propagate dawnward (mainly) and earthward. Strongly enhanced field-aligned beams are observed coincidently with $(B_{x},B_{y})$ fluctuations, while an enhancement of cross-tail currents is associated with the DFs. From the observed pressure imbalance and flux-tube entropy changes between the two regions separated by the DF, we speculate that interchange instability destabilizes the DFs and causes the deformation of the mid-tail magnetic topology. This process generates significant field-aligned currents, and might power the auroral brightening in the ionosphere. However, this event is neither associated with the main substorm auroral breakup nor the poleward expansion, which might indicate that the observed multiple DFs have been dissipated before they reach the inner plasma sheet boundary.

Hwang, Kyoung-Joo↗

DOE OSTI2025

ν-point energy correletors with F AST EEC: Small-x physics from LHC jets

In recent years, energy correlators have emerged as a powerful tool for studying jet substructure, with promising applications such as probing the hadronization transition, analyzing the quark-gluon plasma, and improving the precision of top quark mass measurements. The projected N-point correlator measures correlations between N final-state particles by tracking the largest separation between them, showing a scaling behavior related to DGLAP splitting functions. These correlators can be analytically continued in N, commonly referred to as ν-correlators, allowing access to non-integer moments of the splitting functions. Of particular interest is the ν → 0 limit, where the small momentum fraction behavior of the splitting functions requires resummation. Originally, the computational complexity of evaluating ν-correlators for M particles scaled as 2 2M , making it impractical for real-world analyses. However, by using recursion, we reduce this to M 2M , and through the FastEEC method of dynamically resolving subjets, M is replaced by the number of subjets. This breakthrough enables, for the first time, the computation of ν-correlators for LHC data. In practice, limiting the number of subjets to 16 is sufficient to achieve percent-level precision, which we validate using known integer-ν results and convergence tests for non-integer ν. We have implemented this in an update to FastEEC and conducted an initial study of power-law scaling in the perturbative regime as a function of ν, using CMS Open Data on jets. The results agree with DGLAP evolution, except at small ν, where the anomalous dimension saturates to a value that matches the BFKL anomalous dimension.

Energy correlators↗

NASA NTRS2002

Prediction of Aerodynamic Coefficient using Genetic Algorithm Optimized Neural Network for Sparse Data

Wind tunnels use scale models to characterize aerodynamic coefficients, Wind tunnel testing can be slow and costly due to high personnel overhead and intensive power utilization. Although manual curve fitting can be done, it is highly efficient to use a neural network to define the complex relationship between variables. Numerical simulation of complex vehicles on the wide range of conditions required for flight simulation requires static and dynamic data. Static data at low Mach numbers and angles of attack may be obtained with simpler Euler codes. Static data of stalled vehicles where zones of flow separation are usually present at higher angles of attack require Navier-Stokes simulations which are costly due to the large processing time required to attain convergence. Preliminary dynamic data may be obtained with simpler methods based on correlations and vortex methods; however, accurate prediction of the dynamic coefficients requires complex and costly numerical simulations. A reliable and fast method of predicting complex aerodynamic coefficients for flight simulation I'S presented using a neural network. The training data for the neural network are derived from numerical simulations and wind-tunnel experiments. The aerodynamic coefficients are modeled as functions of the flow characteristics and the control surfaces of the vehicle. The basic coefficients of lift, drag and pitching moment are expressed as functions of angles of attack and Mach number. The modeled and training aerodynamic coefficients show good agreement. This method shows excellent potential for rapid development of aerodynamic models for flight simulation. Genetic Algorithms (GA) are used to optimize a previously built Artificial Neural Network (ANN) that reliably predicts aerodynamic coefficients. Results indicate that the GA provided an efficient method of optimizing the ANN model to predict aerodynamic coefficients. The reliability of the ANN using the GA includes prediction of aerodynamic coefficients to an accuracy of 110% . In our problem, we would like to get an optimized neural network architecture and minimum data set. This has been accomplished within 500 training cycles of a neural network. After removing training pairs (outliers), the GA has produced much better results. The neural network constructed is a feed forward neural network with a back propagation learning mechanism. The main goal has been to free the network design process from constraints of human biases, and to discover better forms of neural network architectures. The automation of the network architecture search by genetic algorithms seems to have been the best way to achieve this goal.

Rajkumar, T.↗

NASA NTRS2011

Cluster Observations of Multiple Dipolarization Fronts

We present Cluster observations of a series of dipolarization fronts (DF 1 to 6) at the central current sheet in Earth's magnetotail. The velocities of fast earthward flow following behind each DF 1.3 are comparable to the Alfven velocity, indicating that the flow bursts might have been generated by bursty reconnection that occurred tailward of the spacecraft. Based on multispacecraft timing analysis, DF normals are found to propagate mainly earthward at 160.335 km/s with a thickness of 900-1500 km, which corresponds to the ion inertial length or gyroradius scale. Each DF is followed by significant fluctuations in the x and y components of the magnetic field whose peaks are found 1.2 min after the DF passage. These (B(sub x), B(sub y)) fluctuations propagate dawnward (mainly) and earthward. Strongly enhanced field-aligned beams are observed coincidently with (B(sub x), B(sub y)) fluctuations, while an enhancement of cross-tail currents is associated with the DFs. From the observed pressure imbalance and flux tube entropy changes between the two regions separated by the DF, we speculate that interchange instability destabilizes the DFs and causes the deformation of the midtail magnetic topology. This process generates significant field-aligned currents and might power the auroral brightening in the ionosphere. However, this event is associated with neither the main substorm auroral breakup nor the poleward expansion, which might indicate that the observed multiple DFs have been dissipated before they reach the inner plasma sheet boundary.

Hwang, K.-J.↗

NASA NTRS2012

Upper-Ocean Heat Balance Processes and the Walker Circulation in CMIP5 Model Projections

Considerable uncertainty remains as to the importance of mechanisms governing decadal and longer variability of the Walker Circulation, its connection to the tropical climate system, and prospects for tropical climate change in the face of anthropogenic forcing. Most contemporary climate models suggest that in response to elevated CO2 and a warmer but more stratified atmosphere, the required upward mass flux in tropical convection will diminish along with the Walker component of the tropical mean circulation as well. Alternatively, there is also evidence to suggest that the shoaling and increased vertical stratification of the thermocline in the eastern Pacific will enable a muted SST increase there-- preserving or even enhancing some of the dynamical forcing for the Walker cell flow. Over the past decade there have been observational indications of an acceleration in near-surface easterlies, a strengthened Pacific zonal SST gradient, and globally-teleconnected dislocations in precipitation. But is this evidence in support of an ocean dynamical thermostat process posited to accompany anthropogenic forcing, or just residual decadal fluctuations associated with variations in warm and cold ENSO events and other stochastic forcing? From a modeling perspective we try to make headway on this question by examining zonal variations in surface energy fluxes and dynamics governing tropical upper ocean heat content evolution in the WCRP CMIP5 model projections. There is some diversity among model simulations; for example, the CCSM4 indicates net ocean warming over the IndoPacific region while the CSIRO model concentrates separate warming responses over the central Pacific and Indian Ocean regions. The models, as with observations, demonstrate strong local coupling between variations in column water vapor, downward surface longwave radiation and SST; but the spatial patterns of changes in the sign of this relationship differ among models and, for models as a whole, with observations. Our analysis focuses initially on probing the inter-model differences in energy fluxes / transports and Walker Circulation response to forcing. We then attempt to identify statistically the El Nino- / La Nina-related ocean heat content variability unique to each model and regress out the associated energy flux, ocean heat transport and Walker response on these shorter time scales for comparison to that of the anthropogenic signals.

Robertson, F. R.↗

NASA NTRS2019

Magnetic Flux Cancelation as the Buildup and Trigger Mechanism for CME-Producing Eruptions in Two Small Active Regions

We follow two small, magnetically isolated CME (Coronal Mass Ejection)-producing solar active regions (ARs) from the time of their emergence until several days later, when their core regions erupt to produce the CMEs. In both cases, magnetograms show: (a) following an initial period where the poles of the emerging regions separate from each other, the poles then reverse direction and start to retract inward; (b) during the retraction period, flux cancelation occurs along the main neutral line of the regions, (c) this cancelation builds the sheared core field/flux rope that eventually erupts to make the CME. In the two cases, respectively 30 percent and 50 percent of the maximum flux of the region cancels prior to the eruption. Recent studies indicate that solar coronal jets frequently result from small-scale filaments eruptions (Sterling et al. 2015), with those “minifilament” eruptions also being built up and triggered by cancelation of magnetic flux (Panesar et al. 2016). Together, the small-AR eruptions here and the coronal jet results suggest that isolated bipolar regions tend to erupt when some threshold fraction, perhaps in the range of 50 percent, of the region’s maximum flux has canceled. Our observed erupting filaments/flux ropes form at sites of flux cancelation, in agreement with previous observations. Thus, the recent finding that minifilaments that erupt to form jets also form via flux cancelation is further evidence that minifilaments are small-scale versions of the long-studied full-sized filaments. (Details are in Sterling et al. 2018, ApJ, 864, 68.) Supported by NASA's Heliophysics Guest Investigators (HGI) Program and the MSFC (Marshall Space Flight Center)/Hinode project.

Sterling, Alphonse C.↗

NASA NTRS2020

Streaming Instability with Multiple Dust Species – I. Favourable Conditions for the Linear Growth

A recent study suggests that the streaming instability, one of the leading mechanisms for driving the formation of planetesimals, may not be as efficient as previously thought. Under some disc conditions, the growth time-scale of the instability can be longer than the disc lifetime when multiple dust species are considered. To further explore this finding, we use both linear analysis and direct numerical simulations with gas fluid and dust particles tomutually validate and study the unstablemodes of the instability in more detail.We extend the previously studied parameter space by one order of magnitude in both the range of the dust-size distribution [Ts, min, Ts, max] and the total solid-to-gas mass ratio ε and introduce a third dimension with the slope q of the size distribution. We find that the fast-growth regime and the slow-growth regime are distinctly separated in the ε–Ts, max space, while this boundary is not appreciably sensitive to q or Ts, min.With a wide range of dust sizes present in the disc (e.g. Ts, min 10−3), the growth rate in the slow-growth regime decreases as more dust species are considered.With a narrow range of dust sizes (e.g. Ts, max/Ts, min = 5), on the other hand, the growth rate in most of the ε–Ts, max space is converged with increasing dust species, but the fast and the slow growth regimes remain clearly separated. Moreover, it is not necessary that the largest dust species dominate the growth of the unstable modes, and the smaller dust species can affect the growth rate in a complicated way. In any case, we find that the fast-growth regime is bounded by ε 1 or Ts, max 1, which may represent the favourable conditions for planetesimal formation.

Zhaohuan Zhu↗

DOE OSTI2024

A thermal-driven graupel generation process to explain dry-season convective vigor over the Amazon

Large-eddy simulations (LESs) are conducted for each day of the intensive observation periods (IOPs) of the Green Ocean Amazon (GoAmazon) field campaign to characterize the updrafts and microphysics within deep convective cores while contrasting those properties between Amazon wet and dry seasons. Mean Doppler velocity (V dop ) values simulated using LESs are compared with 2-year measurements from a radar wind profiler (RWP) as viewed by statistical composites separated according to wet- and dry-season conditions. In the observed RWP and simulated LES V dop composites, we find more intense low-level updraft velocity, vigorous graupel generation, and intense surface rain during the dry periods compared with the wet periods. To investigate coupled updraft–microphysical processes further, single-day golden cases are selected from the wet and dry periods to conduct detailed cumulus thermal tracking analysis. Tracking analysis reveals that simulated dry-season environments generate more droplet-loaded low-level thermals than wet-season environments. This tendency correlates with seasonal contrasts in buoyancy and vertical moisture advection profiles in large-scale forcing. Employing a normalized time series of mean thermal microphysics, the simulated cumulus thermals appear to be the primary generator of cloud droplets. When subsequent thermals penetrate the ice crystal layer, droplets within the thermals interact with entrained ice crystals, which enhances riming in the thermals. This appears to be a production pathway of graupel/hail particles within simulated deep convective cores. In addition, less-diluted dry-case thermals tend to be elevated higher, and graupel grows further during sedimentation after spilling out from thermals. Therefore, greater concentrations of low-level moist thermals likely result in more graupel/hail production and associated dry-season convective vigor.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1976

Characteristics of wake vortex generated by a Boeing 727 jet transport during two-segment and normal ILS approach flight paths

A series of flight tests was conducted to evaluate the vortex wake characteristics of a Boeing 727 (B727-200) aircraft during conventional and two-segment ILS approaches. Twelve flights of the B727, which was equipped with smoke generators for vortex marking, were flown and its vortex wake was intentionally encountered by a Lear Jet model 23 (LR-23) and a Piper Twin Comanche (PA-30). Location of the B727 vortex during landing approach was measured using a system of photo-theodolites. The tests showed that at a given separation distance there were no readily apparent differences in the upsets resulting from deliberate vortex encounters during the two types of approaches. Timed mappings of the position of the landing configuration vortices showed that they tended to descend approximately 91 m(300 ft) below the flight path of the B727. The flaps of the B727 have a dominant effect on the character of the trailed wake vortex. The clean wing produces a strong, concentrated vortex but as the flaps are lowered, the vortex system becomes more diffuse. Pilot opinion and roll acceleration data indicate that 4.5 n.mi. would be a minimum separation distance at which roll control of light aircraft (less than 5,670 kg (12,500 lb) could be maintained during parallel encounters of the B727's landing configuration wake. This minimum separation distance is generally in scale with results determined from previous tests of other aircraft using the small roll control criteria.

Kurkowski, R. L.↗

NASA NTRS1977

Vorticity imbalance and stability in relation to convection

A complete synoptic-scale vorticity budget was related to convection storm development in the eastern two-thirds of the United States. The 3-h sounding interval permitted a study of time changes of the vorticity budget in areas of convective storms. Results of analyses revealed significant changes in values of terms in the vorticity equation at different stages of squall line development. Average budgets for all areas of convection indicate systematic imbalance in the terms in the vorticity equation. This imbalance resulted primarily from sub-grid scale processes. Potential instability in the lower troposphere was analyzed in relation to the development of convective activity. Instability was related to areas of convection; however, instability alone was inadequate for forecast purposes. Combinations of stability and terms in the vorticity equation in the form of indices succeeded in depicting areas of convection better than any one item separately.

Read, W. L.↗

DOE OSTI2026

Block encoding of the three-dimensional heterogeneous Poisson equation with application to fracture flow

Quantum linear system (QLS) algorithms offer the potential to solve large-scale linear systems exponentially faster than classical methods. However, applying QLS algorithms to real-world problems remains challenging due to issues such as state preparation, data loading, and efficient information extraction. In this work, we study the feasibility of applying QLS algorithms to solve discretized three-dimensional (3D) heterogeneous Poisson equations, with specific examples relating to groundwater flow through geologic fracture networks. We explicitly construct a block encoding for the 3D heterogeneous Poisson matrix by leveraging the sparse local structure of the discretized operator. While classical solvers benefit from preconditioning, we show that block encoding the system matrix and preconditioner separately does not improve the effective condition number that dominates the QLS run-time. This differs from classical approaches where the preconditioner and the system matrix can often be implemented independently. Nevertheless, due to the structure of the problem in three dimensions, the quantum algorithm achieves a run-time of 𝑂⁡(𝑁 2/3 polylog 𝑁 ⋅log (1/𝜖)), outperforming the best classical methods (with run times of 𝑂⁡(𝑁⁢log 𝑁 ⋅log (1/𝜖))) and offering exponential memory savings. These results highlight both the promise and limitations of QLS algorithms for practical scientific computing, and point to effective condition-number reduction as a key barrier in achieving quantum advantages.

58 GEOSCIENCES↗