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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

A method for extracting meridonal winds from ionosonde measurements by using ionospheric models

There has been great progress in modelling and measuring the dynamics of the neutral upper atmosphere in recent years. However, future progress will depend on the availability of global measurements of neutral winds. Attention is drawn to a relatively cheap means of supplementing the data base of neutral winds provided by radar and optical measurements with data obtained by ionosondes. Rishbeth in his review of F-region dynamics, derived the relationship between the height of the F2 layer and the component of the neutral wind parallel to the magnetic field of the Earth. The sensitivity is examined of the height and density of the F2 layer over Boulder, Colorado on 30 July 1982 to changes in meridional wind speed using a comprehensive interhemispheric numerical model that solves the continuity and momentum equations for H+ and O+, the energy equations for Te and Ti, and the 2-stream photoelectron equations to obtain electron heating rates. For the neutral atmosphere temperature and densities, the mass spectrometer incoherent scatter radar model of Hedin was used.

Richards, P. G.↗
NASA NTRS1986

Aspects of the determination of winds by means of scatterometry and of the utilization of vector wind data for meteorological forecasts

The present paper provides a description of four aspects of scatterometer winds and their uses. The theory of wave generation by the wind is considered along with an analysis of the properties of superobservations, and studies of intermittent versus continuous data assimilation methods for numerical weather predictions which use remotely sensed data. A comparison of the sum of squares versus the maximum likelihood method for recovering the vector winds is also conducted. Questions regarding wind speed, friction velocity, or normal stress are discussed and synoptic scale fields from Seasat-SASS data are examined.

Pierson, W. J., Jr.↗
NASA NTRS1988

A mathematical model of the middle and high latitude ionosphere

A time-dependent three-dimensional model of the middle and high latitude ionosphere is described. The density distributions of six ion species NO(+), O(2+), N(2+), O(+), N(+), He(+), and the electron and ion temperatures are obtained from a numerical solution of the appropriate continuity, momentum, and energy equations. The equations are solved as a function of height for an inclined magnetic field at E and F region altitudes. The three-dimensional nature of the model is obtained by following flux tubes of plasma as they convect or corotate through a moving neutral atmosphere. The model takes account of field-aligned diffusion, cross-field electrodynamic drifts, thermospheric winds, polar wind escape, energy-dependent chemical reactions, neutral composition changes, ion production due to solar EUV radiation and auroral precipitation, thermal conduction, diffusion-thermal heat flow and local heating and cooling processes. The model also takes account of the offset between the geomagnetic and geographic poles.

Schunk, R. W.↗
NASA NTRS1990

A new class of solutions for interstellar magnetohydrodynamic shock waves

An analysis is presented of the equations of motion for steady MHD shock waves proopagating in interstellar clouds, for boundary conditions that preclude C shocks. In addition to J shocks, in which the neutral fluid component becomes subsonic at an adiabatic jump front, the equations admit a new class of solutions, called C-asterisk shocks, in which the transition to subsonic flow occurs continuously at a sonic point. Numerical methods are developed for computing the structure of J and C-asterisk shocks propagating in diffuse interstellar clouds. The effects of chemical, ionization, and recombination processes are included in this treatment. An alternative numerical method, which uses artificial viscosity to facilitate integration through sonic points, is analyzed and shown to be invalid. A set of exemplary solutions, computed for realistic shock parameters, shows that C-asterisk shocks occur for a broad range of conditions relevant to diffuse interstellar clouds.

Roberge, W. G.↗
NASA NTRS1994

Penetration experiments in aluminum and Teflon targets of widely variable thickness

The morphologies and detailed dimensions of hypervelocity craters and penetration holes on space-exposed surfaces faithfully reflect the initial impact conditions. However, current understanding of this postmortem evidence and its relation to such first-order parameters as impact velocity or projectile size and mass is incomplete. While considerable progress is being made in the numerical simulation of impact events, continued impact simulations in the laboratory are needed to obtain empirical constraints and insights. This contribution summarizes such experiments with Al and Teflon targets that were carried out in order to provide a better understanding of the crater and penetration holes reported from the Solar Maximum Mission (SMM) and the Long Duration Exposure Facility (LDEF) satellites. A 5-mm light gas gun was used to fire spherical soda-lime glass projectiles from 50 to 3175 microns in diameter (D(sub P)), at a nominal 6 km/s, into Al (1100 series; annealed) and Teflon (Teflon(sup TFE)) targets. Targets ranged in thickness (T) from infinite halfspace targets (T approx. equals cm) to ultrathin foils (T approx. equals micron), yielding up to 3 degrees of magnitude variation in absolute and relative (D(sub P)/T) target thickness. This experimental matrix simulates the wide range in D(sub P)/T experienced by a space-exposed membrane of constant T that is being impacted by projectiles of widely varying sizes.

