Search NASA⌕ Search

SEARCH · Search NASA

Results for “splitting methods”

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

Parallel-in-Time Methods for Method-of-Lines Discretizations of Nonlinear Hyperbolic PDEs and Systems (Final Report)

The work for the subcontract is situated in the area of parallel-in-time integration for hyperbolic partial differential equations (PDEs). Parallel-in-time integration is an active area of research due to its ability to enable faster numerical simulations for applications throughout many areas of science. The work in this subcontract builds on a variety of results that were obtained, as part of the work performed for Subcontract No. B648355, for the Multigrid Reduction-in-Time (MGRIT) method from [1] applied to hyperbolic PDEs. This subcontract extends these results further to more efficient methods and to the case of method-of-lines discretizations for nonlinear hyperbolic PDES and systems of PDEs. The following is a summary of the research performed and results achieved during milestone periods 1, 2 and 3 by the PI (Hans De Sterck) and Postdoctoral Research Associate (Oliver Krzysik), for required tasks 1-4 (as listed in the Statement of Work): Research over the previous year has been split into three main projects: (i) solution of acoustic equation system; (ii) solution of nonlinear scalar hyperbolic PDEs; (iii) solution of nonlinear hyperbolic systems of PDEs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Domain Knowledge Guided Bayesian Optimization For Autonomous Alignment Of Complex Scientific Instruments

Bayesian Optimization (BO) is a powerful tool for optimizing complex non-linear systems. However, its performance degrades in high-dimensional problems with tightly coupled parameters and highly asymmetric objective landscapes, where rewards are sparse. In such needle-in-a-haystack scenarios, even advanced methods like trust-region BO (TurBO) often lead to unsatisfactory results. We propose a domain knowledge guided Bayesian Optimization approach, which leverages physical insight to fundamentally simplify the search problem by transforming coordinates to decouple input features and align the active subspaces with the primary search axes. We demonstrate this approach's efficacy on a challenging 12-dimensional, 6-crystal Split-and-Delay optical system, where conventional approaches, including standard BO, TuRBO and multi-objective BO, consistently led to unsatisfactory results. When combined with an reverse annealing exploration strategy, this approach reliably converges to the global optimum. The coordinate transformation itself is the key to this success, significantly accelerating the search by aligning input co-ordinate axes with the problem's active subspaces. As increasingly complex scientific instruments, from large telescopes to new spectrometers at X-ray Free Electron Lasers are deployed, the demand for robust high-dimensional optimization grows. Our results demonstrate a generalizable paradigm: leveraging physical insight to transform high-dimensional, coupled optimization problems into simpler representations can enable rapid and robust automated tuning for consistent high performance while still retaining current optimization algorithms.

FOS: Computer and information sciences↗
NASA NTRS1992

Crew decision making under stress

Flight crews must make decisions and take action when systems fail or emergencies arise during flight. These situations may involve high stress. Full-missiion flight simulation studies have shown that crews differ in how effectively they cope in these circumstances, judged by operational errors and crew coordination. The present study analyzed the problem solving and decision making strategies used by crews led by captains fitting three different personality profiles. Our goal was to identify more and less effective strategies that could serve as the basis for crew selection or training. Methods: Twelve 3-member B-727 crews flew a 5-leg mission simulated flight over 1 1/2 days. Two legs included 4 abnormal events that required decisions during high workload periods. Transcripts of videotapes were analyzed to describe decision making strategies. Crew performance (errors and coordination) was judged on-line and from videotapes by check airmen. Results: Based on a median split of crew performance errors, analyses to date indicate a difference in general strategy between crews who make more or less errors. Higher performance crews showed greater situational awareness - they responded quickly to cues and interpreted them appropriately. They requested more decision relevant information and took into account more constraints. Lower performing crews showed poorer situational awareness, planning, constraint sensitivity, and coordination. The major difference between higher and lower performing crews was that poorer crews made quick decisions and then collected information to confirm their decision. Conclusion: Differences in overall crew performance were associated with differences in situational awareness, information management, and decision strategy. Captain personality profiles were associated with these differences, a finding with implications for crew selection and training.

Orasanu, J.↗
NASA NTRS1995

Sensitivity analysis for aeroacoustic and aeroelastic design of turbomachinery blades

