Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Correlated Response of Protein Side-Chain Fluctuations and Conformational Entropy to Ligand Binding

The heterogeneous fast side-chain dynamics of proteins plays crucial roles in molecular recognition and binding. Site-specific NMR experiments quantify these motions by measuring the model-free order parameter (O axis 2 ) on a scale of 0 (most flexible) to 1 (least flexible) for each methyl-containing residue of proteins. In this work, we have examined ligand-induced variations in the fast side-chain dynamics and conformational entropy of calmodulin (CaM) using five different CaM–peptide complexes. O axis 2 of CaM in the ligand-free (O axis,U 2 ) and ligand-bound (O axis,B 2 ) states are calculated from molecular dynamics trajectories and conformational energy surfaces obtained using the adaptive biasing force (ABF) method. ΔO axis 2 = O axis,B 2 – O axis,U 2 follows a Gaussian-like unimodal distribution whose second moment is a potential indicator of the binding affinity of these complexes. The probability for the binding-induced O axis,U 2 → O axis,B 2 transition decreases with increasing magnitude of ΔO axis 2 , indicating that large flexibility changes are improbable for side chains of CaM after ligand binding. A linear correlation established between ΔOaxis2 and the conformational entropy change of the protein makes possible the determination of the conformational entropy of binding of protein–ligand complexes. The results not only underscore the functional importance of fast side-chain fluctuations but also highlight key motional and thermodynamic correlates of protein–ligand binding.

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

Active learning of ternary alloy structures and energies

Abstract Machine learning models with uncertainty quantification have recently emerged as attractive tools to accelerate the navigation of catalyst design spaces in a data-efficient manner. Here, we combine active learning with a dropout graph convolutional network (dGCN) as a surrogate model to explore the complex materials space of high-entropy alloys (HEAs). We train the dGCN on the formation energies of disordered binary alloy structures in the Pd-Pt-Sn ternary alloy system and improve predictions on ternary structures by performing reduced optimization of the formation free energy, the target property that determines HEA stability, over ensembles of ternary structures constructed based on two coordinate systems: (a) a physics-informed ternary composition space, and (b) data-driven coordinates discovered by the Diffusion Maps manifold learning scheme. Both reduced optimization techniques improve predictions of the formation free energy in the ternary alloy space with a significantly reduced number of DFT calculations compared to a high-fidelity model. The physics-based scheme converges to the target property in a manner akin to a depth-first strategy, whereas the data-driven scheme appears more akin to a breadth-first approach. Both sampling schemes, coupled with our acquisition function, successfully exploit a database of DFT-calculated binary alloy structures and energies, augmented with a relatively small number of ternary alloy calculations, to identify stable ternary HEA compositions and structures. This generalized framework can be extended to incorporate more complex bulk and surface structural motifs, and the results demonstrate that significant dimensionality reduction is possible in thermodynamic sampling problems when suitable active learning schemes are employed.

Chemistry↗
DOE OSTI2020

Graphene-extracted membrane lipids facilitate the activation of integrin α v β 8

Despite the remarkable electrochemical properties of graphene, strong van der Waals attraction between graphene and biomolecules often causes cytotoxicity, which hinders its applications in the biomedical field. Unfortunately, surface passivation of graphene might stimulate undesired immune response as the nanomaterial triggers cytokine production through membrane receptor activation. In this work, we use all-atom Molecular Dynamics (MD) simulations to unravel the underlying mechanism of graphene-induced inside-out activation of integrin αvβ8, a prominent membrane receptor expressed in immune cells. We model the transmembrane (TM) domains of integrin αvβ8 in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer and observe the structural changes in the integrin–membrane complex when interacting with a graphene nanosheet across the membrane. We find that the β8 TM domain interacts with the graphene nanosheet directly or indirectly through extracted lipids, facilitating the pulling of a β8 subunit away from an αv subunit and thus leading to the disruption of the TM domain association by breaking the hydrophobic cluster in the cytoplasmic domains of the αv and β8 subunits. Alanine substitution of two conserved phenylalanine residues on the αv subunit at this hydrophobic cluster further reveals the importance of a stable T-shaped structure in retaining integrin in its inactive state. Our results agree with previous studies on the interactions between other integrin subtypes and their endogenous activators, suggesting an intriguing role that the graphene nanosheet may play in the integrin-related signal transduction during its interaction with the membrane.

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