Hoerz, F.↗
NASA NTRS2000

The Axisymmetric Pulsar Magnetosphere

We present the structure of the axisymmetric force-free magnetosphere of an aligned rotating magnetic dipole, in the case in which there exists a sufficiently large charge density (whose origin we do not question) to satisfy the ideal MHD condition, E (raised dot) B = 0, everywhere. The unique distribution of electric current along the open magnetic field lines which is required for the solution to be continuous and smooth is obtained numerically. We expect that our solution will be useful as the starting point for detailed studies of pulsar magnetospheres under more general conditions, namely when either the force-free and/or the ideal MHD condition E (raised dot) B = 0 are not valid in the entire magnetosphere. Based on our solution, we consider that the most likely positions of such an occurrence are the polar cap, the crossings of the zero space charge surface by open field lines, and the return current boundary, but not the light cylinder.

Contopoulos, J.↗
NASA NTRS2015

Investigating the Effects of Temperature on the Signatures of Shocks Propagated Through Impacts into Minerals Found in Comets and Asteroids

Comets and asteroids are subjected to extremely cold conditions throughout their lifetimes. During their sojourns in the solar system, they are subjected to collisions at speeds that are easily capable of generating shock waves in their constituent materials. In addition to ices, more common silicate minerals such as olivines and pyroxenes are important components of these objects. The collision-induced shocks could affect the spectral signatures of those mineral components, which could in turn be detected telescopically. We have embarked on a project to determine how impact-generated shock might affect the reflectance spectra and structures of select silicates as both impact speed and target temperature are varied systematically. While the effects of impact speed (in the form of shock stress) on numerous materials have been and continue to be studied, the role of target temperature has received comparatively little attention, presumably because of the operational difficulties it can introduce to experimentation. Our experiments were performed with the vertical gun in the Experimental Impact Laboratory of the Johnson Space Center. A liquid-nitrogen system was plumbed to permit cooling of the target container and its contents under vacuum to temperatures as low as -100 C (173 K). Temperatures were monitored by thermocouples mounted on the outside of the target container. Because those sensors were not in contact with the target material at impact, the measured temperatures are treated as lower limits for the actual values. Peridot (Mg-rich olivine) and enstatite (Mg-rich orthopyroxene) were used as targets, which involved the impact of alumina (Al2O3) spheres at speeds of 2.0 - 2.7 km s(exp -1) and temperatures covering 25 C to -100 C (298 K to 173 K). We have begun collecting and analyzing data in the near to mid-IR with a Fourier-transform infrared spectrometer, and preliminary analyses show that notable differences in absorption-band strength and position occur as functions of both impact speed (peak shock stress) and initial temperature.

Lederer, Susan M.↗
NASA NTRS2002

Multidisciplinary Optimization of a Transport Aircraft Wing using Particle Swarm Optimization

The purpose of this paper is to demonstrate the application of particle swarm optimization to a realistic multidisciplinary optimization test problem. The paper's new contributions to multidisciplinary optimization is the application of a new algorithm for dealing with the unique challenges associated with multidisciplinary optimization problems, and recommendations as to the utility of the algorithm in future multidisciplinary optimization applications. The selected example is a bi-level optimization problem that demonstrates severe numerical noise and has a combination of continuous and truly discrete design variables. The use of traditional gradient-based optimization algorithms is thus not practical. The numerical results presented indicate that the particle swarm optimization algorithm is able to reliably find the optimum design for the problem presented here. The algorithm is capable of dealing with the unique challenges posed by multidisciplinary optimization as well as the numerical noise and truly discrete variables present in the current example problem.

Sobieszczanski-Sobieski, Jaroslaw↗
DOE OSTI2026

High-Fidelity, Low-Dissipation/Symmetry-Preserving Numerical Scheme for Solving the Euler Equations with Unstructured, Metric-Based Mesh Adaptation

