Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Three-dimensional numerical simulation of near-surface flows over the Martian north polar cap

Measurements made by Viking Lander VL-2 (48 N) have shown that the near-surface wind and temperature regime on Mars displays striking similarities to terrestrial counterparts. The diurnal radiative cycle is responsible for establishment of a well-defined thermal circulation in which downslope (Katabatic) flows prevail during the nighttime hours and weak upslope (anabatic) conditions prevail during the daytime. Previous work has indicated that the slope flows are much like those found on Earth, particularly the Katabatic winds, which show striking similarities to drainage flows observed over Antarctica. The low-level wind regime appears to be an important factor in the scouring of the martian landscape. The north polar cap shows evidence of eolian features such as dunes, frost streaks, and grooves from Viking imagery. The direction of the prevailing wind can in cases be inferred from the eolian features. We examine the thermally induced flows that result from the radiative heating and cooling of the martian north polar region using a comprehensive three-dimensional atmospheric mesoscale numerical model. The same model has been used previously for simulation of Antarctic Katabatic winds. The model equations are written in terrain-following coordinates to allow for irregular terrain; prognostic equations include the flux forms of the horizontal momentum equations, temperature, continuity. A surface energy budget equation is also incorporated in which the surface temperature is determined. Explicit parameterization of both terrestrial (longwave) and solar (shortwave) radiation is included. Turbulent transfer of heat and momentum in the martian atmosphere is assumed to follow the similarity expressions in the surface boundary layer on Earth. The terrain heights for the martian north polar region have been obtained from the U.S. Geological Survey map and digitized onto a 57x57 grid with a spacing of 75 km. The resulting terrain map is shown in Fig. 1. The vertical grid consists of 15 levels (delta equals 0.998, 0.99, 0.98, 0.97, 0.96, 0.94, 0.92, 0.90, 0.85, 0.775, 0.70, 0.60, 0.50, 0.30, 0.10). The high resolution in the lower atmosphere is necessary to capture details of the boundary layer flows. The lowest level corresponds to a height of approximately 20 m above the ground, the seconds level 100 m.

Parish, Thomas R.↗
NASA NTRS1994

Gravitational instabilities in protostellar disks

The nonaxisymmetric stability of self-gravitating, geometrically thick accretion disks has been studied for protostellar systems having a wide range of disk-to-central object mass ratios. Global eigenmodes with four distinctly different characters were identified using numerical, nonlinear hydrodynamic techniques. The mode that appears most likely to arise in normal star formation settings, however, resembles the 'eccentric instability' that was identified earlier in thin, nearly Keplerian disks: It presents an open, one-armed spiral pattern that sweeps continuously in a trailing direction through more than 2-pi radians, smoothly connecting the inner and outer edges of the disk, and requires cooperative motion of the point mass for effective amplification. This particular instability promotes the development of a single, self-gravitating clump of material in orbit about the point mass, so its routine appearance in our simulations supports the conjecture that the eccentric instability provides a primary route to the formation of short-period binaries in protostellar systems.

Tohline, J. E.↗
NASA NTRS1995

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

Mesoscale model simulations provide insight into the complex jet streak adjustments on 11-12 July 1981 that preceded the first of two significant gravity wave events to have been generated over the Rocky Mountains in Montana. Simulations employing a variety of terrain treatments indicate that prior to wave formation, geostrophic adjustment processes modified the structure of the mid-upper tropospheric jet streak by creating secondary jetlets to the southeast of the polar jet streak in proximity to the gravity wave generation region. This simulated restructuring of the mid-upper tropospheric jet streak is the result of a four stage process. During stage 1, the wind adjusts to the mass field as the jet streak exit region propagates into the inflection point between the upstream trough and downstream ridge in the height field. Stage 2 is initiated as the mass field is forced to adjust to the new ageostrophic wind field created during stage 1. Stage 3 is defined by a second geostrophic adjustment process occurring in a similar manner but to the south and east of the adjustment which occurs during stage 1. A low-level mesoscale jetlet is formed during stage 4 in response to the low-level pressure falls that are established during stage 3. The perturbation of this jetlet, caused by orographically-induced adiabatic and diabatic physical processes, is the likely mechanism responsible for the generation of the first and second episode of observed gravity waves. The dynamics responsible for this wave episode are discussed as differential surface sensible heating inducing an orographically-forced mountain-plains solenoid, resulting in the formation of additional mesoscale jetlets and internal gravity waves. Also discussed is how convective latent heating modifies the numerically simulated terrain-induced internal gravity waves, especially their amplitude and phase velocities, which provide better agreement with those wave characteristics observed in nature. Finally, the three-dimensional linear response of a zonally uniform barotropic flow in a vertically unbounded, continuously stratified, Boussinesq atmosphere which is perturbed from geostrophic equilibrium is investigated.

