Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational tractability”

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 361 records · Page 20

A Theory of Ultrafast Charge Transfer Relaxation with Non-Innocent Solvent Molecules

In this work we revisit the photodynamics of tetracyanoethylene-hexamethylbenzene (TCNE- HMB), the molecular complex studied by Hochstrasser et al. [J. Chem. Phys. 100, 4797–4810, 1994] that has long challenged the applicability of Marcus’s theory of elec- tron transfer for predicting photochemical reactions. Using a novel black-box electronic structure algorithm (time-dependent density functional theory with one double, TD- DFT-1D) to efficiently run molecular dynamics that can treat charge recombination, we run ab initio surface hopping molecular dynamics and confirm that, for a polar solvent, charge recombination rates can be incredibly fast (indeed faster than the sol- vent relaxation time); for non-polar solvents, the rate is much slower. We demonstrate that, although Marcus theory cannot be directly applied, these nonequilibrium (and sometimes incredibly fast) photoexcited dynamics can be effectively explained within a two-state model without any evidence of a transition through a conical intersection. Most importantly, for this paradigmatic model system, we are able to identify two nuclear coordinates of interest (rather than the single coordinate predicted by Marcus or a full set of internal quantum modes studied by Bixon and Jortner): the solvent relaxation in the first shell (that strongly modulates the energies of the charge trans- fer state and differentiates time scales for relaxation) and a nuclear displacement in the TCNE-HMB complex arising from a handful of vibrations that induces non-Born Oppenheimer motion and eventually facilitates an abrupt electronic transition to the ground state. Altogether, these findings suggest a tractable generalization of Marcus theory for future simulations of photochemistry with non-innocent solvent environ- ments in the spirit of a Hamiltonian suggested by Stuchebrukhov (J. Chem. Phys. 107, 3821, 1997).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic, Adaptive, Systems and Materials: Complex, Simple and Emergent Behaviors

