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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

The role of the sea-surface temperature distribution on numerically simulated cyclogenesis during ERICA

The goal was to quantify the extent to which a sea surface temperature (SST) front can influence cyclogenesis. The approach was to use the Drexel Limited-Area Mesoscale Prediction System (LAMPS) dynamical model to simulate cyclogenesis over various SST fields. Research during the past year focused on the development and testing of a four dimensional data assimilation (FDDA) technique within LAMPS. The technique is a continuous dynamical assimilation where forcing terms are added to the governing model equations to gradually nudge the model solution toward a gridded analysis. Here, the nudging is used as a dynamic initialization tool during a 12 hour preforecast to generate model balanced initial conditions for a subsequent 24 hour numerical prediction. Tests were performed to determine which variables to nudge and how to specify the four dimensional weighting function used to scale the nudging terms. To date, optimal results were obtained by nudging the u and v components of the wind along with the potential temperature. The weighting function ranged from 0 to 1 and varies in time as a quadratic polynomial. It was initialized at 0, reached its maximum at 9 hours into the preforecast, and fell back at 0 to 12 hours. The nudging terms are included in the model equations for all grid points except those within the model predicted oceanic boundary layer. This design attempts to confine changes imposed by the specified SST field to the oceanic boundary layer during the preforecast period.

Lapenta, William M.↗
NASA NTRS2021

Plumelets: Dynamic Filamentary Structures in Solar Coronal Plumes

Solar coronal plumes long seemed to possess a simple geometry supporting spatially coherent, stable outflow without significant fine structure. Recent high-resolution observations have challenged this picture by revealing numerous transient, small-scale, collimated out-flows (“jetlets”) at the base of plumes. The dynamic filamentary structure of solar plumes above these outflows, and its relationship with the overall plume structure, have remained largely unexplored. We analyzed the statistics of continuously observed fine structure in-side a single representative bright plume within a mid-latitude coronal hole during 2016 July2-3. By applying advanced edge-enhancement and spatiotemporal analysis techniques to ex-tended series of high-resolution images from the Solar Dynamics Observatory’s Atmospheric Imaging Assembly, we determined that the plume was composed of numerous time-evolving filamentary substructures, referred to as “plumelets” in this paper, that accounted for most of the plume emission. The number of simultaneously identifiable plumelets was positively correlated with plume brightness, peaked in the fully formed plume, and remained saturated thereafter. The plumelets had transverse widths of 10 Mm and intermittently supported upwardly propagating periodic disturbances with phase speeds of 190-260 km s−1and longitudinal wavelengths of 55-65 Mm. The characteristic frequency (3.5 mHz) is commensurate with that of solar p-modes. Oscillations in neighboring plumelets are uncorrelated, indicating that the waves could be driven by p-mode flows at spatial scales smaller than the plumelet separation. Multiple independent sources of outflow within a single coronal plume should impart significant fine structure to the solar wind that may be detectable by Parker Solar Probe and Solar Orbiter.

V M Uritsky↗
NASA NTRS1993

Tetravalent chromium (Cr(4+)) as laser-active ion for tunable solid-state lasers

