Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Implicit functions”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Numerical Method for Fractional Advection-Diffusion Equation with Heredity

We propose a method of construction of difference schemes for fractional partial differential equations with delay in time. For the fractional equation with two-sided diffusion, fractional transfer in time, and a functional aftereffect, we construct an implicit difference scheme. We use the shifted Grünwald–Letnikov formulas for the approximation of fractional derivatives with respect to spatial variables and the L1-algorithm for the approximation of fractional derivatives in time. Also we use piecewise constant interpolation and extrapolation by extending the discrete prehistory of the model in time. The algorithm is a fractional analog of a purely implicit method; on each time step, it is reduced to the solution of linear algebraic systems. We prove the stability of the method and find its order of convergence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Fracture-based shape optimization built upon the topological derivative

In Silva et al. (2011) and Alidoost et al. (2020), the authors developed an approximation of the energy release rate field associated with a small edge or surface crack at any boundary location and with any orientation using the topological derivative. The approximation is computationally attractive because it requires only a single analysis on the non-cracked domain in contrast with conventional boundary-element and finite-element-based methods, which require a separate and costlier analysis for each crack length-location-orientation combination. Here, a shape optimization scheme for fracture-resistant structures is developed using the energy release rate approximation. In the gradient-based optimization scheme, the domain and its boundary are defined implicitly using level-set functions. The level-set functions of arbitrary geometries are constructed using Boolean operations from the level-set functions of simple primitives. This geometrical representation has the dual advantage of (i) allowing shapes to intersect and/or separate during the optimization and (ii) simplifying the computation of the shape sensitivities.

42 ENGINEERING↗
DOE OSTI2022

Quantum Orbital Minimization Method for Excited States Calculation on a Quantum Computer

Herein we propose a quantum-classical hybrid variational algorithm, the quantum orbital minimization method (qOMM), for obtaining the ground state and low-lying excited states of a Hermitian operator. Given parametrized ansatz circuits representing eigenstates, qOMM implements quantum circuits to represent the objective function in the orbital minimization method and adopts a classical optimizer to minimize the objective function with respect to the parameters in ansatz circuits. The objective function has an orthogonality constraint implicitly embedded, which allows qOMM to apply a different ansatz circuit to each input reference state. We carry out numerical simulations that seek to find excited states of H 2 , LiH, and a toy model consisting of four hydrogen atoms arranged in a square lattice in the STO-3G basis with UCCSD ansatz circuits. Comparing the numerical results with existing excited states methods, qOMM is less prone to getting stuck in local minima and can achieve convergence with more shallow ansatz circuits.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Intermolecular Proton Transfer Enabled Reactive CO 2 Capture by the Malononitrile Anion

Task-specific ionic liquids (ILs) employing carbanions represent a new class of ILs for carbon capture. The deprotonated malononitrile carbanion, [CH(CN) 2 ] - , has shown close to equimolar capacity for reactive CO 2 capture. Although the formation of the [C(CN) 2 COOH] - carboxylic acid was found to be the final product, how the hydrogen atom on the [CH(CN) 2 ] - carbanion transfers to the carboxylate group as a proton has not been fully understood. In this work, we employ density functional theory calculations with an implicit solvation model to investigate the proton transfer mechanisms in forming carboxylic acid from the reaction of the [CH(CN) 2 ] - carbanion with CO 2 . We find that the intramolecular proton-transfer pathway in [CH(CN) 2 COO] - to form [C(CN) 2 COOH] - is unlikely due to the high energy barrier of 152 kJ/mol. Instead, the intermolecular proton transfer pathway between two [CH(CN) 2 COO] - anions is more feasible to form two molecules of [C(CN) 2 COOH] - , with a significantly lower activation energy of 50 kJ/mol. Moreover, the [C(CN) 2 COOH] - dimer is further stabilized by the intermolecular hydrogen bonds of the two –COOH groups in the Z-configuration of the π-conjugated planar geometry. This insight of reactive CO 2 capture enabled by intermolecular proton transfer will be useful in designing novel carbanions and ILs for carbon capture and conversion.

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

Risk-averse optimal control of semilinear elliptic PDEs

In this paper, we consider the optimal control of semilinear elliptic PDEs with random inputs. These problems are often nonconvex, infinite-dimensional stochastic optimization problems for which we employ risk measures to quantify the implicit uncertainty in the objective function. In contrast to previous works in uncertainty quantification and stochastic optimization, we provide a rigorous mathematical analysis demonstrating higher solution regularity (in stochastic state space), continuity and differentiability of the control-to-state map, and existence, regularity and continuity properties of the control-to-adjoint map. Our proofs make use of existing techniques from PDE-constrained optimization as well as concepts from the theory of measurable multifunctions. We illustrate our theoretical results with two numerical examples motivated by the optimal doping of semiconductor devices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Iterative derivation of effective potentials to sample the conformational space of proteins at atomistic scale

