Search NASASearch

SEARCH · Search NASA

Results for “SOLUTION”

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 73 records · Page 4

Gravitational Effects on the Morphology and Kinetics of Photodeposition of Polydiacetylene Thin Films From Monomer Solutions

The goal of this proposed work is to study gravitational effects on the photodeposition of polydiacetylene thin films from monomer solutions onto transparent substrates. Polydiacetylenes have been an extensively studied class of organic polymers because they exhibit many unusual and interesting properties, including electrical conductivity and optical nonlinearity. Their long polymeric chains render polydiacetylenes readily conducive to thin film formation, which is necessary for many applications. These applications require thin polydiacetylene films possessing uniform thicknesses, high purity, minimal inhomogeneities and defects (such as scattering centers), etc. Also, understanding and controlling the microstructure and morphology of the films is important for optimizing their electronic and optical properties. The lack of techniques for processing polydiacetylenes into such films has been the primary limitation to their commercial use. We have recently discovered a novel method for the formation of polydiacetylene thin films using photo-deposition from monomer solutions onto transparent substrates with UV light. This technique is very simple to carry out, and can yield films with superior quality to those produced by conventional methods. Furthermore, these films exhibit good third-order properties and are capable of waveguiding. We have been actively studying the chemistry of diacetylene polymerization in solution and the photo-deposition of polydiacetylene thin films from solution. It is well-known that gravitational factors such as buoyancy-driven convection and sedimentation can affect chemical and mass transport processes in solution. One important aspect of polydiacetylene thin film photodeposition in solution, relevant to microgravity science, is that heat generated by absorption of UV radiation induces thermal density gradients that under the influence of gravity, can cause fluid flows (buoyancy-driven convection). Additionally, changes in the chemical composition of the solution during polymerization may cause solutal convection. These fluid flows affect transport of material to and from the film surface and thereby affect the kinetics of the growth process. This manifests itself in the morphology of the resulting films; films grown under the influence of convection tend to have less uniform thicknesses, and can possess greater inhomogeneities and defects. Specifically, polydiacetylene films photodeposited from solution, when viewed under a microscope, exhibit very small particles of solid polymer which get transported by convection from the bulk solution to the surface of the growing film and become embedded. Even when carried out under conditions designed to minimize unstable density gradients (i.e., irradiating the solution from the top), some fluid flow still takes place (particles remain present in the films). It is also possible that defect nucleation may be occurring within the films or on the surface of the substrate; this, too, can be affected by convection (as is the case with crystal growth). Hence films grown in 1-g will, at best, still possess some defects. The objective of this proposal is to investigate, both in 1-g and in low-g, the effects of gravitational factors (primarily convection) on the dynamics of these processes, and on the quality, morphology, and properties of the films obtained.

Paley, Mark S.

A coupled Eulerian/Lagrangian method for the solution of three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of three-dimensional rotational flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method using particle markers is added to the Eulerian time-marching procedure and provides a correction of the Eulerian solution. In turn, the Eulerian solutions is used to integrate the Lagrangian state-vector along the particles trajectories. The Lagrangian correction technique does not require any a-priori information on the structure or position of the vortical regions. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers, used as 'accuracy boosters,' take advantage of the accurate convection description of the Lagrangian solution and enhance the vorticity and entropy capturing capabilities of standard Eulerian finite-volume methods. The combined solution procedures is tested in several applications. The convection of a Lamb vortex in a straight channel is used as an unsteady compressible flow preservation test case. The other test cases concern steady incompressible flow calculations and include the preservation of turbulent inlet velocity profile, the swirling flow in a pipe, and the constant stagnation pressure flow and secondary flow calculations in bends. The last application deals with the external flow past a wing with emphasis on the trailing vortex solution. The improvement due to the addition of the Lagrangian correction technique is measured by comparison with analytical solutions when available or with Eulerian solutions on finer grids. The use of the combined Eulerian/Lagrangian scheme results in substantially lower grid resolution requirements than the standard Eulerian scheme for a given solution accuracy.

Felici, Helene Marie

Application of matched asymptotic expansions to lunar and interplanetary trajectories. Volume 2: Derivations of second-order asymptotic boundary value solutions

