Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solve time”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗
NASA NTRS1980

Prediction and measurement of turbulent aerodynamic trailing edge flows

A viscous-inviscid interaction algorithm is developed for prediction of two-dimensional mean and fluctuating velocity distributions in the wake immediately downstream of an airfoil trailing edge. A composite pressure field is defined, and a Poisson equation solved for transverse pressure variations. A parabolized form of the time-averaged steady Navier-Stokes equations are solved in conjunction with a viscous-augmented two-dimensional inviscid potential flow analysis. A tensor constitutive equation is employed to predict Reynolds stress distributions from solutions of a turbulence kinetic energy two equation closure model. Numerical predictions compared favorably with detailed experimental data for mean and fluctuating velocities, and Reynolds shear stress distributions, in the trailing edge region of a NACA 63-012 airfoil.

Baker, A. J.↗
NASA NTRS1992

Combining constraint satisfaction and local improvement algorithms to construct anaesthetists' rotas

A system is described which was built to compile weekly rotas for the anaesthetists in a large hospital. The rota compilation problem is an optimization problem (the number of tasks which cannot be assigned to an anaesthetist must be minimized) and was formulated as a constraint satisfaction problem (CSP). The forward checking algorithm is used to find a feasible rota, but because of the size of the problem, it cannot find an optimal (or even a good enough) solution in an acceptable time. Instead, an algorithm was devised which makes local improvements to a feasible solution. The algorithm makes use of the constraints as expressed in the CSP to ensure that feasibility is maintained, and produces very good rotas which are being used by the hospital involved in the project. It is argued that formulation as a constraint satisfaction problem may be a good approach to solving discrete optimization problems, even if the resulting CSP is too large to be solved exactly in an acceptable time. A CSP algorithm may be able to produce a feasible solution which can then be improved, giving a good, if not provably optimal, solution.

Smith, Barbara M.↗
NASA NTRS1990

A multiblock multigrid three-dimensional Euler equation solver

Current aerodynamic designs are often quite complex (geometrically). Flexible computational tools are needed for the analysis of a wide range of configurations with both internal and external flows. In the past, geometrically dissimilar configurations required different analysis codes with different grid topologies in each. The duplicity of codes can be avoided with the use of a general multiblock formulation which can handle any grid topology. Rather than hard wiring the grid topology into the program, it is instead dictated by input to the program. In this work, the compressible Euler equations, written in a body-fitted finite-volume formulation, are solved using a pseudo-time-marching approach. Two upwind methods (van Leer's flux-vector-splitting and Roe's flux-differencing) were investigated. Two types of explicit solvers (a two-step predictor-corrector and a modified multistage Runge-Kutta) were used with multigrid acceleration to enhance convergence. A multiblock strategy is used to allow greater geometric flexibility. A report on simple explicit upwind schemes for solving compressible flows is included.

Cannizzaro, Frank E.↗
DOE OSTI2022

On the Convergence of Overlapping Schwarz Decomposition for Nonlinear Optimal Control

Here, we study the convergence properties of an overlapping Schwarz decomposition algorithm for solving nonlinear optimal control problems (OCPs). The algorithm decomposes the time domain into a set of overlapping subdomains, and solves all subproblems defined over subdomains in parallel. The convergence is attained by updating primal-dual information at the boundaries of overlapping subdomains. We show that the algorithm exhibits local linear convergence, and that the convergence rate improves exponentially with the overlap size. We also establish global convergence results for a general quadratic programming, which enables the application of the Schwarz scheme inside second-order optimization algorithms (e.g., sequential quadratic programming). The theoretical foundation of our convergence analysis is a sensitivity result of nonlinear OCPs, which we call "exponential decay of sensitivity" (EDS). Intuitively, EDS states that the impact of perturbations at domain boundaries (i.e., initial and terminal time) on the solution decays exponentially as one moves into the domain. Here, we expand a previous analysis available in the literature by showing that EDS holds for both primal and dual solutions of nonlinear OCPs, under uniform second-order sufficient condition, controllability condition, and boundedness condition. We conduct experiments with a quadrotor motion planning problem and a partial differential equations (PDE) control problem to validate our theory, and show that the approach is significantly more efficient than alternating direction method of multipliers and as efficient as the centralized interior-point solver.