This work presents an overview of a high-fidelity compressible Euler solver that utilizes the continuous Galerkin (CG) method with added artificial numerical diffusion for stabilization to solve a variety of unsteady and steady benchmark inviscid flow problems. This work shows that discretizing the Euler equations with this CG approach and first order basis functions produces a cost-effective stencil as well as simple well-posed boundary conditions. We show through convergence testing with manufactured solutions that the reduced stencil of CG, combined with the low amount of artificial diffusion required when using the stabilization method outlined in this work, leads to stable and highly accurate results for a variety of unsteady and steady applications. When combined with the adaptive mesh refinement approach used for many of the cases in this work, our results show that the flow solver achieves even more accurate results. A variety of inviscid flow cases are presented in this work, including transient 2D cases with complex shock structures and several steady 3D airfoils sections with a constant span.

Doetsch, Kevin [ORNL] (ORCID:0000000267051705)↗
NASA NTRS1995

Solar Thermal Concept Evaluation

Concentrated solar thermal energy can be utilized in a variety of high temperature applications for both terrestrial and space environments. In each application, knowledge of the collector and absorber's heat exchange interaction is required. To understand this coupled mechanism, various concentrator types and geometries, as well as, their relationship to the physical absorber mechanics were investigated. To conduct experimental tests various parts of a 5,000 watt, thermal concentrator, facility were made and evaluated. This was in anticipation at a larger NASA facility proposed for construction. Although much of the work centered on solar thermal propulsion for an upper stage (less than one pound thrust range), the information generated and the facility's capabilities are applicable to material processing, power generation and similar uses. The numerical calculations used to design the laboratory mirror and the procedure for evaluating other solar collectors are presented here. The mirror design is based on a hexagonal faceted system, which uses a spherical approximation to the parabolic surface. The work began with a few two dimensional estimates and continued with a full, three dimensional, numerical algorithm written in FORTRAN code. This was compared to a full geometry, ray trace program, BEAM 4, which optimizes the curvatures, based on purely optical considerations. Founded on numerical results, the characteristics of a faceted concentrator were construed. The numerical methodologies themselves were evaluated and categorized. As a result, the three-dimensional FORTRAN code was the method chosen to construct the mirrors, due to its overall accuracy and superior results to the ray trace program. This information is being used to fabricate and subsequently, laser map the actual mirror surfaces. Evaluation of concentrator mirrors, thermal applications and scaling the results of the 10 foot diameter mirror to a much larger concentrator, were studied. Evaluations, recommendations and pit falls regarding the structure, materials and facility design are presented.

Hawk, Clark W.↗
DOE OSTI2022

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

42 ENGINEERING↗
NASA NTRS2003

Recent Advances In Science Support For Isolated Droplet Combustion Experiments

In a joint program involving Prof. F.A. Williams of the University of California, San Diego and Dr. V. Nayagam of the National Center for Microgravity Research, the combustion characteristics of isolated liquid fuel droplets of n-heptane, n-decane, methanol, methanol-water, ethanol and ethanol-water having initial diameters between about 1 mm and 6 mm continues to be investigated. The objectives of the work are to improve fundamental knowledge of droplet combustion dynamics for pure fuels and fuel-water mixtures through microgravity experiments and theoretical analyses. The Princeton contributions support the engineering design, data analysis, and data interpretation requirements for the study of initially single component, spherically symmetric, isolated droplet combustion studies through experiments and numerical modeling. UCSD contributions are described in a companion communication in this conference. The Princeton effort also addresses the analyses of Fiber Supported Droplet Combustion (FSDC) experiments conducted with the above fuels and collaborative work with others who are investigating droplet combustion in the presence of steady convection. A thorough interpretation of droplet burning behavior for n-heptane and n-decane over a relatively wide range of conditions also involves the influences of sooting on the combustion behavior, and this particular aspect on isolated burning of droplets is under consideration in a collaborative program underway with Drexel University. This collaboration is addressed in another communication at this conference. The one-dimensional, time-dependent, numerical modeling approach that we have continued to evolve for analyzing isolated, quiescent droplet combustion data has been further applied to investigate several facets of isolated droplet burning of simple alcohols, n-heptane, and n-decane. Some of the new results are described below.

Dryer, F. L.↗
NASA NTRS1979

Effects of sounding temperature assimilation on weather forecasting - Model dependence studies