Previously published asymptotic solutions for lunar and interplanetery trajectories have been modified and combined to formulate a general analytical solution to the problem of N-bodies. The earlier first-order solutions, derived by the method of matched asymptotic expansions, have been extended to second order for the purpose of obtaining increased accuracy. The complete derivation of the second-order solution, including the application of a regorous matching principle, is given. It is shown that the outer and inner expansions can be matched in a region of order mu to the alpha power, where 2/5 alpha 1/2, and mu (the moon/earth or planet/sun mass ratio) is much less than one. The second-order asymptotic solution has been used as a basis for formulating a number of analytical two-point boundary value solutions. These include earth-to-moon, one- and two-impulse moon-to-Earth, and interplanetary solutions. Each is presented as an explicit analytical solution which does not require iterative steps to satisfy the boundary conditions. The complete derivation of each solution is shown, as well as instructions for numerical evaluation. For Vol. 1, see N73-27738.

Lancaster, J. E.

A similarity solution of the Navier-Stokes equations with wall catalysis and slip for hypersonic, low Reynolds number flow over spheres

The slip conditions for a multicomponent mixture with diffusion, wall-catalyzed atom recombination and thermal radiation are derived. The more realistic multicomponent species slip conditions are shown to be necessary for accurate merged shock layer solutions on a sphere. These slip conditions are used in a first-order similarity solution of the Navier-Stokes equations with nonequilibrium chemistry for the merged shock layer. Results of this quick numerical solution are compared with a time dependent solution around the sphere and with measured arc jet results at low Reynolds numbers. The similarity solution, unlike the time dependent solution, shows smooth radial profiles of the pressure and smooth variations of velocity slip, skin friction, temperature slip and heat transfer around the body. The present first-order similarity solution is valid up to 25 deg from the stagnation point and takes less than 1% of the computer time to run a time dependent scheme. The smaller stand-off distance obtained from the similarity solution is supported by experimental data. The measured heat flux is closer to the similarity solution than the time dependent method at the stagnation point and shows the proper variation with circumferential angle up to at least 40 deg.

Hendricks, W. L.

Solution of the minimum time-to-climb problem by matched asymptotic expansions

Application of singular perturbation techniques to trajectory optimization problems of flight mechanics is discussed. The method of matched asymptotic expansions is used to obtain an approximate solution to the aircraft minimum time-to-climb problem. Outer, boundary-layer, and composite solutions are obtained to zeroth and first orders. A stability criterion is derived for the zeroth-order boundary-layer solutions (the theory requires a form of boundary-layer stability). A numerical example is considered for which it is shown that the stability criterion is satisfied and a useful numerical solution is obtained. The zeroth-order solution proves to be a poor approximation, but the first-order solution gives a good approximation for both the trajectory and the minimum time-to-climb. The computational cost of the singular perturbation solution is considerably less than that of a steepest descent solution. Thus singular perturbation methods appear to be promising for the solution of optimal control problems.

Ardema, M. D.

Numerical solution of random singular integral equation appearing in crack problems

The solution of several elasticity problems, and particularly crack problems, can be reduced to the solution of one-dimensional singular integral equations with a Cauchy-type kernel or to a system of uncoupled singular integral equations. Here a method for the numerical solution of random singular integral equations of Cauchy type is presented. The solution technique involves a Chebyshev series approximation, the coefficients of which are the solutions of a system of random linear equations. This method is applied to the problem of periodic array of straight cracks inside an infinite isotropic elastic medium and subjected to a nonuniform pressure distribution along the crack edges. The statistical properties of the random solution are evaluated numerically, and the random solution is used to determine the values of the stress-intensity factors at the crack tips. The error, expressed as the difference between the mean of the random solution and the deterministic solution, is established. Values of stress-intensity factors at the crack tip for different random input functions are presented.

Sambandham, M.

Radiation boundary condition and anisotropy correction for finite difference solutions of the Helmholtz equation