The current capacity of computers makes it possible to perform simulations of small systems with portable, explicit-solvent potentials achieving high degree of accuracy. However, simplified models must be employed to exploit the behavior of large systems or to perform systematic scans of smaller systems. While powerful algorithms are available to facilitate the sampling of the conformational space, successful applications of such models are hindered by the availability of simple enough potentials able to satisfactorily reproduce known properties of the system. We develop an interatomic potential to account for a number of properties of proteins in a computationally economic way. The potential is defined within an all-atom, implicit solvent model by contact functions between the different atom types. The associated numerical values can be optimized by an iterative Monte Carlo scheme on any available experimental data, provided that they are expressible as thermal averages of some conformational properties. We test this model on three different proteins, for which we also perform a scan of all possible point mutations with explicit conformational sampling. The resulting models, optimized solely on a subset of native distances, not only reproduce the native conformations within a few Angstroms from the experimental ones, but show the cooperative transition between native and denatured state and correctly predict the measured free-energy changes associated with point mutations. Moreover, differently from other structure-based models, our method leaves a residual degree of frustration, which is known to be present in protein molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1975

Starting vortex, separation bubbles and stall - A numerical study of laminar unsteady flow around an airfoil

The stalling characteristics of an airfoil in a laminar viscous incompressible fluid are investigated. The governing equations in terms of the vorticity and stream function are solved using an implicit finite-difference scheme and point successive relaxation procedure. The development of the impulsively started flow, the initial generation of circulation, and the behavior of the forces at large times are studied. Following the impulsive start, the lift is at first very large and then it rapidly drops. The subsequent growth of circulation and lift is associated with the starting vortex. After incipient separation, the lift increases owing to enlargement of the separation bubble and intensification of the flow rotation in it. The extension of this bubble beyond the trailing edge causes it to rupture and brings about the stalling characteristics of the airfoil. Subsequently, new bubbles are formed near the upper surface of the airfoil and are swept away. The behavior of the lift acting on the airfoil is explained in terms of the strength and sense of these bubbles.

Mehta, U. B.↗
NASA NTRS1996

A New Procedure for Simulating Unsteady Flows Through Turbomachinery Blade Passages

The development of two new unsteady wake-blade row aerodynamic interaction models and of a rotor-stator unsteady aerodynamic interaction model are outlined. The solutions of Adamczyk's average-passage flow model were used. The responses to the potential disturbances through a blade row were calculated using the MSUTC code. This code can run with and without the use of wall functions. The solver is an implicit finite volume method with flux Jacobians which are evaluated by the flux-vector splitting and the residual fluxes by the Roe's flux-difference splitting.

Chen, Jen Ping↗
NASA NTRS1983

A direct method for the solution of unsteady two-dimensional incompressible Navier-Stokes equations

The unsteady incompressible Navier-Stokes equations are formulated in terms of vorticity and stream function in generalized curvilinear orthogonal coordinates to facilitiate analysis of flow configurations with general geometries. The numerical method developed solves the conservative form of the transport equation using the alternating-direction implicit method, whereas the stream-function equation is solved by direct block Gaussian elimination. The method is applied to a model problem of flow over a back-step in a doubly infinite channel, using clustered conformal coordinates. One-dimensional stretching functions, dependent on the Reynolds number and the asymptotic behavior of the flow, are used to provide suitable grid distribution in the separation and reattachment regions, as well as in the inflow and outflow regions. The optimum grid distribution selected attempts to honor the multiple length scales of the separated-flow model problem. The asymptotic behavior of the finite-differenced transport equation near infinity is examined and the numerical method is carefully developed so as to lead to spatially second-order accurate wiggle-free solutions, i.e., with minimum dispersive error. Results have been obtained in the entire laminar range for the backstep channel and are in good agreement with the available experimental data for this flow problem.

Ghia, K. N.↗
DOE OSTI2023

A Level Set Model for Structured Mesh Representation of Fibrous Preforms used in Composites Manufacturing

Here, a model for implicit representation of fibrous and woven preforms used for composites manufacturing is described. The method is based on a level set function defined on a structured mesh to implicitly capture the complex fiber and weave geometries. Since most software packages used to model woven preforms are based on unstructured mesh, a numerical model for transforming a discretized surface to a structured mesh representation and vice-a-versa is presented. Specific emphasis is on scalable computation of a level set function for large and complex 3D woven structures from a triangulated surface of the smallest unit (yarn). The model is applied to capture a periodic 3D layered 5-harness satin (5HS) weave geometry on a structured mesh. Geometric transformations needed to efficiently compute the level set function for the 3D structure on a high resolution mesh are described. Effects of mesh resolution on the geometric features such as number of resolved fibers, effective fiber size, and the surface area to volume ratio of the iso-surface are analyzed. Using the distance property of the level set function, an approach to evaluate the structure function resulting from transient evolving topology representing densification due to matrix phase depositions is also presented. The approach captures transition of the geometric characteristic from fiber governed at early stages of densification to yarn governed at the terminal stage.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Formal optimization of hovering performance using free wake lifting surface theory

