Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient calculation”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Development of new flux splitting schemes

Maximizing both accuracy and efficiency has been the primary objective in designing a numerical algorithm for computational fluid dynamics (CFD). This is especially important for solutions of complex three dimensional systems of Navier-Stokes equations which often include turbulence modeling and chemistry effects. Recently, upwind schemes have been well received for their capability in resolving discontinuities. With this in mind, presented are two new flux splitting techniques for upwind differencing. The first method is based on High-Order Polynomial Expansions (HOPE) of the mass flux vector. The second new flux splitting is based on the Advection Upwind Splitting Method (AUSM). The calculation of the hypersonic conical flow demonstrates the accuracy of the splitting in resolving the flow in the presence of strong gradients. A second series of tests involving the two dimensional inviscid flow over a NACA 0012 airfoil demonstrates the ability of the AUSM to resolve the shock discontinuity at transonic speed. A third case calculates a series of supersonic flows over a circular cylinder. Finally, the fourth case deals with tests of a two dimensional shock wave/boundary layer interaction.

Liou, Meng-Sing↗
NASA NTRS1999

What is the Best Alternative to Diagonalization of the Hamiltonian in Large Scale Semiempirical Calculations?

Recently, several linear scaling approaches have been introduced which replace the time dominating diagonalization step in semiempirical methods, enabling practical calculations to be performed, on very large molecules. This paper compares the accuracy and performance of pseudodiagonalization (PD), conjugate gradient density matrix search (CG-DMS), the Chebyshev polynomial expansion method (CEM), and purification of the density matrix (PDM) as linear scaling substitutions for diagonalization. The scaling, speed, and reliability of these methods are compared for AMI single point energy calculations on polyglycine chains (up to 20,000 atoms), water clusters (up to 12,300 atoms), and nucleic acids (up to 6300 atoms).

Daniels, Andrew D.↗
DOE OSTI2015

Chiral drag force

Here, we provide a holographic evaluation of novel contributions to the drag force acting on a heavy quark moving through strongly interacting plasma. The new contributions are chiral in the sense that they act in opposite directions in plasmas containing an excess of left- or right-handed quarks. The new contributions are proportional to the coefficient of the axial anomaly, and in this sense also are chiral. These new contributions to the drag force act either parallel to or antiparallel to an external magnetic field or to the vorticity of the fluid plasma. In all these respects, these contributions to the drag force felt by a heavy quark are analogous to the chiral magnetic effect (CME) on light quarks. However, the new contribution to the drag force is independent of the electric charge of the heavy quark and is the same for heavy quarks and antiquarks, meaning that these novel effects do not in fact contribute to the CME current. We show that although the chiral drag force can be non-vanishing for heavy quarks that are at rest in the local fluid rest frame, it does vanish for heavy quarks that are at rest in a suitably chosen frame. In this frame, the heavy quark at rest sees counterpropagating momentum and charge currents, both proportional to the axial anomaly coefficient, but feels no drag force. This provides strong concrete evidence for the absence of dissipation in chiral transport, something that has been predicted previously via consideration of symmetries. Along the way to our principal results, we provide a general calculation of the corrections to the drag force due to the presence of gradients in the flowing fluid in the presence of a nonzero chemical potential. We close with a consequence of our result that is at least in principle observable in heavy ion collisions, namely an anticorrelation between the direction of the CME current for light quarks in a given event and the direction of the kick given to the momentum of all the heavy quarks and antiquarks in that event.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Electronic structure and thermoelectric properties of CoAsSb with post-DFT approaches