This program has been funded by DoE/BES for twenty years. It has moved into and out of various subjects as it has developed, but it has retained its focus on complexity and complex systems. The project has evolved in the following way: Self-Assembly and Biomimetic Self-Assembly: All self-assembling systems depend upon a minimum of two types of interaction: a repulsion and an attraction. For the familiar molecular systems, attractive interactions are typically hydrogen bonds and electrostatic interactions. Repulsive interactions include steric effects, hydrophobic effects (in biological systems), and charge-charge repulsion. We have expanded this repertoire to include surface interactions, magnetic interactions, and others. I list these systems in the order in which we have explored them: i) A key emphasis in current work is in understanding how the movement of ions in a magnetic field (the Lorentz effect) interacts with catalytic systems. We have demonstrated that an acceleration in rate of reduction of CO 2 to CO can be accomplished by applying an external magnetic field. This acceleration is largely due to the application of the Lorentz effect on mass transport at the catalyst’s surface. ii) We have also extensively explored the influence of electrostatics, as exhibited in self-assembling systems, by tribocharging. iii) Another key system involves surface tension effects; examples include interactions between heavy particles floating at a liquid-air interface, and interacting by changes in surface area; interactions of bubbles and bubble rafts, behaviors of bubble trains in microfluidic networks, and behaviors of microorganisms in constraining environments. iv) This work has intentionally de-emphasized biological systems; but it does include some work on protein-ligand interactions and interactions among microorganisms. v) We have also explored applications of some of these effects, these explorations include bubble rafts as diffraction gratings, exploration of the structures that can be obtained by tribocharging and uses of these structures in exploring nucleation and melting of crystals. vi) Although not a major focus of this work, several other topics have emerged and offer opportunities for future work. These include the behavior of bubble trains and bubble rafts in microfluidic systems. A particularly interesting example is the formation of bubble trains that repeat in the alteration of large and small bubbles according to rules we do not presently understand, but are uniquely large-period oscillating systems. These systems offer a new route into understanding the instabilities of the type represented by oscillations. vii) We have also begun exploratory projects on magnetic levitation (especially to determine molecular density), and information storage (in molecules). Magnetic Levitation: Self-assembling and biomimetic systems require both attraction and repulsion. We have used electrostatics (tribocharging), interfacial free-energies (surface tension and related forces) and others. Potential uses include reconfigurable diffraction gratings and liquid lenses; exploration of mechanisms in tribocharging; tunneling in EGaIn junctions; and bubble trains (especially in micro-fluidic systems). Examples of systems representing these topics is included in the following papers: Complexity: Disks rotating at a water-air interface; Benard-Marangoni effects; Vortex-Crystals from spinning magnetic disks (Marangoni effects); EGaIn Electrode to study quantum tunneling; Self-Assembly of electrostatically-charged metallic spheres (electrets); Dynamically reconfigurable lens; Using computational designs of ligands for enzymes; Electrostatic self-assembly by tribocharging; Monodisperse bubble trains in microchannel systems; Inverted dripping faucet; Flames; Printing of micro-organisms to regenerate the “ink” of printing device; Using micro-organisms to move loads (“microoxen”); Motion of bacterial swarms near surfaces; Making monodisperse particles in microfluidic systems; Coding/decoding of information stored in droplet trains in microfluidic networks; Magnetic levitation; and Information storage. i) Tribocharging. The change in focus of this work on electrets from the fundamentals of charging to applications of these materials in studying self-assembly using electrostatic interactions. ii) Bubbles in Microchannels. The realization that systems of bubbles in microchannels represented a major opportunity to study complexity in a very tractable system, and the development of a semi-quantitative theory of this subject. iii) Flames. The growth of “flames” remains an exploratory subject for the research, although their currently relatively little active work involving it ongoing. iv) Systems with Microorganisms. The removal of work in biological systems from this project. Based on work supported in this program, we now have a significant project on the development of microfluidic tools for studying C. elegans (a nematode), but this work was not appropriate for a program focused on complexity, and we developed separate support for it. (It is, however, an example of successful seeding of a new area by BES.) The work on electrets has gone through a period in which a part of the program was the subject of a MURI; the focus of this work was to develop materials that did not charge electrostatically on friction or contact. The MURI is now over, and the work on dynamic self-assembly (supported by BES) is the major focus. “Flames” has also enjoyed synergistic support, in terms of a project supported by DARPA on flame suppression (in the absence of extinguishing agents, using acoustic and electrostatic interactions). This work was helpful in understanding some of the basics of flames, but is entirely distinct from the BES focus in complexity. A growing interest is in the Lorentz effect. The Lorentz effect is the force exerted on charged particles (electrons, ions, charged molecules) when they move through a perpendicular magnetic field. The Lorentz effect is almost ubiquitous in modern technology: examples of applications include electric motors, dynamos, cathode ray tubes, many batteries, and most systems that control electrical currents with magnetic forces. We have begun to explore the Lorentz effect in electrochemical systems and heterogeneous catalytic systems involving charged organic species and inorganic ions. This work is still at an early stage, but initial studies that Lorentz effects can be large when ions move through magnetic fields, or magnetic fields move in the presence of ions.

36 MATERIALS SCIENCE↗

Multidisciplinary optimization applied to a transport aircraft

Decomposition of a large optimization problem into several smaller subproblems has been proposed as an approach to making large-scale optimization problems tractable. To date, the characteristics of this approach have been tested on problems of limited complexity. The objective of the effort is to demonstrate the application of this multilevel optimization method on a large-scale design study using analytical models comparable to those currently being used in the aircraft industry. The purpose of the design study which is underway to provide this demonstration is to generate a wing design for a transport aircraft which will perform a specified mission with minimum block fuel. A definition of the problem; a discussion of the multilevel composition which is used for an aircraft wing; descriptions of analysis and optimization procedures used at each level; and numerical results obtained to date are included. Computational times required to perform various steps in the process are also given. Finally, a summary of the current status and plans for continuation of this development effort are given.

Giles, G. L.↗

Hypothesis testing via AI: Generating physically interpretable models of scientific data with machine learning (Full Technical Report)