In this paper finite-difference solutions of the Helmholtz equation in an open domain are considered. By using a second-order central difference scheme and the Bayliss-Turkel radiation boundary condition, reasonably accurate solutions can be obtained when the number of grid points per acoustic wavelength used is large. However, when a smaller number of grid points per wavelength is used excessive reflections occur which tend to overwhelm the computed solutions. Excessive reflections are due to the incompability between the governing finite difference equation and the Bayliss-Turkel radiation boundary condition. The Bayliss-Turkel radiation boundary condition was developed from the asymptotic solution of the partial differential equation. To obtain compatibility, the radiation boundary condition should be constructed from the asymptotic solution of the finite difference equation instead. Examples are provided using the improved radiation boundary condition based on the asymptotic solution of the governing finite difference equation. The computed results are free of reflections even when only five grid points per wavelength are used. The improved radiation boundary condition has also been tested for problems with complex acoustic sources and sources embedded in a uniform mean flow. The present method of developing a radiation boundary condition is also applicable to higher order finite difference schemes. In all these cases no reflected waves could be detected. The use of finite difference approximation inevita bly introduces anisotropy into the governing field equation. The effect of anisotropy is to distort the directional distribution of the amplitude and phase of the computed solution. It can be quite large when the number of grid points per wavelength used in the computation is small. A way to correct this effect is proposed. The correction factor developed from the asymptotic solutions is source independent and, hence, can be determined once and for all. The effectiveness of the correction factor in providing improvements to the computed solution is demonstrated in this paper.

Tam, Christopher K. W.

Ca-Rich Carbonate Melts: A Regular-Solution Model, with Applications to Carbonatite Magma + Vapor Equilibria and Carbonate Lavas on Venus

A thermochemical model of the activities of species in carbonate-rich melts would be useful in quantifying chemical equilibria between carbonatite magmas and vapors and in extrapolating liquidus equilibria to unexplored PTX. A regular-solution model of Ca-rich carbonate melts is developed here, using the fact that they are ionic liquids, and can be treated (to a first approximation) as interpenetrating regular solutions of cations and of anions. Thermochemical data on systems of alkali metal cations with carbonate and other anions are drawn from the literature; data on systems with alkaline earth (and other) cations and carbonate (and other) anions are derived here from liquidus phase equilibria. The model is validated in that all available data (at 1 kbar) are consistent with single values for the melting temperature and heat of fusion for calcite, and all liquidi are consistent with the liquids acting as regular solutions. At 1 kbar, the metastable congruent melting temperature of calcite (CaCO3) is inferred to be 1596 K, with (Delta)bar-H(sub fus)(calcite) = 31.5 +/- 1 kJ/mol. Regular solution interaction parameters (W) for Ca(2+) and alkali metal cations are in the range -3 to -12 kJ/sq mol; W for Ca(2+)-Ba(2+) is approximately -11 kJ/sq mol; W for Ca(2+)-Mg(2+) is approximately -40 kJ/sq mol, and W for Ca(2+)-La(3+) is approximately +85 kJ/sq mol. Solutions of carbonate and most anions (including OH(-), F(-), and SO4(2-)) are nearly ideal, with W between 0(ideal) and -2.5 kJ/sq mol. The interaction of carbonate and phosphate ions is strongly nonideal, which is consistent with the suggestion of carbonate-phosphate liquid immiscibility. Interaction of carbonate and sulfide ions is also nonideal and suggestive of carbonate-sulfide liquid immiscibility. Solution of H2O, for all but the most H2O-rich compositions, can be modeled as a disproportionation to hydronium (H3O(+)) and hydroxyl (OH(-)) ions with W for Ca(2+)-H3O(+) (approximately) equals 33 kJ/sq mol. The regular-solution model of carbonate melts can be applied to problems of carbonatite magma + vapor equilibria and of extrapolating liquidus equilibria to unstudied systems. Calculations on one carbonatite (the Husereau dike, Oka complex, Quebec, Canada) show that the anion solution of its magma contained an OH mole fraction of (approximately) 0.07, although the vapor in equilibrium with the magma had P(H2O) = 8.5 x P(CO2). F in carbonatite systems is calculated to be strongly partitioned into the magma (as F(-)) relative to coexisting vapor. In the Husereau carbonatite magma, the anion solution contained an F(-) mole fraction of (approximately) 6 x 10(exp -5).