Here, we study the electronic structure and thermoelectric properties of recently synthesized CoAsSb. The calculated bandgap becomes more accurate for increasingly-complex electronic structure methods: generalized gradient approximation, hybrid functionals, selfconsistent linearized quasiparticle GW method (LQSGW), and LQSGW with simplified vertex corrections. By equating the bandgaps of each method from a rigid shift of the bands, we evaluate the contributions made to thermoelectric properties beyond the bandgap. In doing so, we evaluate the efficacy of a commonpractice: a rigid shift applied to less-costly electronic structure methods in order to gain some accuracy of the more-costly methods. We find that while the shift made the Seebeck coefficients much closer to one another than from the original electronic structures, there remain differences in the Goldsmid-Sharp (thermoelectric) bandgap between the methods and from the intended electronic bandgap. Additionally, some lasting differences in temperature dependence remain between the methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Effect of Naturally Deposited Film on the Sub-rapid Solidification of Medium Manganese Steel by Using Droplet Solidification Technique

In this study, the effect of naturally deposited film on the interfacial heat transfer during the sub-rapid solidification of medium manganese steel was studied by means of a droplet solidification technique. The heat transfer rates were calculated by the Inverse Heat Conduction Program (IHCP), through computing the responding temperatures’ gradient inside the cooling mold. The results showed that the contact angle between the mold and solidified droplet becomes lower (from 117.539 to 82.217 deg) with the repeat of the dropping tests, which means the wettability improves after the deposition of the natural film. Meanwhile the liquid time (the time from the start of ejection to the initiation of solidification) decreases from 1.76 to 1.28 seconds during the tests, suggesting that the steel keeps in the liquid state shorter. Therefore, the initiation of solidification is advanced, such that there would be less time for matching the scale of the shell and the mold. Consequently, the maximum heat flux increases from 2.90 to 6.35 MW/m{sup 2} with the deposition of the natural film, as the wettability is getting better leading to the reduction of the interfacial thermal resistance. However, the final maximum heat flux was decreasing to 4.98 MW/m{sup 2} due to the increased thickness of the film and corresponding interfacial thermal resistance. Moreover, it was found that the secondary dendrite arm spacing decreased from 10.84 to 5.01 μm at first, and then it increased to 8.11 μm with the variations of heat flux and cooling rate.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Film and membrane-model thermodynamics of free thin liquid films

In spite of over 7 decades of effort, the thermodynamics of thin free liquid films (as in emulsions and foams) lacks clarity. Following a brief review of the meaning and measurement of thin-film forces (i.e., conjoining/disjoining pressures), we offer a consistent analysis of thin-film thermodynamics. By carefully defining film reversible work, two distinct thermodynamic formalisms emerge: a film model with two zero-volume membranes each of film tension γ ƒ and a membrane model with a single zero-volume membrane of membrane tension 2γ m . In both models, detailed thermodynamic analysis gives rise to thin-film Gibbs adsorption equations that allow calculation of film and membrane tensions from measurements of disjoining-pressure isotherms. A modified Young–Laplace equation arises in the film model to calculate film-thickness profiles from the film center to the surrounding bulk meniscus. No corresponding relation exists in the membrane model. Illustrative calculations of disjoining-pressure isotherms for water are presented using square-gradient theory. In conclusion, we report considerable deviations from Hamaker theory for films less than about 3 nm in thickness. Such thin films are considerably more attractive than in classical Hamaker theory. Available molecular simulations reinforce this finding.

Conjoining/disjoining pressure↗
DOE OSTI2021

Analytic evaluation of energy first derivatives for spin–orbit coupled-cluster singles and doubles augmented with noniterative triples method: General formulation and an implementation for first-order properties

The paper reports a formulation of analytic energy first derivatives for the coupled-cluster singles and doubles augmented with noniterative triples [CCSD(T)] method with spin–orbit coupling included at the orbital level and an implementation for evaluation of first-order properties. Here, the new development extends the applicability of relativistic coupled cluster methods for heavy-element chemistry and spectroscopy. Benchmark calculations of first-order properties including dipole moments and electric-field gradients using the corresponding exact two-component property integrals are presented for heavy-element containing molecules to demonstrate the usefulness of the present analytic scheme.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Analytical gradients for molecular-orbital-based machine learning