Deep learning has demonstrated an exceptional ability to solve complex tasks (an engineering success); however, it has done so at the expense of the ability to generate new knowledge (a scientific failure). We propose an alternative framework—entitled Deep Symbolic Regression (DSR)—in which artificial neural networks (NNs) rapidly generate hypotheses about physical relationships among inputs. This framework bypasses the need to interpret an NN altogether, while still leveraging the representational power of deep learning. The resulting models are tractable mathematical expressions, which are inherently and readily human interpretable and can provide insights into underlying physical phenomena. Further, we fold this methodology into the scientific process by allowing the scientist to directly integrate a priori knowledge and beliefs to accelerate learning. We demonstrate this methodology on symbolic regression—the problem of rediscovering underlying expressions describing a dataset—and achieve state-of-the-art performance across a wide variety of symbolic regression problems. Further, we generalize our DSR framework to apply to the more general class of symbolic optimization problems, in which one seeks to optimize a sequence of symbols or “tokens” under a black-box reward function. Examples of other symbolic optimization problems include neural architecture search and computational antibody design. Our generalized tool, Deep Symbolic Optimization (DSO), has been demonstrated on the task of learning symbolic control policies for reinforcement learning environments, and has been adopted as an enabling capability for computational antibody design.

97 MATHEMATICS AND COMPUTING↗

Mean field model parameterization to recover coverage-dependent kinetics

Lateral interactions between adsorbates introduce coverage dependence into adsorption energies and activation barriers of surface reactions. Lattice-based kinetic Monte Carlo (kMC) simulations can capture these interactions quantitatively but are laborious to parameterize and solve. Mean field models are more tractable, but protocols to construct and parameterize them are unclear. Here we explore the ability of a coverage-aware mean-field model to map to a lattice-kMC model of a generic two-step reaction network, including quasi-equilibrated adsorption and rate-limiting dissociation steps. We derive expressions for mean-field and coverage-dependent adsorption energies and dissociation barriers and parameterize against lattice-kMC predictions. We show that the parameterized mean-field rates correlate with ground truth lattice-kMC results across a wide range of reaction conditions and identify regions where the mean field fails. The mean field model similarly captures kMC-derived rate-order, Arrhenius and Sabatier plots at a greatly reduced computational cost. Further, the results provide guidance for parameterizing mean-field models, benchmarked against explicit lattice-based approaches for incorporating the influence of coverage effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Bilevel Approach for Identifying the Worst Contingencies for Nonconvex Alternating Current Power Systems

We address the bilevel optimization problem of identifying the most critical attacks to an alternating current (AC) power flow network. The upper-level binary maximization problem consists of choosing an attack that is treated as a parameter in the lower-level defender minimization problem. Instances of the lower-level global minimization problem by themselves are NP-hard due to the nonconvex AC power flow constraints, and bilevel solution approaches commonly apply a convex relaxation or approximation to allow for tractable bilevel reformulations at the cost of underestimating some power system vulnerabilities. Our main contribution is to provide an alternative branch-and-bound algorithm whose upper bounding mechanism (in a maximization context) is based on a reformulation that avoids relaxation of the AC power flow constraints in the lower-level defender problem. Lower bounding is provided with semidefinite programming (SDP) relaxed solutions to the lower-level problem. We establish finite termination with guarantees of either a globally optimal solution to the original bilevel problem, or a globally optimal solution to the SDP-relaxed bilevel problem which is included in a vetted list of upper-level attack solutions, at least one of which is a globally optimal solution to the bilevel problem. We demonstrate through computational experiments applied to IEEE case instances both the relevance of our contribution, and the effectiveness of our contributed algorithm for identifying power system vulnerabilities without resorting to convex relaxations of the lower-level problem. We conclude with a discussion of future extensions and improvements.

97 MATHEMATICS AND COMPUTING↗

Lighting up yeast: overview of optogenetics in yeast and their applications to yeast biotechnology

Abstract Optogenetics is an empowering technology that uses light-responsive proteins to control biological processes. Because of its genetic tractability, abundance of genetic tools, and robust culturing conditions, Saccharomyces cerevisiae has served for many years as an ideal platform in which to study, develop, and apply a wide range of optogenetic systems. In many instances, yeast has been used as a steppingstone in which to characterize and optimize optogenetic tools to later be deployed in higher eukaryotes. More recently, however, optogenetic tools have been developed and deployed in yeast specifically for biotechnological applications, including in nonconventional yeasts. In this review, we summarize various optogenetic systems responding to different wavelengths of light that have been demonstrated in diverse yeast species. We then describe various applications of these optogenetic tools in yeast, particularly in metabolic engineering and recombinant protein production. Finally, we discuss emerging applications in yeast cybergenetics—the interfacing of yeast and computers for closed-loop controls of yeast bioprocesses—and the potential impact of optogenetics in other future biotechnological applications.