A new method for computing the effect that small changes in the airfoil shape and cascade geometry have on the aeroacoustic and aeroelastic behavior of turbomachinery cascades is presented. The nonlinear unsteady flow is assumed to be composed of a nonlinear steady flow plus a small perturbation unsteady flow that is harmonic in time. First, the full potential equation is used to describe the behavior of the nonlinear mean (steady) flow through a two-dimensional cascade. The small disturbance unsteady flow through the cascade is described by the linearized Euler equations. Using rapid distortion theory, the unsteady velocity is split into a rotational part that contains the vorticity and an irrotational part described by a scalar potential. The unsteady vorticity transport is described analytically in terms of the drift and stream functions computed from the steady flow. Hence, the solution of the linearized Euler equations may be reduced to a single inhomogeneous equation for the unsteady potential. The steady flow and small disturbance unsteady flow equations are discretized using bilinear quadrilateral isoparametric finite elements. The nonlinear mean flow solution and streamline computational grid are computed simultaneously using Newton iteration. At each step of the Newton iteration, LU decomposition is used to solve the resulting set of linear equations. The unsteady flow problem is linear, and is also solved using LU decomposition. Next, a sensitivity analysis is performed to determine the effect small changes in cascade and airfoil geometry have on the mean and unsteady flow fields. The sensitivity analysis makes use of the nominal steady and unsteady flow LU decompositions so that no additional matrices need to be factored. Hence, the present method is computationally very efficient. To demonstrate how the sensitivity analysis may be used to redesign cascades, a compressor is redesigned for improved aeroelastic stability and two different fan exit guide vanes are redesigned for reduced downstream radiated noise. In addition, a framework detailing how the two-dimensional version of the method may be used to redesign three-dimensional geometries is presented.

Lorence, Christopher B.↗
DOE OSTI2017

Local nematic susceptibility in stressed BaFe 2 As 2 from NMR electric field gradient measurements

The electric field gradient (EFG) tensor at the 75 As site couples to the orbital occupations of the As p orbitals and is a sensitive probe of local nematicity in BaFe 2 As 2 . We use nuclear magnetic resonance to measure the nuclear quadrupolar splittings and find that the EFG asymmetry responds linearly to the presence of a strain field in the paramagnetic phase. We extract the nematic susceptibility from the slope of this linear response as a function of temperature and find that it diverges near the structural transition, in agreement with other measures of the bulk nematic susceptibility. In conclusion, our work establishes an alternative method to extract the nematic susceptibility which, in contrast to transport methods, can be extended inside the superconducting state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1994

Dynamic analysis of a magnetic bearing system with flux control

Using measured values of two-dimensional forces in a magnetic actuator, equations of motion for an active magnetic bearing are presented. The presence of geometric coupling between coordinate directions causes the equations of motion to be nonlinear. Two methods are used to examine the unbalance response of the system: simulation by direct integration in time; and determination of approximate steady state solutions by harmonic balance. For relatively large values of the derivative control coefficient, the system behaves in an essentially linear manner, but for lower values of this parameter, or for higher values of the coupling coefficient, the response shows a split of amplitudes in the two principal directions. This bifurcation is sensitive to initial conditions. The harmonic balance solution shows that the separation of amplitudes actually corresponds to a change in stability of multiple coexisting solutions.

Knight, Josiah↗
DOE OSTI2023

A Stabilizer Framework for the Contextual Subspace Variational Quantum Eigensolver and the Noncontextual Projection Ansatz

Quantum chemistry is a promising application for noisy intermediate-scale quantum (NISQ) devices. However, quantum computers have thus far not succeeded in providing solutions to problems of real scientific significance, with algorithmic advances being necessary to fully utilize even the modest NISQ machines available today. We discuss a method of ground state energy estimation predicated on a partitioning of the molecular Hamiltonian into two parts: one that is noncontextual and can be solved classically, supplemented by a contextual component that yields quantum corrections obtained via a Variational Quantum Eigensolver (VQE) routine. This approach has been termed Contextual Subspace VQE (CS-VQE); however, there are obstacles to overcome before it can be deployed on NISQ devices. The problem we address here is that of the ansatz, a parametrized quantum state over which we optimize during VQE; it is not initially clear how a splitting of the Hamiltonian should be reflected in the CS-VQE ansätze. We propose a “noncontextual projection” approach that is illuminated by a reformulation of CS-VQE in the stabilizer formalism. This defines an ansatz restriction from the full electronic structure problem to the contextual subspace and facilitates an implementation of CS-VQE that may be deployed on NISQ devices. We validate the noncontextual projection ansatz using a quantum simulator and demonstrate chemically precise ground state energy calculations for a suite of small molecules at a significant reduction in the required qubit count and circuit depth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1981