Detection of multi-reference character imbalances enables a transfer learning approach for virtual high throughput screening with coupled cluster accuracy at DFT cost

Appropriately identifying and treating molecules and materials with significant multi-reference (MR) character is crucial for achieving high data fidelity in virtual high-throughput screening (VHTS). Despite development of numerous MR diagnostics, the extent to which a single value of such a diagnostic indicates the MR effect on a chemical property prediction is not well established. We evaluate MR diagnostics for over 10 000 transition-metal complexes (TMCs) and compare to those for organic molecules. We observe that only some MR diagnostics are transferable from one chemical space to another. By studying the influence of MR character on chemical properties (i.e., MR effect) that involve multiple potential energy surfaces (i.e., adiabatic spin splitting, ΔE H–L , and ionization potential, IP), we show that differences in MR character are more important than the cumulative degree of MR character in predicting the magnitude of an MR effect. Motivated by this observation, we build transfer learning models to predict CCSD(T)-level adiabatic ΔE H–L and IP from lower levels of theory. By combining these models with uncertainty quantification and multi-level modeling, we introduce a multi-pronged strategy that accelerates data acquisition by at least a factor of three while achieving coupled cluster accuracy (i.e., to within 1 kcal mol –1 MAE) for robust VHTS.

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

Impacts of phase state and water content on secondary organic aerosol formation and partitioning (Final Technical Report)

Secondary organic aerosols (SOA) account for a major fraction of particulate matter in the atmosphere, affecting climate, air quality, and public health. SOA formation and evolution are highly complex processes involving both chemical reactions and transport of molecules in air, at particle surfaces, and within particles. The state of the matter, or “phase state,” comprising SOA can vary from a liquid, over an amorphous semi-solid, to a glassy solid, depending on chemical composition, water content, relative humidity (RH), and temperature. The occurrence of glassy and amorphous semi-solid states can pose limitations on the rate of transport of molecules, affecting gas-particle interactions and challenging the treatment of SOA in atmospheric models. The objective of this project was to improve fundamental understanding of the interplay of the phase state of particles and water content on the evolution of SOA formation. We have developed a method to estimate the glass transition temperature (T g ) of organic compounds to predict the phase state and viscosity of SOA. This method has been applied to high resolution mass spectrometry data of various types of SOA, achieving a good agreement with viscosity measurements. In addition, the viscosity estimation method was implemented into a regional model CMAQ to simulate diel and seasonal variations of SOA phase state over the U.S. Applying kinetic multilayer modeling, we analyzed laboratory experiments and ARM field measurements (e.g., HI-SCALE) on kinetic limitations of amine uptake. We have also conducted numerical simulations to estimate equilibration timescales of SOA partitioning. Overall, we have evaluated the impacts of phase state and water content on SOA lifecycle that should contribute to reducing the uncertainty of SOA representation in regional climate and air quality predictions. The project has yielded 17 publications, in which DE-SC0018349 is explicitly acknowledged. Please see the list of publications as below, followed by summary of research activities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Zentropy Theory for Transformative Functionalities of Magnetic and Superconducting Materials

The proposed research developed the zentropy theory through applications to complex magnetic materials and superconductors under the hypothesis that the emergent properties of complex magnetic materials and superconductors can be predicted by statistical mechanics of ergodic microstates with their partition functions computed from DFT-predicted free energies. The key objective is to develop approaches to systematically determine the types and number of microstates and the supercell size in DFT-based calculations through convergency of macroscopic functionalities, with the incorporation of our mixed-space approach accounting for the interactions between periodic supercells. In addition to use scientific intuitions to guide the design of important microstates, the key innovation of the proposed research is to integrate the domain knowledge and the material-property-descriptor database (MPDD) with 4 million microstates, which is supported by our deep neural network machine learning models (SIPFENN: structure-informed prediction of formation energy using neural networks) and integrated with our high throughput DFT Tool Kit (DFTTK). For complex magnetic materials, one of the objectives is to develop approaches to calculate short-range ordering from the statistical distribution of each microstate. For superconductors, the divergency of quasiparticle effective mass at a quantum critical point will be investigated, and the superconducting and non-superconducting microstates will be delineated through analysis of electronic band structure, density of states, charge density, and Fermi surface.