42 ENGINEERING↗
DOE OSTI2023

A pseudospectral implicit particle-in-cell method with exact energy and charge conservation

The standard particle-in-cell (PIC) method employs explicit finite-difference (FD) methods (e.g. the leap-frog scheme) for both spatial and temporal integrations. Here, we employ a pseudospectral method for solving the Poisson equation and a fully implicit time integration to achieve exact energy conservation. The advantage of a pseudospectral field solver is its spectral accuracy in solving field solutions. Earlier studies of implicit time integration of PIC FD equations can enforce exact energy exchange between field and particles, resulting in exact energy-conserving schemes. Here, we prove that the exact energy conservation property can be carried over to the pseudospectral scheme. Simultaneously, we provide a solution to ensure a pseudospectral charge continuity equation. We demonstrate the new scheme in a 2D electrostatic PIC code. In conclusion, theoretical results are confirmed via numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Influence of polycarboxylate superplasticizer, citric acid and their combination on the hydration and workability of calcium sulfoaluminate cement

Highlights: • The dispersion efficiency of polycarboxylate superplasticizers drops quickly over time. • Citric acid maintains the dispersion efficiency of PCE over time. • Competitive adsorption between citric acid and PCE occurs. • The combination of citric acid with PCE decreases the initial dispersion efficiency of PCE. PCEs are well known to improve the initial fluidity of CSA. However, their dispersion efficiency drops quickly over time. This issue can be solved by incorporating retarders. In this context, this paper deals with the influence of citric acid, used as a retarder, PCE and their combination on the hydration and workability of CSA. Isothermal calorimetry, XRD and TG analysis were used to describe the hydration process, while workability was characterized with the mini-cone test. Adsorption behavior was investigated using total organic carbon analyzer coupled with ion chromatography. Results show that the introduction of citric acid retained the dispersion efficiency of PCE over time. However, the initial dispersion efficiency of PCE was decreased by citric acid as the latter tend to adsorb first on the surface of cement grains, inhibiting the adsorption of PCE. A dispersion model was proposed to describe the acting mechanism of these admixtures on CSA.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Countermeasures to Neurobehavioral Deficits from Cumulative Partial Sleep Deprivation During Space Flight

This project is concerned with identifying ways to prevent neurobehavioral and physical deterioration due to inadequate sleep in astronauts during long-duration manned space flight. The performance capability of astronauts during extended-duration space flight depends heavily on achieving recovery through adequate sleep. Even with appropriate circadian alignment, sleep loss can erode fundamental elements of human performance capability including vigilance, cognitive speed and accuracy, working memory, reaction time, and physiological alertness. Adequate sleep is essential during manned space flight not only to ensure high levels of safe and effective human performance, but also as a basic regulatory biology critical to healthy human functioning. There is now extensive objective evidence that astronaut sleep is frequently restricted in space flight to averages between 4 hr and 6.5 hr/day. Chronic sleep restriction during manned space flight can occur in response to endogenous disturbances of sleep (motion sickness, stress, circadian rhythms), environmental disruptions of sleep (noise, temperature, light), and curtailment of sleep due to the work demands and other activities that accompany extended space flight operations. The mechanism through which this risk emerges is the development of cumulative homeostatic pressure for sleep across consecutive days of inadequate sleep. Research has shown that the physiological sleepiness and performance deficits engendered by sleep debt can progressively worsen (i.e., accumulate) over consecutive days of sleep restriction, and that sleep limited to levels commonly experienced by astronauts (i.e., 4 - 6 hr per night) for as little as 1 week, can result in increased lapses of attention, degradation of response times, deficits in complex problem solving, reduced learning, mood disturbance, disruption of essential neuroendocrine, metabolic, and neuroimmune responses, and in some vulnerable persons, the emergence of uncontrolled sleep attacks. The prevention of cumulative performance deficits and neuroendocrine disruption from sleep restriction during extended duration space flight involves finding the most effective ways to obtain sleep in order to maintain the high-level cognitive and physical performance functions required for manned space flight. There is currently a critical deficiency in knowledge of the effects of how variations in sleep duration and timing relate to the most efficient return of performance per unit time invested in sleep during long-duration missions, and how the nature of sleep physiology (i.e., sleep stages, sleep electroencephalographic [EEG] power spectral analyses) change as a function of sleep restriction and performance degradation. The primary aim of this project is to meet these critical deficiencies through utilization of a response surface experimental paradigm, testing in a dose-response manner, varying combinations of sleep duration and timing, for the purpose of establishing how to most effectively limit the cumulative adverse effects on human performance and physiology of chronic sleep restriction in space operations.

