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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Semi-Analytical Solution Approach for Solving Constant-Coefficient First-Order Partial Differential Equations

Simulation and control of many dynamic systems involve solving partial differential equations (PDE). This letter proposes a semi-analytical solution (SAS) approach for fast and high-quality solution of first-order PDEs. The region of interest of the studied PDE is divided into a grid, and an SAS is derived for each grid cell in the form of the multivariate polynomials, of which the coefficients are identified using initial value and boundary value conditions. The solutions are solved in a “time-stepping” manner, i.e. within one time step, the coefficients of the SAS are identified and the initial value of the next time step is evaluated. This approach achieves a significantly larger grid cell than the widely used finite difference method, and thus enhances the computational efficiency significantly. Furthermore, the simulation result on the natural gas pipeline model demonstrates the advantages of SAS in accuracy and computational efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Quantum Electromagnetic Transients Program

This work devises a quantum electromagnetic transients program (QEMTP) which is the first attempt to tackle the computational challenges in solving EMTP through quantum computing. The main contributions lie in: (1) A quantum-enabled EMTP formulation with Dommel's model encoded in quantum states; 2) An Harrow-Hassidim-Lloyd (HHL)-based QEMTP to solve the discrete-time nodal equations; 3) An iteration-based HHL revision for mitigating temporal errors to achieve high accuracy of QEMTP with limited quantum resources. Case studies verify the correctness and efficacy of QEMTP in simulating both fast and slow dynamics.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Design of a Nanophotonic Hybrid Coupler

This project aims to design and simulate a nanophotonic hybrid coupler for detecting vacuum ultraviolet (VUV) light with higher efficiency. VUV light is of significant interest because of its applications ranging from high energy physics to space science and electronic industry. However, current VUV photodetectors have very low efficiencies because most materials absorb in the VUV region (100-200 nm). To address this challenge, I have designed a plasmonic-dielectric coupler with wavelength shifting properties, enabling high Purcell enhancement and far-field confinement. The coupler consists of a cylindrical cSi resonator on top of a thin ZnO layer, and Ag and SiO2 are used as substrates. Simulations are done using the finite-difference time domain method to solve Maxwell’s equations, and specifically with Meep and Lumerical. The effectiveness of the coupler is characterized by the Purcell factor and far-field radiation. The Purcell factor reaches a maximum near 390 nm and is in the o rder of 10^4. Furthermore, the coupler also allows partial control of the directionality of the far-field radiation. With the presence of the resonator, the radiation is directed towards the coupled metalens. Overall, our design combines the advantages of plasmonic and dielectric couplers, along with the intrinsic wavelength shifting properties of materials to enhance VUV detection significantly, and such a device can be used in, for example, the DUNE project.

Zhou, Jiayang↗
DOE OSTI2024

On a Simplified Approach to Achieve Parallel Performance and Portability Across CPU and GPU Architectures

This paper presents software advances to easily exploit computer architectures consisting of a multi-core CPU and CPU+GPU to accelerate diverse types of high-performance computing (HPC) applications using a single code implementation. The paper describes and demonstrates the performance of the open-source C++ matrix and array (MATAR) library that uniquely offers: (1) a straightforward syntax for programming productivity, (2) usable data structures for data-oriented programming (DOP) for performance, and (3) a simple interface to the open-source C++ Kokkos library for portability and memory management across CPUs and GPUs. The portability across architectures with a single code implementation is achieved by automatically switching between diverse fine-grained parallelism backends (e.g., CUDA, HIP, OpenMP, pthreads, etc.) at compile time. The MATAR library solves many longstanding challenges associated with easily writing software that can run in parallel on any computer architecture. This work benefits projects seeking to write new C++ codes while also addressing the challenges of quickly making existing Fortran codes performant and portable over modern computer architectures with minimal syntactical changes from Fortran to C++. We demonstrate the feasibility of readily writing new C++ codes and modernizing existing codes with MATAR to be performant, parallel, and portable across diverse computer architectures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Subphotospheric Emission from Short Gamma-Ray Bursts: Protons Mold the Multimessenger Signals