36 MATERIALS SCIENCE↗
NASA NTRS1977

On the origin of stable remanence in pseudo-single domain grains

A critique is presented of the quantitative model for the magnetic moment of pseudo-single domain grains (hypothetical magnetite grains larger than the critical size threshold for single domain behavior, yet also difficult to demagnetize), derived by Stacey and Banerjee (1974). Evidence from theoretical studies in micromagnetics demonstrates that the spin orientations in such grains are too complex to permit ready prediction of the magnetic moments. However, this limitation in the theory may be overcome by experiments involving rare earth-cobalt alloys and yttrium iron garnet crystals; these studies have suggested that surface anisotropy is the predominant cause of the high coercivity of pseudo-single domain grains.

Banerjee, S. K.↗
NASA NTRS1985

Impurity Extraction by Droplets

The goals are to model and to measure the phase equilibrium properties of a finely divided fluid containing a large number of chemically similar species. The objective is to develop an accurate, usable model for such phenomena as pollutant extraction of rain clouds, industrial separation in spray towers, and separation in emulsions. The project was designed as a hierarchy of complementary theoretical and experimental steps. A theory was developed to describe the segregation of complex impurities at the interface of a solvent. This phenomenon is important in phase behavior when a large fraction of molecules in a material are near an interface, the situation in a finely divided material. The theory will be modified to account for the effect of surface curvature on the surface tension. The study of mixtures differs from pure fluids not only because of the surface effects but also because composition differences between the droplet and the surrounding vapor can stabilize a droplet with respect to a bulk phase.

Morrison, G.↗
NASA NTRS1989

Advances in fluid dynamics

Various papers on advances in fluid dynamics are presented. Individual topics addressed include: nonunique solutions of the Euler equations; numerical solution of the biharmonic equation using a spectral multigrid method; finite-volume shape and placement for efficiency, accuracy, and physics; NND schemes and their application to numerical simulation of complex plume flow problems; symbolic analysis and computational algorithm construction; and velocity and mass flux distribution measurements of spherical glass beads in air flow in a 90-deg vertical-to-horizontal bend. Also discussed are: hyperbolic solution adaptive mesh generation; progress and future directions in CFD; numerical simulation of shock focusing; microbubble formation and splitting in a turbulent boundary layer for turbulence reduction; numerical study of spherical Taylor-Couette flow, trajectories of vortex lines beneath a free surface or above a plane; and a mathematical model for response of the coronary circulation to high sustained gravitational force fields.

Ballhaus, W. F., Jr.↗
NASA NTRS1992

Hypersonic aerodynamic/aerothermodynamic testing capabilities at Langley Research Center

Capabilities of the Langley Hypersonic Facilities Complex that consists of nine blowdown-to-vacuum wind tunnels are reviewed. These tunnels complement one another to provide ranges of Mach number from 6 to 20, unit Reynolds number from 0.03 to 40 million per ft, and normal shock density ratio from 2 to 12. Topics under consideration include the basic components of each facility, status, and upgrades recently performed; instrumentation and testing techniques routinely used to measure forces and moments, surface pressures, surface temperature-time histories, flow properties within the shock layer about the model via survey probes, and flow visualization; and data acquisition systems and upgrades performed via a NASA revitalization program.

Miller, C. G., III↗
NASA NTRS1999

Bubble Generation in a Flowing Liquid Medium and Resulting Two-Phase Flow in Microgravity

