Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “system partition”

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 433 records · Page 24

Calculating Free Energies Using Scaled-Force Molecular Dynamics Algorithm

One common objective of molecular simulations in chemistry and biology is to calculate the free energy difference between different states of the system of interest. Examples of problems that have such an objective are calculations of receptor-ligand or protein-drug interactions, associations of molecules in response to hydrophobic, and electrostatic interactions or partition of molecules between immiscible liquids. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), 'native' state. Perhaps the best example of such a problem is folding of proteins or short RNA molecules. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to 'quasi non-ergodicity', whereby a part of phase space is inaccessible on timescales of the simulation. A host of strategies has been developed to improve efficiency of sampling the phase space. For example, some Monte Carlo techniques involve large steps which move the system between low-energy regions in phase space without the need for sampling the configurations corresponding to energy barriers (J-walking). Most strategies, however, rely on modifying probabilities of sampling low and high-energy regions in phase space such that transitions between states of interest are encouraged. Perhaps the simplest implementation of this strategy is to increase the temperature of the system. This approach was successfully used to identify denaturation pathways in several proteins, but it is clearly not applicable to protein folding. It is also not a successful method for determining free energy differences. Finally, the approach is likely to fail for systems with co-existing phases, such as water-membrane systems, because it may lead to spontaneous mixing. A similar difficulty may be encountered in any method relying on global modifications of phase space.

Darve, Eric↗

Local-Field Effects in Linear Response Properties within a Polarizable Frozen Density Embedding Method

In this work, we present a polarizable frozen density embedding (FDE) method for calculating polarizabilities of coupled subsystems. The method (FDE-pol) combines a FDE method with an explicit polarization model such that the expensive freeze/thaw cycles can be bypassed, and approximate nonadditive kinetic potentials are avoided by enforcing external orthogonality between the subsystems. To describe the polarization of the frozen environment, we introduce a Hirshfeld partition-based density-dependent method for calculating the atomic polarizabilities of atoms in molecules, which alleviates the need to fit the atomic parameters to a specific system of interest or to a larger general set of molecules. Further, we show that the Hirshfeld partition-based method predicts molecular polarizabilities close to the basis set limit, and thus, a single basis set-dependent scaling parameter can be introduced to improve the agreement against the reference polarizability data. To test the model, we characterized the uncoupled and coupled response of small interacting molecular complexes. Here, the coupled response properties include the perturbation of the frozen system due to the external perturbation which is ignored in the uncoupled response. We show that FDE-pol can accurately reproduce both the exact uncoupled polarizability and the coupled polarizabilities of the supermolecular systems. Using damped response theory, we also demonstrate that the coupled frequency-dependent polarizability can be described by including local field effects. The results emphasize the necessity of including local-field effects for describing the response properties of coupled subsystems, as well as the importance of accurate atomic polarizability models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Distributed state-space generation of discrete-state stochastic models

High-level formalisms such as stochastic Petri nets can be used to model complex systems. Analysis of logical and numerical properties of these models of ten requires the generation and storage of the entire underlying state space. This imposes practical limitations on the types of systems which can be modeled. Because of the vast amount of memory consumed, we investigate distributed algorithms for the generation of state space graphs. The distributed construction allows us to take advantage of the combined memory readily available on a network of workstations. The key technical problem is to find effective methods for on-the-fly partitioning, so that the state space is evenly distributed among processors. In this paper we report on the implementation of a distributed state-space generator that may be linked to a number of existing system modeling tools. We discuss partitioning strategies in the context of Petri net models, and report on performance observed on a network of workstations, as well as on a distributed memory multi-computer.

Ciardo, Gianfranco↗

Spatial heterogeneity effects on land surface modeling of water and energy partitioning