Lin, Yuh-Lang↗
NASA NTRS1997

Electromagnetic Scattering from Arbitrarily Shaped Aperture Backed by Rectangular Cavity Recessed in Infinite Ground Plane

The electromagnetic scattering from an arbitrarily shaped aperture backed by a rectangular cavity recessed in an infinite ground plane is analyzed by the integral equation approach. In this approach, the problem is split into two parts: exterior and interior. The electromagnetic fields in the exterior part are obtained from an equivalent magnetic surface current density assumed to be flowing over the aperture and backed by an infinite ground plane. The electromagnetic fields in the interior part are obtained in terms of rectangular cavity modal expansion functions. The modal amplitudes of cavity modes are determined by enforcing the continuity of the electric field across the aperture. The integral equation with the aperture magnetic current density as an unknown is obtained by enforcing the continuity of magnetic fields across the aperture. The integral equation is then solved for the magnetic current density by the method of moments. The electromagnetic scattering properties of an aperture backed by a rectangular cavity are determined from the magnetic current density. Numerical results on the backscatter radar cross-section (RCS) patterns of rectangular apertures backed by rectangular cavities are compared with earlier published results. Also numerical results on the backscatter RCS patterns of a circular aperture backed by a rectangular cavity are presented.

Cockrell, C. R.↗
NASA NTRS2017

Finite Element Modeling and Analysis of Mars Entry Aeroshell Baseline Concept

The structure that is developed and analyzed in this project must be able to survive all the various load conditions that it will encounter along its course to Mars with the minimal amount of weight and material. At this stage, the goal is to study the capability of the structure using a finite element model (FEM). This FEM is created using a python script, and is numerically solved in Nastran. The purpose of the model is to achieve an optimization of mass given specific constraints on launch and entry. The generation and analysis of the baseline Rigid Mid-Range Lift to Drag Ratio Aeroshell model is a continuation and an improvement on previous work done for the FEM. The model is generated using Python programming with the axisymmetric placement of nodes for beam and shell elements. The shells are assigned a honeycomb sandwich material with an aluminum honeycomb core and composite face sheets, and the beams are assigned the same material as the shell face sheets. There are two load cases assigned to the model: Earth launch and Mars entry. The Earth launch case consists of pressure, gravity, and vibration loads, and the Mars entry case consists of just pressure and gravity loads. The Earth launch case was determined to be the driving case, though the analyses are performed for both cases to ensure the constraints are satisfied. The types of analysis performed with the model are design optimization, statics, buckling, normal modes, and frequency response, the last of which is only for the Earth launch load case. The final results indicated that all of the requirements are satisfied except the thermal limits, which could not yet be tested, and the normal modes for the Mars entry. However, the frequency limits during Mars entry are expected to be much higher than the lower frequency limits set for the analysis. In addition, there are still improvements that can be made in order to reduce the weight while still meeting all requirements.

Ahmed, Samee W.↗
NASA NTRSDate not supplied

Updates and Correlation of EMU System-Level Model (SINDA EMU)

During United States Extravehicular Activity 80 (US EVA 80), water was observed in the helmet of an Extravehicular Mobility Unit (EMU) during cabin repressurization. One of the primary methods of determining the likely cause of this failure was through a comparison of EVA 80 to other historical EVAs using an analytical approach. The Systems Improved Numerical Differential Analysis EMU model (SINDA EMU) is a system-level model of the EMU that was used in this investigation. SINDA EMU was initially developed and correlated to test data in the 1980s. Since its conception, SINDA EMU has been continually adjusted based on new test data and changes to the EMU design. To support the water in the helmet investigation, SINDA EMU needed to be further updated and recorrelated to ensure accurate results. These changes included changing the primary carbon dioxide (CO2) removal technology, implementing logic to allow for re-evaporation of sweat runoff from the liquid cooling and ventilation garment (LCVG), and improving the transient modeling capabilities. To validate the implementation of these adjustments, SINDA EMU was correlated to test data from the 1990s and human-in-the-loop (HITL) testing from 2014. These updates and correlation efforts proved that SINDA EMU is an effective tool for investigating the EVA 80 water in the helmet failure event.