Free wake techniques for performance prediction and optimization of hovering rotor are discussed. The influence functions due to vortex ring, vortex cylinder, and source or vortex sheets are presented. The vortex core sizes of rotor wake vortices are calculated and their importance is discussed. Lifting body theory for finite thickness body is developed for pressure calculation, and hence performance prediction of hovering rotors. Numerical optimization technique based on free wake lifting line theory is presented and discussed. It is demonstrated that formal optimization can be used with the implicit and nonlinear objective or cost function such as the performance of hovering rotors as used in this report.

Chung, S. Y.↗
DOE OSTI2024

Open Circuit Potential of a Au Catalyst during Selective Oxidation of Glycerol

It has been of great interest to understand relationships between the potential of a metal catalyst and its thermochemical catalytic activity, especially in an aqueous phase environment. In the literature, there are correlations of open circuit potential with reaction rates or surface concentrations of reaction intermediates. In this study, we measured the open circuit potential (OCP) of a Au gauze during selective oxidation of glycerol to glyceric acid and hydrogen peroxide in a basic solution, as a function of hydroxide and alcoholate concentrations. It is found that applying a potential to the Au catalyst has no influence on the reaction rate. Although a rough correlation appears to exist between OCP and reaction rate, the data are better fit to an equation which assumes that the potential of the metal (i.e. the OCP) is in equilibrium with the electrochemical potential of the solution, which is defined by the thermodynamic activities of the oxidizing and reducing species. The equation: OCP = constant + Σ ox [(RT/FZ ox )ln[Ox]] - Σ red [(RT/FZ red )ln[Red]] implicitly assumes that the Au metal functions as a probing electrode. It is further found that this equation also applies to the Au-catalyzed H 2 O 2 oxidation/decomposition reaction in a basic medium, and possibly to formic acid dehydrogenation. Here, we postulate that the apparent correlation between OCP and reaction rate is due to the fact that the reaction rate is proportional to the concentrations of the reaction products and/or reactants that define the electrochemical potential of the solution.

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

Regulation of Viable and Optimal Cohorts

This study deals with the evolution of (scalar) attributes (resources or income in evolutionary demography or economics, position in traffic management, etc.) of a population of “mobiles” (economic agents, vehicles, etc.). The set of mobiles sharing the same attributes is regarded as an instantaneous cohort described by the number of its elements. The union of instantaneous cohorts during a mobile window between two attributes is a cohort. Given a measure defining the number of instantaneous cohorts, the accumulation of the mobile attributes on a evolving mobile window is the measure of the cohort on this temporal mobile window. Imposing accumulation constraints and departure conditions, this study is devoted to the regulation of the evolutions of the attributes which are1.viable in the sense that the accumulations constraints are satisfied at each instant;2.and, among them, optimal, in the sense that both the duration of the temporal mobile window is maximum and that the accumulation on this temporal mobile window is the largest viable one. This value is the “accumulation valuation” function. Viable and optimal evolutions under accumulation constraints are regulated by an “implicit Volterra integro-differential inclusion” built from the accumulation valuation function, solution to an Hamilton–Jacobi–Bellman partial differential equation under constraints which is constructed for this purpose.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1988

High-order low-pass implicit tangent filters for use in finite area calculations

High-order implicit tangent filters are developed. The implicit tangent filters possess a highly selective amplitude response function and they can be applied relatively close to a boundary. Comparisons are made between the implicit and the traditional or explicit filters. Numerical simulations are performed to test the response in a limited-area model. An algorithm to maximize computing efficiency is presented. All tests indicate the desirability and utility of the implicit tangent filter.

Raymond, William H.↗
DOE OSTI2015

Fully-Implicit Orthogonal Reconstructed Discontinuous Galerkin for Fluid Dynamics with Phase Change