The present investigation reports a study of bubble generation under reduced gravity conditions, using both a co-flow and a cross-flow configuration. This study may be used in the conceptual design of a space-based thermal management system. Ensuing two-phase flow void fraction can be accurately monitored using a single nozzle gas injection system within a continuous liquid flow conduit, as utilized in the present investigation. Accurate monitoring of void fraction leads to precise control of heat and mass transfer coefficients related to a thermal management system; hence providing an efficient and highly effective means of removing heat aboard spacecraft or space stations. Our experiments are performed in parabolic flight aboard the modified DC-9 Reduced Gravity Research Aircraft at NASA Lewis Research Center, using an air-water system. For the purpose of bubble dispersion in a flowing liquid, we use both a co-flow and a cross-flow configuration. In the co-flow geometry, air is introduced through a nozzle in the same direction with the liquid flow. On the other hand, in the cross-flow configuration, air is injected perpendicular to the direction of water flow, via a nozzle protruding inside the two-phase flow conduit. Three different flow conduit (pipe) diameters are used, namely, 1.27 cm, 1.9 cm and 2.54 cm. Two different ratios of nozzle to pipe diameter (D(sub N))sup * are considered, namely (D(sub N))sup * = 0.1 and 0.2, while superficial liquid velocities are varied from 8 to 70 cm/s depending on flow conduit diameter. It is experimentally observed that by holding all other flow conditions and geometry constant, generated bubbles decrease in size with increase in superficial liquid velocity. Detached bubble diameter is shown to increase with air injection nozzle diameter. Likewise, generated bubbles grow in size with increasing pipe diameter. Along the same lines, it is shown that bubble frequency of formation increases and hence the time to detachment of a forming bubble decreases, as the superficial liquid velocity is in-creased. Furthermore, it is shown that the void fraction of the resulting two-phase flow increases with volumetric gas flow rate Q(sub d), pipe diameter and gas injection nozzle diameter, while they decrease with surrounding liquid flow. The important role played by flowing liquid in detaching bubbles in a reduced gravity environment is thus emphasized. We observe that the void fraction can be accurately controlled by using single nozzle gas injection, rather than by employing multiple port injection, since the later system gives rise to unpredictable coalescence of adjacent bubbles. It is of interest to note that empirical bubble size and corresponding void fraction are somewhat smaller for the co-flow geometry than the cross-flow configuration at similar flow conditions with similar pipe and nozzle diameters. In order to supplement the empirical data, a theoretical model is employed to study single bubble generation in the dynamic (Q(sub d) = 1 - 1000 cu cm/s) and bubbly flow regime within the framework of the co-flow configuration. This theoretical model is based on an overall force balance acting on the bubble during the two stages of generation, namely the expansion and the detachment stage. Two sets of forces, one aiding and the other inhibiting bubble detachment are identified. Under conditions of reduced gravity, gas momentum flux enhances, while the surface tension force at the air injection nozzle tip inhibits bubble detachment. In parallel, liquid drag and inertia can act as both attaching and detaching forces, depending on the relative velocity of the bubble with respect to the surrounding liquid. Predictions of the theoretical model compare well with our experimental results. However, at higher superficial liquid velocities, as the bubble loses its spherical form, empirical bubble size no longer matches the theoretical predictions. In summary, we have developed a combined experimental and theoretical work, which describes the complex process of bubble generation and resulting two-phase flow in a microgravity environment. Results of the present study can be used in a wide range of space-based applications, such as thermal energy and power generation, propulsion, cryogenic storage and long duration life support systems, necessary for programs such as NASA's Human Exploration for the Development of Space (HEDS).

Pais, S. C.↗
NASA NTRS2002

A Mechanistic Study of Nucleate Boiling Under Microgravity Conditions