Major accomplishments under NASA grant NAG-1-1346 are summarized. (1) numerical modeling of the four mirror astigmatically compensated, Z-fold cavity was performed and several design parameters to be used for the construction of a femtosecond forsterite laser were revealed by simulation. (2) femtosecond pulses from a continuous wave mode-locked chromium doped forsterite laser were generated. The forsterite laser was actively mode-locked using an acousto-optic modulator operating at 78 MHz with two Brewster high dispersion glass prisms for intra-cavity chirp compensation. Transform-limited sub-100-fs pulses were routinely generated in the TEM(sub 00) mode with 85 mW of continuous power tunable over 1230-1280 nm. The shortest pulses of 60-fs pulsewidth were measured. (3) Self-mode-locked operation of the Cr:forsterite laser was achieved. Synchronous pumping was used to mode lock the forsterite laser resulting in picosecond pulses, which in turn provided the starting mechanism for self-mode-locking. The pulses generated had an FWHM of 105 fs and were tunable between 1230-1270 nm. (4) Numerical calculations indicated that the pair of SF 14 prisms used in the cavity compensated for quadratic phase but introduced a large cubic phase term. Further calculations of other optical glasses indicated that a pair of SFN 64 prisms can introduce the same amount of quadratic phase as SF 14 prisms but introduce a smaller cubic phase. When the SF 14 prisms were replaced by SFN 64 prisms the pulsewidth was reduced to 50 fs. Great improvement was observed in the stability of the self mode-locked forsterite laser and in the ease of achieving mode locking. Using the same experimental arrangement and a new forsterite crystal with improved FOM the pulse width was reduced to 36 fs.

Seas, A.↗
NASA NTRS1971

Non-asymptotic effects in the approach to the far field sonic boom.

An extension is described for the Whitham theory (1952, 1954) which indicates small deviations from the ultimate N-wave for bodies with a continuous or discontinuous tangent. Some simple formulas are derived which yield a good reduction in the sonic boom intensity. The theory is elucidated by some numerical examples for a class of simple theoretically relevant shapes.

Henderson, L. F.↗
NASA NTRS1985

Putting the 'super' in supercomputers

Computers used for numerical simulations of physical phenomena, e.g., flowfields, meteorology, structural analysis, etc., replace physical experiments that are too expensive or impossible to perform. The problems considered continually become increasingly more complex and thus demand faster processing times to do all necessary computations. The effects components technologies have on computer speed are leveling off, leaving new architectures and programming as the only currently viable means to upgrade speed. Parallel computations, either in the form of array processors, assembly line processing or multiprocessors are being explored using existing microprocessor technologies. Slower hardware configurations can also be made equivalent to faster supercomputers by economic programming. The availability of rudimentary parallel architecture supercomputers for general industrial use is increasing. Scientific applications continue to drive the development of more sophisticated parallel machines.

Schulbach, C.↗
NASA NTRS1987

Lanczos modes for reduced-order control of flexible structures

Lanczos mode models represent low-frequency forced response better than do normal mode models and can be developed for both continuous and finite element structural representations. It was recommended that Lanczos mode models for systems with multiple input and/or rigid body modes should be developed; numerical stability of the Lanczos algorithm should be assessed; and control system designs employing the Lanczos mode models should be attempted.

Craig, Roy R., Jr.↗
NASA NTRS1987

The dynamics and control of large flexible space structures X, part 1

The effect of delay in the control system input on the stability of a continuously acting controller which is designed without considering the delay is studied. The stability analysis of a second order plant is studied analytically and verified numerically. For this example it is found that the system becomes unstable for a delay which is equivalent to only 16 percent of its natural period of motion. It is also observed that even a small amount of natural damping in the system can increase the amount of delay that can be tolerated before the onset of instability. The delay problem is formulated in the discrete time domain and an analysis procedure suggested. The maximum principle from optimal control theory is applied to minimize the time required for the slewing of a general rigid spacecraft. The slewing motion need not be restricted to a single axis maneuver. The minimum slewing time is calculated based on a quasi-linearization algorithm for the resulting two point boundary value problem. Numerical examples based on the rigidized in-orbit model of the SCOLE also include the more general reflector line-of-sight slewing maneuvers.

Bainum, Peter M.↗
NASA NTRS1986

Geomagnetic temporal change: 1903-1982 - A spline representation

The secular variation of the earth's magnetic field is itself subject to temporal variations. These are investigated with the aid of the coefficients of a series of spherical harmonic models of secular variation deduced from data for the interval 1903-1982 from the worldwide network of magnetic observatories. For some studies it is convenient to approximate the time variation of the spherical harmonic coefficients with a smooth, continuous, function; for this a spline fitting is used. The phenomena that are investigated include periodicities, discontinuities, and correlation with the length of day. The numerical data presented will be of use for further investigations and for the synthesis of secular variation at any place and at any time within the interval of the data - they are not appropriate for temporal extrapolations.