Treiman, Allan H.

Protein Crystal Growth Under Forced Solution Flow: Experimental Setup and General Response of Lysozyme

We have experimentally studied the effects of solution flow on the growth kinetics of the protein lysozyme. To this end, we have expanded our interferometry setup by a novel crystallization cell and solution recirculation system. This combination permits monitoring of interface morphology and kinetics with a depth resolution of 200 A at bulk flow rates of up to 2000 micron/s. Particular attention was paid to the prevention of protein denaturation that is often associated with the pumping of protein solutions. We found that at bulk flow rates it less than 250 microns/s the average growth rate and step velocity, R(sub avg) and upsilon(sub avg) increase with increasing it. This can be quantitatively understood in terms of the enhanced, convective solute supply to the interface. With high-purity solutions, it u greater than 250 microns/s lead to growth deceleration, and, at low supersaturations sigma, to growth cessation. When solutions containing approx. 1% of other protein impurities were used, growth deceleration occurred at any u greater than 0 and cessation in the low sigma experiments was reached at about half the it causing cessation with pure solution. The flow-induced changes in R(sub avg) and upsilon(sub avg) including growth cessation, were reversible and reproducible, independent of the direction of the u-changes and solution purity. Hence, we attribute the deceleration to the convection-enhanced supply of impurities to the interface, which at higher flow rates overpowers the effects of enhanced interfacial solute concentration. Most importantly, we found that convective transport leads to a significant reduction in kinetics fluctuations, in agreement with our earlier expectations for the lysozyme system. This supports our hypothesis that these long-term fluctuations represent an intrinsic response feature of the coupled bulk transport-interfacial kinetics system in the mixed growth control regime.

Vekilov, P. G.

Modelling and Estimation of Dynamic Ocean Topography Within Global Geopotential Solutions

In this poster, we will present a report on the ongoing investigation "Improved Gravitational and Dynamic Height Models Through the use of Oceanographic Data." we have used a two year mean (1993-1994) of the Dynamic ocean Topography (DOT) field implied by the Semtner and Chervin POCM_4B model, and developed normal equations in surface spherical harmonics to degree and order 30. These normal equations, were combined with normal equations derived from TOPEX and ERS-1 altimeter data over the same time period. Combination solutions (based on satellite tracking data, altimeter data, surface gravity data and OCM data) were were developed Test solutions were obtained estimating the DOT field to 2Ox2O and 3ox3O. These solutions were tested with independent DOT values computed over 38 WOCE hydrographic sections, which contained a total of 3072 stations and represented 216000 km of travelled lines. The weighted standard deviation of the differences between the DOT obtained from the hydrographic data and and the field estimated from the joint combination solutions was computed for each of the test models, the weighted standard deviation for the baseline combination solution excluding the POCM4_B data was 9.7 cm for a DOT solution to 30x30, where introducing the POCM4_B data into the combination model reduced the standard deviation to 9.2 cm, indicating the the introduction of oceanographic information benefits the solution for the dynamic ocean topography. We will discuss the weighting schemes applied and the method of solution. Another aspect of our investigation involves testing alternate parametric representations of the dynamic height field. We looked at the alternative representations in terms of the Proudman functions (PF), and compared these to the use of spherical harmonics (SH) to represent the dynamic ocean topography, using once again the 1993 and 1994 output of the POCM4_B OCM as the reference model. A significant advantage of PF's over SH's is that the former require no "fill-in" values over areas where the DOT is undefined. We show that the PF and SH results using equal number of parameters agree quite well in overall content. In addition, PFs appear better suited for representation of high frequency signals close to the signal boundary. The development of the PF solutions for the DOT will be described and the statistics of the comparisons will be presented.

Lemoine, Frank G.

Advanced in Visualization of 3D Time-Dependent CFD Solutions