A new reconstructed Discontinuous Galerkin (rDG) method, based on orthogonal basis/test functions, is developed for fluid flows on unstructured meshes. Orthogonality of basis functions is essential for enabling robust and efficient fully-implicit Newton-Krylov based time integration. The method is designed for generic partial differential equations, including transient, hyperbolic, parabolic or elliptic operators, which are attributed to many multiphysics problems. We demonstrate the method’s capabilities for solving compressible fluid-solid systems (in the low Mach number limit), with phase change (melting/solidification), as motivated by applications in Additive Manufacturing (AM). We focus on the method’s accuracy (in both space and time), as well as robustness and solvability of the system of linear equations involved in the linearization steps of Newton-based methods. The performance of the developed method is investigated for highly-stiff problems with melting/solidification, emphasizing the advantages from tight coupling of mass, momentum and energy conservation equations, as well as orthogonality of basis functions, which leads to better conditioning of the underlying (approximate) Jacobian matrices, and rapid convergence of the Krylov-based linear solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Self-consistent field theory simulations of polymers on arbitrary domains

We introduce a framework for simulating the mesoscale self-assembly of block copolymers in arbitrary confined geometries subject to Neumann boundary conditions. We employ a hybrid finite difference/volume approach to discretize the mean-field equations on an irregular domain represented implicitly by a level-set function. The numerical treatment of the Neumann boundary conditions is sharp, i.e. it avoids an artificial smearing in the irregular domain boundary. This strategy enables the study of self-assembly in confined domains and enables the computation of physically meaningful quantities at the domain interface. In addition, we employ adaptive grids encoded with Quad-/Oc-trees in parallel to automatically refine the grid where the statistical fields vary rapidly as well as at the boundary of the confined domain. This approach results in a significant reduction in the number of degrees of freedom and makes the simulations in arbitrary domains using effective boundary conditions computationally efficient in terms of both speed and memory requirement. Finally, in the case of regular periodic domains, where pseudo-spectral approaches are superior to finite differences in terms of CPU time and accuracy, we use the adaptive strategy to store chain propagators, reducing the memory footprint without loss of accuracy in computed physical observables.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Machine Learning-Driven Solvent Screening for Biobased 2,3-Butanediol Extraction

Biobased 2,3-butanediol (2,3-BDO) is a valuable biomass-derived chemical due to its versatility in being transformed into a wide variety of products. However, the separation and purification of 2,3-BDO from fermentation broth remain a significant challenge owing to its high boiling point and hydrophilic nature. Herein, we developed a machine learning (ML)-based screening workflow that uses molecular calculations as training data and requires only a small number of experimental measurements for validation to identify alternative solvent candidates for the liquid–liquid extraction (LLE) of 2,3-BDO from aqueous solution. In particular, 130 density functional theory (DFT) calculations with the implicit solvation method not only built a correlation between the computational partition coefficient and the experimental distribution coefficient of 2,3-BDO but also parameterized an Extra-Trees ML model to screen the distribution coefficient for a wider range of 6717 organic solvents. The experimental measurements of only 24 solvents were needed to validate the computational results. A list of 50 prioritized solvents was proposed for 2,3-BDO LLE, and seven additional experimental measurements were conducted to further verify our selected solvents. The impact of the extraction temperature and solvent-to-feed ratio was also investigated for selected solvents in experiments. Furthermore, this work suggested alternative solvents for 2,3-BDO LLE and proposed a versatile workflow that requires fewer experiments and can be applied to a broader range of LLE studies.

Extraction↗
DOE OSTI2021

Directing the far-from-equilibrium assembly of nanoparticles in confined liquid crystals by hydrodynamic fields

The assembly of nematic colloids relies on long-range elastic interactions that can be manipulated through external stimuli. Confinement and the presence of a hydrodynamic field alter the defect structures and the energetic interactions between the particles. Here, the assembly landscape of nanoparticles embedded in a nematic liquid crystal confined in a nanochannel under a pressure-driven flow is determined. The dynamics of the liquid crystal tensor alignment field is determined through a Poisson-Bracket framework, namely the Stark–Lubensky equations, coupled with the zero-Reynolds momentum equations and the liquid crystal Landau-de Gennes free energy functional. A second order semi-implicit time integration and a three-dimensional Galerkin finite element method are used to resolve flow and nematic fields under several conditions. In general, the zero Reynolds flow displaces the defects around the particles in the upstream direction and renders the surface anchoring ineffective when the flow strength dominates over the nematic elasticity. More importantly, the potential of mean force for particle assembly is non-monotonic independent of surface anchoring. Our results show that the confinement length scale determines the repulsion/attraction transition between colloids, while the flow strength modifies the static defect structure surrounding the particles and determines the magnitude of the energetic barrier for successful assembly. In the attractive regime, the particles move at different rates through the nematic until one particle eventually catches up with the other. This process occurs against or along the direction of flow depending on the flow strength. Ultimately, these results provide a template for engineering and controlling the transport and assembly of nanoparticles under far-from equilibrium conditions in anisotropic media.

36 MATERIALS SCIENCE↗