Langel, R. A.↗
NASA NTRS1990

Integrated system identification and modal state estimation for control of flexible space structures

A novel approach of integrated system identification and modal state estimation is proposed for control of linear dynamical systems including flexible space structures. There are four steps involved in this approach. First, the relation between a stochastic state space model of a dynamical system and the coefficients of its autoregressive model with exogenous input is derived. Second, an adaptive least-squares transversal predictor is used to estimate the coefficients of the model. Third, a state space model and a steady state Kalman filter gain of the dynamical system are then identified from the coefficients of the model by using the eigensystem realization algorithm. Fourth, a modal state estimator is constructed using the modal parameters of the identified model. On-line implementation of this algorithm can continually improve the modal parameters and the filter gain. It can also gradually update the system model when the system characteristics are slowly changing. A numerical example is used to illustrate the feasibility of the new approach.

Chen, Chung-Wen↗
NASA NTRS1994

Towards enhancing and delaying disturbances in free shear flows

The family of shear flows comprising the jet, wake, and the mixing layer are subjected to perturbations in an inviscid incompressible fluid. By modeling the basic mean flows as parallel with piecewise linear variations for the velocities, complete and general solutions to the linearized equations of motion can be obtained in closed form as functions of all space variables and time when posed as an initial value problem. The results show that there is a continuous as well as the discrete spectrum that is more familiar in stability theory and therefore there can be both algebraic and exponential growth of disturbances in time. These bases make it feasible to consider control of such flows. To this end, the possibility of enhancing the disturbances in the mixing layer and delaying the onset in the jet and wake is investigated. It is found that growth of perturbations can be delayed to a considerable degree for the jet and the wake but, by comparison, cannot be enhanced in the mixing layer. By using moving coordinates, a method for demonstrating the predominant early and long time behavior of disturbances in these flows is given for continuous velocity profiles. It is shown that the early time transients are always algebraic whereas the asymptotic limit is that of an exponential normal mode. Numerical treatment of the new governing equations confirm the conclusions reached by use of the piecewise linear basic models. Although not pursued here, feedback mechanisms designed for control of the flow could be devised using the results of this work.

Criminale, W. O.↗
NASA NTRS1998

Prediction of Transitional Flows in the Low Pressure Turbine