Biotechnology & Applied Microbiology↗

Human Factors Directions for Civil Aviation

Despite considerable progress in understanding human capabilities and limitations, incorporating human factors into aircraft design, operation, and certification, and the emergence of new technologies designed to reduce workload and enhance human performance in the system, most aviation accidents still involve human errors. Such errors occur as a direct or indirect result of untimely, inappropriate, or erroneous actions (or inactions) by apparently well-trained and experienced pilots, controllers, and maintainers. The field of human factors has solved many of the more tractable problems related to simple ergonomics, cockpit layout, symbology, and so on. We have learned much about the relationships between people and machines, but know less about how to form successful partnerships between humans and the information technologies that are beginning to play a central role in aviation. Significant changes envisioned in the structure of the airspace, pilots and controllers' roles and responsibilities, and air/ground technologies will require a similarly significant investment in human factors during the next few decades to ensure the effective integration of pilots, controllers, dispatchers, and maintainers into the new system. Many of the topics that will be addressed are not new because progress in crucial areas, such as eliminating human error, has been slow. A multidisciplinary approach that capitalizes upon human studies and new classes of information, computational models, intelligent analytical tools, and close collaborations with organizations that build, operate, and regulate aviation technology will ensure that the field of human factors meets the challenge.

Hart, Sandra G.↗

A Crop Yield Change Emulator for Use in GCAM and Similar Models: Persephone v1.0

Future changes in Earth system state will impact agricultural yields and, through these changed yields, can have profound impacts on the global economy. Global gridded crop models estimate the influence of these Earth system changes on future crop yields but are often too computationally intensive to dynamically couple into global multisector economic models, such as the Global Change Assessment Model (GCAM) and other similar-in-scale models. Yet, generalizing a faster site-specific crop model’s results to be used globally will introduce inaccuracies, and the question of which model to use is unclear given the wide variation in yield response across crop models. To examine the feedback loop among socioeconomics, Earth system changes, and crop yield changes, rapidly generated yield responses with some quantification of crop response uncertainty are desirable. The Persephone v1.0 response functions presented in this work are based on the Agricultural Model Intercomparison and Improvement Project (AgMIP) Coordinated Climate-Crop Modeling Project (C3MP) sensitivity test data set and are focused on providing GCAM and similar models with a tractable number of rapid to evaluate dynamic yield response functions corresponding to a range of the yield response sensitivities seen in the C3MP data set. With the Persephone response functions, a new variety of agricultural impact experiments will be open to GCAM and other economic models: for example, examining the economic impacts of a multi-year drought in a key agricultural region and how economic changes in response to the drought can, in turn, impact the drought.

Snyder, Abigail↗

Formal Requirements-Based Programming for Complex Systems

Computer science as a field has not yet produced a general method to mechanically transform complex computer system requirements into a provably equivalent implementation. Such a method would be one major step towards dealing with complexity in computing, yet it remains the elusive holy grail of system development. Currently available tools and methods that start with a formal model of a system and mechanically produce a provably equivalent implementation are valuable but not sufficient. The gap that such tools and methods leave unfilled is that the formal models cannot be proven to be equivalent to the system requirements as originated by the customer For the classes of complex systems whose behavior can be described as a finite (but significant) set of scenarios, we offer a method for mechanically transforming requirements (expressed in restricted natural language, or appropriate graphical notations) into a provably equivalent formal model that can be used as the basis for code generation and other transformations. While other techniques are available, this method is unique in offering full mathematical tractability while using notations and techniques that are well known and well trusted. We illustrate the application of the method to an example procedure from the Hubble Robotic Servicing Mission currently under study and preliminary formulation at NASA Goddard Space Flight Center.

Rash, James L.↗

MINE: a new way to design genetics experiments for discovery

Abstract The Maximally Informative Next Experiment or MINE is a new experimental design approach for experiments, such as those in omics, in which the number of effects or parameters p greatly exceeds the number of samples n (p > n). Classical experimental design presumes n > p for inference about parameters and its application to p > n can lead to over-fitting. To overcome p > n, MINE is an ensemble method, which makes predictions about future experiments from an existing ensemble of models consistent with available data in order to select the most informative next experiment. Its advantages are in exploration of the data for new relationships with n < p and being able to integrate smaller and more tractable experiments to replace adaptively one large classic experiment as discoveries are made. Thus, using MINE is model-guided and adaptive over time in a large omics study. Here, MINE is illustrated in two distinct multiyear experiments, one involving genetic networks in Neurospora crassa and a second one involving a genome-wide association study in Sorghum bicolor as a comparison to classic experimental design in an agricultural setting.

