Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical differentiation”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Self-Similar Compressible Free Vortices

Lie group methods are used to find both exact and numerical similarity solutions for compressible perturbations to all incompressible, two-dimensional, axisymmetric vortex reference flow. The reference flow vorticity satisfies an eigenvalue problem for which the solutions are a set of two-dimensional, self-similar, incompressible vortices. These solutions are augmented by deriving a conserved quantity for each eigenvalue, and identifying a Lie group which leaves the reference flow equations invariant. The partial differential equations governing the compressible perturbations to these reference flows are also invariant under the action of the same group. The similarity variables found with this group are used to determine the decay rates of the velocities and thermodynamic variables in the self-similar flows, and to reduce the governing partial differential equations to a set of ordinary differential equations. The ODE's are solved analytically and numerically for a Taylor vortex reference flow, and numerically for an Oseen vortex reference flow. The solutions are used to examine the dependencies of the temperature, density, entropy, dissipation and radial velocity on the Prandtl number. Also, experimental data on compressible free vortex flow are compared to the analytical results, the evolution of vortices from initial states which are not self-similar is discussed, and the energy transfer in a slightly-compressible vortex is considered.

vonEllenrieder, Karl↗
NASA NTRS1995

Multithreaded Model for Dynamic Load Balancing Parallel Adaptive PDE Computations