In comparing various methods for the assimilation of remote sounding information into numerical weather prediction (NWP) models, the problem of model dependence for the different results obtained becomes important. The paper investigates two aspects of the model dependence question: (1) the effect of increasing horizontal resolution within a given model on the assimilation of sounding data, and (2) the effect of using two entirely different models with the same assimilation method and sounding data. Tentative conclusions reached are: first, that model improvement as exemplified by increased resolution, can act in the same direction as judicious 4-D assimilation of remote sounding information, to improve 2-3 day numerical weather forecasts. Second, that the time continuous 4-D methods developed at GLAS have similar beneficial effects when used in the assimilation of remote sounding information into NWP models with very different numerical and physical characteristics.

Ghil, M.↗
NASA NTRS1988

A theoretical framework for convergence and continuous dependence of estimates in inverse problems for distributed parameter systems

Numerical techniques for parameter identification in distributed-parameter systems are developed analytically. A general convergence and stability framework (for continuous dependence on observations) is derived for first-order systems on the basis of (1) a weak formulation in terms of sesquilinear forms and (2) the resolvent convergence form of the Trotter-Kato approximation. The extension of this framework to second-order systems is considered.

Banks, H. T.↗
NASA NTRS2001

Dual-Mode Combustion of Hydrogen in a Mach 5, Continuous-Flow Facility

Results of an experimental and numerical study of a dual-mode scramjet combustor are reported. The experiment consisted of a direct-connect test of a Mach 2 hydrogen-air combustor with a single unswept-ramp fuel injector. The flow stagnation enthalpy simulated a flight Mach number of 5. Measurements were obtained using conventional wall instrumentation and a particle-imaging laser diagnostic technique. The particle imaging was enabled through the development of a new apparatus for seeding fine silicon dioxide particles into the combustor fuel stream. Numerical simulations of the combustor were performed using the GASP code. The modeling, and much of the experimental work, focused on the supersonic combustion mode. Reasonable agreement was observed between experimental and numerical wall pressure distributions. However, the numerical model was unable to predict accurately the effects of combustion on the fuel plume size, penetration, shape, and axial growth.

Goyne, C. P.↗
NASA NTRS1996

Mixing and Transition Control Studied

Considerable progress in understanding nonlinear phenomena in both unbounded and wallbounded shear flow transition has been made through the use of a combination of high- Reynolds-number asymptotic and numerical methods. The objective of this continuing work is to fully understand the nonlinear dynamics so that ultimately (1) an effective means of mixing and transition control can be developed and (2) the source terms in the aeroacoustic noise problem can be modeled more accurately. Two important aspects of the work are that (1) the disturbances evolve from strictly linear instability waves on weakly nonparallel mean flows so that the proper upstream conditions are applied in the nonlinear or wave-interaction streamwise region and (2) the asymptotic formulations lead to parabolic problems so that the question of proper out-flow boundary conditions--still a research issue for direct numerical simulations of convectively unstable shear flows--does not arise. Composite expansion techniques are used to obtain solutions that account for both mean-flow-evolution and nonlinear effects. A previously derived theory for the amplitude evolution of a two-dimensional instability wave in an incompressible mixing layer (which is in quantitative agreement with available experimental data for the first nonlinear saturation stage for a plane-jet shear layer, a circular-jet shear layer, and a mixing layer behind a splitter plate) have been extended to include a wave-interaction stage with a three-dimensional subharmonic. The ultimate wave interaction effects can either give rise to explosive growth or an equilibrium solution, both of which are intimately associated with the nonlinear self-interaction of the three dimensional component. The extended theory is being evaluated numerically. In contrast to the mixing-layer situation, earlier comparisons of theoretical predictions based on asymptotic methods and experiments in wall-bounded shear-flow transition have been somewhat lacking in one aspect or another. The current work strongly suggests that the main weakness is the underlying asymptotic representation of the linear "part" of the problem and not the explicit modeling of the nonlinear/wave-interaction effects. Consequently, the long-wave-length/high-Reynolds-number asymptotic limit for the Blasius boundary-layer stability problem was reexamined, and a new dispersion relationship for the instability waves that is uniformly valid for both the upper- and lower branch regions to the required order of approximation was obtained. A comparison with numerical results, obtained by solving the Orr-Sommerfeld stability problem, shows that the asymptotic formula provides surprisingly good results, even for values of the frequency parameter usually encountered in experimental investigations. This is particularly evident in the dynamically important upper-branch region, where much of the nonlinear interactions in transition experiments are believed to take place. The result is important in that it can be used to greatly improve the accuracy of weakly nonlinear critical-layer-based theories, and a consistent nonlinear theory is currently under evaluation.

Source record↗