Small-signal current transients in insulators with traps.
Time-dependent current flow equations solved for insulator with traps if injected space charge produces negligible perturbation on external electric field
SPACE CHARGE↗
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.
Solidification of a binary mixture
The time dependent concentration and temperature profiles of a finite layer of a binary mixture are investigated during solidification. The coupled time dependent Stefan problem is solved numerically using an implicit finite differencing algorithm with the method of lines. Specifically, the temporal operator is approximated via an implicit finite difference operator resulting in a coupled set of ordinary differential equations for the spatial distribution of the temperature and concentration for each time. Since the resulting differential equations set form a boundary value problem with matching conditions at an unknown spatial point, the method of invariant imbedding is used for its solution.
Antar, B. N.↗
Dynamic Restructuring Of Problems In Artificial Intelligence
"Dynamic tradeoff evaluation" (DTE) denotes proposed method and procedure for restructuring problem-solving strategies in artificial intelligence to satisfy need for timely responses to changing conditions. Detects situations in which optimal problem-solving strategies cannot be pursued because of real-time constraints, and effects tradeoffs among nonoptimal strategies in such way to minimize adverse effects upon performance of system.
Schwuttke, Ursula M.↗
Pseudo-time methods for constrained optimization problems governed by PDE
In this paper we present a novel method for solving optimization problems governed by partial differential equations. Existing methods are gradient information in marching toward the minimum, where the constrained PDE is solved once (sometimes only approximately) per each optimization step. Such methods can be viewed as a marching techniques on the intersection of the state and costate hypersurfaces while improving the residuals of the design equations per each iteration. In contrast, the method presented here march on the design hypersurface and at each iteration improve the residuals of the state and costate equations. The new method is usually much less expensive per iteration step since, in most problems of practical interest, the design equation involves much less unknowns that that of either the state or costate equations. Convergence is shown using energy estimates for the evolution equations governing the iterative process. Numerical tests show that the new method allows the solution of the optimization problem in a cost of solving the analysis problems just a few times, independent of the number of design parameters. The method can be applied using single grid iterations as well as with multigrid solvers.
Taasan, Shlomo↗
Computation of a Single Airfoil Gust Response and Gust-Cascade Interaction Using the CE/SE Method
The problems 1 and 2 in Category 3 are solved using the space-time conservation element and solution element (CE/SE) method. Problem 1 concerns the acoustic field generated by the interaction of a harmonic vertical gust with a single isolated airfoil. Problem 2 models rotor-stator interaction in a 2D cascade. Both problems involve complex geometries and flow physics including vortex shielding and acoustic radiation. An unstructured triangular mesh is used to solve both problems. For problem 2, the Giles approach is incorporated with the CE/SE method to handle non-equal pitches of the rotor and stator. Numerical solution of both near and far fields of problem 1 are presented and compared with a frequency-domain solver GUST3D and a time-domain high-order Discontinuous Spectra Element Method (DSEM) solutions. For problem 2, numerical solutions on the blade surface, inlet and outlet planes are presented.
Wang, Xiao-Yen↗
A combined algorithm for minimum time slewing of flexible spacecraft
The use of Pontryagin's Maximum Principle for the large-angle slewing of large flexible structures usually results in the so-called two-point boundary-value problem (TPBVP), in which many requirements (e.g., minimum time, small flexible amplitude, and limited control powers, etc.) must be satisfied simultaneously. The successful solution of this problem depends largely on the use of an efficient numerical computational algorithm. There are many candidate algorithms available for his problem (e.g., quasilinearization, gradient, and shooting, etc.). In this paper, a proposed algorithm, which combines the quasilinearization method with a time shortening technique and a shooting method, is applied to the minimum-time, three-dimensional, and large-angle maneuver of flexible spacecraft, particularly the orbiting Spacecraft Control Laboratory Experiment (SCOLE) configuration. Theoretically, the nonlinear TPBVP can be solved only through the shooting method to find the 'exact' switching times for the bang-bang controls. However, computationally, a suitable guess for the missing initial costates is crucial because the convergence range of the unknown initial costates is usually narrow, especially for systems with high dimensions and when a multi-bang-bang control strategy is needed. On the other hand, the problems of near minimum time attitude maneuver of general rigid spacecraft and fast slewing of flexible spacecraft have been examined by the authors through a numerical approach based on the quasilinearization algorithm with a time shortening technique. Computational results have demonstrated its broad convergence range and insensitivity to initial costate choices. Consequently, a combined approach is naturally suggested here to solve the minimum time slewing problem. That is, in the computational process, the quasilinearization method is used first to obtain a near minimum time solution. Then, the acquired converged initial costates from the quasilinearization approach are transformed (tailored) to and used as the initial costate guess for starting the shooting method. Finally, the shooting method takes over the remaining calculations until the minimum-time solution converges. The nonlinear equations of motion of the SCOLE are formulated by using Lagrange's equations, with the mast modeled as a continuous beam subject to three-dimensional deformations. The numerical results will be presented and some related computational issues will also be discussed.
Bainum, P. M.↗
Descending Layer Variability Over Arecibo
Descending layers of ionization over Arecibo exhibit very diverse behavior from night to night that does not appear to be strongly correlated to geomagnetic activity, solar forcing, or average semidiurnal tidal winds. On some nights, three or more distinct layers are observed to form near 170 km over timescales of approx. 2 hours. Rather than descending smoothly over periods of several hours, these layers stall, abruptly disappear, or even reverse direction in the midst of their descent. The time scales for their disappearance are examined and compared to loss rates arising from diffusion and recombination. Diffusion alone is found to be too slow to account for the observations, but recombination is fast enough provided that the convergent wind shear that forms the layer is relatively weak coincident with their disappearance. The continuity equation is solved in conjunction with a time sequence of radar profiles to estimate the vertical drift and horizontal neutral wind consistent with the observed behavior. The resultant wind field is northward, has an average speed of approx. 80 m/ s, and varies significantly near the altitude where the layers are observed. These inferred winds are consistent with the presence of the observed layers, and their magnitudes as obtained from the classical continuity and momentum equations are reasonable for this altitude range.
Earle, G. D.↗
Living Color Frame System: PC graphics tool for data visualization
Living Color Frame System (LCFS) is a personal computer software tool for generating real-time graphics applications. It is highly applicable for a wide range of data visualization in virtual environment applications. Engineers often use computer graphics to enhance the interpretation of data under observation. These graphics become more complicated when 'run time' animations are required, such as found in many typical modern artificial intelligence and expert systems. Living Color Frame System solves many of these real-time graphics problems.
Truong, Long V.↗
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↗
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↗
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.↗
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.↗
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.↗
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.↗
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.↗
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.↗
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.↗
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↗