Noah Lial Andersen↗
DOE OSTI2024

Quantum spatial search with multiple excitations

Spatial search is the problem of finding a marked vertex in a graph. A continuous-time quantum walk in the single-excitation subspace of an $n$ spin system solves the problem of spatial search by finding the marked vertex in $O(\sqrt{n})$ time. Here, we investigate a natural extension of the spatial search problem, marking multiple vertices of a graph, which are still marked with local fields. We prove that a continuous-time quantum walk in the $k$-excitation subspace of $n$ spins can determine the binary string of $k$ marked vertices with an asymptotic fidelity in time $O(\sqrt{n})$, despite the size of the state space growing as $O(n^k)$. Numerically, we show that this algorithm can be implemented with interactions that decay as $1/r^\alpha$, where $r$ is the distance between spins, and an $\alpha$ that is readily available in current ion trap systems.

Lewis, Dylan↗
NASA NTRS2001

Global Discrete Artificial Boundary Conditions for Time-Dependent Wave Propagation

We construct global artificial boundary conditions (ABCs) for the numerical simulation of wave processes on unbounded domains using a special non-deteriorating algorithm that has been developed previously for the long-term computation of wave-radiation solutions. The ABCs are obtained directly for the discrete formulation of the problem; in so doing, neither a rational approximation of 'non-reflecting kernels,' nor discretization of the continuous boundary conditions is required. The extent of temporal nonlocality of the new ABCs appears fixed and limited; in addition, the ABCs can handle artificial boundaries of irregular shape on regular grids with no fitting/adaptation needed and no accuracy loss induced. The non-deteriorating algorithm, which is the core of the new ABCs is inherently three-dimensional, it guarantees temporally uniform grid convergence of the solution driven by a continuously operating source on arbitrarily long time intervals, and provides unimprovable linear computational complexity with respect to the grid dimension. The algorithm is based on the presence of lacunae, i.e., aft fronts of the waves, in wave-type solutions in odd-dimension spaces, It can, in fact, be built as a modification on top of any consistent and stable finite-difference scheme, making its grid convergence uniform in time and at the same time keeping the rate of convergence the same as that of the non-modified scheme. In the paper, we delineate the construction of the global lacunae-based ABCs in the framework of a discretized wave equation. The ABCs are obtained for the most general formulation of the problem that involves radiation of waves by moving sources (e.g., radiation of acoustic waves by a maneuvering aircraft). We also present systematic numerical results that corroborate the theoretical design properties of the ABCs' algorithm.

Ryabenkii, V. S.↗
DOE OSTI2020

Quantifying the generalization error in deep learning in terms of data distribution and neural network smoothness

We report the accuracy of deep learning, i.e., deep neural networks, can be characterized by dividing the total error into three main types: approximation error, optimization error, and generalization error. Whereas there are some satisfactory answers to the problems of approximation and optimization, much less is known about the theory of generalization. Most existing theoretical works for generalization fail to explain the performance of neural networks in practice. To derive a meaningful bound, we study the generalization error of neural networks for classification problems in terms of data distribution and neural network smoothness. We introduce the cover complexity (CC) to measure the difficulty of learning a data set and the inverse of the modulus of continuity to quantify neural network smoothness. A quantitative bound for expected accuracy/error is derived by considering both the CC and neural network smoothness. Although most of the analysis is general and not specific to neural networks, we validate our theoretical assumptions and results numerically for neural networks by several data sets of images. The numerical results confirm that the expected error of trained networks scaled with the square root of the number of classes has a linear relationship with respect to the CC. We also observe a clear consistency between test loss and neural network smoothness during the training process. In addition, we demonstrate empirically that the neural network smoothness decreases when the network size increases whereas the smoothness is insensitive to training dataset size.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Stabilizing a strongly nonlinear structure through shaker dynamics in fixed frequency voltage control tests