Biochemistry & Molecular Biology↗

Study and modeling of finite rate chemistry effects in turbulent non-premixed flames

The development of numerical models that reflect some of the most important features of turbulent reacting flows requires information about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between turbulent and chemical processes is so strong that it is extremely difficult to isolate the role played by one individual physical phenomenon. Direct numerical simulation (hereafter DNS) allows us to study in detail the turbulence-chemistry interaction in some restricted but completely defined situations. Globally, non-premixed flames are controlled by two limiting regimes: the fast chemistry case, where the turbulent flame can be pictured as a random distribution of local chemical equilibrium problems; and the slow chemistry case, where the chemistry integrates in time the turbulent fluctuations. The Damkoehler number, ratio of a mechanical time scale to chemical time scale, is used to distinguish between these regimes. Today most of the industrial computer codes are able to perform predictions in the hypothesis of local equilibrium chemistry using a presumed shape for the probability density function (pdt) of the conserved scalar. However, the finite rate chemistry situation is of great interest because industrial burners usually generate regimes in which, at some points, the flame is undergoing local extinction or at least non-equilibrium situations. Moreover, this variety of situations strongly influences the production of pollutants. To quantify finite rate chemistry effect, the interaction between a non-premixed flame and a free decaying turbulence is studied using DNS. The attention is focused on the dynamic of extinction, and an attempt is made to quantify the effect of the reaction on the small scale mixing process. The unequal diffusivity effect is also addressed. Finally, a simple turbulent combustion model based on the DNS observations and tractable in real flow configurations is proposed.

Vervisch, Luc↗

Bayesian Estimation of Precipitation from Satellite Passive Microwave Observations Using Combined Radar-Radiometer Retrievals

Precipitation estimation from satellite passive microwave radiometer observations is a problem that does not have a unique solution that is insensitive to errors in the input data. Traditionally, to make this problem well posed, a priori information derived from physical models or independent, high-quality observations is incorporated into the solution. In the present study, a database of precipitation profiles and associated brightness temperatures is constructed to serve as a priori information in a passive microwave radiometer algorithm. The precipitation profiles are derived from a Tropical Rainfall Measuring Mission (TRMM) combined radar radiometer algorithm, and the brightness temperatures are TRMM Microwave Imager (TMI) observed. Because the observed brightness temperatures are consistent with those derived from a radiative transfer model embedded in the combined algorithm, the precipitation brightness temperature database is considered to be physically consistent. The database examined here is derived from the analysis of a month-long record of TRMM data that yields more than a million profiles of precipitation and associated brightness temperatures. These profiles are clustered into a tractable number of classes based on the local sea surface temperature, a radiometer-based estimate of the echo-top height (the height beyond which the reflectivity drops below 17 dBZ), and brightness temperature principal components. For each class, the mean precipitation profile, brightness temperature principal components, and probability of occurrence are determined. The precipitation brightness temperature database supports a radiometer-only algorithm that incorporates a Bayesian estimation methodology. In the Bayesian framework, precipitation estimates are weighted averages of the mean precipitation values corresponding to the classes in the database, with the weights being determined according to the similarity between the observed brightness temperature principal components and the brightness temperature principal components of the classes. Because the classes are stratified by the sea surface temperature and the echo-top-height estimator, the number of classes that are considered for retrieval is significantly smaller than the total number of classes, making the algorithm computationally efficient. The radiometer-only algorithm is applied to TMI observations, and precipitation estimates are compared with combined TRMM precipitation radar (PR) TMI reference estimates. The TMI-only algorithm, supported by the empirically derived database, produces estimates that are more consistent with the reference values than the precipitation estimates from the version-6 TRMM facility TMI algorithm. Cloud-resolving model simulations are used to assign a latent heating profile to each precipitation profile in the empirically derived database, making it possible to estimate latent heating using the radiometer-only algorithm. Although the evaluation of latent heating estimates in this study is preliminary, because realistic conditional probability distribution functions are attached to latent heating structures in the algorithm s database, a generally positive impact on latent heating estimation from passive microwave observations is expected.