Proposed dynamic phase difference method for the detection of tile debonding from the space shuttle orbiter

A noncontracting, semi-global, dynamic technique was developed for detecting loose tiles on the space shuttle orbiter. In laboratory tests on a single tile, the substrate was excited into lateral motion at a constant frequency and amplitude of 2g. The phase relationship between the motions of tile and substrate was examined by noncontacting probes in order to relate the dynamic properties of the tile SIP system to its fatigue history; by a visual technique using a stroboscope and split screen video monitor for practical application in the field. When the substrate is excited at an appropriate frequency (between 30 and 60 Hz), a good tile moves in phase and a loose tile out of phase with the substrate. The out of phase motion is readily observable in the form of a "beat" between the tile and a reference marker on the substrate.

Zuckerwar, A. J.↗
DOE OSTI2024

Cover crops and poultry litter impact on soil structural stability in dryland soybean production in southeastern United States

Abstract This study explored the efficacy of soil aggregate indices in quantifying soil structural development, utilizing 5‐year field experiment data from the Southeastern United States. The experiment utilized a split‐plot design with cover crops (native vegetation as control, cereal rye (Secale cerealeL.), winter wheat (Triticum aestivum), hairy vetch (Vicia villosa), and mustard (Brassica rapa) plus cereal rye as the main factor and fertilizer source (no fertilizer as control, inorganic fertilizer with phosphorus, potassium, and elemental sulfur, and poultry litter) as the secondary factor. Aggregate size fractions were determined using the wet‐sieving method, and aggregate stability index (ASI), mean weight diameter (MWD), geometric mean diameter (GMD), and fractal dimension (FD) were calculated to assess soil structural stability. Main effects results indicated that cereal rye (55.11%) and poultry litter (50.97%) exhibited the highest ASI values. The highest MWD, GMD, and FD were observed under mustard plus cereal rye (1.187 mm), cereal rye (0.462 mm), and hairy vetch (2.573), respectively. Principal component analysis revealed that cover crops significantly improved soil aggregate structure and stability, overcoming limitations of sole fertilization practices. Regression analysis suggested that ASI, MWD, and GWD positively correlated with soil organic carbon, whereas FD negatively correlated with MWD, GMD, and ASI. Principal component analysis exhibited that FD decreased with increasing soil organic carbon, ASI, MWD, and GMD, demonstrating that lower FD values indicate enhanced soil aggregation and structure. Assessed indices, FD included, effectively gauged soil structural stability. These metrics should be prioritized in managerial decisions to support soil productivity and health in agricultural systems.

Agriculture↗
DOE OSTI2020

Formation of monomeric Sn(ii) and Sn(iv) perfluoropinacolate complexes and their characterization by 119 Sn Mössbauer and 119 Sn NMR spectroscopies

The synthesis and characterization of a series of Sn(II) and Sn(IV) complexes supported by the highly electron-withdrawing dianionic perfluoropinacolate (pinF) ligand are reported herein. Three analogs of [Sn IV (pin F ) 3 ] 2- with NEt 3 H + (1), K + (2), and {K(18C6)} + (3) counter cations and two analogs of [SnII(pinF)2]2- with K+ (4) and {K(15C5) 2 } + (5) counter cations were prepared and characterized by standard analytical methods, single-crystal X-ray diffraction, and 119 Sn Mössbauer and NMR spectroscopies. The six-coordinate SnIV(pinF) complexes display 119 Sn NMR resonances and 119 Sn Mössbauer spectra similar to SnO 2 (cassiterite). In contrast, the four-coordinate SnII(pinF) complexes, featuring a stereochemically-active lone pair, possess low 119 Sn NMR chemical shifts and relatively high quadrupolar splitting. Furthermore, the Sn(II) complexes are unreactive towards both Lewis bases (pyridine, NEt 3 ) and acids (BX 3 , Et 3 NH + ). Calculations confirm that the Sn(II) lone pair is localized within the 5s orbital and reveal that the Sn 5p x LUMO is energetically inaccessible, which effectively abates reactivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Narrow band laser source for efficient 3 He polarization

Polarized 3 He is widely in polarized neutron scattering as neutron polarizers or analyzers. Increasing 3 He polarization can improve both neutron polarization and transmission which are important for any polarized neutron experiments. Here at the US, 3 Hee is usually polarized using spin-exchange optical pumping (SEOP), where 3 He gains polarization via spin-exchange collisions with laser pumped alkali metals, typically rubidium (Rb) and potassium (K). Over the last two decades, 3 He polarization has improved from around 50% to around 80% due to the application of spectral narrowed diode lasers which narrows the laser linewidth from 2 nm to between 0.2 nm and 0.3 nm. However, the highest 3 He polarization achieved so far is limited to about 85%. The mechanism behind this is still not well understood, but it is suspected that a laser with an ultra-narrow linewidth and lower broadband background may further improve the 3 He polarization. In addition, if we define the spin-exchange efficiency as the ratio of the rate at which the polarized 3 He is produced to the rate at which the alkali atom polarization is lost through various collisions, K has much higher spin-exchange efficiency than Rb does. In principle, direct pumping of K should improve both the pumping efficiency and the production rate of polarized 3 He. However, because of the small fine structure splitting in K (D1 770 nm, D2 766.7 nm), the laser used to pump K needs to have a ultra-narrow spectrum tuned at 770 nm and minimal broadband background as the pumping of the K D2 transition will reduce the pumping efficiency. Such a laser is not readily available. Therefore, the current practice of polarizing 3 He is through a hybrid SEOP method in which the laser tuned at Rb D1 transition (795 nm) polarizes Rb first and then Rb transfers its polarization to K. This hybrid method takes advantages of the more accessible 795 nm laser and gains some pumping efficiency. But to realize the full potential in SEOP, direct pumping of K would be necessary. The new Optigrate laser provides a tool to investigate the impact on the pumping efficiency and achievable 3 He polarization under the condition of ultra-narrow spectrum and very low background.

3He polarizing filters↗
DOE OSTI2022

IM3/HyperFACETS Thermodynamic Global Warming (TGW) Simulation Datasets

Publication For a thorough description of the methods, see the peer-reviewed paper: Jones, A.D., Rastogi, D., Vahmani, P. et al. Continental United States climate projections based on thermodynamic modification of historical weather. Sci Data 10, 664 (2023). https://doi.org/10.1038/s41597-023-02485-5 Overview The IM3 / HyperFACETS climate simulations provide 40-year historical (1980-2019) as well as four 80-year future simulations (2020-2099) over the U.S. The future simulations are split into near (2020-2059) and far future (2060-2099) segments. The future scenarios span a range of plausible changes in future climate (both Global Circulation Model (GCM) and Representative Concentration Pathways/Shared Socioeconomic Pathway (RCP/SSP) dimensions). The simulations provide climate variables with high spatiotemporal resolution (25 hourly variables and 207 3-hourly variables at 12 km2). The datasets are generated using dynamical downscaling with the WRF (Weather Research and Forecasting) model (version 4.2.1) and therefore preserve physical consistency across variables. WRF is a state-of-the-art, fully compressible, non-hydrostatic, mesoscale numerical weather prediction model. WRF is coupled with an urban canopy model (UCM), which resolves urban surfaces. The future scenarios were developed using a thermodynamic global warming approach where past events are replayed under a range of future warming conditions. These scenarios therefore provide a perspective on potential increases in extreme event intensity, geographic scope, and duration, with previously non-extreme conditions potentially crossing new thresholds to be considered extreme by today's standards. This approach is not intended to estimate future changes in extreme event frequency that might result from changes in large-scale atmospheric dynamics. This dataset has NOT been bias corrected. A bias corrected version of selected variables is under development and will be released here when available. Scenarios Files Data for each scenario is provided in weekly NetCDF files. 25 variables are available at hourly resolution, and 207 variables are available at three-hourly resolution. Spatial resolution is 12km and spans the conterminous United States (CONUS), including some areas of Canada and Mexico, resulting in a grid of 424 by 299 cells. The spatial projection is a Lambert Conformal Conic with the following proj-string: "+proj=lcc +lat_0=40.0000076293945 +lon_0=-97 +lat_1=30 +lat_2=45 +x_0=0 +y_0=0 +R=6370000 +units=m +no_defs". The available scenarios and simulation periods are listed below: historical | 1980 - 2019 rcp45cooler | 2020 - 2059 rcp45cooler | 2060 - 2099 rcp45hotter | 2020 - 2059 rcp45hotter | 2060 - 2099 rcp85cooler | 2020 - 2059 rcp85cooler | 2060 - 2099 rcp85hotter | 2020 - 2059 rcp85hotter | 2060 - 2099 * The first year (1979, 2019, and 2059) of data within each scenario represents a model warmup period and should not be used. These are located in the `spinup_files` directory. Historical year 2020 is considered an extra year of data beyond the simulation period and can be found in the `additional_files` directory. For information on specific variables and a more in-depth discussion of methodology, please refer to the data landing page at https://tgw-data.msdlive.org. Delta Warming Files The global and CONUS warming deltas for each scenario are provided in degrees Celsius annually and monthly. Restart Files Yearly restart files are provided for each scenario which can be used to restart the WRF model at a particular point in time. Spinup Files The first year of data within each simulation period represents a model warmup period and should not be used. The files are provided here for the sake of reproducibility. Additional Files Additional years of data are provided as an extension of the historic simulation.

Jones, Andrew D.↗
NASA NTRS1991

Interlaminar fracture characterization: A current review

Interlaminar fracture characterization has been investigated for several years. Only now is it well enough understood for standardization organizations to attempt to write standard test methods. A review of the current philosophies in characterizing interlaminar fracture is presented. First, the mode 1, double cantilever beam test (DCB) for measuring G(sub 1c) and the end notched flexure (ENF) test for measuring G(sub 2c) are discussed. These tests have undergone the most extensive research throughout the years and are furthest towards standardization. In addition, the mode 2, end loaded split (ELS) specimen is discussed. Mixed mode fracture is also discussed and the recently developed mixed mode bending (MMB) test is detailed. Then, tests for evaluating mode 3 fracture toughness, including the split cantilever beam (SCB), are reviewed. Finally, the work done on interlaminar fracture characterization in fatigue is reviewed.

Martin, Roderick H.↗
NASA NTRS2006

High-frequency ECG

The standard ECG is by convention limited to 0.05-150 Hz, but higher frequencies are also present in the ECG signal. With high-resolution technology, it is possible to record and analyze these higher frequencies. The highest amplitudes of the high-frequency components are found within the QRS complex. In past years, the term "high frequency", "high fidelity", and "wideband electrocardiography" have been used by several investigators to refer to the process of recording ECGs with an extended bandwidth of up to 1000 Hz. Several investigators have tried to analyze HF-QRS with the hope that additional features seen in the QRS complex would provide information enhancing the diagnostic value of the ECG. The development of computerized ECG-recording devices that made it possible to record ECG signals with high resolution in both time and amplitude, as well as better possibilities to store and process the signals digitally, offered new methods for analysis. Different techniques to extract the HF-QRS have been described. Several bandwidths and filter types have been applied for the extraction as well as different signal-averaging techniques for noise reduction. There is no standard method for acquiring and quantifying HF-QRS. The physiological mechanisms underlying HF-QRS are still not fully understood. One theory is that HF-QRS are related to the conduction velocity and the fragmentation of the depolarization wave in the myocardium. In a three-dimensional model of the ventricles with a fractal conduction system it was shown that high numbers of splitting branches are associated with HF-QRS. In this experiment, it was also shown that the changes seen in HF-QRS in patients with myocardial ischemia might be due to the slowing of the conduction velocity in the region of ischemia. This mechanism has been tested by Watanabe et al by infusing sodium channel blockers into the left anterior descending artery in dogs. In their study, 60 unipolar ECGs were recorded from the entire ventricular surface and were signal-averaged and filtered in the 30-250 Hz frequency range. The results showed that the decrease noted in the HF-QRS correlated linearly with the local conduction delay. The results suggest that HF-QRS is a potent indicator of disturbed local conduction. An alternative theory is that HF-QRS reflect the shape of the original electrocardiographic signal. Bennhagen et al showed that root mean square (RMS) voltage values of the depolarization signal correlate poorly with the signal amplitude but highly with the first and second derivatives, i.e. the velocity and the acceleration of the signal. It has also been suggested that the autonomic nervous system affects HF-QRS. For example, sitting up causes significant changes in HF-QRS in some leads compared to the supine position [Douglas et al., 2006]. Unpublished results indicate that familial dysautonomic patients (both vagal and sympathetic degeneration) have very little Reduced Amplitude Zones (RAZ) formation . Athletic individuals, especially elite athletes, who have vagally-mediated changes on the conventional ECG (i.e. early repolarization, bradycardia) have increased RAZ formation. Further electrophysiological studies are needed, however, to better understand the underlying mechanisms of HF-QRS. Several investigators have studied HF-QRS in different cardiac conditions, including acute myocardial ischemia and myocardial infarction (MI). However, in order for clinicians to confidently use HF-QRS as an adjunct to standard ECG, more knowledge about the characteristics of HF-QRS is needed.

Tragardh, Elin↗
DOE OSTI2024

Metal-facilitated, sustainable nitroarene hydrogenation under ambient conditions

Hydrogenation is a critical reaction in the chemical industry, yielding a range of important compounds such as fine chemicals, pharmachemicals and agrochemicals. However, conventional hydrogenation typically requires pressurized hydrogen, high temperatures and involves noble metal catalysts. Here, we proposed a two-step hydrogenation process, utilizing water as the hydrogen source for the industrially important reduction of nitroarenes to anilines. A metal or reduced metal oxide, which can be obtained from solar thermal or electrochemical reduction, acts as the active site for nitrobenzene adsorption, H 2 O dissociation and in-situ hydrogen generation. Among the 15 metal and reduced metal oxides investigated, Zn and Sn emerged as highly efficient catalysts for the reduction of a broad range of organic nitro compounds under mild conditions, with H 2 utilization efficiency 1-2 orders of magnitude above the state-of-the-art. The presented protocol provides extra dimensions for designing and optimizing conventional hydrogenation process with an alternative pathway. The reactive hydrogen atoms generated in-situ effectively overcome the barriers associated with hydrogen gas dissolution and its subsequent dissociation on the catalyst surface, thereby greatly enhancing the overall effectiveness for the hydrogenation reaction. This research potentially establishes a sustainable, generally applicable alternative to conventional hydrogenation methods, simultaneously presenting a viable solution for renewable energy storage.

25 ENERGY STORAGE↗
DOE OSTI2017

Atomic and electronic structure of the CdTe(111)B–(2√3 × 4) orthogonal surface

The atomic and electronic structure of four variants of Te-terminated CdTe(111)B–(2√3 × 4) orthogonal polar surface (ideal, relaxed, reconstructed, and reconstructed with subsequent relaxation) are calculated ab initio for the first time. The surface is modeled by a film composed of 12 atomic layers with a vacuum gap of ~16 Å in the layered superlattice approximation. To close Cd dangling bonds on the opposite side of the film, 24 fictitious hydrogen atoms with a charge of 1.5 electrons each are added. Ab initio calculations are performed using the Quantum Espresso program based on density functional theory. It is demonstrated that relaxation leads to splitting of the four upper layers. The band energy structures and total and layer-by-layer densities of electronic states for the four surface variants are calculated and analyzed.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A new relaxed splitting preconditioner for the generalized saddle point problems from the incompressible Navier–Stokes equations

Based on the modified relaxed splitting (MRS) preconditioner proposed by Fan and Zhu (Appl Math Lett 55:18–26, 2016), a new relaxed splitting preconditioner is proposed for the generalized saddle point problems arising from the incompressible Navier–Stokes equations. The eigenvalue distribution and an upper bound of the degree of the minimal polynomial of the preconditioned matrix are analyzed. Numerical results show that the proposed preconditioner is superior to the MRS preconditioner when they are used as preconditioners to accelerate the convergence rate of the Krylov subspace methods, such as GMRES.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Electrochemical pH regulation in droplet microfluidics

Here we report a method for electrochemical pH regulation in microdroplets generated in a microfluidic device. The key finding is that controlled quantities of reagents can be generated electrochemically in moving microdroplets confined within a microfluidic channel. Additionally, products generated at the anode and cathode can be isolated within descendant microdroplets. Specifically, ~5 nL water-in-oil microdroplets are produced at a T-junction and then later split into two descendant droplets. During splitting, floor-patterned microelectrodes drive water electrolysis within the aqueous microdroplets to produce H + and OH - . This results in a change in the pHs of the descendant droplets. The droplet pH can be regulated over a range of 5.9 to 7.7 by injecting controlled amounts of charge into the droplets. When the injected charge is between -6.3 and 54.5 nC nL -1 , the measured pH of the resulting droplets is within ±0.1 pH units of that predicted based on the magnitude of the injected charge. This technique can likely be adapted to electrogeneration of other reagents within microdroplets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