Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytic gradients”

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 775 records · Page 43

On the Electron Diffusion Region in Asymmetric Reconnection with a Guide Magnetic Field

Particle-in-cell simulations in a 2.5-D geometry and analytical theory are employed to study the electron diffusion region in asymmetric reconnection with a guide magnetic field. The analysis presented here demonstrates that similar to the case without guide field, in-plane flow stagnation and null of the in-plane magnetic field are well separated. In addition, it is shown that the electric field at the local magnetic X point is again dominated by inertial effects, whereas it remains dominated by nongyrotropic pressure effects at the in-plane flow stagnation point. A comparison between local electron Larmor radii and the magnetic gradient scale lengths predicts that distribution should become nongyrotropic in a region enveloping both field reversal and flow stagnation points. This prediction is verified by an analysis of modeled electron distributions, which show clear evidence of mixing in the critical region.

reconnection↗

Micrometeorite dynamic pyrometamorphism: Observation of a thermal gradient in iron-nickel sulfide

Dynamic pyrometamorphism describes the mineralogical changes in a micrometeorite that occur in response to flash-heating during atmospheric entry. Mineral reactions during this event will be endothermic and act as local sinks for thermal energy which delays a uniform temperature distribution in decelerating micrometeorites. The most common pyrometamorphic reaction is formation of magnetite nanocrystal decorations on Mg,Fe-silicate grains and magnetite or (partial) maghemite rims. Constituent silicates also respond to this thermal event and show formation of olivine-maghemite symplectites and growth of garnet and partial laihunitisation. In continued effort to constrain dynamic pyrometamorphism of unmelted interplanetary dust particles (IDP's), I determined the mineralogical composition, and Fe,Ni-sulfide chemistry, in the sulfide-rich particles L2005C39, L2005E40, and L2006A28 using a JEOL 2000FX analytical electron microscope equipped with a TN5500 energy dispersive spectrometer for in situ microanalysis. Sulfide compositions in these three IDP's are distinctly different from those in silicate-rich particle L2005T13 which might suggest a delicate balance of the fO2/fS2 ratios during dynamic pyrometamorphism.

Rietmeijer, Frans J. M.↗

Accumulation of Antifreeze Proteins on Ice Is Determined by Adsorption

Antifreeze proteins (AFPs) enable the survival of organisms in freezing environments by binding to ice crystals and suppressing their growth. The rate at which AFPs accumulate on ice is determined by an interplay between their diffusion from the aqueous solution and their adsorption to the ice crystal. To interrogate the relative importance of these two processes, we combined fluorescence experiments with reaction–diffusion modeling. We found that as AFPs accumulate on ice, their concentration in the aqueous solution does not develop a gradient but remains equal to its bulk concentration throughout our experiments. Furthermore, these findings led us to conclude that AFP accumulation on ice crystals, which are smaller than 100 μm in radius, is not limited by the diffusion of AFPs, but by the kinetics of AFP adsorption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Harnessing deep reinforcement learning to construct time-dependent optimal fields for quantum control dynamics

Here, we present an efficient deep reinforcement learning (DRL) approach to automatically construct time-dependent optimal control fields that enable desired transitions in dynamical chemical systems. Our DRL approach gives impressive performance in constructing optimal control fields, even for cases that are difficult to converge with existing gradient-based approaches. We provide a detailed description of the algorithms and hyperparameters as well as performance metrics for our DRL-based approach. Our results demonstrate that DRL can be employed as an effective artificial intelligence approach to efficiently and autonomously design control fields in quantum dynamical chemical systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

PyFLOSIC: Python-based Fermi–Löwdin orbital self-interaction correction

We present pyflosic, an open-source, general-purpose python implementation of the Fermi–Löwdin orbital self-interaction correction (FLO-SIC), which is based on the python simulation of chemistry framework (pyscf) electronic structure and quantum chemistry code. Thanks to pyscf, pyflosic can be used with any kind of Gaussian-type basis set, various kinds of radial and angular quadrature grids, and all exchange-correlation functionals within the local density approximation, generalized-gradient approximation (GGA), and meta-GGA provided in the libxc and xcfun libraries. A central aspect of FLO-SIC is the Fermi-orbital descriptors, which are used to estimate the self-interaction correction. Importantly, they can be initialized automatically within pyflosic; they can also be optimized within pyflosic with an interface to the atomic simulation environment, a python library that provides a variety of powerful gradient-based algorithms for geometry optimization. Although pyflosic has already facilitated applications of FLO-SIC to chemical studies, it offers an excellent starting point for further developments in FLO-SIC approaches, thanks to its use of a high-level programming language and pronounced modularity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Permutation invariant polynomial neural network based diabatic ansatz for the ( E + A ) × ( e + a ) Jahn–Teller and Pseudo-Jahn–Teller systems

In this work, the permutation invariant polynomial neural network (PIP-NN) approach is employed to construct a quasi-diabatic Hamiltonian for system with non-Abelian symmetries. It provides a flexible and compact NN-based diabatic ansatz from the related approach of Williams, Eisfeld, and co-workers. The example of $H^3_+$ is studied, which is an (E + A) × (e + a) Jahn–Teller and Pseudo-Jahn–Teller system. The PIP-NN diabatic ansatz is based on the symmetric polynomial expansion of Viel and Eisfeld, the coefficients of which are expressed with neural network functions that take permutation-invariant polynomials as input. This PIP-NN-based diabatic ansatz not only preserves the correct symmetry but also provides functional flexibility to accurately reproduce ab initio electronic structure data, thus resulting in excellent fits. Further, the adiabatic energies, energy gradients, and derivative couplings are well reproduced. A good description of the local topology of the conical intersection seam is also achieved. Therefore, this diabatic ansatz completes the PIP-NN based representation of DPEM with correct symmetries and will enable us to diabatize even more complicated systems with complex symmetries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Non-local large-strain FFT-based formulation and its application to interface-dominated plasticity of nano-metallic laminates

Here, this paper presents a novel formulation and its robust numerical implementation of strain-gradient (SG) crystal plasticity within a large-strain (LS) elasto-viscoplastic (EVP) fast Fourier transform (FFT)-based micromechanical model. The resulting non-local SG-LS-EVPFFT formulation is used to model and understand the process of kink band formation during layer-parallel compression of nano-metallic laminates (NMLs). NMLs are layered composites with nanoscale thicknesses, thus requiring consideration of the interaction between dislocations and interfaces within the micromechanical model. The length-scale parameter of the SG model is calibrated by simulating a double pile-up and comparing predictions to analytical solution. This required new expressions for the defect energy, resulting in more accurate double pile-up predictions. The calibrated SG-LS-EVPFFT model is then used to simulate layer-parallel compression of copper-niobium NML. Formation of kink bands is predicted, and the model is used to rationalize the microscopic mechanisms enabling the formation process. It is found that accumulation of dislocations at interfaces leads to activation of layer-parallel slip, which in turn leads to kink band formation.

36 MATERIALS SCIENCE↗

Analytical and experimental investigation of circulation control by means of a turbulent Coanda jet

An analytical and experimental investigation of circulation control on a circular cylinder by means of tangential blowing (Coanda effect) is presented. The analytical method developed has also been used to estimate the blowing coefficients required for achieving potential flow on airfoils with flaps. The analysis is presented for conditions for which the flow in the boundary layer ahead of the jet exit is turbulent. The turbulent boundary layer and the jet layer on the upper surface, and the turbulent boundary layer on the lower surface are computed by a multi-strip integral method. The region of integration is between the correponding transition and separation points on each surface. Longitudinal curvature effects, which give rise to a radial pressure gradient across the jet layer and to an additional adverse tangential pressure gradient just upstream of the separation point, are included in the jet layer analysis in an approximate manner. The longitudinal curvature effect is found to have a pronounced influence on the separation of the jet layer.

Levinsky, E. S.↗

Study of casing treatment stall margin improvement phenomena

The results of a program of experimental and analytical research in casing treatments over axial compressor rotor blade tips are presented. Circumferential groove, axial-skewed slot, and blade angle slot treatments were tested. These yielded, for reduction in stalling flow and loss in peak efficiency, 5.8% and 0 points, 15.3% and 2.0 points, and 15.0% and 1.2 points, respectively. These values are consistent with other experience. The favorable stalling flow situations correlated well with observations of higher-then-normal surface pressures on the rotor blade pressure surfaces in the tip region, and with increased maximum diffusions on the suction surfaces. Annular wall pressure gradients, especially in the 50-75% chord region, are also increased and blade surface pressure loadings are shifted toward the trailing edge for treated configurations. Rotor blade wakes may be somewhat thinner in the presence of good treatments, particularly under operating conditions close to the baseline stall.

Prince, D. C., Jr.↗

Enriching Cations Using Electric Field Gradients Generated by Bipolar Electrodes in the Absence of Buffer

We discuss that ionic current gradients and corresponding electric field gradients formed at bipolar electrodes are useful for enriching cations up to ~350-fold at specific locations in microfluidic devices. The enrichment mechanism is determined by collecting current data, measuring solution conductivity, and by performing pH indicator experiments. These data reveal the dynamic nature of the enrichment mechanism. Specifically, they show the complex relationship between the ionic current gradients formed at the bipolar electrode and the extent of cation enrichment. Finally, we demonstrate cation enrichment is possible in buffer-free solutions containing up to 50.0 mM KCl. This finding is significant because it broadens the range of solutions in which membrane-free, electrochemical separations using bipolar electrodes can be performed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

One-pot ternary sequential reactions for photopatterned gradient multimaterials

Seamless multimaterial construction is a common motif in animal physiology. Such continuous mechanical gradients remain challenging to reproduce in engineered systems, as current resin chemistries typically result in a single fixed set of properties. As an alternative to single-property materials, we introduce a thiol-ene-epoxy-based photothermal reaction scheme that produces multimaterials by altering the polymer microstructure within a single resin. In this system, the photodosage during the first stage of processing dictates the extent of conversion for each subsequent reaction. As a result, our photochemistry can exhibit a diverse range of soft (Young’s modulus, E ~ 400 kPa; elongation, dL/L0 ~ 300%) and stiff (E ~ 1.6 GPa; dL/L 0 ~ 3%) mechanical properties. Furthermore, we pattern photostable and mechanically robust modulus gradients (d[E r, stiff /E r , soft] /dx > 1,000 mm -1 ) that exceed those found in squid beaks and human knee entheses. We demonstrate the ability to build intricate multimaterial architectures including a soft, wearable braille display.

36 MATERIALS SCIENCE↗

Electronic Structure Calculations of the Fundamental Interactions in Solvent Extraction Desalination

Solvent extraction desalination (SED) is a promising membrane-free desalination technology, but new solvents are still needed that can improve the molecular-level water selectivity, while minimizing solvent contamination in the aqueous phase. Improvements in solvent design for SED require a fundamental understanding of the molecular interactions governing solvation behavior. Herein, quantum chemical calculations are used to delineate different fundamental interactions in water-solvent systems involving 38 different solvents, with a particular focus on the role of the electrostatic potential characteristics. Furthermore, these interactions are rigorously analyzed using the independent gradient model based on Hirshfeld partition (IGMH), Symmetry-Adapted Perturbation Theory (SAPT), as well as the energy decomposition analysis based on force field (EDA-FF) to compare against a standard force field model. The solvation free energies and partition coefficients are also calculated to estimate water extraction performance in a bulk solvent system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analytically differentiable metrics for phase stability

Here, in this work, a long-established but sparsely documented method of obtaining semi-analytic derivatives of thermodynamic properties with respect to equilibrium conditions is briefly reviewed and rigorously derived. This procedure is then leveraged to construct general forms of derivatives of the residual driving force, a metric for measuring phase stability used in CALPHAD model optimization, with respect to overall system and individual phase compositions. Applied examples – calculating heat capacity in the Al-Fe system, thermodynamic factors in the Nb-V-W system, and residual driving force derivatives in the Ni-Ti system – demonstrate the versatility, accuracy, and extensibility of this method. Using the developed method, residual driving force gradients can be applied directly in CALPHAD model optimizers, as well as in materials design frameworks, to identify regions of phase stability with an efficient, gradient-based approach.

36 MATERIALS SCIENCE↗

The role of gas flow distributions on CO 2 mineralization within monolithic cemented composites: coupled CFD-factorial design approach

The carbonation kinetics of monolithic cementing composites are strongly affected by gas transport which is, in turn, influenced by microstructural resistances and the presence of liquid water within pore networks. The non-uniform gas flow distribution within the CO 2 mineralization reactor can impart mass transfer resistance in the monolith microstructure, which affects the uptake of CO 2 (“carbonation”) of the cementing composites. This paper demonstrates how the gas spatial distribution (velocity and flow rate; quantified by CFD analysis) and processing conditions (temperature, relative humidity, and flow rate; quantified by factorial design) affect drying and carbonation, and in turn, the engineering properties of a representative ‘monolithic’ carbonate-cemented concrete component (i.e., herein concrete masonry unit: CMUs, also known as concrete block). It is shown that the gas flow distribution affects drying front penetration and results in moisture and carbonation gradients within the monolith. Particularly, variations in drying kinetics caused by non-uniformity of the contacting gas velocity impose gradients in moisture saturation, which results in increasing microstructural resistance to CO 2 transport. The resultant non-uniform carbonate-mineral formation (i.e., carbonate cementation), if not controlled, can produce gradients in mechanical properties and may alter failure patterns upon loading. Finally, these insights inform the optimal design of gas flow distribution systems and processing conditions within CO 2 mineralization reactors for the manufacturing of low-CO 2 concrete components using CO 2 -dilute industrial flue gas streams.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design of Structurally Efficient Tapered Struts (SETS)

A study was conducted to develop mass efficient composite struts. A closed-form design methodology for composite struts was developed using well established analyses to predict Euler buckling, local wall buckling; compression strength, damage tolerance, and interlaminar shear at geometric gradients. The methodology was coded in a spreadsheet suitable for convenient and rapid sizing of tapered composite struts. This spreadsheet analysis was used to determine the influence of several variables such as material stiffness, strut diameter, and material allowables on strut weight for given loading conditions. The comparison showed that, while the Park Aerospace design method was well suited to preliminary sizing for a conservative design, the closed-form-analyses-based spreadsheet accounts for all possible failure modes and is a good optimum strut design tool. The report concludes with a set of recommendations for future work in analytical design and analysis methodology enhancements.

Ravi Deo↗

Self-consistent theory of the Darrieus–Landau and Rayleigh–Taylor instabilities with self-generated magnetic fields

The Rayleigh–Taylor (RT) and Darrieus–Landau (DL) instabilities are studied in an inertial confinement fusion context within the framework of small critical-to-shell density ratio D R and weak acceleration regime, i.e., large Froude number Fr. The quasi-isobaric analysis in Sanz et al. is completed with the inclusion of non-isobaric and self-generated magnetic-field effects. The analysis is restricted to perturbation wavelengths k~1 larger than the conduction length scale at the ablation front, yet its validity ranges from wavelengths shorter and larger than the conduction layer width (distance between the ablation front and the critical surface). The use of a sharp boundary model leads to a single analytical expression of the dispersion relation encompassing both instabilities. The two new effects come into play by modifying the perturbed mass and momentum fluxes at the ablation front. The momentum flux (perturbed pressure at the spike) is the predominant stabilizing mechanism in the RT instability (overpressure) and the driving mechanism in the DL instability (underpressure). The non-isobaric effects notably modify the scaling laws in the DL limit, leading to an underpressure scaling as ~k –11/15 rather than ~k –2/5 obtained in the quasi-isobaric model. The magnetic fields are generated due to misalignment between pressure and density gradients (Biermann battery effect). They affect the hydrodynamics by bending the heat flux lines. Within the framework of this paper, they enhance ablation, resulting in a stabilizing effect that peaks for perturbation wavelengths comparable to the conduction layer width. Here, the combination of parameters D R Fr 2/3 defines the region of predominance of each instability in the dispersion relation. It is proven that the DL region falls outside of the parameter range in inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three-dimensional inviscid analysis of radial turbine flow and a limited comparison with experimental data

The three-dimensional inviscid DENTON code is used to analyze flow through a radial-inflow turbine rotor. Experimental data from the rotor are compared with analytical results obtained by using the code. The experimental data available for comparison are the radial distributions of circumferentially averaged values of absolute flow angle and total pressure downstream of the rotor exit. The computed rotor-exit flow angles are generally underturned relative to the experimental values, which reflect the boundary-layer separation at the trailing edge and the development of wakes downstream of the rotor. The experimental rotor is designed for a higher-than-optimum work factor of 1.126 resulting in a nonoptimum positive incidence and causing a region of rapid flow adjustment and large velocity gradients. For this experimental rotor, the computed radial distribution of rotor-exit to turbine-inlet total pressure ratios are underpredicted due to the errors in the finite-difference approximations in the regions of rapid flow adjustment, and due to using the relatively coarser grids in the middle of the blade region where the flow passage is highly three-dimensional. Additional results obtained from the three-dimensional inviscid computation are also presented, but without comparison due to the lack of experimental data. These include quasi-secondary velocity vectors on cross-channel surfaces, velocity components on the meridional and blade-to-blade surfaces, and blade surface loading diagrams. Computed results show the evolution of a passage vortex and large streamline deviations from the computational streamwise grid lines. Experience gained from applying the code to a radial turbine geometry is also discussed.

Choo, Y. K.↗

Three-dimensional inviscid analysis of radial-turbine flow and a limited comparison with experimental data

The three-dimensional inviscid DENTON code is used to analyze flow through a radial-inflow turbine rotor. Experimental data from the rotor are compared with analytical results obtained by using the code. The experimental data available for comparison are the radial distributions of circumferentially averaged values of absolute flow angle and total pressure downstream of the rotor exit. The computed rotor-exit flow angles are generally underturned relative to the experimental values, which reflect the boundary-layer separation at the trailing edge and the development of wakes downstream of the rotor. The experimental rotor is designed for a higher-than-optimum work factor of 1.126 resulting in a nonoptimum positive incidence and causing a region of rapid flow adjustment and large velocity gradients. For this experimental rotor, the computed radial distribution of rotor-exit to turbine-inlet total pressure ratios are underpredicted due to the errors in the finite-difference approximations in the regions of rapid flow adjustment, and due to using the relatively coarser grids in the middle of the blade region where the flow passage is highly three-dimensional. Additional results obtained from the three-dimensional inviscid computation are also presented, but without comparison due to the lack of experimental data. These include quasi-secondary velocity vectors on cross-channel surfaces, velocity components on the meridional and blade-to-blade surfaces, and blade surface loading diagrams. Computed results show the evolution of a passage vortex and large streamline deviations from the computational streamwise grid lines. Experience gained from applying the code to a radial turbine geometry is also discussed.

Choo, Y. K.↗