The overall objective of this work is to study nucleate boiling heat transfer under microgravity conditions in such a way that while providing basic knowledge of the phenomena, it also leads to development of simulation models and correlations that can be used as design tools for a wide range of gravity levels. In the study a building block type of approach is used and both pool and low velocity flow boiling are investigated. Starting with experiments using a single bubble, the complexity of the experiments is increased to two or three inline bubbles, to five bubbles placed on a two-dimensional grid. Finally, experiments are conducted where a large number of prescribed cavities nucleate on the heater and when a commercial surface is used. So far experiments have been conducted at earth normal gravity and in the reduced gravity environment of the KC-135 aircraft whereas experiments on the space station are planned. Modeling/complete numerical simulation of the boiling process is an integral part of the total effort. Experiments conducted with single bubbles formed on a nucleation site microfabricated on a polished silicon wafer show that for gravity levels (g) varying from 1.5g(sub e) to 0.01g(sub e), the bubble diameter at departure varies approximately as (g(sub e)/g)(exp 1/2) and the growth period as (g(sub e)/g). When bubbles merge either inline or in a plane, the bubble diameter at departure is found to be smaller than that obtained for a single bubble and shows a weaker dependence on the level of gravity. The possible reason is that as the bubbles merge they create fluid circulation around the bubbles, which in turn induces a lift force that is responsible for the earlier departure of the bubbles. The verification of this proposition is being sought through numerical simulations. There is a merger of two inline, three inline, and several bubbles in a plane in the low gravity environment of the KC-135 aircraft. After merger and before departure, a mushroom type of bubble with several stems attached to the heater surface is clearly evident. Local heat fluxes during growth and departure of a single bubble were also measured. It was found that during most of the growth period of the bubble, generally the wall heat flux decreased with time because of the increased dry area under the bubble. However, the heat flux increased rapidly just prior to departure of the bubble because of the transient conduction into the cold liquid rushing to fill the space vacated by the bubble as the bubble base shrinks. The measured heat fluxes at various radial locations are found to be in qualitative agreement with the numerical predictions. Single bubble studies at earth normal gravity have also been performed on surfaces oriented at different angles to the gravitational acceleration with flow parallel to the surface. It is found that in all cases the bubbles slide along the surface before lift-off from the surface. The lift force generated as a result of the relative motion between the sliding bubbles and the imposed flow is found to play an important role when the normal force due to buoyancy is reduced. An experimental apparatus for the study of the bubble behavior with imposed flow under reduced gravity conditions has been developed and will soon be employed for experiments in the KC-135 aircraft.

Dhir, V. K.↗
NASA NTRS2010

Introduction to Vector Field Visualization

Vector field visualization techniques are essential to help us understand the complex dynamics of flow fields. These can be found in a wide range of applications such as study of flows around an aircraft, the blood flow in our heart chambers, ocean circulation models, and severe weather predictions. The vector fields from these various applications can be visually depicted using a number of techniques such as particle traces and advecting textures. In this tutorial, we present several fundamental algorithms in flow visualization including particle integration, particle tracking in time-dependent flows, and seeding strategies. For flows near surfaces, a wide variety of synthetic texture-based algorithms have been developed to depict near-body flow features. The most common approach is based on the Line Integral Convolution (LIC) algorithm. There also exist extensions of LIC to support more flexible texture generations for 3D flow data. This tutorial reviews these algorithms. Tensor fields are found in several real-world applications and also require the aid of visualization to help users understand their data sets. Examples where one can find tensor fields include mechanics to see how material respond to external forces, civil engineering and geomechanics of roads and bridges, and the study of neural pathway via diffusion tensor imaging. This tutorial will provide an overview of the different tensor field visualization techniques, discuss basic tensor decompositions, and go into detail on glyph based methods, deformation based methods, and streamline based methods. Practical examples will be used when presenting the methods; and applications from some case studies will be used as part of the motivation.

Kao, David↗
NASA NTRS2016

Effect of Molecular Flexibility upon Ice Adhesion Shear Strength

Ice formation on aircraft surfaces effects aircraft performance by increasing weight and drag leading to loss of lift. Current active alleviation strategies involve pneumatic boots, heated surfaces, and usage of glycol based de-icing fluids. Mitigation or reduction of in-flight icing by means of a passive approach may enable retention of aircraft capabilities, i.e., no reduction in lift, while reducing the aircraft weight and mechanical complexity. Under a NASA Aeronautics Research Institute Seedling activity, the effect of end group functionality and chain length upon ice adhesion shear strength (IASS) was evaluated with the results indicating that chemical functionality and chain length (i.e. molecular flexibility) affected IASS. Based on experimental and modeling results, diamine monomers incorporating molecular flexibility as either a side chain or in between diamine functionalities were prepared, incorporated into epoxy resins that were subsequently used to fabricate coatings on aluminum substrates, and tested in a simulated icing environment. The IASS was found to be lower when molecular flexibility was incorporated in the polymer chain as opposed to a side chain.