Abstract The origin of the observed Band-like photon spectrum in short gamma-ray bursts (sGRBs) is a long-standing mystery. We carry out the first general relativistic magnetohydrodynamic simulation of an sGRB jet with initial magnetizationσ 0 = 150 in dynamical ejecta from a binary merger. From this simulation, we identify regions along the jet of efficient energy dissipation due to magnetic reconnection and collisionless subshocks. Taking into account electron and proton acceleration processes, we solve for the first time the coupled transport equations for photons, electrons, protons, neutrinos, and intermediate particle species up to close to the photosphere (i.e., up to 1 × 10 12 cm), accounting for all relevant radiative and cooling processes. We find that the subphotospheric multimessenger signals carry strong signatures of the hadronic interactions and their resulting particle cascades. Importantly, the spectral energy distribution of photons is significantly distorted with respect to the Wien one, commonly assumed below the photosphere. Our findings suggest that the bulk of the nonthermal photon spectrum observed in sGRBs can stem from hadronic processes occurring below the photosphere and previously neglected, with an accompanying energy flux of neutrinos peaking in the GeV energy range.

Astronomy & Astrophysics↗
NASA NTRS1971

The pressure and deformation profile between two colliding lubricated cylinders

The pressure and deformation profiles between two colliding lubricated cylinders are obtained by solving the coupled, time-dependent elastohydrodynamic equations with an iterative procedure. The analysis includes effects which were not considered in a previous solution, namely, the effect of the lubricant compressibility and the effect of a lubricant with composite pressure-viscosity coefficients. It is found that the local approach velocity plays an important role during final stages of normal approach. It causes the lubricant to be entrapped within the contact region; neither the pressure nor the deformation profile converges to the Hertzian profile for a dry contact. The use of a smaller pressure-viscosity coefficient at high pressures reduces the sharp pressure gradient at the center of the contact and produces a much milder variation of load with respect to the film thickness. The effect of compressibility of the lubricant is found to be relatively small.

Lee, K.↗
NASA NTRS1974

Significant applications of ERTS-1 data to resource management activities at the state level in Ohio

Described are techniques utilized and the progress made in applying ERTS-1 data to (1) detecting, inventorying, and monitoring surface mining activities, particularly in relation to recently passed strip mine legislation in Ohio; (2) updating current land use maps at various scales for multiagency usage, and (3) solving other real-time problems existing throughout the various Ohio governmental agencies. General conclusions regarding current user views as to the opportunities and limitations of operationally using ERTS-1 data at the state level are also noted.

Sweet, D. C.↗
NASA NTRS1974

Viscous reacting flows with wall slip and catalysis applied to spheres in arc jets and flight

The influence of wall slip and catalytic atom-recombination on the flow field and wall heat flux are calculated for high altitude flight and arc jet flow conditions. Boundary equations, which include velocity slip, temperature jump, and wall catalytic atom recombination, are coupled to the viscous reacting multicomponent Navier-Stokes equations. These equations are solved using a time-dependent finite difference technique applied to spheres in an arc jet flow (Reynolds number of 550) and a high altitude flight case representative of the Space Shuttle Orbiter (Reynolds number of 450). The results indicate that catalysis strongly influences the temperature jump, but not the velocity slip. Slip increases the atom fraction and temperature at both the wall and the flow field. Likewise, the shock stand-off distance, the wall heat flux, and friction coefficient are increased over the nonslip cases. The reacting gas calculations confirm the chemically frozen nature of the shock layer in arc jet flows.

Scott, C. D.↗
NASA NTRS1974

Quasi-one-dimensional solution for the power of CO2 gasdynamic lasers

A theory has been developed for calculating the laser power of CO2 gasdynamic lasers. A three-mode vibrational gas model was assumed, and the coupled rate and one-dimensional conservation equations were solved by the time-dependent technique. The results were compared with experimental data, and the agreement was good.

Lee, G.↗
NASA NTRS1975

Modal analysis and dynamic stresses for acoustically excited shuttle insulation tiles

Improvements and extensions to the RESIST computer program developed for determining the normalized modal stress response of shuttle insulation tiles are described. The new version of RESIST can accommodate primary structure panels with closed-cell stringers, in addition to the capability for treating open-cell stringers. In addition, the present version of RESIST numerically solves vibration problems several times faster than its predecessor. A new digital computer program, titled ARREST (Acoustic Response of Reusable Shuttle Tiles) is also described. Starting with modal information contained on output tapes from RESIST computer runs, ARREST determines RMS stresses, deflections and accelerations of shuttle panels with reusable surface insulation tiles. Both programs are applicable to stringer stiffened structural panels with or without reusable surface insulation titles.

Ojalvo, I. U.↗
NASA NTRS1976

Numerical simulation of supersonic and hypersonic turbulent compression corner flows using relaxation models

