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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Non-cancellation of the parity anomaly in the strong-field regime of QED 2+1

Quantum fluctuations lead to an anomalous violation of parity symmetry in quantum electrodynamics for an even number of spatial dimensions. While the leading parity-odd electric current vanishes in vacuum, we uncover a non-cancellation of the anomaly for strong electric fields with distinct macroscopic signatures. We perform real-time lattice simulations with fully dynamical gauge fields and Wilson fermions in 2 + 1 space-time dimensions. In the static field limit, relevant at early times, we solve the problem analytically. Our results point out the fundamental role of quantum anomalies for strong-field phenomena, relevant for a wide range of condensed matter and high-energy applications, but also for the next generation of gauge theory quantum simulators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Optimizing Individualized Treatment Planning for Parkinson’s Disease Using Deep Reinforcement Learning

More than one million people currently live with Parkinson's Disease (PD) in the U.S. alone. Medications, such as levodopa, can help manage PD symptoms. However, medication treatment planning is generally based on patient history and limited interaction between physicians and patients during office visits. This limits the extent of benefit that may be derived from the treatment as disease/patient characteristics are generally non-stationary. Wearable sensors that provide continuous monitoring of various symptoms, such as bradykinesia and dyskinesia, can enhance symptom management. However, using such data to overhaul the current static medication treatment planning approach and prescribe personalized medication timing and dosage that accounts for patient/care-giver/physician feedback/preferences remains an open question. We develop a model to prescribe timing and dosage of medications, given the motor fluctuation data collected using wearable sensors in real-time. We solve the resulting model using deep reinforcement learning (DRL). The prescribed policy determines the optimal treatment plan that minimizes patient's symptoms. Our results show that the model-prescribed policy outperforms the static a priori treatment plan in improving patients' symptoms, providing a proof-of-concept that DRL can augment medical decision making for treatment planning of chronic disease patients.

Watts, Jeremy↗

DOE OSTI2023

DT-HYDRO

The software solves the time dependent, one dimensional (1D) coupled mass and momentum balance equations governing the elastic flow of water through the penstock, turbine and draft tube in a hydroelectric facility using a high order finite volume based method. The numerical method is based on the Kurganov-Tadmor central method paired with the Monotonic Upstream-centered Scheme for Conservation Laws (MUSCL). This solution method accurately resolves the fast transient behavior of the flow, including water hammer. Additionally, the software estimates the full 3D flow field within the turbine chamber in real time, a feat that is made possible by leveraging pre-computed CFD results by utilizing a reduced order modeling method based on an efficient randomized singular value decomposition (SVD) driven proper orthogonal decomposition (POD) with POD-mode weight regression. The reduced order model of the 3D flow is directly coupled to the 1D elastic flow model so the entire flow field through the penstock and turbine system is resolved quickly and with high fidelity.