Abstract. Understanding the influence of land surface heterogeneity on surface water and energy fluxes is crucial for modeling earth system variability and change. This study investigates the effects of four dominant heterogeneity sources on land surface modeling, including atmospheric forcing (ATM), soil properties (SOIL), land use and land cover (LULC), and topography (TOPO). Our analysis focused on their impacts on the partitioning of precipitation (P) into evapotranspiration (ET) and runoff (R), partitioning of net radiation into sensible heat and latent heat, and corresponding water and energy fluxes. An initial set of 16 experiments were performed over the continental US (CONUS) using the E3SM land model (ELMv1) with different combinations of heterogeneous and homogeneous datasets. The Sobol' total and first-order sensitivity indices were utilized to quantify the relative importance of the four heterogeneity sources. Sobol' total sensitivity index measures the total heterogeneity effects induced by a given heterogeneity source, consisting of the contribution from its own heterogeneity (i.e., the first-order index) and its interactions with other heterogeneity sources. ATM and LULC are the most dominant heterogeneity sources in determining spatial variability of water and energy partitioning, mainly contributed by their own heterogeneity and slightly contributed by their interactions with other heterogeneity sources. Their heterogeneity effects are complementary, both spatially and temporally. The overall impacts of SOIL and TOPO are negligible, except TOPO dominates the spatial variability of R/P across the transitional climate zone between the arid western and humid eastern CONUS. Accounting for more heterogeneity sources improves the simulated spatial variability of water and energy fluxes when compared with ERA5-Land reanalysis dataset. An additional set of 13 experiments identified the most critical components within each heterogeneity source, which are precipitation, temperature, and longwave radiation for ATM, soil texture, and soil color for SOIL and maximum fractional saturated area parameter for TOPO.

54 ENVIRONMENTAL SCIENCES↗

Three-dimensional, conservative, Euler computations using patched grid systems and explicit methods

The method of 'zonal approach' (in which the flow field is partitioned into regions with independent grids) for computation of flow over complex geometries, such as aircraft configuration, requires application of a grid-interfacing procedure. A three-dimensional conservative, boundary scheme for patched grids, applicable in generalized coordinates for arbitrary point distributions on a planar zonal surface is presented. The computation technique is derived within the framework of the Osher (1983) upwind scheme, using Euler equations. The three-dimensional interfacing method is applied to the computation of flow about a wing-canard combination using a two-zone, patched grid.

Hessenius, K. A.↗

Global Influence of Organic Aerosol Volatility on Aerosol Microphysical Processes: Composition and Number