Numerical simulations of complex 3D time-dependent (unsteady) flows are becoming increasingly feasible because of the progress in computing systems. Unfortunately, many existing flow visualization systems were developed for time-independent (steady) solutions and do not adequately depict solutions from unsteady flow simulations. Furthermore, most systems only handle one time step of the solutions individually and do not consider the time-dependent nature of the solutions. For example, instantaneous streamlines are computed by tracking the particles using one time step of the solution. However, for streaklines and timelines, particles need to be tracked through all time steps. Streaklines can reveal quite different information about the flow than those revealed by instantaneous streamlines. Comparisons of instantaneous streamlines with dynamic streaklines are shown. For a complex 3D flow simulation, it is common to generate a grid system with several millions of grid points and to have tens of thousands of time steps. The disk requirement for storing the flow data can easily be tens of gigabytes. Visualizing solutions of this magnitude is a challenging problem with today's computer hardware technology. Even interactive visualization of one time step of the flow data can be a problem for some existing flow visualization systems because of the size of the grid. Current approaches for visualizing complex 3D time-dependent CFD solutions are described. The flow visualization system developed at NASA Ames Research Center to compute time-dependent particle traces from unsteady CFD solutions is described. The system computes particle traces (streaklines) by integrating through the time steps. This system has been used by several NASA scientists to visualize their CFD time-dependent solutions. The flow visualization capabilities of this system are described, and visualization results are shown.

Lane, David A.

Lunar Dust and Lunar Simulant Activation, Monitoring, Solution and Cellular Toxicity Properties

During the Apollo missions, many undesirable situations were encountered that must be mitigated prior to returning humans to the moon. Lunar dust (that part of the lunar regolith less than 20 m in diameter) was found to produce several problems with astronaut s suits and helmets, mechanical seals and equipment, and could have conceivably produced harmful physiological effects for the astronauts. For instance, the abrasive nature of the dust was found to cause malfunctions of various joints and seals of the spacecraft and suits. Additionally, though efforts were made to exclude lunar dust from the cabin of the lunar module, a significant amount of material nonetheless found its way inside. With the loss of gravity correlated with ascent of the lunar module from the lunar surface to rendezvous with the command module, much of the major portions of the contaminating soil and dust began to float, irritating the astronaut s eyes and being inhaled into their lungs. Our goal has been to understand some of the properties of lunar dust that could lead to possible hazards for humans. Due to the lack of an atmosphere, there is nothing to protect the lunar soil from ultraviolet radiation, solar wind, and meteorite impacts. These processes could all serve to activate the soil, or produce reactive surface species. In order to understand the possible toxic effects of the reactive dust, it is necessary to reactivate the dust, as samples returned during the Apollo missions were exposed to the atmosphere of the Earth. We have used grinding and UV exposure to mimic some of the processes occurring on the Moon. The level of activation has been monitored using two methods: fluorescence spectroscopy and electron paramagnetic resonance spectroscopy (EPR). These techniques allow the monitoring of hydroxyl radical production in solution. We have found that grinding of lunar dust produces 2-3 times the concentration of hydroxyl radicals as lunar simulant and 10 times that of quartz. Exposure of the lunar dust to UV radiation under vacuum was also found to lead to hydroxyl radical production. After grinding, we have also monitored loss of reactivity of the dusts by exposing them to conditions of known humidity and temperature. From these tests, it was found that the reactivity half-life of lunar simulant is approximately 3 hours, while that of quartz is approximately 2 hours. Placing lunar dust in solution could lead to effects on mechanical and physiological systems, as well as other biological systems. For instance, while it is known that lunar dust is highly abrasive and caused a variety of problems with suits and equipment during Apollo, it is unknown as to how these properties might be affected in the presence of water or other liquids. It is possible that the dust may release minerals (e.g., metallic nanophase Fe) into solution that could speed corrosion or rust. Also, as lunar dust produces hydroxyl radicals (and possibly other reactive oxygen species) in solution, these radicals could also lead to the breakdown of suit or habitat materials. In the body (i.e., in lung solution), the effects could be two-fold. First, if the lunar dust dissolves, it may release an excess of elements (such as zero-valence metallic Fe) that are necessary for bodily functions but only in certain concentration ranges. For lunar dust, the presence of nanophase iron being released into the body is a concern. Secondly, the hydroxyl radicals or other reactive oxygen species produced by the dust in solution could conceivably interact with cells, leading to various problems. We have studied the dissolution of both ground and unground lunar simulant in buffer solutions of different pH. The concentration of a number of species was determined using mass spectrometry. These studies showed that lowering the pH of the solution causes a dramatic increase in the amount of each element released into solution and that grinding also produces higher concentrations. Finally, we have perfmed initial tests aimed at understanding the effects of lunar simulant on cellular systems. Alveolar epithelial cells were cultured and exposed to different concentrations of dust suspended in cell culture media. After predetermined amounts of time, the media was removed and the concentrations of important inflammatory cytokines (IL6, IL8, and TNF-alpha ) were measured. The results of these tests are being used to develop the correct protocols for tests to be performed using lunar dust samples.