We report molecular-orbital-based machine learning (MOB-ML) enables the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. Here, we present the derivation, implementation, and numerical demonstration of MOB-ML analytical nuclear gradients, which are formulated in a general Lagrangian framework to enforce orthogonality, localization, and Brillouin constraints on the molecular orbitals. The MOB-ML gradient framework is general with respect to the regression technique (e.g., Gaussian process regression or neural networks) and the MOB feature design. We show that MOB-ML gradients are highly accurate compared to other ML methods on the ISO17 dataset while only being trained on energies for hundreds of molecules compared to energies and gradients for hundreds of thousands of molecules for the other ML methods. The MOB-ML gradients are also shown to yield accurate optimized structures at a computational cost for the gradient evaluation that is comparable to a density-corrected density functional theory calculation.

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

Instabilities in mimetic matter perturbations

We study cosmological perturbations in mimetic matter scenario with a general higher derivative function. We calculate the quadratic action and show that both the kinetic term and the gradient term have the wrong sings. We perform the analysis in both comoving and Newtonian gauges and confirm that the Hamiltonians and the associated instabilities are consistent with each other in both gauges. The existence of instabilities is independent of the specific form of higher derivative function which generates gradients for mimetic field perturbations. It is verified that the ghost instability in mimetic perturbations is not associated with the higher derivative instabilities such as the Ostrogradsky ghost.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Flexoinduced ferroelectricity in low-dimensional transition metal dichalcogenides

We developed a Landau-type theory for the description of polar phenomena in low-dimensional transition metal dichalcogenides (TMDs), specifically exploring the flexoelectric origin of the polarization induced by a spontaneous bending and by inversion symmetry breaking due to the interactions with substrate. We consider the appearance of the spontaneous out-of-plane polarization due to the flexoelectric coupling with the strain gradient of the spontaneous surface rippling and surface-induced piezoelectricity. Additionally, performed calculations proved that the out-of-plane spontaneous polarization, originating from the flexoelectric effect in a rippled TMD, is bistable and reversible by a nonuniform electric field. In contrast, the spontaneous polarization induced by a misfit strain and symmetry-sensitive surface-induced piezoelectric coupling, cannot be reversed by an external electric field. Special attention is paid to the spectral analysis of the linear dielectric susceptibility and gain factor, which enhancement is critically important for the observation of the polar phenomena in low-dimensional TMDs by the surface-enhanced vibrational spectroscopy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
DOE OSTI2017

Enhancement of the Accelerating Gradient in Superconducting Microwave Resonators

The accelerating gradient of superconducting resonators can be enhanced by engineering the thickness of a dirty layer grown at the cavity's rf surface. In this paper the description of the physics behind the accelerating gradient enhancement by meaning of the dirty layer is carried out by solving numerically the the Ginzburg-Landau (GL) equations for the layered system. The calculation shows that the presence of the dirty layer stabilizes the Meissner state up to the lower critical field of the bulk, increasing the maximum accelerating gradient.

43 PARTICLE ACCELERATORS↗
NASA NTRS1963

Application of Dorodnitsyn's Integral Method to Nonequilibrium Flows over Pointed Bodies

Dorodnitsyn's integral method is used to obtain an approximate solution to the supersonic nonequilibrium flow over pointed bodies with attached shock waves. The partial differential equations governing the flow are converted to an approximate set of ordinary equations, which are solved by numerical integration starting at the body tip. Detailed analytical and numerical results for the first approximation are presented, considering the vibrational relaxation of a diatomic gas over a wedge or cone. It is shown that the first approximation yields: (1) The exact flow-variable gradients at the wedge tip (2) Expressions for the flow-variable gradients at the cone tip which are in agreement with extrapolations of characteristics calculations (3) A good approximate algebraic solution for frozen or equilibrium conical flow (4) An approximate expression for the nonequilibrium-flow stream function which affords a means of obtaining variations across the shock layer of the temperature and vibrational energy. Numerical results for both the wedge and cone compare favorably with identical cases computed by the method of characteristics.

