Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differentiable physics”

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,261 records · Page 70

Cooling histories of Apollo 15 quartz-normative basalts

Controlled cooling rate studies of vitrophyre 15597 are used to calculate the physical properties of a quartz-normative magma as well as cooling histories of Apollo 15 quartz-normative basalts (QNB). Pyroxene phenocrysts of a QNB magma do not settle appreciably during the cooling of the flow unit; this behavior is consistent with the small degree of chemical differentiation that occurred between members of the suite. Two types of cooling histories for QNB are discerned. The vitrophyres experienced initial slow cooling followed by rapid cooling, and the coarse-grained microgabbros underwent slow nearly linear cooling histories. Textural evidence suggests that lunar vitrophyres crystallized near the top of the flow. It is suggested that the pattern of cooling may have been caused by the chilling of liquids that inherited early crystals in the interior of the flow.

Grove, T. L.↗

Lubrication of Machine Elements

The understanding of hydrodynamic lubrication began with the classical experiments of Tower and Petrov. Reynolds used a reduced form of the Navier-Stokes equations and the continuity equation to generate a second order differential equation for the pressure in the narrow, converging gap of a bearing contact. Such a pressure enables a load to be transmitted between the surfaces with very low friction since the surfaces are completely separated by a film of fluid. In such a situation it is the physical properties of the lubricant, notably the dynamic viscosity, that dictate the behavior of the contact. The understanding of boundary lubrication is normally attributed to Hardy and Doubleday. In boundary lubrication it is the physical and chemical properties of thin films of molecular proportions and the surfaces to which they are attached that determine contact behavior. The lubricant viscosity is not an influential parameter. Research is devoted to a better understanding and more precise definition of other lubrication regimes between these extremes. One such regime, elastohydrodynamic lubrication, occurs in nonconformal contacts, where the pressures are high and the bearing surfaces deform elastically. In this situation the viscosity of the lubricant may raise considerably, and this further assists the formation of an effective fluid film. The science of these three lubrication regimes (hydrodynamic, elastohydrodynamic, and boundary) are described and the manner in which this science is used in the design of machine elements is examined.

Hamrock, Bernard J.↗

Exotic containers for capillary surfaces

This paper discusses 'exotic' rotationally symmetric containers that admit an entire continuum of distinct equilibrium capillary free surfaces. The paper extends earlier work to a larger class of parameters and clarifies and simplifies the governing differential equations, while expressing them in a parametric form appropriate for numerical integration. A unified presentation suitable for both zero and nonzero gravity is given. Solutions for the container shapes are depicted graphically along with members of the free-surface continuum, and comments are given concerning possible physical experiments.

Concus, Paul↗

Computation of unsteady shock-induced combustion using logarithmic species conservation equations

Numerical simulations are used to investigate periodic combustion instabilities observed in ballistic-range experiments of blunt bodies flying at supersonic speeds through hydrogen-air mixtures. The computations are validated by comparing experimental shadowgraphs with shadowgraphs created from the computed flowfields and by comparing the experimentally measured instability frequencies with computed frequencies. The numerical simulations use a logarithmic transformation of the species conservation equations as a way to reduce the grid requirements for computing shock-induced combustion. The transformation is applied to the Euler equations coupled to a detailed hydrogen-air chemical reaction mechanism with 13 species and 33 reactions. The resulting differential equations are solved using a finite volume formulation and a two-step predictor-corrector scheme to advance the solution in time. Results are presented and compared for both a flux-vector splitting scheme and an upwind TVD scheme. The computations add insight to the physical processes observed in the experiments and the numerical methods needed to simulate them. The usefulness of the ballistic-range experiments for the validation of numerical techniques and chemical kinetic models is also demonstrated.

Wilson, Gregory J.↗

The combustion program at CTR