Grecu, Mircea↗

3D Cloud Tomography and Droplet Size Retrieval from Multi-Angle Polarimetric Imaging of Scattered Sunlight from Above

Tomography aims to recover a three-dimensional (3D) density map of a medium or an object. In medical imaging,it is extensively used for diagnostics via X-ray computed tomography (CT). We define and derive a tomographyof cloud droplet distributions via passive remote sensing. We use multi-view polarimetric images to fit a 3Dpolarized radiative transfer (RT) forward model. Our motivation is 3D volumetric probing of vertically-developedconvectively-driven clouds that are ill-served by current methods in operational passive remote sensing. Currenttechniques are indeed based on strictly 1D RT modeling and applied to a single cloudy pixel, where cloud geometrydefaults to that of a plane-parallel slab. Incident unpolarized sunlight, once scattered by cloud droplets, changesits polarization state according to droplet size. Therefore, polarimetric measurements in the rainbow and gloryangular regions can be used to infer the droplet size distribution. This work defines and derives a framework for afull 3D tomography of cloud droplets for both their mass concentration in space and their distribution across arange of sizes. This gridded 3D retrieval of key microphysical properties is made tractable by our novel approachthat involves a restructuring and partial linearization of an open-source polarized 3D RT code to accommodate aspecial two-step iterative optimization technique. Physically-realistic synthetic clouds are used to demonstrate themethodology with rigorous uncertainty quantification, while a real-world cloud imaged by AirMSPI is processedto illustrate the new remote sensing capability

Schechner, Yoav Y.↗

Combining Force Fields and Neural Networks for an Accurate Representation of Chemically Diverse Molecular Interactions

A key goal of molecular modeling is the accurate reproduction of the true quantum mechanical potential energy of arbitrary molecular ensembles with a tractable classical approximation. The challenges are that analytical expressions found in general purpose force fields struggle to faithfully represent the intermolecular quantum potential energy surface at close distances and in strong interaction regimes; that the more accurate neural network approximations do not capture crucial physics concepts, e.g., nonadditive inductive contributions and application of electric fields; and that the ultra-accurate narrowly targeted models have difficulty generalizing to the entire chemical space. We therefore designed a hybrid wide-coverage intermolecular interaction model consisting of an analytically polarizable force field combined with a short-range neural network correction for the total intermolecular interaction energy. Here, we describe the methodology and apply the model to accurately determine the properties of water, the free energy of solvation of neutral and charged molecules, and the binding free energy of ligands to proteins. The correction is subtyped for distinct chemical species to match the underlying force field, to segment and reduce the amount of quantum training data, and to increase accuracy and computational speed. For the systems considered, the hybrid ab initio parametrized Hamiltonian reproduces the two-body dimer quantum mechanics (QM) energies to within 0.03 kcal/mol and the nonadditive many-molecule contributions to within 2%. Simulations of molecular systems using this interaction model run at speeds of several nanoseconds per day.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gyrokinetic simulation of turbulence and transport in the SPARC tokamak

The turbulence and transport expected in the SPARC tokamak Primary Reference Discharge (PRD) [P. Rodriguez-Fernandez et al., J. Plasma Phys. 86, 865860503 (2020)] have been investigated with the gyrokinetic code CGYRO [J. Candy et al., J. Comput. Phys. 324, 73–93 (2016)]. Linear and nonlinear simulations that focus on ion (kθρs<1.0) and electron-scale (kθρs>1.0) turbulence were used to probe the nature of the turbulence and the resulting transport in the fusion core. It is found that in the SPARC PRD, ion temperature gradient (ITG) turbulence is expected to dominate transport over most of the profile with some potential trapped electron mode impact in the near edge. Stiff turbulence is observed over a part of the plasma core such that SPARC's ion temperature profile will likely be pinned to just above the critical gradient for ITG. The role of electromagnetic turbulence, rotation, and electron-scale turbulence was investigated to provide some insight into the physics required to accurately predict SPARC performance via gyrokinetics. Additionally, predictions of impurity peaking for potential low- and high-Z SPARC first-wall materials are probed using ion-scale simulation. The dominance of low-k turbulence in SPARC provides a potential opportunity for more tractable prediction of plasma profiles using nonlinear gyrokinetics. This work is the first step toward full gyrokinetic profile prediction of SPARC kinetic profiles and the resulting fusion power and plasma gain.