Current turbulence models tend to give too early and too short a length of flow transition to turbulence, and hence fail to predict flow separation induced by the adverse pressure gradients and streamline flow curvatures. Our discussion will focus on the development and validation of transition models. The baseline data for model comparisons are the T3 series, which include a range of free-stream turbulence intensity and cover zero-pressure gradient to aft-loaded turbine pressure gradient flows. The method will be based on the conditioned N-S equations and a transport equation for the intermittency factor. First, several of the most popular 2-equation models in predicting flow transition are examined: k-e [Launder-Sharina], k-w [Wilcox], Lien-Leschiziner and SST [Menter] models. All models fail to predict the onset and the length of transition, even for the simplest flat plate with zero-pressure gradient(T3A). Although the predicted onset position of transition can be varied by providing different inlet turbulent energy dissipation rates, the appropriate inlet conditions for turbulence quantities should be adjusted to match the decay of the free-stream turbulence. Arguably, one may adjust the low-Reynolds-number part of the model to predict transition. This approach has so far not been very successful. However, we have found that the low-Reynolds-number model of Launder and Sharma [1974], which is an improved version of Jones and Launder [1972] gave the best overall performance. The Launder and Sharma model was designed to capture flow re-laminarization (a reverse of flow transition), but tends to give rise to a too early and too fast transition in comparison with the physical transition. The three test cases were for flows with zero pressure gradient but with different free-stream turbulent intensities. The same can be said about the model when considering flows subject to pressure gradient(T3C1). To capture the effects of transition using existing turbulence models, one approach is to make use of the concept of the intermittency to predict the flow transition. It was originally based on the intermittency distribution of Narasimha [1957], and then gradually evolved into a transport equation for the intermittency factor. Gostelow and associates [1994,1995] have made some improvements to Narasimha's method in an attempt to account for both favorable and adverse pressure gradients. Their approach is based on a linear, explicit combination of laminar and turbulent solutions. This approach fails to predict the overshoot of the skin friction on a flat plate near the end of transition zone, even though the length of transition is well predicted. The major flaw of Gostelow's approach is that it assumes the non-turbulent part being the laminar solution and the turbulent part being the turbulent solution and they do not interact across the transitional region. The technique in condition averaging the flow equations in intermittent flows was first introduced by Libby [1975] and Dopazo [1977] and further refined by Dick and associates [1988, 1996]. This approach employs two sets of transport equations for the non-turbulent part and the other for the turbulent part. The advantage of this approach is that it allows the interaction of non-turbulent and turbulent velocities through the introduction of additional source terms in the continuity and momentum equations for the non-turbulent and turbulent velocities. However, the strong coupling of the two sets of equations has caused some numerical difficulties, which requires special attention. The prediction of the skin friction can be improved by this approach via the implicit coupling of non-turbulent and turbulent velocity flelds. Another improvement of the interrmittency model can be further made by allowing the intermittency to vary in the cross-stream direction. This is one step prior to testing any proposal for the transport equation for the intermittency factor. Instead of solving the transport equation for the intermittency factor, the distribution for the intermittency factor is prescribed by Klebanoff's empirical formula [1955]. The skin friction is very well predicted by this new modification, including the overshoot of the profile near the end of the transition zone. The outcome of this study is very encouraging since it indicates that the proper description of the intermittency distribution is the key to the success of the model prediction. This study will be used to guide us on the modelling of the intermittency transport equation.

Huang, George↗
NASA NTRS2004

Characteristics of Electromagnetic Pulse Propagation in Metal

It is well known that the solution of the diffusion equation for an electromagnetic field with a time harmonic term, e(sup iwt), is in the form of a traveling wave whose amplitude attenuates over distance into a conducting medium. As the attenuation is an increasing function of frequency, the high frequency components attenuate more rapidly than those of low ones upon entering a well conducting object. At the same time, the phase velocity of an individual component is also an increasing function of frequency causing a broadening of the pulse traveling inside a conductor. In the results of our previous study of numerical simulations, the problem of using a gaussian input pulse was immediately clear. First, having the dominant frequency components distributed around zero, the movement of the peak was not well defined. Second, with the amplitude of fourier components varying slowly over a wide range, the dispersion-induced blurring of the peak position was seen to be severe. For the present study, we have used a gaussian modulated single frequency sinusoidal wave, i. e., the carrier, as an input pulse in an effort to improve the issues related to the unclear movement of peak and dispersion as described above. This was based on the following two anticipated advantages: First, the packet moves in a conductor at the group velocity calculated at the carrier frequency, which means it is well controllable. Second, the amplitude of frequency components other than that of the carrier can be almost negligible, such that the effect of dispersion can be significantly reduced. A series of experiments of transmitting electromagnetic pulses through aluminum plates of various thickness was performed to test the validity of the above points. The results of numerical simulation based on wave propagation are discussed with respect to the experimental results. Finally, a simple simulation was performed based on diffusion of a continuous sine wave input and the results are compared with those of a single frequency sinusoidal wave observed over time at difference locations inside a conductor.

Namkung, M.↗
NASA NTRS2008

LETS: Lunar Environments Test System