Understanding and modeling of turbulent combustion are key problems in the computation of numerous practical systems. Because of the lack of analytical theories in this field and of the difficulty of performing precise experiments, direct numerical simulation (DNS) appears to be one of the most attractive tools to use in addressing this problem. The general objective of DNS of reacting flows is to improve our knowledge of turbulent combustion but also to use this information for turbulent combustion models. For the foreseeable future, numerical simulation of the full three-dimensional governing partial differential equations with variable density and transport properties as well as complex chemistry will remain intractable; thus, various levels of simplification will remain necessary. On one hand, the requirement to simplify is not necessarily a handicap: numerical simulations allow the researcher a degree of control in isolating specific physical phenomena that is inaccessible in experiments. CTR has pursued an intensive research program in the field of DNS for turbulent reacting flows since 1987. DNS of reacting flows is quite different from DNS of non-reacting flows: without reaction, the equations to solve are clearly the five conservation equations of the Navier Stokes system for compressible situations (four for incompressible cases), and the limitation of the approach is the Reynolds number (or in other words the number of points in the computation). For reacting flows, the choice of the equations, the species (each species will require one additional conservation equation), the chemical scheme, and the configuration itself is more complex.

Poinsot, Thierry J.↗

Predicting nonequilibrium Green’s function dynamics and photoemission spectra via nonlinear integral operator learning

Understanding the dynamics of nonequilibrium quantum many-body systems is an important research topic in a wide range of fields across condensed matter physics, quantum optics, and high-energy physics. However, numerical studies of large-scale nonequilibrium phenomena in realistic materials face serious challenges due to intrinsic high-dimensionality of quantum many-body problems and the absence of time-invariance. The nonequilibrium properties of many-body systems can be described by the dynamics of the correlator, or the Green's function of the system, whose time evolution is given by a high-dimensional system of integro-differential equations, known as the Kadanoff–Baym equations (KBEs). The time-convolution term in KBEs, which needs to be recalculated at each time step, makes it difficult to perform long-time numerical simulation. In this paper, we develop an operator-learning framework based on recurrent neural networks (RNNs) to address this challenge. We utilize RNNs to learn the nonlinear mapping between Green's functions and convolution integrals in KBEs. By using the learned operators as a surrogate model in the KBE solver, we obtain a general machine-learning scheme for predicting the dynamics of nonequilibrium Green's functions. Besides significant savings per each time step, the new methodology reduces the temporal computational complexity from $O(N_t^3)$ to $O(N_t)$ where N t is the number of steps taken in a simulation, thereby making it possible to study large many-body problems which are currently infeasible with conventional KBE solvers. Through various numerical examples, we demonstrate the effectiveness of the operator-learning based approach in providing accurate predictions of physical observables such as the reduced density matrix and time-resolved photoemission spectra. Moreover, our framework exhibits clear numerical convergence and can be easily parallelized, thereby facilitating many possible further developments and applications.

97 MATHEMATICS AND COMPUTING↗

Multiamorphous Phases in Diketopyrrolopyrrole-Based Conjugated Polymers: From Bulk to Ultrathin Films

The glass transition temperature (T g ) of conjugated polymers is a crucial physical parameter that governs their mechanical and electrical properties for applications from sensor technology to skin-like electronics. Despite some well-developed methods that exist for detecting T g of conjugated polymers, there is still a need for the combination of advanced calorimetric methodologies to characterize T g in a broad range of temperatures and scales from bulk to nanometer-thin films. In this work, alternating current (ac) chip calorimetry and flash differential scanning calorimetry (flash DSC) are applied to successfully identify the Tgs and segmental relaxations of conjugated polymers. The detected relaxations were interpreted, in agreement with existing observations in similar systems, as a side-chain relaxation (low temperature), mobile amorphous fractions (MAF), and rigid amorphous fractions (RAF). T g s of MAF and RAF increase and decrease with increasing size of fused thiophene units in the backbone, respectively. The T g s of diketopyrrolopyrrole-based conjugated polymer thin films are also identified by the same method, with both T g s showing a decreased value compared to those of bulk samples. Furthermore, the weight fractions of MAF and RAF with various film thicknesses are qualitatively analyzed based on the heat capacity step of each amorphous fraction during the glass transition. This facile method gives insights into the physical nature of amorphous phases in conjugated polymers for better revealing their structure-property relationships, which will guide the designs of new materials with optimized electrical and mechanical properties in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A combinatorially complete epistatic fitness landscape in an enzyme active site