Dinges, David F.↗
DOE OSTI2024

Quantum mechanical model of crossing and anti-crossing points in 3D full-band Monte Carlo simulations

This work presents a 3D quantum mechanics based model to address the physics at band structure crossing/anti-crossing points in full band Monte Carlo (FBMC) simulations. The model solves the Krieger and Iafrate (KI) equations in real time using pre-computed coefficients at k-points spatially sampled within the first Brillouin zone. Solving the KI equations in real time makes this model applicable for all electric fields, which enables its use in FBMC device simulations. In this work, a two-level refinement scheme is used to aggressively sample regions in proximity to band crossings for accurate solutions to the KI equations and coarsely sample everywhere else to limit the number of k-points used. The presented sampling method is demonstrated on the band structure of silicon but is effective for the band structure of any semiconductor material. Next, the adaptation of the fully quantum KI model into the semi-classical FBMC method is discussed. Finally, FBMC simulations of hole transport in 4H silicon carbide with and without the KI model are performed. Results along different crystallographic directions for a wide range of electric fields are compared to previously published simulation and experimental values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Feedforward equilibrium trajectory optimization with GSPulse

One of the common tasks required for designing new plasma scenarios or evaluating capabilities of a tokamak is to design the desired equilibria using a Grad-Shafranov (GS) equilibrium solver. However, most standard equilibrium solvers are time-independent and do not include dynamic effects such as plasma current flux consumption, induced vessel currents, or voltage constraints. Another class of tools, plasma equilibrium evolution simulators, do include time-dependent effects. These are generally structured to solve the forward problem of evolving the plasma equilibrium given feedback-controlled voltages. In this work, we introduce GSPulse, a novel algorithm for equilibrium trajectory optimization, that is more akin to a pulse planner than a pulse simulator. GSPulse includes time-dependent effects and solves the inverse problem: given a user-specified set of target equilibrium shapes, as well as limits on the coil currents and voltages, the optimizer returns trajectories of the voltages, currents, and achievable equilibria. This task is useful for scoping performance of a tokamak and exploring the space of achievable pulses. The computed equilibria satisfy both Grad-Shafranov force balance and axisymmetric circuit dynamics. The optimization is performed by restructuring the free-boundary equilibrium evolution equations into a form where it is computationally efficient to optimize the entire dynamic sequence. GSPulse can solve for hundreds of equilibria simultaneously within a few minutes. GSPulse has been validated against NSTX-U and MAST-U experiments and against SPARC feedback control simulations, and is being used to perform scenario design for SPARC. The computed trajectories can be used as feedforward inputs that are connected to the feedback controller to inform and improve feedback performance. The code for GSPulse is available open-source at github.com/jwai-cfs/GSPulse_public.

equilibrium↗
NASA NTRS1994

Hypersonic flow past open cavities