Bifurcations are commonly encountered during force controlled swept and stepped sine testing of nonlinear structures, which generally leads to the so-called jump-down or jump-up phenomena between stable solutions. There are various experimental closed-loop control algorithms, such as control-based continuation and phase-locked loop, to stabilize dynamical systems through these bifurcations, but they generally rely on specialized control algorithms that are not readily available with many commercial data acquisition software packages. A recent method was developed to experimentally apply sequential continuation using the shaker voltage that can be readily deployed using commercially available software. By utilizing the stabilizing effects of electrodynamic shakers and the force dropout phenomena in fixed frequency voltage control sine tests, this approach has been demonstrated to stabilize the unstable branch of a nonlinear system with three branches, allowing for three multivalued solutions to be identified within a specific frequency bandwidth near resonance. Recent testing on a strongly nonlinear system with vibro-impact nonlinearity has revealed jumping behavior when performing sequential continuation along the voltage parameter, like the jump phenomena seen during more traditional force controlled swept and stepped sine testing. Here, this paper investigates the stabilizing effects of an electrodynamic shaker on strongly nonlinear structures in fixed frequency voltage control tests using both numerical and experimental methods. The harmonic balance method is applied to the coupled shaker-structure system with an electromechanical model to simulate the fixed voltage control tests and predict the stabilization for different parameters of the model. The simulated results are leveraged to inform the design of a set of experiments to demonstrate the stabilization characteristics on a fixture-pylon assembly with a vibro-impact nonlinearity. Through numerical simulation and experimental testing on two different strongly nonlinear systems, the various parameters that influence the stability of the coupled shaker-structure are revealed to better understand the performance of fixed frequency voltage control tests.

42 ENGINEERING↗
NASA NTRS1973

Alternating method applied to edge and surface crack problems.

The alternating method, which intimately combines analytical results with numerical calculations, as applied to edge crack problems in two dimensions and surface crack problems in three dimensions, is treated. The case of a crack perpendicular to the edge of a semiinfinite material is considered. One of the crack geometries that has received continual interest in fracture mechanics is that of a semielliptical crack whose major axis lies on a stress free surface. In order to demonstrate the sensitivity of the solution to the influence of the free surface the semicircular crack problem is again treated by the alternating method.

Hartranft, R. J.↗
NASA NTRS1979

Dynamic properties of elastomer cartridge specimens under a rotating load

This paper presents the results of a program of analysis and test to determine the dynamic properties of elastomer cartridges operating under a rotating load. These measured properties were compared to predictions based on results of unidirectional tests with the same elastomer material. The test method for the dynamic stiffness and damping measurements was essentially the same as the Base Excitation Resonant Mass Method. The primary difference is that the exciting force used for these most recent tests was exerted by rotating unbalance in a rotational test rig rather than a shake table. The specimens tested were: two rectangular cross-section, continuous ring cartridges of different cross-section and three cylindrical button cartridges of different button thickness. Tests were performed for strains from about 0.0001 to about 0.01 (double amplitude). Material properties and prediction equations determined from reciprocating tests were used to make numerical predictions of stiffness, damping, and loss coefficient for the test elements, with encouraging results. Strain was shown to be an important parameter in determining these dynamic properties, particularly damping and loss coefficient.

Darlow, M. S.↗
NASA NTRS1980

Generation of C-type cascade grids for viscous flow computation

A rapid procedure for generating C-type cascade grids suitable for viscous flow computations in turbomachinery blade rows is presented. The resulting mesh is periodic from one blade passage to the next, nearly orthogonal, and continuous across the wake downstream of a blade. The procedure employs a pair of conformal mappings that take the exterior of the cascade into the interior of an infinite strip with curved boundaries. The final transformation to a rectangular computational domain is accomplished numerically. The boundary values are obtained from a panel solution of an integral equation and the interior values by a rapid ADI solution of Laplace's equation. Examples of C-type grids are presented for both compressor and turbine blades and the extension of the procedure to three dimensions is briefly outlined.

Sockol, P. M.↗
NASA NTRS1982

A viscous-inviscid interaction algorithm for three-dimensional turbulent subsonic aerodynamic juncture region flow