Protein engineering often targets binding pockets or active sites which are enriched in epistasis—nonadditive interactions between amino acid substitutions—and where the combined effects of multiple single substitutions are difficult to predict. Few existing sequence-fitness datasets capture epistasis at large scale, especially for enzyme catalysis, limiting the development and assessment of model-guided enzyme engineering approaches. We present here a combinatorially complete, 160,000-variant fitness landscape across four residues in the active site of an enzyme. Assaying the native reaction of a thermostable β-subunit of tryptophan synthase (TrpB) in a nonnative environment yielded a landscape characterized by significant epistasis and many local optima. These effects prevent simulated directed evolution approaches from efficiently reaching the global optimum. There is nonetheless wide variability in the effectiveness of different directed evolution approaches, which together provide experimental benchmarks for computational and machine learning workflows. The most-fit TrpB variants contain a substitution that is nearly absent in natural TrpB sequences—a result that conservation-based predictions would not capture. Thus, although fitness prediction using evolutionary data can enrich in more-active variants, these approaches struggle to identify and differentiate among the most-active variants, even for this near-native function. Overall, this work presents a large-scale testing ground for model-guided enzyme engineering and suggests that efficient navigation of epistatic fitness landscapes can be improved by advances in both machine learning and physical modeling.

biocatalysis↗

TNet: A Model-Constrained Tikhonov Network Approach for Inverse Problems

Deep learning (DL), in particular deep neural networks, by default is purely data-driven and in general does not require physics. This is the strength of DL but also one of its key limitations when applied to science and engineering problems in which underlying physical properties—such as stability, conservation, and positivity—and accuracy are required. DL methods in their original forms are often not capable of respecting the underlying mathematical models or achieving desired accuracy even in big-data regimes. On the other hand, many data-driven science and engineering problems, such as inverse problems, typically have limited experimental or observational data, and DL would overfit the data in this case. Leveraging information encoded in the underlying mathematical models, we argue, not only compensates for missing information in low data regimes but also provides opportunities to equip DL methods with the underlying physics, hence promoting better generalization. This paper develops a model-constrained DL approach and its variant TNet—a Tikhonov neural network—which are capable of learning not only information hidden in the training data but also in the underlying mathematical models to solve inverse problems governed by partial differential equations in low data regimes. We provide the constructions and some theoretical results for the proposed approaches for both linear and nonlinear inverse problems. Since TNet is designed to learn inverse solutions with Tikhonov regularization, it is interpretable: in fact it recovers Tikhonov solutions for linear cases while potentially approximating Tikhonov solutions for nonlinear inverse problems. We also prove that data randomization can enhance not only the smoothness of the networks but also their generalizations. Comprehensive numerical results confirm the theoretical findings and show that with even as little as 1 training data sample for one-dimensional (1D) deconvolution, 5 for an inverse 2D heat conductivity problem, 100 for inverse initial conditions for a time-dependent 2D Burgers’s equation, and 50 for inverse initial conditions for 2D Navier–Stokes equations, TNet solutions can be as accurate as Tikhonov solutions while being several orders of magnitude faster. Furthermore, this is possible owing to the model-constrained term, replications, and randomization.

97 MATHEMATICS AND COMPUTING↗

Characterization of Soil Thermal and Electrical Properties along Multiple Hillslope Transects at Teller Road Site, Seward Peninsula, Alaska, 2017

