Method for calculation of compressible laminar boundary layer with axial pressure gradient and heat transfer
Explore the source record for details and available documents.
Libby, Paul A↗
SEARCH · Engineering Papers
Results for “gradient methods”
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.
An autoencoder compression approach for accelerating large-scale inverse problems
Abstract Partial differential equation (PDE)-constrained inverse problems are some of the most challenging and computationally demanding problems in computational science today. Fine meshes required to accurately compute the PDE solution introduce an enormous number of parameters and require large-scale computing resources such as more processors and more memory to solve such systems in a reasonable time. For inverse problems constrained by time-dependent PDEs, the adjoint method often employed to compute gradients and higher order derivatives efficiently requires solving a time-reversed, so-called adjoint PDE that depends on the forward PDE solution at each timestep. This necessitates the storage of a high-dimensional forward solution vector at every timestep. Such a procedure quickly exhausts the available memory resources. Several approaches that trade additional computation for reduced memory footprint have been proposed to mitigate the memory bottleneck, including checkpointing and compression strategies. In this work, we propose a close-to-ideal scalable compression approach using autoencoders to eliminate the need for checkpointing and substantial memory storage, thereby reducing the time-to-solution and memory requirements. We compare our approach with checkpointing and an off-the-shelf compression approach on an earth-scale ill-posed seismic inverse problem. The results verify the expected close-to-ideal speedup for the gradient and Hessian-vector product using the proposed autoencoder compression approach. To highlight the usefulness of the proposed approach, we combine the autoencoder compression with the data-informed active subspace (DIAS) prior showing how the DIAS method can be affordably extended to large-scale problems without the need for checkpointing and large memory.
Mathematics↗
Locally adaptive artificial viscosity strategies for Lagrangian hydrodynamics
To accurately model inviscid flow with shock waves using staggered-grid Lagrangian hydrodynamics, artificial viscosity is introduced to convert kinetic energy into internal energy, thereby providing a mechanism to generate the required entropy increase across shocks. In this paper, we propose a new method for constructing an adaptive, artificial viscosity in the context of one-dimensional, staggered-grid Lagrangian hydrodynamics. Our adaptive, artificial viscosity is defined in terms of two parameters that depend on density locally in a neighborhood around a shock, and hence, vary cell-by-cell. Our methodology is based on building a reference set of pre-computed optimal, globally constant artificial viscosity coefficients for a family of isolated shock test problems. For arbitrary flows, the evaluation of the unknown coefficients is automated by first estimating shock intensity locally, and second computing the corresponding adaptive parameter value via interpolation over the data from the pre-computed reference set of optimal, globally constant coefficient values. To illustrate the performance of our new approach, we compare our results against two existing methods. The first method is a limiter-based approach, which relies on estimating velocity gradients of the flow, and the second method is utilized in several commercial codes. Finally, we demonstrate that our new adaptive methodology produces more accurate results for a variety of tests with propagating shock waves, as well as for the aforementioned family of isolated shock problems.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Three dimensional numerical scattering from ocean-like surfaces
A new method of calculating the electric field scattered from two-dimensional conducting rough surfaces is presented. This method uses the operator conjugate gradients technique and the extended boundary condition equation. This avoids the build-up of round-off errors encountered in the inversion of large matrices. It also avoids the singularities of the integral kernels usually encountered in the method of moments. Techniques for improving the rate of convergence by applying physical constraints are discussed. This method is presently being used to calculate the polarization signatures from simulated ocean-like surfaces.
Rodriguez, Ernesto↗
Gradient-based optimization of 3D MHD equilibria
Using recently developed adjoint methods for computing the shape derivatives of functions that depend on magnetohydrodynamic (MHD) equilibria (Antonsen et al. , J. Plasma Phys. , vol. 85, issue 2, 2019; Paul et al. , J. Plasma Phys. , vol. 86, issue 1, 2020), we present the first example of analytic gradient-based optimization of fixed-boundary stellarator equilibria. We take advantage of gradient information to optimize figures of merit of relevance for stellarator design, including the rotational transform, magnetic well and quasi-symmetry near the axis. With the application of the adjoint method, we reduce the number of equilibrium evaluations by the dimension of the optimization space ( ${\sim }50\text {--}500$ ) in comparison with a finite-difference gradient-based method. We discuss regularization objectives of relevance for fixed-boundary optimization, including a novel method that prevents self-intersection of the plasma boundary. We present several optimized equilibria, including a vacuum field with very low magnetic shear throughout the volume.
Physics↗
Mapping structures and dynamics with frequency-correlated diffusion exchange
Understanding molecular motion in diffusion-driven complex environments is critical for designing sustainable materials and improving chemical processes. Here, we introduce a multidimensional nuclear magnetic resonance (NMR) method that captures how molecular populations exchange across different dynamic regimes. By extending the modulated gradient spin-echo technique to include frequency-frequency correlations, our approach reveals diffusion pathways that are otherwise obscured in heterogeneous systems. Implemented on a unilateral NMR magnet, the method eliminates gradient pulsing constraints and accesses dynamics in the kilohertz regime. We apply this technique to swelling and acid-catalyzed deconstruction of cross-linked and linear polymers to observe how structural heterogeneity evolves over time. By linking molecular motion to topology and chemical state, we extract physical metrics such as fractal surface dimensionality and reaction wavefront velocity, properties inaccessible with standard diffusion measurements. This work expands the capabilities of NMR for probing soft matter, with implications for polymer recycling and materials design.
Fricke, Sophia N. [University of California, Berke↗
Hybrid analysis of deformation using interferometric measurements
Two hybrid methods of deformation analysis combining moire interferometry and finite element techniques are introduced. The hybrid methods are applied to determine the strain fields in the test section of a composite cruciform specimen in order to evaluate the uniformity of the shear strain field. Both methods show that there are significant shear strain gradients in the specimen test section. One of the hybrid methods - the full-field analysis using the random-node mesh generator - detects the presence of material property nonuniformity which is assumed uniform in the boundary-node method.
Lin, S.-Y.↗
Imaging extended single crystal lattice distortion fields with multi-peak Bragg ptychography
Recent advances in phase-retrieval-based x-ray imaging methods have demonstrated the ability to reconstruct 3D distortion vector fields within a nanocrystal by using coherent diffraction information from multiple crystal Bragg reflections. However, these works do not provide a solution to the challenges encountered in imaging lattice distortions in crystals with significant defect content that result in phase wrapping. Moreover, these methods only apply to isolated crystals smaller than the x-ray illumination, and therefore cannot be used for imaging of distortions in extended crystals. We introduce multi-peak Bragg ptychography which addresses both challenges via an optimization framework that combines stochastic gradient descent and phase unwrapping methods for robust image reconstruction of lattice distortions and defects in extended crystals. Our work uses modern automatic differentiation toolsets so that the method is easy to extend to other settings and easy to implement in high-performance computers. This work is particularly timely given the broad interest in using the increased coherent flux in fourth-generation synchrotrons for innovative material research.
36 MATERIALS SCIENCE↗
Adjoint DSMC for nonlinear Boltzmann equation constrained optimization
Applications for kinetic equations such as optimal design and inverse problems often involve finding unknown parameters through gradient-based optimization algorithms. Based on the adjoint-state method, we derive two different frameworks for approximating the gradient of an objective functional constrained by the nonlinear Boltzmann equation. While the forward problem can be solved by the DSMC method, it is difficult to efficiently solve the high-dimensional continuous adjoint equation obtained by the “optimize-then-discretize” approach. This challenge motivates us to propose an adjoint DSMC method following the “discretize-then-optimize” approach for Boltzmann-constrained optimization. We also analyze the properties of the two frameworks and their connections. Here, several numerical examples are presented to demonstrate their accuracy and efficiency.
97 MATHEMATICS AND COMPUTING↗
HPLC/ESI-quadrupole ion trap mass spectrometry for characterization and direct quantification of amphoteric and nonionic surfactants in aqueous samples
An amphoteric (cocamidopropylbetaine, CAPB) and a nonionic (alcohol polyethoxylate, AE) surfactant were characterized by electrospray ionization quadrupole ion trap mass spectrometry (ESI-MS) as to their homologue distribution and ionization/fragmentation chemistry. Quantitative methods involving reversed-phase gradient HPLC and (+)ESI-MSn were developed to directly determine these surfactants in hydroponic plant growth medium that received simulated graywater. The predominant homologues, 12 C alkyl CAPB and 9 EO AE, were monitored to represent the total amount of the respective surfactants. The methods demonstrated dynamic linear ranges of 0.5-250 ng (r2 > 0.996) for CAPB and 8-560 ng (r2 > 0.998) for AE homologue mixture, corresponding to minimum quantification limits of 25 ppb CAPB and 0.4 ppm AE with 20-microL injections. This translated into an even lower limit for individual components due to the polydispersive nature of the surfactants. The procedure was successfully employed for the assessment of CAPB and AE biodegradation in a hydroponic plant growth system used as a graywater bioreactor.
NASA Discipline Life Support Systems↗
Research on diabatic initialization
The objective of this research project is to contribute to improvement in the synoptic analyses in the tropics for numerical weather prediction and climate research. In addition to a prediction model, four dimensional data assimilation systems have two principal components. One is objective analysis and the other is initialization. Various methods of objective analysis are designed primarily to analyze the mass and rotational wind fields. Methods of initialization are developed to obtain the irrotational wind and its associated vertical velocity field which are balanced with the mass field and free from meteorological noise. There are essentially three approaches to the problem of initialization: quasi-geostrophic theory, bounded derivative method and nonlinear normal mode method. In the midlatitudes, these approaches generally produce satisfactory results even without diabatic effects for large-scale motions. In the tropics, the situation is quite different from that in the midlatitudes. Because of a small magnitude of the Coriolis parameter and a weak horizontal temperature gradient in the tropics, any method of initialization must incorporate diabatic effects. In fact, it can be said that understanding the problem of diabatic initialization is the key to improving the analysis and weather forecasting in the tropics.
Kasahara, Akira↗
Dataset from: "Adjoint methods for quasisymmetry of vacuum fields on a surface"
Adjoint methods can speed up stellarator optimisation by providing gradient information more efficiently compared to finite-difference evaluations. Adjoint methods are herein applied to vacuum magnetic fields, with objective functions targeting quasisymmetry and a rotational transform value on a surface. To measure quasisymmetry, a novel way of evaluating approximate flux coordinates on a single flux surface without the assumption of a neighbourhood of flux surfaces is proposed. The shape gradients obtained from the adjoint formalism are evaluated numerically and verified against finite-difference evaluations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Derivative-free stochastic optimization via adaptive sampling strategies
In this paper, we present a novel derivative-free framework for solving unconstrained stochastic optimization problems. Many problems in fields ranging from simulation optimization to reinforcement learning to quantum computing involve settings where only stochastic function values are obtained via a zeroth-order oracle, which has no available gradient information and necessitates the usage of derivative-free optimization methodologies. Our approach includes estimating gradients using stochastic function evaluations and integrating adaptive sampling techniques to control the accuracy in these stochastic approximations. Our framework encapsulates several gradient estimation techniques, including standard finite-difference, Gaussian smoothing, sphere smoothing, randomized coordinate finite-difference, and randomized subspace finite-difference methods. We provide theoretical convergence guarantees for our framework and analyze the worst-case iteration and sample complexities associated with each gradient estimation method. Finally, we demonstrate the empirical performance of the methods on logistic regression and nonlinear least squares problems.
Adaptive sampling↗
Adjoint methods for quasi-symmetry of vacuum fields on a surface
Adjoint methods can speed up stellarator optimisation by providing gradient information more efficiently compared with finite-difference evaluations. Adjoint methods are herein applied to vacuum magnetic fields, with objective functions targeting quasi-symmetry and a rotational transform value on a surface. Here to measure quasi-symmetry, a novel way of evaluating approximate flux coordinates on a single flux surface without the assumption of a neighbourhood of flux surfaces is proposed. The shape gradients obtained from the adjoint formalism are evaluated numerically and verified against finite-difference evaluations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Function space quasi-Newton algorithms for optimal control problems with bounded controls and singular arcs
Two existing function space algorithms, Davidon and projected gradient, are modified so that they may handle directly control variable inequality constraints. A third quasi-Newton type algorithm developed by C. G. Broyden is extended to optimal control problems. The Broyden algorithm is further modified so that it also may handle directly control variable inequality constraints. These methods along with a pure gradient and two conjugate gradient algorithms are simulated on three relatively simple yet representative bounded control problems, two of which have singular subarcs. Overall the Broyden algorithm was found to be superior. The most notable result of the study was the clear superiority of the Broyden and Davidon algorithms in producing a sharp interior control subarc.
Edge, E. R.↗
Novel Relativistic Electronic Structure Theories for Actinide-Containing Compounds
Actinides of importance to basic energy sciences contain electrons moving at speed comparable to the speed of light. Reliable computational simulation of these electrons and hence actinide chemistry requires accurate description of relativistic effects. The present project advances computational actinide chemistry with development of new methodologies, algorithms, and computer programs in relativistic quantum chemistry, as well as applications to actinide chemistry and spectroscopy. A new “electrons-only” exact two-component approach has been developed to provide efficient treatments of relativistic effects, while maintaining chemical accuracy. New computational algorithms developed here extend the applicability of relativistic electron-correlation methods to larger molecules. The method-development work in this project also features the first implementation of analytic gradient technique for relativistic electron-correlation methods, which provides significantly enhanced ability to compute properties for molecules containing actinides. The applicability and usefulness of these new methods and computer programs have been demonstrated in calculations of actinide-containing molecules to facilitate understanding of actinide chemistry and spectroscopy.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry
Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A proposed method for wind velocity measurement from space
An investigation was made of the feasibility of making wind velocity measurements from space by monitoring the apparent change in the refractive index of the atmosphere induced by motion of the air. The physical principle is the same as that resulting in the phase changes measured in the Fizeau experiment. It is proposed that this phase change could be measured using a three cornered arrangement of satellite borne source and reflectors, around which two laser beams propagate in opposite directions. It is shown that even though the velocity of the satellites is much larger than the wind velocity, factors such as change in satellite position and Doppler shifts can be taken into account in a reasonable manner and the Fizeau phase measured. This phase measurement yields an average wind velocity along the ray path through the atmosphere. The method requires neither high accuracy for satellite position or velocity, nor precise knowledge of the refractive index or its gradient in the atmosphere. However, the method intrinsically yields wind velocity integrated along the ray path; hence to obtain higher spatial resolution, inversion techniques are required.
Censor, D.↗