The Environmental Effects Branch (EM50) at the Marshall Space Flight Center has developed a unique capability within the agency, namely the Lunar Environment Test System (LETS). LETS is a cryo-pumped vacuum chamber facility capable of high vacuum (10-7 Torr). LETS is a cylindrical chamber, 30 in. (0.8 m) diameter by 48 in. (1.2 m) long thermally controlled vacuum system. The chamber is equipped with a full array of radiation sources including vacuum ultraviolet, electron, and proton radiation. The unique feature of LETS is that it contains a large lunar simulant bed (18 in. x 40 in. x 6 in.) holding 75 kg of JSC-1a simulant while operating at a vacuum of 10-7 Torr. This facility allows three applications: 1) to study the charging, levitation and migration of dust particles, 2) to simulate the radiation environment on the lunar surface, and 3) to electrically charge the lunar simulant enhancing the attraction and adhesion of dust particles to test articles more closely simulating the lunar surface dust environment. LETS has numerous diagnostic instruments including TREK electrostatic probes, residual gas analyzer (RGA), temperature controlled quartz crystal microbalance (TQCM), and particle imaging velocimeter (PIV). Finally, LETS uses continuous Labview data acquisition for computer monitoring and system control.

Vaughn, Jason A.↗
NASA NTRSDate not supplied

Development of a C-ELS Specimen-Based Numerical Benchmark for Mode II Delamination and Assessment of Two VCCT-Based Propagation Strategies

A finite element (FE) benchmark example inspired by the calibrated end-loaded split (C-ELS) specimen is developed and used to assess the performance of delamination propagation capabilities based on linear elastic fracture mechanics (LEFM). The C-ELS specimen has the advantage of a longer region of stable delamination propagation compared to the existing mode II benchmark case. The new benchmark example may therefore provide a better assessment tool by enabling more stable crack growth in regions further away from the boundary conditions or load application. First, a benchmark result is created manually using two-dimensional finite element models of the C-ELS specimen with different delamination lengths. Second, the performance of the virtual crack closure technique (VCCT) delamination propagation capabilities in the Abaqus/Standard®1 FE code and the recently developed Progressive Release eXplicit-VCCT (PRX-VCCT) method are assessed by comparing the results to the benchmark case. Two examples with different starter delamination lengths are studied. A shorter starter length is chosen to create a scenario with unstable delamination propagation. A longer delamination encourages stable delamination propagation. Detailed results from three-dimensional analyses with aligned and misaligned meshes and two levels of mesh refinement are provided. In general, good agreement can be achieved between the results obtained from the quasi-static propagation analysis and the benchmark analysis. Numerical artifacts including anomalous unreleased nodes in the crack wake and zig-zag crack fronts occur for propagation analyses using Abaqus/Standard VCCT. In comparison, continuous, smooth, delamination fronts are observed for PRX-VCCT. The use of the benchmark case to assess different VCCT-based propagation strategies illustrates the value of establishing benchmark cases.

Composite Materials↗
NASA NTRS2023

Universal Two-Phase Convection Heat Transfer Correlations for Cryogenic Pipe Chilldown

This paper presents a set of universal two-phase convection heat transfer correlations for modeling boiling heat transfer during pipe chilldown fit over the widest available range of cryogenic fluids and thermodynamic conditions. The correlations improve upon prior correlations that were developed separately for liquid nitrogen (LN2) and liquid hydrogen (LH2) pipe quenching datasets. The new correlations include equations to calculate the single-phase vapor heat transfer, film boiling heat transfer, transition boiling heat transfer, nucleate boiling heat transfer, single-phase liquid heat transfer, bulk vapor temperature during high quality film boiling, Leidenfrost temperature, critical heat flux, critical heat flux temperature, and the onset of nucleate boiling temperature. The correlations were validated against LH2, LN2, liquid methane, liquid oxygen, and liquid argon pipe quenching datasets. The eight datasets cover the following parameter ranges: pipe lengths of 0.1 to 6.5m; outer pipe diameters of 12.7 to 25.4 mm; pipe wall thicknesses from 0.51 to 1.64 mm; flow directions of upward, downward, and horizontal; gravity levels of 1g and 0g±0.01g. A numerical model assumes homogeneous mixing between the vapor and liquid and implicitly integrates the coupled energy equations for the pipe and fluid, as well as continuity for the fluid. The model was used to estimate the equilibrium quality and fluid mass accumulation along the pipes during chilldown from estimates of the fluid-pipe heat transfer extracted from temperature measurements. The correlations can be implemented in lumped parameter codes such as SINDA/FLUINT and the Generalized Fluid System Simulation Program (GFSSP) to improve the accuracy of chilldown time and chilldown boiloff mass predictions. Such predictions are useful for designing ground or in-space cryogenic liquid transfer systems.