This dataset has been acquired along five-119 m long transects located on the bottom part of the watershed hillslope at the NGEE Arctic Teller Road site at mile marker 27 (TL_MM27) on the Seward Peninsula, Alaska in July and September 2017. The Distributed Temperature Profiling (DTP) system dataset consist in vertically-resolved profile of soil temperature covering the top 0.8 m of soil with 8 cm interval. In addition to DPT data, electrical resistivity tomography (ERT) data, soil moisture, depth to rock or thaw layer thickness (no differentiation) and ground elevations data have been acquired along each of the transects. A UAV-based geotiff mosaic of the investigated site is also provided. The four data types provided with this dataset of 37 files (*.csv, *.tif, *.DATA): (1) soil temperature profiles, (2) ERT data, (3) the physical measurements of the thaw layer, and (4) an orthomosaic GeoTIFF of the transect study area.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Efficient stream distributions in radiative transfer theory

This paper discusses a new, computationally-efficient method for approximating the integro-differential equation of radiative transfer with a finite set of coupled differential equations for discrete streams. The method uses recommended spatial distributions of streams that are quite different from those typically used in that they are based on the symmetry of several regular Platonic solids. To facilitate the use of such distributions, an explicit, one-parameter relationship between the physical radiance and the abstract stream is formulated. The parameter is used to determine the minimum number of streams required in the radiative transfer model. Accuracy and computational efficiency are shown to be served best by choosing a stream distribution that is invariant to a large number of three space rotations. For various values of the above-mentioned parameter, the resulting recommended stream distribution is shown to be more computationally efficient than more conventional stream distributions. Finally, the incorporation of polarization in the stream definition is described.

Whitney, C.↗

Axial-flow compressor turning angle and loss by inviscid-viscous interaction blade-to-blade computation

A method for computation of the flow field around an arbitrary airfoil cascade on an axially symmetric blade-to-blade surface was developed which takes into account the development and separation of the blade surface boundary layers and mixing in the wake. The method predicts the overall fluid turning and total pressure loss in the context of an inviscid-viscous interaction scheme. The inviscid flow solution is obtained from a compressible flow matrix method. The viscous flow is obtained from a differential boundary layer method which calculates laminar, transitional and turbulent boundary layers. Provisions for the calculation of laminar and turbulent separation regions were added to the viscous scheme. The combined inviscid-viscous interaction scheme described yields results which are quantitatively consistent with experimental data. This suggests that the physical basis for the interactive system is correct and justifies continued exploration and use of the method.

Hansen, E. C.↗

An efficient and robust algorithm for two dimensional time dependent incompressible Navier-Stokes equations: High Reynolds number flows

An algorithm is presented for unsteady two-dimensional incompressible Navier-Stokes calculations. This algorithm is based on the fourth order partial differential equation for incompressible fluid flow which uses the streamfunction as the only dependent variable. The algorithm is second order accurate in both time and space. It uses a multigrid solver at each time step. It is extremely efficient with respect to the use of both CPU time and physical memory. It is extremely robust with respect to Reynolds number.

Goodrich, John W.↗

An efficient and robust algorithm for two dimensional time dependent incompressible Navier-Stokes equations - High Reynolds number flows

An algorithm is presented for unsteady two-dimensional incompressible Navier-Stokes calculations. This algorithm is based on the fourth order partial differential equation for incompressible fluid flow which uses the streamfunction as the only dependent variable. The algorithm is second order accurate in both time and space. It uses a multigrid solver at each time step. It is extremely efficient with respect to the use of both CPU time and physical memory. It is extremely robust with respect to Reynolds number.

Goodrich, John W.↗

Performance and durability of anion exchange membrane water electrolyzers using down-selected polymer electrolytes