The hypersonic flow over a cavity is investigated. The time-dependent compressible Navier-Stokes equations are numerically solved. An implicit algorithm, with a subiteration procedure to recover time accuracy, is used to perform the time-accurate computations. The objective of the study is to investigate the effects of Reynolds number and cavity dimensions. The comparsion of the computations with available experimental data, in terms of time mean static pressure, heat transfer, and Mach number, show good agreement. In the computations large vortex structures, which adversely affect the cavity flow characteristics, are observed at the rear of the cavity. A self-sustained oscillatory motion occurs within the cavity over a range of Reynolds number and cavity dimensions. The frequency spectra of the oscillations show good agreement with a modified semiempirical relation.

Morgenstern, Algacyr, Jr.↗
DOE OSTI2023

Demonstrating Computational Equivalence Between Continuous and Discrete Adjoint Methods by Calculating Time-Dependent Adjoint Solutions with Neutron Diffusion Models

The continuous adjoint method and the discrete adjoint method are two alternative approaches used to calculate adjoint solutions for adjoint systems. The continuous adjoint method derives adjoint equations analytically from continuous forward equations and then solves the adjoint equations either analytically or numerically in a discretized form whereas the discrete adjoint method calculates the adjoint solutions directly from the discretized forward equations. With regard to the methodology development and calculation procedure, distinct differences are well recognized between the two methods. For certain reasons, both methods are exclusively preferred and commonly used by different computational communities, but limited studies clarify the connections between the two adjoint methods from either of the communities. Herein, this paper demonstrates the computational equivalence between the continuous and discrete adjoint methods by investigating time-dependent adjoint solutions to the two-group neutron diffusion model in nuclear reactor analysis problems using both methods. Adjoint solutions can be used to estimate system parameters for reactor safety analysis. Appropriate final state conditions for the adjoint systems are specified in both of the methods, and the conditions are clarified with proper physical explanations. With the help of an event-based case study on neutron diffusion models, the accuracy of the time-dependent adjoint fluxes obtained from both methods is verified, and the pros and cons of both adjoint methods are examined. More importantly, the computational equivalence of both methods is demonstrated when they are applied to multigroup neutron diffusion systems. The advantage of calculating time-dependent adjoint fluxes by directly solving time-dependent adjoint systems rather than taking steady-state approximations as in common practice is also demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

Incompressible spectral-element method: Derivation of equations

A fractional-step splitting scheme breaks the full Navier-Stokes equations into explicit and implicit portions amenable to the calculus of variations. Beginning with the functional forms of the Poisson and Helmholtz equations, we substitute finite expansion series for the dependent variables and derive the matrix equations for the unknown expansion coefficients. This method employs a new splitting scheme which differs from conventional three-step (nonlinear, pressure, viscous) schemes. The nonlinear step appears in the conventional, explicit manner, the difference occurs in the pressure step. Instead of solving for the pressure gradient using the nonlinear velocity, we add the viscous portion of the Navier-Stokes equation from the previous time step to the velocity before solving for the pressure gradient. By combining this 'predicted' pressure gradient with the nonlinear velocity in an explicit term, and the Crank-Nicholson method for the viscous terms, we develop a Helmholtz equation for the final velocity.

Deanna, Russell G.↗
DOE OSTI2020

A low-rank solver for the stochastic unsteady Navier–Stokes problem

Here we study a low-rank iterative solver for the unsteady Navier–Stokes equations for incompressible flows with a stochastic viscosity. The equations are discretized using the stochastic Galerkin method, and we consider an all-at-once formulation where the algebraic systems at all the time steps are collected and solved simultaneously. The problem is linearized with Picard’s method. To efficiently solve the linear systems at each step, we use low-rank tensor representations within the Krylov subspace method, which leads to significant reductions in storage requirements and computational costs. Combined with effective mean-based preconditioners and the idea of inexact solve, we show that only a small number of linear iterations are needed at each Picard step. The proposed algorithm is tested with a model of flow in a two-dimensional symmetric step domain with different settings to demonstrate the computational efficiency.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1980