NONEQUILIBRIUM FLOW↗
NASA NTRS1973

CONMIN: A FORTRAN program for constrained function minimization: User's manual

CONMIN is a FORTRAN program, in subroutine form, for the solution of linear or nonlinear constrained optimization problems. The basic optimization algorithm is the Method of Feasible Directions. The user must provide a main calling program and an external routine to evaluate the objective and constraint functions and to provide gradient information. If analytic gradients of the objective or constraint functions are not available, this information is calculated by finite difference. While the program is intended primarily for efficient solution of constrained problems, unconstrained function minimization problems may also be solved, and the conjugate direction method of Fletcher and Reeves is used for this purpose. This manual describes the use of CONMIN and defines all necessary parameters. Sufficient information is provided so that the program can be used without special knowledge of optimization techniques. Sample problems are included to help the user become familiar with CONMIN and to make the program operational.

Vanderplaats, G. N.↗
NASA NTRS1976

Dynamic response of membranes with finite elements

A general set of finite-element equations is developed for treating the nonlinear dynamic response of membranes consisting of complex nonlinear materials. The approach proposed by Key (1974) for a most complete treatment of the finite-element method to the nonlinear dynamic analysis of solids is used in developing the nonlinear equations of motion for membranes with cubic quadrilateral finite elements. The computations are kept general with respect to material modeling by computing the deformation gradient and requiring a constitutive subroutine to supply the Cauchy stress. An example calculation of the vibration of a square membrane demonstrates the accuracy of the formulation for small deformation theory and shows the nonlinear effect of large deformations.

Benzley, S. E.↗
NASA NTRS1976

Gravitational dynamics of biosystems - Some speculations

The response of organisms to gravity is generally discussed in terms of hypotheses involving sedimentation and other static effects. This paper considers several complex, inhomogeneous fluid-containing systems that are intended to model some possible dynamic effects of gravity on biosystems. It is shown that the presence of gravity may result in modified long range transport, concentration oscillations, and broken symmetries. The magnitude of density-gradient-driven convective transport times, and their ratios to diffusive transport times, are calculated for cell dimensions of six different plant varieties. The results indicate that further investigation of gravitational convection effects may be realistic in some cases and is definitely not in others. The results of this paper should aid in the planning of 'zero-gravity' experiments concerning plant geotropism and bio-materials processing.

Kessler, J. O.↗
NASA NTRS1976

Similarity solution of the boundary-layer equations for laser heated flows

The laminar boundary-layer equations with local similarity approximation are solved in order to estimate the heat transfer rate to the wall in a laser-heated rocket thruster using pure hydrogen, where the average temperature of the hot plasma core is about 14,000 K and the core Reynolds number based on nozzle throat diameter is about 2000. Under these conditions, the density-viscosity product at the wall can be 10 times the freestream value, and the hydrogen is completely dissociated. Hence the boundary layer equations with variable transport properties are solved by a quasi-linearization technique, and the equilibrium properties of hydrogen are used in the calculations. Velocity profiles are obtained, and the wall shear and wall stagnation enthalpy gradient are plotted against the pressure gradient parameter.

Wu, P. K.↗
NASA NTRS1977

Some design considerations for large space structures

Physical characteristics of large skeletal frameworks for space applications are investigated by analyzing one concept: the tetrahedral truss, which is idealized as a sandwich plate with isotropic faces. Appropriate analytical relations are presented in terms of the truss column-element properties, which for calculations were taken as slender graphite/epoxy tubes. Column loads, resulting from gravity-gradient control and orbital transfer, are found to be small for the class structure investigated. Fundamental frequencies of large truss structures are shown to be an order of magnitude lower than large earth-based structures. Permissible loads are shown to result in small lateral deflections of the truss due to low strain at Euler buckling of the slender graphite/epoxy truss column elements. Lateral thermal deflections are found to be a fraction of the truss depth using graphite/epoxy columns.

Bush, H. G.↗