Jeevarajan, A.S.

A General Closed-Form Solution for the Lunar Reconnaissance Orbiter (LRO) Antenna Pointing System

The National Aeronautics and Space Administration s (NASA) Lunar Reconnaissance Orbiter (LRO) launched on June 18, 2009 from the Cape Canaveral Air Force Station aboard an Atlas V launch vehicle into a direct insertion trajectory to the Moon LRO, designed, built, and operated by the NASA Goddard Space Flight Center in Greenbelt, MD, is gathering crucial data on the lunar environment that will help astronauts prepare for long-duration lunar expeditions. During the mission s nominal life of one year its six instruments and one technology demonstrator will find safe landing site, locate potential resources, characterize the radiation environment and test new technology. To date, LRO has been operating well within the bounds of its requirements and has been collecting excellent science data images taken from the LRO Camera Narrow Angle Camera (LROC NAC) of the Apollo landing sites have appeared on cable news networks. A significant amount of information on LRO s science instruments is provided at the LRO mission webpage. LRO s Attitude Control System (ACS), in addition to controlling the orientation of the spacecraft is also responsible for pointing the High Gain Antenna (HGA). A dual-axis (or double-gimbaled) antenna, deployed on a meter-long boom, is required to point at a selected Earth ground station. Due to signal loss over the distance from the Moon to Earth, pointing precision for the antenna system is very tight. Since the HGA has to be deployed in spaceflight, its exact geometry relative to the spacecraft body is uncertain. In addition, thermal distortions and mechanical errors/tolerances must be characterized and removed to realize the greatest gain from the antenna system. These reasons necessitate the need for an in-flight calibration. Once in orbit around the moon, a series of attitude maneuvers was conducted to provide data needed to determine optimal parameters to load onboard, which would account for the environmental and mechanical errors at any antenna orientation. The nominal geometry for the HGA involves an outer gimbal axis that is exactly perpendicular to the inner gimbal axis, and a target direction that is exactly perpendicular to the outer gimbal axis. For this nominal geometry, closed-form solutions of the desired gimbal angles are simple to get for a desired target direction specified in the spacecraft body fame. If the gimbal axes and the antenna boresight are slightly misaligned, the nominal closed-form solution is not sufficiently accurate for computing the gimbal angles needed to point at a target. In this situation, either a general closed-form solution has to be developed for a mechanism with general geometries, or a correction scheme has to be applied to the nominal closed-form solutions. The latter has been adopted for Solar Dynamics Observatory (SDO) as can be seen in Reference 1, and the former has been used for LRO. The advantage of the general closed-form solution is the use of a small number of parameters for the correction of nominal solutions, especially in the regions near singularities. Singularities here refer to cases when the nominal closed-form solutions have two or more solutions. Algorithm complexity, however, is the disadvantage of the general closed-form solution.

Shah, Neerav

On the Solution of the Continuity Equation for Precipitating Electrons in Solar Flares