Steady state solutions to dynamically loaded periodic structures

The general problem of solving for the steady state (time domain) dynamic response (i.e., NASTRAN rigid format-8) of a general elastic periodic structure subject to a phase difference loading of the type encountered in traveling wave propagation problems was studied. Two types of structural configurations were considered; in the first type, the structure has a repeating pattern over a span that is long enough to be considered, for all practical purposes, as infinite; in the second type, the structure has structural rotational symmetry in the circumferential direction. The theory and a corresponding set of DMAP instructions which permits the NASTRAN user to automatically alter the rigid format-8 sequence to solve the intended class of problems are presented. Final results are recovered as with any ordinary rigid format-8 solution, except that the results are only printed for the typical periodic segment of the structure. A simple demonstration problem having a known exact solution is used to illustrate the implementation of the procedure.

Kalinowski, A. J.↗
NASA NTRS1996

Multigrid Acceleration of Time-Accurate Navier-Stokes Calculations

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is fully implicit in time and is unconditionally stable (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability. This is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. A multigrid-multiblock, steady-state, three-dimensional Navier-Stokes solver, TLNS3D, was modified to iteratively invert the equations at each physical time step. The implementation of this procedure in TLNS3D is discussed. The implications of applying several popular turbulence models to unsteady flow are also considered. Numerical results are presented to show the application of the scheme to various two-dimensional turbulent flows. The results of a three-dimensional laminar flow calculation are also given.

Melson, N. Duane↗
DOE OSTI2024

Train small, model big: Scalable physics simulators via reduced order modeling and domain decomposition

Numerous cutting-edge scientific technologies originate at the laboratory scale, but transitioning them to practical industry applications is a formidable challenge. Traditional pilot projects at intermediate scales are costly and time-consuming. An alternative, the pilot-scale model, relies on high-fidelity numerical simulations, but even these simulations can be computationally prohibitive at larger scales. To overcome these limitations, we propose a scalable, physics-constrained reduced order model (ROM) method. The ROM identifies critical physics modes from small-scale unit components, projecting governing equations onto these modes to create a reduced model that retains essential physics details. We also employ Discontinuous Galerkin Domain Decomposition (DG-DD) to apply ROM to unit components and interfaces, enabling the construction of large-scale global systems without data at such large scales. Here this method is demonstrated on the Poisson and Stokes flow equations, showing that it can solve equations about 15–40 times faster with only ~1% relative error. Furthermore, ROM takes one order of magnitude less memory than the full order model, enabling larger scale predictions at a given memory limitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Joint Modeling of Quasar Variability and Accretion Disk Reprocessing Using Latent Stochastic Differential Equations

Quasars are bright active galactic nuclei powered by the accretion of matter around supermassive black holes at the center of galaxies. Their stochastic brightness variability depends on the physical properties of the accretion disk and black hole. The upcoming Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to observe tens of millions of quasars, so there is a need for efficient techniques like machine learning that can handle the large volume of data. Quasar variability is believed to be driven by an X-ray corona, which is reprocessed by the accretion disk and emitted as UV/optical variability. We are the first to introduce an auto-differentiable simulation of the accretion disk and reprocessing. We use the simulation as a direct component of our neural network to jointly model the driving variability and reprocessing, trained with supervised learning on simulated LSST-like 10 yr quasar light curves. We encode the light curves using a transformer encoder, and the driving variability is reconstructed using latent stochastic differential equations, a physically motivated generative deep learning method that can model continuous-time stochastic dynamics. By embedding the physical processes of the driving signal and reprocessing into our network, we achieve a model that is more robust and interpretable. We demonstrate that our model outperforms a Gaussian process regression baseline and can infer accretion disk parameters and time delays between wave bands, even for out-of-distribution driving signals. Our approach provides a powerful framework that can be adapted to solve other inverse problems in multivariate time series.

Fagin, Joshua [City Univ. of New York (CUNY), NY (↗