Gurecky, William [Oak Ridge National Laboratory (O↗

DOE OSTI2023

Superior discretizations and AMG solvers for extremely anisotropic diffusion via hyperbolic operators [Slides]

Diffusion in magnetic confinement fusion is extremely anisotropic in the direction of field lines. Rewrote diffusion system based on directional gradients, apply discretization and solver techniques developed for advection. Orders of magnitude decrease in error and solve wallclock time vs. traditional methods. The next steps include: (1) incorporate into larger MHD simulations, (2) better solvers for closed field lines or mixed regimes, and (3) possibly other extremely anisotropic equations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Learning to Branch with Interpretable Machine Learning Models

This presentation describes an algorithm for applying machine learning to branching to speed up the solution of integer optimization problems. These problems are challenging and solved multiple times a day by power systems operators. We show that our approach speeds up a widely used open-source optimization solver.

Bayramoglu, Selin↗

NASA NTRS1974

Skylab thruster attitude control system

Preflight activities and the Skylab mission support effort for the thruster attitude control system (TACS) are documented. The preflight activities include a description of problems and their solutions encountered in the development, qualification, and flight checkout test programs. Mission support effort is presented as it relates to system performance assessment, real-time problem solving, flight anomalies, and the daily system evaluation. Finally, the detailed flight evaluation is presented for each phase of the mission using system telemetry data. Data assert that the TACS met or exceeded design requirements and fulfilled its assigned mission objectives.

Wilmer, G. E., Jr.↗

NASA NTRS1974

Lateral conduction effects on heat-transfer data obtained with the phase-change paint technique

A computerized tool, CAPE, (Conduction Analysis Program using Eigenvalues) has been developed to account for lateral heat conduction in wind tunnel models in the data reduction of the phase-change paint technique. The tool also accounts for the effects of finite thickness (thin wings) and surface curvature. A special reduction procedure using just one time of melt is also possible on leading edges. A novel iterative numerical scheme was used, with discretized spatial coordinates but analytic integration in time, to solve the inverse conduction problem involved in the data reduction. A yes-no chart is provided which tells the test engineer when various corrections are large enough so that CAPE should be used. The accuracy of the phase-change paint technique in the presence of finite thickness and lateral conduction is also investigated.

Maise, G.↗

NASA NTRS1975

Multiple-compartment venting program

Computer program solves time-dependent energy and state equations for gas reservoirs, using the solutions of conductor conservation equations as mass and energy rate changes to reservoirs.

Leblanc, L. P.↗

NASA NTRS1976

A block iterative finite element algorithm for numerical solution of the steady-state, compressible Navier-Stokes equations

An iterative method for numerically solving the time independent Navier-Stokes equations for viscous compressible flows is presented. The method is based upon partial application of the Gauss-Seidel principle in block form to the systems of nonlinear algebraic equations which arise in construction of finite element (Galerkin) models approximating solutions of fluid dynamic problems. The C deg-cubic element on triangles is employed for function approximation. Computational results for a free shear flow at Re = 1,000 indicate significant achievement of economy in iterative convergence rate over finite element and finite difference models which employ the customary time dependent equations and asymptotic time marching procedure to steady solution. Numerical results are in excellent agreement with those obtained for the same test problem employing time marching finite element and finite difference solution techniques.

Cooke, C. H.↗

NASA NTRS1978

Critical issues in viscous flow computations

In developing computer programs to numerically solve the Navier-Stokes equations, the purpose of the computation must be clearly kept in mind. In the Air Force, the purpose is to provide design information on non-linear aerodynamic phenomenon for aircraft that perform throughout the flight corridor. This translates into the requirement for a computer program which can solve the time averaged compressible Navier-Stokes equations (with a turbulence model) in three dimensions for generalized geometries. The intended application of the results then controls the priorities in addressing critical issues. Recurrent problem areas encountered in the study of viscous flow include: (1) grid generation for arbitrary geometry; (2) numerical difficulties; (3) turbulence models; (4) accuracy and efficiency; and (5) smearing of discontinuities.

Hankey, W. L., Jr.↗

NASA NTRS1977

Nitramine smokeless propellant research

A transient ballistics and combustion model was derived to represent the closed vessel experiment that is widely used to characterize propellants. The model incorporates the nitramine combustion mechanisms. A computer program was developed to solve the time dependent equations, and was applied to explain aspects of closed vessel behavior. It is found that the rate of pressurization in the closed vessel is insufficient at pressures of interest to augment the burning rate by time dependent processes. Series of T-burner experiments were performed to compare the combustion instability characteristics of nitramine (HMX) containing propellants and ammonium perchlorate (AP) propellants. It is found that the inclusion of HMX consistently renders the propellant more stable.

Source record↗

NASA NTRS1977

A block iterative finite element algorithm for numerical solution of the steady-state, compressible Navier-Stokes equations

An iterative method for numerically solving the time independent Navier-Stokes equations for viscous compressible flows is presented. The method is based upon partial application of the Gauss-Seidel principle in block form to the systems of the nonlinear algebraic equations which arise in construction of finite element (Galerkin) models approximating solutions of fluid dynamic problems. The continuous cubic element on triangles is employed for function approximation. Computational results for a free shear flow at Re = 1000 indicate significant achievement of economy in iterative convergence rate over finite element and finite difference models which employ the customary time dependent equations and symptotic time marching procedure to steady solution. Numerical results are in excellent agreement with those obtained for the same test problem employing time marching finite element and finite difference solution techniques.

Cooke, C. H.↗

NASA NTRS1980

Viscous compressible flow about blunt bodies using a numerically generated orthogonal coordinate system

A numerical solution to the Navier-Stokes equations was obtained for blunt axisymmetric entry bodies of arbitrary shape in supersonic flow. These equations are solved on a finite difference mesh obtained from a simple numerical technique which generates orthogonal coordinates between arbitrary boundaries. The governing equations are solved in time dependent form using Stetter's improved stability three step predictor corrector method. For the present application, the metric coefficients were obtained numerically using fourth order accurate, finite difference relations and proved to be totally reliable for the highly stretched mesh used to resolve the thin viscous boundary layer. Solutions are obtained for a range of blunt body nose shapes including concavities.

Graves, R. A., Jr.↗

NASA NTRS1981

Tropospheric chemistry - A global perspective

Processes that affect the concentration of tropospheric OH are considered, taking into account linkages between the chemistry of OH and the chemistry of H, HO2, and H2O2. A review is presented of observational data of special relevance to OH, notably the distributions of CO, CH4, O3, H2O, HNO3, NO, and NO2. Most of the results presented in connection with the description of the model were obtained by solving time-dependent continuity equations. Rates for photolytic processes were allowed to vary diurnally with insolation. It is found that a well-calibrated and reliable model for OH places important constraints on global budgets for a variety of gases including CH4, H2, and CH3Cl in addition to CO. Accurate measurements of species such as CH3CCl3 can provide valuable checks on global models that must integrate over a variety of atmospheric conditions. However, emission rates for the relevant gases must be adequately quantified.

Logan, J. A.↗

NASA NTRS1983

Magnetohydrodynamic modeling of coronal structure and expansion

The presence of a magnetic field in the corona adds structure to the solar wind and almost certainly plays an important role in the energetics of the flow. Analytical and numerical modeling of gas-magnetic field interactions as used to compute steady, global flow are discussed. The approach used in, and results from a recent global model (Steinolfson, Suess and Wu, 1982) are discussed. Ideas on the most effective ways to improve the physical content and numerical efficiency of these models are outlined. Solutions of the MHD equations are discussed only in order to find steady-state flows, even though this often entails solving time-dependent equations.

Suess, S. T.↗

NASA NTRS1985

A computational study of saddle point separation and horseshoe vortex system

Incompressible flow around a cylinder-end wall junction has been simulated by solving the incompressible Navier-Stokes equations in three dimensions. The equations, cast in generalized curvilinear coordinates, are solved in time as a hyperbolic system by adding a pressure term in the continuity equation and are marched to a steady state. Various physical quantities associated with the saddle point of separation and the horseshoe vortex system are calculated. Computational and experimental results are generally consistent. The skin friction and the pressure distribution on the end wall are consistent with the physics of the problem. Secondary flows both in front of the cylinder and behind it are predicted that are in qualitative agreement with flow visualization results. The calculations also indicate a strongly nonuniform pressure loading along the length of the cylinder. A new mechanism for the existence of the recirculation bubbles behind the cylinder-end wall with relatively low ratio of cylinder height to the approaching boundary layer thickness is observed which is markedly different from its two-dimensional counterpart.

Kaul, U. K.↗

NASA NTRS1987

Discontinuous solutions to hyperbolic systems under operator splitting

Two-dimensional systems of linear hyperbolic equations are studied with regard to their behavior under a solution strategy that in alternate time-steps solves exactly the component one-dimensional operators. The initial data is a step function across an oblique discontinuity. The manner in which this discontinuity breaks up under repeated applications of the split operator is analyzed, and it is shown that the split solution will fail to match the true solution in any case where the two operators do not share all their eigenvectors. The special case of the fluid flow equations is analyzed in more detail, and it is shown that arbitrary initial data gives rise to pseudo acoustic waves and a non-physical stationary wave. The implications of these findings for the design of high-resolution computing schemes are discussed.

Roe, P. L.↗

NASA NTRS1988

Newton solution of inviscid and viscous problems

The application of Newton iteration to inviscid and viscous airfoil calculations is examined. Spatial discretization is performed using upwind differences with split fluxes. The system of linear equations which arises as a result of linearization in time is solved directly using either a banded matrix solver or a sparse matrix solver. In the latter case, the solver is used in conjunction with the nested dissection strategy, whose implementation for airfoil calculations is discussed. The boundary conditions are also implemented in a fully implicit manner, thus yielding quadratic convergence. Complexities such as the ordering of cell nodes and the use of a far field vortex to correct freestream for a lifting airfoil are addressed. Various methods to accelerate convergence and improve computational efficiency while using Newton iteration are discussed. Results are presented for inviscid, transonic nonlifting and lifting airfoils and also for laminar viscous cases.

Venkatakrishnan, V.↗