Smith, Joseph G.↗
NASA NTRS1985

The distortion of a supersonic turbulent boundary layer by bulk compression and surface curvature

An experimental investigation is conducted of the effects of bulk compression and streamline curvature on a turbulent supersonic boundary layer whose incoming Mach number is 2.87. Two different, constant-radius curved wall models were used; the total turning angle in both cases was fixed at 8 deg. While the amplification of the normal stress was approximately the same in both models, the amplification of the shear stress was greater in the one with delta(o)/R value 0.1 than in the model with a value of 0.02. The flow over the latter model appears to be especially complex.

Smits, A. J.↗
DOE OSTI2026

Neural network reconstruction of the DIII-D tokamak plasma boundary using a reduced set of diagnostics

This study investigates the feasibility of reconstructing the last closed flux surface in the DIII-D tokamak using neural network models trained on reduced input feature sets, addressing an ill-posed task. Two models are compared: one trained solely on coil currents and another incorporating coil currents, plasma current and loop voltage. The model trained exclusively on coil currents achieved a mean point displacement of $0.04$ m on a held-out test set, while the inclusion of plasma current and loop voltage reduced the error to $0.03$ m. This comparison highlights the trade-offs between input feature complexity and reconstruction accuracy, demonstrating the potential of machine learning algorithms to perform effectively in data-limited environments, such as those expected in fusion power plants due to diagnostic constraints imposed by the presence of blankets and shielding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Metal Independent Correlations for Site-specific Binding Energies of Relevant Catalytic Intermediates

Establishing energy correlations among different metals can accelerate the discovery of efficient and cost-effective catalysts for complex reactions. Using a recently introduced coordination-based model, we can predict site specific metal binding energies (∆EM) that can be used as a descriptor for chemical reactions. In this study, we have examined a range of metals including Ag, Au, Co, Cu, Ir, Ni, Os, Pd, Pt, Rh, and Ru and found linear correlations between predicted ∆EM and adsorption energies of CH and O (∆ECH and ∆EO) at various coordination environment for all the considered metals. Interestingly, all the metals correlate with one another under specific surface site coordination, indicating that different metals are interrelated in a particular coordination environment. Furthermore, we have tested and verified for PtPd and PtIr-based alloys that they follow a similar behavior. Moreover, we have expanded the metal space by taking some early transition metals along with few s-block metals and shown a cyclic behavior of the adsorbate binding energy (∆EA) versus ∆EM. Therefore, ∆ECH and ∆EO can be efficiently interpolated between metals, alloys and intermetallics based on information related to one metal only. This simplifies the process of screening new metal catalyst formulations and their reaction energies.

Density function theory↗
DOE OSTI2021

Analysis of the Reactive CO 2 Surface Flux in Electrocatalytic Aqueous Flow Reactors

We study how the mass transfer of a chemical species in shear flow is suppressed by the production of a second species that partially reacts away the first species. The second species is produced at the surface with the first species as a reactant. Our work is directly motivated by electrochemical CO 2 reduction in aqueous flow reactors, where OH– molecules generated by the CO evolution reaction react away CO 2 molecules, ultimately inhibiting the mass transfer of CO 2 to the cathode surface. We derive a simple approximation of the Sherwood number, a dimensionless measure of the mass flux into the surface, as a function of the Péclet number, surface Damköhler number, and bulk Damköhler number, and we compare our approximation to numerical solutions of the governing equations. We find that in the diffusion-limited regime, the Sherwood number is well described by the classical Graetz-Lévêque result reduced by a reaction factor due to the competing bulk reaction; with the stoichiometry relevant to electrochemical CO 2 reduction, this reaction factor is 1/2. While the model problem we solve provides insight into how OH– production affects CO 2 mass transfer, experimental systems often involve more complex chemistry. We thus also show how a common buffered electrolyte, KHCO 3 , affects our results.

30 DIRECT ENERGY CONVERSION↗