Universal Correlations↗
NASA NTRS1988

The development of flux-split algorithms for flows with non-equilibrium thermodynamics and chemical reactions

A finite-volume method for the numerical computation of flows with nonequilibrium thermodynamics and chemistry is presented. A thermodynamic model is described which simplifies the coupling between the chemistry and thermodynamics and also results in the retention of the homogeneity property of the Euler equations (including all the species continuity and vibrational energy conservation equations). Flux-splitting procedures are developed for the fully coupled equations involving fluid dynamics, chemical production and thermodynamic relaxation processes. New forms of flux-vector split and flux-difference split algorithms are embodied in a fully coupled, implicit, large-block structure, including all the species conservation and energy production equations. Several numerical examples are presented, including high-temperature shock tube and nozzle flows. The methodology is compared to other existing techniques, including spectral and central-differenced procedures, and favorable comparisons are shown regarding accuracy, shock-capturing and convergence rates.

Grossman, B.↗
NASA NTRS1995

Implementation of a state of the art automated system for the production of cloud/water vapor motion winds from geostationary satellites

The research objectives in this proposal were part of a continuing program at UW-CIMSS to develop and refine an automated geostationary satellite winds processing system which can be utilized in both research and operational environments. The majority of the originally proposed tasks were successfully accomplished, and in some cases the progress exceeded the original goals. Much of the research and development supported by this grant resulted in upgrades and modifications to the existing automated satellite winds tracking algorithm. These modifications were put to the test through case study demonstrations and numerical model impact studies. After being successfully demonstrated, the modifications and upgrades were implemented into the NESDIS algorithms in Washington DC, and have become part of the operational support. A major focus of the research supported under this grant attended to the continued development of water vapor tracked winds from geostationary observations. The fully automated UW-CIMSS tracking algorithm has been tuned to provide complete upper-tropospheric coverage from this data source, with data set quality close to that of operational cloud motion winds. Multispectral water vapor observations were collected and processed from several different geostationary satellites. The tracking and quality control algorithms were tuned and refined based on ground-truth comparisons and case studies involving impact on numerical model analyses and forecasts. The results have shown the water vapor motion winds are of good quality, complement the cloud motion wind data, and can have a positive impact in NWP on many meteorological scales.

Velden, Christopher↗
NASA NTRS1984

Numerical modeling of laser thermal propulsion flows

An review of the problems associated with modeling laser thermal propulsion flows, a synopsis of the status of such models, and the attributes of a successful model are presented. The continuous gaseous hydrogen laser-supported combustion wave (LSCW) thruster, for which a high-energy laser system (preferably space-based) should exist by the time the propulsion technology is developed, is considered in particular. The model proposed by Raizer (1970) is based on the assumptions of one-dimensional flow at constant pressure with heat conduction as the principal heat transfer mechanism. Consideration is given to subsequent models which account for radiative transfer into the ambient gas; provide a two-dimensional generalization of Raizer's analysis for the subsonic propagation of laser sparks in air; include the effect of forward plasma radiation in a one-dimensional model; and attempt a time-dependent (elliptic) solution of the full Navier-Stokes equations for the flow in a simple axisymmetric thruster. Attention is also given to thruster and nozzle flow models and thermodynamic and transport properties.

Mccay, T. D.↗