An order-of magnitude analysis of the subsonic, three-dimensional, steady time-averaged Navier-Stokes equations for semibounded aerodynamic juncture geometries yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three-dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement first-order continuity effects is accomplished using a penalty differential constraint concept within a finite element solution algorithm. A Reynolds stress constitutive equation, with low turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. The algorithm is employed to predict the three-dimensional turbulent flowfield in the juncture region formed by two intersecting parabolic arcs at a free-stream Mach-number of 0.08, and Re/C = 600,000/m.

Baker, A. J.↗
NASA NTRS1984

Color and Grey Scale in Sonar Displays

In spite of numerous publications 1 it is still rather unclear, whether color is of any help in sonar displays. The work presented here deals with a particular type of sonar data, i.e., LOFAR-grams (low frequency analysing and recording) where acoustic sensor data are continuously written as a time-frequency plot. The question to be answered quantitatively is, whether color coding does improve target detection when compared with a grey scale code. The data show significant differences in receiver-operating characteristics performance for the selected codes. In addition it turned out, that the background noise level affects the performance dramatically for some color codes, while others remain stable or even improve. Generally valid rules are presented on how to generate useful color scales for this particular application.

Kraiss, K. F.↗
NASA NTRS1986

Theoretical study on the effect of the design of small (milli-Newton) thruster jets on molecular contamination for the space station

The self-induced molecular contamination around the space station could have adverse effects on space station components (for example solar panels) as well as scientific experiments that might be done on or near the space station. Aerospace engineers need to design a space station (SS) propulsion system that keeps the SS in a stable orbit and at the same time does not allow the propellant gases to interfere with the experiments of the user. One scenario that might accomplish the above requirements is to use an electrothermal propulsion system, resistojet, that will thrust continuously in the hundreds of milli-Newton range which will provide a constant altitude for the SS with a low g environment. As a first attempt to understand the contamination from such a propulsion system, a point source model was developed. The numerical results of the point source model are given. Number column densities for CO2 are presented as a function of direction of observation (line of sight), temperature of the exit gas, and mean exit velocity. All the results are for a constant exhaust rate of 5,000 kg/year. In addition, a mathematical model to study the effect of nozzle design on the induced molecular environment around the space station produced by simple gas propellants is described. The mathematical model would allow one to follow the expansion of the gas from the throat of a nozzle to the nozzle exit plane and then into the space external to the nozzle.

Riley, B. R.↗
NASA NTRS1988

Numerical simulations of thermal-chemical instabilities at the core-mantle boundary

Numerical simulations of thermal-chemical instabilities in the D-double-prime layer at the core-mantle boundary are presented which show that strong lateral heterogeneities in the composition and density fields can be initiated and maintained dynamically if there is continuous replenishment of material from subduced slabs coming from the upper mantle. These chemical instabilities have a tendency to migrate laterally and may help to support core-mantle boundary topography with short and long wavelengths. The thermal-chemical flows produce a relatively stagnant D-double-prime layer with strong lateral and temporal variations in basal heat flux, which gives rise to thermal core-mantle interactions influencing the geodynamo.

Hansen, Ulrich↗
NASA NTRS1991

Time-domain finite elements in optimal control with application to launch-vehicle guidance

A time-domain finite element method is developed for optimal control problems. The theory derived is general enough to handle a large class of problems including optimal control problems that are continuous in the states and controls, problems with discontinuities in the states and/or system equations, problems with control inequality constraints, problems with state inequality constraints, or problems involving any combination of the above. The theory is developed in such a way that no numerical quadrature is necessary regardless of the degree of nonlinearity in the equations. Also, the same shape functions may be employed for every problem because all strong boundary conditions are transformed into natural or weak boundary conditions. In addition, the resulting nonlinear algebraic equations are very sparse. Use of sparse matrix solvers allows for the rapid and accurate solution of very difficult optimization problems. The formulation is applied to launch-vehicle trajectory optimization problems, and results show that real-time optimal guidance is realizable with this method. Finally, a general problem solving environment is created for solving a large class of optimal control problems. The algorithm uses both FORTRAN and a symbolic computation program to solve problems with a minimum of user interaction. The use of symbolic computation eliminates the need for user-written subroutines which greatly reduces the setup time for solving problems.

Bless, Robert R.↗