Howard, N. T. (ORCID:0000000287876309)↗

Development of a New Membrane Casting Apparatus for Studying Macrovoid Defects in Low-G

A new membrane-casting apparatus is developed for studying macrovoid defects in polymeric membranes made by the wet- and dry-casting process in low-gravity. Macrovoids are large (10-50 micron), open cavities interspersed among the smaller pores in the substructure under the gelled skin surface layer of the cast membrane. Although their occurrence is considered endemic to the wet- and dry-casting process since they can lead to compaction or skin rupture in the membrane process, recent studies suggest several useful applications such as transdermal and osmotic drug delivery systems, miniature bioreactors, etc. However, lack of knowledge about the macrovoid formation mechanism is an obstacle to further development of applications using them. An on-going debate is the role of the surface-tension-driven solutocapillary convection during macrovoid formation. The rapid growth of macrovoids within 1-5 seconds and the high polymer concentration in and near macrovoids make it difficult to explain the mechanism of macrovoid growth by diffusion alone, which is the widely accepted hypothesis proposed by Reuvers et al. The hypothesis advanced by our research group can explain this rapid growth via a mechanism that involves diffusion from the casting solution in the meta-stable region to the macrovoid enhanced by solutocapillary convection induced by the steep nonsolvent concentration gradient in the vicinity of the macrovoid. Since macrovoid growth is hypothesized to be the interplay of a solutocapillary-induced driving force counteracted by viscous drag and buoyancy, eliminate the latter provides a means for testing this hypothesis. Moreover, free convection mass transfer in the nonsolvent immersion bath used to cause phase-separation in membrane casting complicates developing a model for both the wet-casting process and macrovoid growth. The low-g environment minimizes gravitationally induced free convection thereby permitting a tractable solution to the ternary diffusion equations that characterize membrane formation. NASA's Parabolic Flight Research Aircraft provides a small window of low-g (approximately 25 s) that can be used to study macrovoid development in both wet- and dry-cast membranes if an appropriate casting apparatus is used. This casting apparatus should be able to cast the membrane in both low- and high-g in a manner so that essential one-dimensional mass transfer conditions are achieved to insure lateral uniformity in the membrane. The apparatus used in previous research on membrane casting in low-gravity was operated with the plunger driven mechanism. The spring-loaded plunger pushes the bottom block containing the polymer casting solution well directly under the absorbent chamber located in the upper stationary block. However, membranes made via this casting apparatus often displayed lateral nonuniformities that precluded obtaining quantitative information on the macrovoid growth process. Thus, it was necessary to determine the reason for these structural irregularities observed in the low-g casting apparatus. Both experimental as well as computer simulation studies of the low-g casting apparatus established that the impulsive action of the plunger caused the undesired structural nonuniformities. The simulation results showed that the width-to-depth aspect ratio of the shallow well that contains the casting solution in this apparatus was not an important factor in minimizing this problem. Even for a 40:1 (width : depth) aspect ratio, any convection induced by the horizontal motion of the interface of the casting solution will be damped out within 6.25x10(exp 4) seconds. However, the experimental studies revealed that the impulsive motion of the plunger caused a 'sloshing' of the casting solution that had to be eliminated. Therefore, the plungerdriven mechanism was changed to a cam-driven mechanism that did not cause any impulsive motion of the casting solution. Other refinements to this new membrane-casting apparatus include provision for removing the membranes from the casting wells in a less destructive manner. This was accomplished by using a slit geometry for the casting well that permitted disassembly for removal of the cast membrane. The materials used in the construction of this casting apparatus were chosen to insure wetting at the side walls and to maintain precise control of the thickness of the polymer solution in the casting well. An additional provision in this new casting apparatus is the ability to carry out both wet- as well as dry-casting. As such, this apparatus permitted the first studies of the wet-casting of polymeric membranes in low-g. Both wet- and dry-casting experiments on NASA's KC-135 research aircraft employing this new membrane-casting apparatus are scheduled in July 2002. The morphology of the resulting membranes will be characterized using an environmental scanning electron microscope (ESEM). The results of these low-g studies will be reported later.

Lee, Hanyong↗