Over the last decade, several stable anion exchange polymer electrolytes have been developed for electrochemical devices. Herein, we investigate how chemical structure and physical properties of polymer electrolytes affect performance and durability of anion exchange membrane water electrolyzers (AEMWEs). We select polymer electrolytes with high alkaline stability and consider their polymer properties including conductivity, mechanical/chemical stability, and material interactions to interpret the performance and durability of AEMWEs. Here, the AEMWE with a poly(phenylene) membrane and a poly(fluorene) ionomeric binder exhibited the best performance among those tested in this study; the AEMWE showed ~1 A cm –2 at 2 V under 1 wt% K 2 CO 3 -fed conditions. The voltage degradation rate was 270–550 μV h –1 for several hundred operating hours at a constant current density of 750 mA cm –2 and a differential pressure of 100 pounds per square inch gauge. Based on these results, we discuss research needs of polymer electrolytes for practical AEMWEs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Crystal Structure Evolution of UCl 3 from Room Temperature to Melting

We report uranium trichloride (UCl 3 ) is actively researched to develop and improve applications ranging from molten salt reactors to actinide processing, including spent fuel reprocessing. Here, we report for the first time the crystal structure evolution between room temperature and melting point from in situ high-temperature neutron diffraction to quantify, for example, the thermal expansion of the hexagonal a and c lattice parameters. The results are compared with density functional theory calculations. The melting point of UCl 3 is determined by differential scanning calorimetry to be 1108.2 ± 0.2 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First-Principles Approach to Extracting Chemical Information from X-ray Absorption Near-Edge Spectra of Ga-Containing Materials

The X-ray absorption near-edge structure (XANES) can provide uniquely detailed information on the coordination environments of important Ga-containing materials with unknown structures, including catalytically active materials. In this work, the Ga K-edge XANES was simulated using first-principles-based methods for seven molecular Ga complexes as well β-Ga 2 O 3 , in order to explore the chemical origins of the experimentally observed features. The theoretical spectra were computed using FEFF, CASTEP, and StoBe, to assess the sensitivity of the results to the computational approach. While the XANES features depend on the Ga coordination environment, they are also sensitive to the electronegativity of the ligands and the symmetry at Ga. The white line position responds to changes in both the core state (due to differential screening) and the valence “p” states (arising from differences in ligand coordination).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improved Self-Assembly of P3HT with Pyrene-Functionalized Methacrylates

A block copolymer with discotic liquid crystalline behavior was synthesized using Grignard metathesis polymerization (GRIM) and initiators for continuous activator regeneration atom transfer radical polymerization (ICAR-ATRP). A novel discotic liquid crystalline mesogen, 6-(pyren-1-yloxy)hexyl methacrylate (PyMA), comprises a block that is attached to regioregular poly(3- hexylthiophene) (rr-P 3 HT) generated by GRIM and subjected to end-group modification. Due to the continuous regeneration of Cu + in the reaction mixture in ICAR-ATRP compared to conventional methods, the synthesis was successfully performed with less catalyst. The purity and yield of the final product are increased by eliminating rigorous post-synthesis purification. Stacked pyrene units have contributed to the enhanced long-range π-π interactions and aligning of the P 3 HT block as observed in thin-film X-ray diffraction (XRD). Furthermore, field-effect mobilities in the order of 10 -2 cm 2 V -1 s -1 in bottom-gate, top-contact organic field-effect transistors (OFETs) suggest an enhancement in charge transport due to the discotic electron-rich pyrene units that help mitigate the insulating effect of the methacrylate backbone. The formation of uniform microdomains of P 3 HT-bpoly(PyMA) observed with tapping mode atomic force microscopy (TMAFM) on the channel regions of OFETs indicates the unique packing of the block copolymer in comparison to pristine P 3 HT. Thermotropic properties of the novel discotic mesogen in the presence and absence of P3HT were observed with both the poly(3-hexylthiophene)-b-poly(6-(pyren-1-yloxy)hexyl methacrylate) (P3HT-b-poly(PyMA)) block copolymer and poly(6-(pyren-1-yloxy)hexyl methacrylate) (poly(PyMA)) homopolymer using polarized optical microscopy (POM) and differential scanning calorimetry (DSC).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