We present MATRIX-VBS, a new aerosol scheme that simulates organic partitioning in an aerosol microphysics model, as part of the NASA GISS ModelE Earth System Model. MATRIX-VBS builds on its predecessor aerosol microphysics model MATRIX (Bauer et al., 2008, https://doi.org/10.5194/acp-8-6003-2008) and was developed in the box model framework (Gao et al., 2017, https://doi.org/10.5194/gmd-10-751-2017). The scheme features the inclusion of organic partitioning between the gas and particle phases and the photochemical aging process using the volatility-basis set (Donahue et al., 2006, https://doi.org/10.1021/es052297c). To assess the performance of the new model, we compared its mass concentration, number concentration, and activated number concentration to the original scheme MATRIX, and evaluated its mass concentrations in four seasons against data from the NASA Atmospheric Tomography Mission (ATom) aircraft campaign. Results from MATRIX-VBS show that organics are transported further away from their source, and their mass concentration increases aloft and decreases at the surface compared to those in MATRIX. The mass concentration of organics at the surface agrees well with measurements, and there are discrepancies for vertical profiles aloft. In the new scheme, the global mass load of organic aerosols increased by 50%, there is also an increased number of particles at the surface and fewer activated ones in most regions. The new scheme presents advanced and more comprehensive capability in simulating aerosol processes.

aerosol modeling↗

Tunable-diode laser absorption spectrometry

Tunable-diode laser absorption spectrometry (TDLAS) affords a number of advantages for atmospheric measurements. It is a universal method, applicable, in principle, to all gases of atmospheric interest. Because of its extremely high spectral resolution it provides unequivocal identification of the target species, with no interferences from other gases. It provides real-time, in situ measurements with time resolutions better than 1 minute. The sensitivity of the current TDLAS system is marginally capable of measuring HO2. This species exists in the troposphere at concentrations which are up to 2 orders of magnitude higher than those of HO and, in addition, is much less susceptible to removal by the surfaces of the instrument and its sampling system. HO2 is an important HO sub x species in its own right but can also give direct information on the HO concentration by virtue of the rapid partitioning between these two species. The addition of the high-frequency modulation technique to the TDLAS system would ensure its ability to measure HO2 under most atmospheric conditions. The ability of the TDLAS to measure hydrogen peroxide (H2O2) in the ambient atmosphere was also demonstrated. H2O2 measurements give a clear indication of HO sub x mixing ratios and are also important as a photolytic source of HO and as an important oxidant for other atmospheric consitituents such as SO2.

Schiff, H. I.↗

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↗

Integrated flight/propulsion control - Subsystem specifications for performance

A procedure is presented for calculating multiple subsystem specifications given a number of performance requirements on the integrated system. This procedure applies to problems where the control design must be performed in a partitioned manner. It is based on a structured singular value analysis, and generates specifications as magnitude bounds on subsystem uncertainties. The performance requirements should be provided in the form of bounds on transfer functions of the integrated system. This form allows the expression of model following, command tracking, and disturbance rejection requirements. The procedure is demonstrated on a STOVL aircraft design.

Neighbors, W. K.↗

Dynamics of Quantum Adiabatic Evolution Algorithm for Number Partitioning

We have developed a general technique to study the dynamics of the quantum adiabatic evolution algorithm applied to random combinatorial optimization problems in the asymptotic limit of large problem size n. We use as an example the NP-complete Number Partitioning problem and map the algorithm dynamics to that of an auxiliary quantum spin glass system with the slowly varying Hamiltonian. We use a Green function method to obtain the adiabatic eigenstates and the minimum exitation gap, gmin = O(n2(sup -n/2)), corresponding to the exponential complexity of the algorithm for Number Partitioning. The key element of the analysis is the conditional energy distribution computed for the set of all spin configurations generated from a given (ancestor) configuration by simultaneous flipping of a fixed number of spins. For the problem in question this distribution is shown to depend on the ancestor spin configuration only via a certain parameter related to the energy of the configuration. As the result, the algorithm dynamics can be described in terms of one-dimensional quantum diffusion in the energy space. This effect provides a general limitation of a quantum adiabatic computation in random optimization problems. Analytical results are in agreement with the numerical simulation of the algorithm.

Smelyanskiy, Vadius↗

A model of asynchronous iterative algorithms for solving large, sparse, linear systems

Solving large, sparse, linear systems of equations is one of the fundamental problems in large scale scientific and engineering computation. A model of a general class of asynchronous, iterative solution methods for linear systems is developed. In the model, the system is solved by creating several cooperating tasks that each compute a portion of the solution vector. This model is then analyzed to determine the expected intertask data transfer and task computational complexity as functions of the number of tasks. Based on the analysis, recommendations for task partitioning are made. These recommendations are a function of the sparseness of the linear system, its structure (i.e., randomly sparse or banded), and dimension.

Reed, D. A.↗

Progress in Parallel Schur Complement Preconditioning for Computational Fluid Dynamics

We consider preconditioning methods for nonself-adjoint advective-diffusive systems based on a non-overlapping Schur complement procedure for arbitrary triangulated domains. The ultimate goal of this research is to develop scalable preconditioning algorithms for fluid flow discretizations on parallel computing architectures. In our implementation of the Schur complement preconditioning technique, the triangulation is first partitioned into a number of subdomains using the METIS multi-level k-way partitioning code. This partitioning induces a natural 2X2 partitioning of the p.d.e. discretization matrix. By considering various inverse approximations of the 2X2 system, we have developed a family of robust preconditioning techniques. A computer code based on these ideas has been developed and tested on the IBM SP2 and the SGI Power Challenge array using MPI message passing protocol. A number of example CFD calculations will be presented to illustrate and assess various Schur complement approximations.

Barth, Timothy J.↗

Dynamics of Quantum Adiabatic Evolution Algorithm for Number Partitioning

We have developed a general technique to study the dynamics of the quantum adiabatic evolution algorithm applied to random combinatorial optimization problems in the asymptotic limit of large problem size n. We use as an example the NP-complete Number Partitioning problem and map the algorithm dynamics to that of an auxiliary quantum spin glass system with the slowly varying Hamiltonian. We use a Green function method to obtain the adiabatic eigenstates and the minimum excitation gap. g min, = O(n 2(exp -n/2), corresponding to the exponential complexity of the algorithm for Number Partitioning. The key element of the analysis is the conditional energy distribution computed for the set of all spin configurations generated from a given (ancestor) configuration by simultaneous flipping of a fixed number of spins. For the problem in question this distribution is shown to depend on the ancestor spin configuration only via a certain parameter related to 'the energy of the configuration. As the result, the algorithm dynamics can be described in terms of one-dimensional quantum diffusion in the energy space. This effect provides a general limitation of a quantum adiabatic computation in random optimization problems. Analytical results are in agreement with the numerical simulation of the algorithm.

Smelyanskiy, V. N.↗

Seasat synthetic-aperture radar data reduction using parallel programmable array processors

This paper presents a digital processing system that produces the Seasat synthetic-aperture radar (SAR) imagery. The system consists of a SEL 32/77 host minicomputer and three AP-120B array processors. The partitioning of the SAR processing functions and the design of software modules is described. The rationale for selecting the parallel array processor architecture and the methodology for developing the parallel processing scheme on this system is described. This system attains a Seasat SAR data reduction speed of 2.5 h per 25-m resolution 4-look and 100 km x 100 km image frame. A preliminary performance evaluation of this parallel processing system and potential future applications for remote sensing data reduction are described.

Wu, C.↗

Aqueous two-phase partition applied to the isolation of plasma membranes and Golgi apparatus from cultured mammalian cells

Partitioning in dextran-poly(ethylene)glycol (PEG) aqueous-aqueous phase systems represents a mature technology with many applications to separations of cells and to the preparation of membranes from mammalian cells. Most applications to membrane isolation and purification have focused on plasma membranes, plasma membrane domains and separation of right side-out and inside-out plasma membrane vesicles. The method exploits a combination of membrane properties, including charge and hydrophobicity. Purification is based upon differential distributions of the constituents in a sample between the two principal compartments of the two phases (upper and lower) and at the interface. The order of affinity of animal cell membranes for the upper phase is: endoplasmic reticulum <mitochondria<Golgi apparatus<lysosomes and endosomes< plasma membranes. Salt concentrations and temperature affect partitioning behavior and must be precisely standardized. In some cases, it is more fortuitous to combine aqueous two-phase partition with other procedures to obtain a more highly purified preparation. A procedure is described for preparation of Golgi apparatus from transformed mammalian cells that combines aqueous two-phase partition and centrifugation. Also described is a periodic NADH oxidase, a new enzyme marker for right side-out plasma membrane vesicles not requiring detergent disruptions for measurement of activity.

Non-NASA Center↗

Demixing of aqueous polymer two-phase systems in low gravity

When polymers such as dextran and poly(ethylene glycol) are mixed in aqueous solution biphasic systems often form. On Earth the emulsion formed by mixing the phases rapidly demixes because of phase density differences. Biological materials can be purified by selective partitioning between the phases. In the case of cells and other particulates the efficiency of these separations appears to be somewhat compromised by the demixing process. To modify this process and to evaluate the potential of two-phase partitioning in space, experiments on the effects of gravity on phase emulsion demixing were undertaken. The behavior of phase systems with essentially identical phase densities was studied at one-g and during low-g parabolic aircraft maneuvers. The results indicate the demixing can occur rather rapidly in space, although more slowly than on Earth. The demixing process was examined from a theoretical standpoint by applying the theory of Ostwald ripening. This theory predicts demizing rates many orders of magnitude lower than observed. Other possible demixing mechanisms are considered.

Bamberger, S.↗

Thermodynamic analysis of hydrogel swelling in aqueous sodium chloride solutions

Here, a comprehensive thermodynamic model is presented for swelling and salt partitioning of poly (n-isopropyl acrylamide) hydrogels in aqueous sodium chloride (NaCl) solutions at 298 K. Together with a Helmholtz energy expression for network elastic energy, a modified electrolyte Nonradom Two-Liquid activity coefficient model is used to represent both the long-range electrostatic interactions and the short-range van der Waals interactions present in the gel phase. The model parameters include one network parameter per hydrogel system plus two binary parameters per interaction pair in the polymer–solvent-salt systems. With semiquantitave agreement with experimental data for both hydrogel swelling and salt partitioning up to salt saturation, the analysis suggests the ion hydration between NaCl and water is the root cause for hydrogel deswelling in aqueous NaCl solution at high NaCl concentrations, i.e., NaCl weight fraction >0.03. On the other hand, the strong repulsive interaction between NaCl and hydrogel polymer is the dominant factor for salt partitioning.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