Electrons accelerated in solar flares are injected into the surrounding plasma, where they are subjected to the influence of collisional (Coulomb) energy losses. Their evolution is modeled by a partial differential equation describing continuity of electron number. In a recent paper, Dobranskis & Zharkova claim to have found an "updated exact analytical solution" to this continuity equation. Their solution contains an additional term that drives an exponential decrease in electron density with depth, leading them to assert that the well-known solution derived by Brown, Syrovatskii & Shmeleva, and many others is invalid. We show that the solution of Dobranskis & Zharkova results from a fundamental error in the application of the method of characteristics and is hence incorrect. Further, their comparison of the "new" analytical solution with numerical solutions of the Fokker-Planck equation fails to lend support to their result.We conclude that Dobranskis & Zharkova's solution of the universally accepted and well-established continuity equation is incorrect, and that their criticism of the correct solution is unfounded. We also demonstrate the formal equivalence of the approaches of Syrovatskii & Shmeleva and Brown, with particular reference to the evolution of the electron flux and number density (both differential in energy) in a collisional thick target. We strongly urge use of these long-established, correct solutions in future works.

Sun: flares

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES

Pushing the limits of size selectivity in nanoscale solute separations

Transport of a spherical solute through a cylindrical pore has been modelled for decades using well-established hindered transport theory, predicting solutes with a size smaller than the pore to be rejected nonetheless because of convective and diffusive hindrance; this rejection mechanism prevents extremely sharp solute separations by a membrane. Whereas the model has been historically verified, solute transport through near-perfect isoporous membranes may finally overcome this limitation. Here, in this study, encouraging solute rejections are achieved using nanofabricated, defect-free silicon nitride isoporous membranes. The membrane is challenged by a recirculated feed to increase the opportunity for interactions between solutes and the pore array. Results show the membrane completely reject solutes with greater size than the pore size while effectively allowing smaller solutes to permeate through. With effectively increasing the number of interactions, we propose that a steeper size-selective rejection curve may be achieved. With this traditional hurdle overcome, there is new promise for unprecedented membrane separations through judicious process design and extremely tight pore-size distributions.

Gao, Feng

Isolating solvent–solute hydrogen bonding interactions via 2D IR solvation shell spectroscopy

The solvation shell around a solute is a fundamental feature of liquid-phase solutions, determining the behavior and properties of both the solute and the overall solution. Direct experimental measurements of the solvation shell properties are challenging due to the strong signals generated from the bulk solvent, which overwhelm the small contribution of the solvation shell. Here, we use ultrafast two dimensional infrared (2D IR) spectroscopy and intermolecular cross-peaks to isolate the IR absorption spectrum of methanol molecules in the solvation shell surrounding the solute N-methylacetamide. We demonstrate that the intermolecular coupling between the solvent and solute vibrations is indirectly mediated by a low-frequency hydrogen-bonding mode, suggesting an important mechanism for anharmonic coupling induced by hydrogen bonds. From the relative frequency shifts and cross-peak anisotropy, we find that methanol molecules surrounding N-methylacetamide form stronger and distinctly oriented hydrogen bonds than those in the bulk solvent. Here, we also compare these results with the solvent spectra of the solute N,N-dimethylacetamide to investigate how solute structural changes alter the solvation shell and the contribution of N–H hydrogen bond donation. Our results are supported by molecular dynamics simulations, which provide detailed insights into the hydrogen-bonding distributions. Through these results, we demonstrate 2D solvation shell spectroscopy to be a valuable method for investigating solvation structures and dynamics without interference from the bulk solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Hierarchical Network Partitioning for Solution of Potential-Driven, Steady-State Nonlinear Network Flow Equations

The solution of potential-driven steady-state flow in large networks is a task which manifests in various engineering applications, such as transport of natural gas or water through pipeline networks. The resultant system of nonlinear equations depends on the network topology, and in general, there is no numerical algorithm that offers guaranteed convergence to the solution (assuming a solution exists). Some methods offer guarantees in cases where the network topology satisfies certain assumptions, but these methods fail for larger networks. On the other hand, the Newton-Raphson algorithm offers a convergence guarantee if the starting point lies close to the (unknown) solution. It would be advantageous to compute the solution of the large nonlinear system through the solution of smaller nonlinear sub-systems wherein the solution algorithms (Newton-Raphson or otherwise) are more likely to succeed. Here, this letter proposes and describes such a procedure, a hierarchical network partitioning algorithm that enables the solution of large nonlinear systems corresponding to potential-driven steady-state network flow equations.

42 ENGINEERING