Relaxation turbulence eddy viscosity models are incorporated to solve the complete Navier-Stokes equations for supersonic and hypersonic flows over a two-dimensional compression corner. The system of equations is solved by a time-split, second-order accurate numerical scheme. Details of relaxation process are studied, and several relaxation models are tested and compared. Good improvement in the prediction of upstream pressure propagation is obtained for a Mach number of 2.96, and Reynolds number of 10,000,000, with wedge angle of 25 deg. However, the application of the relaxation models to hypersonic flow at Mach number 8.66 and Reynolds number of 22,000,000, with highly cooled wall, shows unfavorable effects on heat transfer and skin friction in the reattachment and recompression regions.

Hung, C. M.↗
NASA NTRS1976

Moisture transport in composites

The moisture distribution and moisture content of composite materials exposed to humid air or to water are discussed. Expressions are presented for estimating the weight gain of both single and multilayered composites when the moisture content and the temperature of the environment are constant. The parameters needed for solving problems in the time varying environmental conditions are specified, and a numerical approach for the solution of such problems is indicated.

Springer, G. S.↗
NASA NTRS1977

Implicit approximate-factorization schemes for the efficient solution of steady transonic flow problems

Implicit approximate-factorization algorithms (AF) are developed for the solution of steady-state transonic flow problems. The performance of the AF solution method is evaluated relative to that of the standard solution method for transonic flow problems, successive line over-relaxation (SLOR). Both methods are applied to the solution of the nonlinear, two-dimensional transonic small-disturbance equation. Results indicate that the AF method requires substantially less computer time than SLOR to solve the nonlinear finite-difference matrix equation for a transonic flow field. This increase in computational efficiency is achieved with no appreciable increase in computer storage or coding complexity.

Ballhaus, W. F.↗
NASA NTRS1977

Combustion response modeling for composite solid propellants

A computerized mathematical model of the combustion response function of composite solid propellants was developed with particular attention to the contributions of the solid phase heterogeneity. The one-dimensional model treats the solid phase as alternating layers of ammonium perchlorate and binder, with an exothermic melt layer at the surface. Solution of the Fourier heat equation in the solid provides temperature and heat flux distributions with space and time. The problem is solved by conserving the heat flux at the surface from that produced by a suitable model of the gas phase. An approximation of the BDP flame model is utilized to represent the gas phase. By the use of several reasonable assumptions, it is found that a significant portion of the problem can be solved in closed form. A method is presented by which the model can be applied to tetramodal particle size distributions. A computerized steady-state version of the model was completed, which served to validate the various approximations and lay a foundation for the combustion response modeling. The combustion response modeling was completed in a form which does not require an iterative solution, and some preliminary results were acquired.

Source record↗
NASA NTRS1979

High-altitude gas releases - Transition from collisionless flow to diffusive flow in a nonuniform atmosphere

The paper is motivated by a need for describing high-altitude gas releases from rockets or the Space Shuttle. In the tenuous upper atmosphere, the injected gases expand from a collisionless to a collision-dominated state. The analysis presented extends the efforts of Baum (1973-1974) to include the effects of a nonuniform background atmosphere and a velocity-dependent collision frequency. The unsteady expansion of gas releases is analyzed using gas kinetic theory. The Boltzmann equation with the Krook collision integral is solved numerically. At late times (after many collisions with the background atmosphere) the solution is identical to the one given by diffusion in a nonuniform atmosphere (Bernhardt, 1979). The theoretical model predicts elongation and heating of the vapor trails due to collisions.

Bernhardt, P. A.↗
NASA NTRS1980

A new method for optimizing multiple flyby trajectories

A new procedure has been developed which minimizes total impulsive Delta V for multiple flyby trajectories with constraints on flyby parameters and maneuver times. The method involves solving a bounds-constrained parameter optimization problem with a Newton algorithm utilizing analytic first and second derivatives. Each trajectory segment connecting consecutive maneuver points is found by first targeting from the preceding maneuver point to the parameters of the upcoming flyby and then propagating the resulting trajectory to the next maneuver point. Multi-conic techniques are used for trajectory propagation and for computation of the state transition matrix. This procedure has successfully optimized Galileo satellite tours containing up to 11 flybys.

Damario, L. A.↗
NASA NTRS1981

An implicit method for the calculation of inlet flow fields

Inlet flow fields are calculated by an implicit, time marching procedure to solve the thin layer Navier-Stokes equations formulated in body fitted coordinates. Because the method can be used for a flow field with both subsonic and supersonic regions, it is applicable to subcritical as well as supercritical inlet operation. Results are presented and discussed for an inlet of current design practice. Results include inviscid calculations performed for supercritical inlet operation with uniform and nonuniform inflow boundary conditions as well as for subcritical inlet operation with uniform inflow boundary conditions. Results for viscous calculations performed for supercritical inlet operation with uniform inflow boundary conditions are also discussed.

Biringen, S.↗