We present a multithreaded model for the dynamic load-balancing of numerical, adaptive computations required for the solution of Partial Differential Equations (PDE's) on multiprocessors. Multithreading is used as a means of exploring concurrency in the processor level in order to tolerate synchronization costs inherent to traditional (non-threaded) parallel adaptive PDE solvers. Our preliminary analysis for parallel, adaptive PDE solvers indicates that multithreading can be used an a mechanism to mask overheads required for the dynamic balancing of processor workloads with computations required for the actual numerical solution of the PDE's. Also, multithreading can simplify the implementation of dynamic load-balancing algorithms, a task that is very difficult for traditional data parallel adaptive PDE computations. Unfortunately, multithreading does not always simplify program complexity, often makes code re-usability not an easy task, and increases software complexity.

Chrisochoides, Nikos↗
DOE OSTI2024

Semi-inclusive single-jet production in DIS at next-to-leading order in the Color Glass Condensate

Within the Color Glass Condensate (CGC) effective field theory, we derive the next-to-leading order (NLO) cross-section for the single-jet semi-inclusive cross-section in deep inelastic scattering (DIS) at small x, for both longitudinally and transversely polarized virtual photons. We provide analytic expressions, valid at finite N c and suitable for numerical evaluation, for both the cross-section differential in rapidity and transverse momentum and the cross-section differential in rapidity only. Our NLO formulae demonstrate that the very forward rapidity regime is plagued by large double logarithmic corrections coming from phase space constraints on soft gluons close to the kinematic threshold for jet production. A joint resummation of small-x and threshold logarithms at single logarithmic accuracy is proposed to remedy the instability of the cross-section in this regime. By integrating over the single-jet phase space, we recover known results for the NLO DIS structure functions at small x, previously obtained using the optical theorem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Hierarchical Model Reduction Driven by Machine Learning for Parametric Advection-Diffusion-Reaction Problems in the Presence of Noisy Data

Abstract We propose a new approach to generate a reliable reduced model for a parametric elliptic problem, in the presence of noisy data. The reference model reduction procedure is the directional HiPOD method, which combines Hierarchical Model reduction with a standard Proper Orthogonal Decomposition, according to an offline/online paradigm. In this paper we show that directional HiPOD looses in terms of accuracy when problem data are affected by noise. This is due to the interpolation driving the online phase, since it replicates, by definition, the noise trend. To overcome this limit, we replace interpolation with Machine Learning fitting models which better discriminate relevant physical features in the data from irrelevant unstructured noise. The numerical assessment, although preliminary, confirms the potentialities of the new approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Parallel-in-Time Solution of Allen-Cahn Equations by Integrating Operator Learning into the Parareal Method

While recent advances in deep learning have shown promising efficiency gains in solving time-dependent partial differential equations (PDEs), matching the accuracy of conventional numerical solvers still remains a challenge. One strategy to improve the accuracy of deep learning-based solutions for time-dependent PDEs is to use the learned model as the coarse propagator in the Parareal method and a traditional numerical method as the fine solver. However, successful integration of deep learning into the Parareal method requires consistency between the coarse and fine solvers, particularly for PDEs exhibiting rapid changes such as sharp transitions. Here, to ensure this consistency, we propose using convolutional neural networks (CNNs) to learn the fully discrete time-stepping operator defined by the same numerical scheme employed as the fine solver. We demonstrate the effectiveness of the proposed method in solving the classical and mass-conservative Allen–Cahn (AC) equations. Through iterative updates in the Parareal algorithm, our approach achieves a significant computational speedup compared to traditional fine solvers while converging to high-accuracy solutions. Our results highlight that the proposed hybrid Parareal algorithm effectively accelerates simulations, particularly when implemented on multiple GPUs, and converges to the desired accuracy in only a few iterations. Another advantage of our method is that the CNN model is trained on trajectory-based data generated from random initial conditions, such that the trained model can be used to solve the AC equations with various initial conditions without retraining. This work demonstrates the potential of integrating neural network methods into parallel-in-time frameworks for efficient and accurate simulations of time-dependent PDEs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Approximating the linear response of physical chaos

Abstract Parametric derivatives of statistics are highly desired quantities in prediction, design optimization and uncertainty quantification. In the presence of chaos, the rigorous computation of these quantities is certainly possible, but mathematically complicated and computationally expensive. Based on Ruelle’s formalism, this paper shows that the sophisticated linear response algorithm can be dramatically simplified in higher-dimensional systems featuring statistical homogeneity in the physical space. We argue that the contribution of the SRB (Sinai–Ruelle–Bowen) measure gradient, which is an integral yet the most cumbersome part of the full algorithm, is negligible if the objective function is appropriately aligned with unstable manifolds. This abstract condition could potentially be satisfied by a vast family of real-world chaotic systems, regardless of the physical meaning and mathematical form of the objective function and perturbed parameter. We demonstrate several numerical examples that support these conclusions and that present the use and performance of a simplified linear response algorithm. In the numerical experiments, we consider physical models described by differential equations, including Lorenz 96 and Kuramoto–Sivashinsky.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

hp -VPINNs: Variational physics-informed neural networks with domain decomposition

We formulate a general framework for hp-variational physics-informed neural networks (hp-VPINNs) based on the nonlinear approximation of shallow and deep neural networks and hp-refinement via domain decomposition and projection onto the space of high-order polynomials. The trial space is the space of neural network, which is defined globally over the entire computational domain, while the test space contains piecewise polynomials. Specifically in this study, the hp-refinement corresponds to a global approximation with a local learning algorithm that can efficiently localize the network parameter optimization. Here, we demonstrate the advantages of hp-VPINNs in both accuracy and training cost for several numerical examples of function approximation and in solving differential equations.

42 ENGINEERING↗
DOE OSTI2025

Accelerating phase field simulations through a hybrid adaptive Fourier neural operator with U-net backbone

Prolonged contact between a corrosive liquid and metal alloys can cause progressive dealloying. For one such process as liquid-metal dealloying (LMD), phase field models have been developed to understand the mechanisms leading to complex morphologies. However, the LMD governing equations in these models often involve coupled non-linear partial differential equations (PDE), which are challenging to solve numerically. In particular, numerical stiffness in the PDEs requires an extremely refined time step size (on the order of 10 -12 s or smaller). This computational bottleneck is especially problematic when running LMD simulation until a late time horizon is required. This motivates the development of surrogate models capable of leaping forward in time, by skipping several consecutive time steps at-once. In this paper, we propose a U-shaped adaptive Fourier neural operator (U-AFNO), a machine learning (ML) based model inspired by recent advances in neural operator learning. U-AFNO employs U-Nets for extracting and reconstructing local features within the physical fields, and passes the latent space through a vision transformer (ViT) implemented in the Fourier space (AFNO). We use U-AFNOs to learn the dynamics of mapping the field at a current time step into a later time step. We also identify global quantities of interest (QoI) describing the corrosion process (e.g., the deformation of the liquid-metal interface, lost metal, etc.) and show that our proposed U-AFNO model is able to accurately predict the field dynamics, in spite of the chaotic nature of LMD. Most notably, our model reproduces the key microstructure statistics and QoIs with a level of accuracy on par with the high-fidelity numerical solver, while achieving a significant 11, 200 × speed-up on a high-resolution grid when comparing the computational expense per time step. Finally, we also investigate the opportunity of using hybrid simulations, in which we alternate forward leaps in time using the U-AFNO with high-fidelity time stepping. We demonstrate that while advantageous for some surrogate model design choices, our proposed U-AFNO model in fully auto-regressive settings consistently outperforms hybrid schemes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Juqbox.jl

The Juqbox.jl package implements functionality for solving the quantum optimal control problem for realizing logical gates in closed quantum systems. The dynamics of the quantum system is modeled by Schroedinger's equation, which takes to form of a linear system of ordinary differential equations (ODE). Juqbox.jl solves this ODE by numerical time stepping and applies a gradient-based optimization technique to determine control pulses for driving an initial state to a final state, according to the desired logical gate transformation. To evaluate the gradient, Juqbox.jl applies the ``first discretize, then optimize'' approach based on a discrete adjoint time stepping technique. The actual optimization is performed by the open source Ipopt libaray. Juqbox.jl is written in the Julia programming language which, among many other features, provides a convenient interface to the Ipopt library.

PETERSSON, NILSA.↗
DOE OSTI2022

Double-power-law Feature of Energetic Particles Accelerated at Coronal Shocks

Recent observations have shown that in many large solar energetic particle (SEP) events the event-integrated differential spectra resemble double power laws. We perform numerical modeling of particle acceleration at coronal shocks propagating through a streamer-like magnetic field by solving the Parker transport equation, including protons and heavier ions. We find that for all ion species the energy spectra integrated over the simulation domain can be described by a double power law, and the break energy depends on the ion charge-to-mass ratio as E B ~ (Q/A) α , with α varying from 0.16 to 1.2 by considering different turbulence spectral indices. We suggest that the double-power-law distribution may emerge as a result of the superposition of energetic particles from different source regions where the acceleration rates differ significantly due to particle diffusion. The diffusion and mixing of energetic particles could also provide an explanation for the increase of Fe/O at high energies as observed in some SEP events. Although further mixing processes may occur, our simulations indicate that either a power-law break or rollover can occur near the Sun and predict that the spectral forms vary significantly along the shock front, which may be examined by upcoming near-Sun SEP measurements from the Parker Solar Probe and Solar Orbiter.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1971

Transient-forced convection film boiling on an isothermal flat plate

An approach for the solution of transient-forced convection film boiling on an isothermal flat plate using the boundary layer model is developed. The similarity variables are used to convert the governing partial differential equations to ordinary ones. The results of numerical solutions of these ordinary equations indicate that the transient process can be classified as one-dimensional conduction, intermediate, and the steady state regions. The time required for the one-dimensional conduction and the time necessary to attain a steady state condition are obtained. The influence of interfacial shear is seen to be negligible while the Prandtl Number and the ratio (C sub p delta T divided by h sub fg times Pr) have major influence. The use of local similarity approximations for the intermediate regime facilitates prediction of complete boundary layer growth. Using the ratio of time at any instant to the steady state time as abscissa, the curves representing the boundary layer growth can be merged into a single mean curve within 5 percent. Further, the analysis shows that the average rate of heat transfer during transient is 50 to 100 percent higher than those at steady state. The average rate of vapor convected away is 10 to 15 percent lower than at steady state while the average rate of accumulation to form the vapor layer is 1 to 14 times larger. Further, the total heat transferred during transient increases and the evaporation decreases for increasing values of C sub p delta T divided by h sub fg times Pr

Nagendra, H. R.↗
NASA NTRS1971

Concepts for a theoretical and experimental study of lifting rotor random loads and vibrations. Phase 5B: Analysis of gust alleviation methods and rotor dynamic stability

The effects of various gust alleviation methods on the random blade response in flapping were studied analytically, assuming a rigid rotor support. The analytical model assumes rigid flapping blades with elastic root restraints. Linearized equations which are approximately valid at low lift conditions were used. Because of the interblade coupling from the feedback devices, the method of multiblade generalized coordinates was most convenient and was extended to include coning, differential coning, and warping of the rotor. The numerical examples cover the dynamic stability characteristics as affected by feedback gains of three- to six-bladed rotors. The number of blades had large effects on stability limits and modal time functions at these limits. The random flapping response of the blades to atmospheric turbulence was determined at 1.6 rotor advance ratio using feedback gains below the stability limit. The most effective reduction of the flapping response per unit gain was achieved with a rotor coning angle feedback.

Hohenemser, K. H.↗
NASA NTRS1972

Transient-forced convection film boiling on an isothermal flat plate.

Development of a new approach for the solution of transient-forced convection film boiling on an isothermal flat plate using the boundary layer model. The similarity variables are used to convert the governing partial differential equations to ordinary ones. The results of numerical solutions of these ordinary equations indicate that the transient process can be classified as one-dimensional conduction, intermediate, and the steady-state regions. The time required for the one-dimensional conduction and the time necessary to attain a steady-state condition are obtained. The use of local similarity approximations for the intermediate regime facilitates prediction of complete boundary layer growth. Using the ratio of time at any instant to the steady-state time as abscissa, the curves representing the boundary layer growth can be merged into a single mean curve within 5%. Further, the analysis shows that the average rate of heat transfer during transient is 50 to 100% higher than those at steady state. The average rate of vapor convected away is 10 to 15% lower than at steady state while the average rate of accumulation to form the vapor layer is 1 to 14 times larger.

Nagendra, H. R.↗
NASA NTRS1972

Space shuttle guidance, navigation and control equation document no. 4: Precision state and filter weighting matrix extrapolation

The Precision State and Filter Weighting Matrix Extrapolation Routine is described which provides the capability to extrapolate any spacecraft geocentric state vector either backwards or forwards in time through a force field consisting of the earth's primary central-force gravitational attraction and a superimposed perturbing acceleration. The routine also provides the capability of extrapolating the filter-weighting matrix along the precision trajectory. This matrix is a square root form of the error covariance matrix and contains statistical information relative to the accuracies of the state vectors and certain other optionally estimated quantities. The routine is a cooled algorithm for the numerical solution of modified forms of the basic differential equations which are satisfied by the geocentric state vector of the spacecraft's center of mass and by the filter-weighting matrix.

Robertson, W. M.↗
NASA NTRS1972

Precision state and filter weighting matrix extrapolation, revision 3

The matrix extrapolation routine presented is a coded algorithm for the numerical solution of modified forms of the basic differential equations which are satisfied by the geocentric state vector of a spacecraft's center of mass and by the filter-weighting matrix. The equations are described in detail and the various input and output variables designated. The capability is thus established to extrapolate any spacecraft's geocentric state vector, either backwards or forwards in time, through a force field consisting of the earth's primary central-force gravitational attraction, and a superimposed perturbing acceleration.

Robertson, W. M.↗
NASA NTRS1975

Stabilization by modification of the Lagrangian

In order to reduce the error growth during a numerical integration, a method of stabilization of the differential equations of the Keplerian motion is offered. It is characterized by the use of the eccentric anomaly as an independent variable in such a way that the time transformation is given by a generalized Lagrange formalism. The control terms in the equations of motion obtained by this modified Lagrangian give immediately a completely Liapunov-stable set of differential equations. In contrast to other publications, here the equation of time integration is modified by a control term which leads to an integral which defined the time element for the perturbed Keplerian motion.

Baumgarte, J. W.↗