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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

BLSTA: A boundary layer code for stability analysis

A computer program is developed to solve the compressible, laminar boundary-layer equations for two-dimensional flow, axisymmetric flow, and quasi-three-dimensional flows including the flow along the plane of symmetry, flow along the leading-edge attachment line, and swept-wing flows with a conical flow approximation. The finite-difference numerical procedure used to solve the governing equations is second-order accurate. The flow over a wide range of speed, from subsonic to hypersonic speed with perfect gas assumption, can be calculated. Various wall boundary conditions, such as wall suction or blowing and hot or cold walls, can be applied. The results indicate that this boundary-layer code gives velocity and temperature profiles which are accurate, smooth, and continuous through the first and second normal derivatives. The code presented herein can be coupled with a stability analysis code and used to predict the onset of the boundary-layer transition which enables the assessment of the laminar flow control techniques. A user's manual is also included.

Wie, Yong-Sun↗
NASA NTRS1995

Transport phenomena in the micropores of plug-type phase separators

This study numerically investigates the transport phenomena within and across a porous-plug phase separator. The effect of temperature differential across a single pore and of the sidewall boundary conditions, i.e., isothermal or linear thermal gradient, are presented and discussed. The effects are quantified in terms of the evaporation mass flux across the boundary and the mean surface temperature. A two-dimensional finite element model is used to solve the continuity, momentum, and energy equations for the liquid. Temperature differentials across the pore interface of 1.0, and 1.5 K are examined and their effect on evaporation flux and mean surface temperature is shown. For isothermal side boundary conditions, the evaporation flux across the pore is directly proportional and linear with Delta T. For the case of an imposed linear thermal gradient on the side boundaries, Biot numbers of 0.0, 0.15, and 0.5 are examined. The most significant effect of Biot number is to lower the overall surface temperature and evaporation flux.

Fazah, M. M.↗
NASA NTRS2004

On the Possibilities of Predicting Geomagnetic Secular Variation with Geodynamo Modeling

We use our MoSST core dynamics model and geomagnetic field at the core-mantle boundary (CMB) continued downward from surface observations to investigate possibilities of geomagnetic data assimilation, so that model results and current geomagnetic observations can be used to predict geomagnetic secular variation in future. As the first attempt, we apply data insertion technique to examine evolution of the model solution that is modified by geomagnetic input. Our study demonstrate that, with a single data insertion, large-scale poloidal magnetic field obtained from subsequent numerical simulation evolves similarly to the observed geomagnetic variation, regardless of the initial choice of the model solution (so long it is a well developed numerical solution). The model solution diverges on the time scales on the order of 60 years, similar to the time scales of the torsional oscillations in the Earth's core. Our numerical test shows that geomagnetic data assimilation is promising with our MoSST model.

Kuang, Wei-Jia↗
NASA NTRS2006

SiC Field Effect Transistor Technology Demonstrating Prolonged Stable Operation at 500 C

While there have been numerous reports of short-term transistor operation at 500 degree C or above, these devices have previously not demonstrated sufficient long-term operational durability at 500 degree C to be considered viable for most envisioned applications. This paper reports the development of Silicone Carbi field effect transistors capable of long-term electrical operation at 500 degree C. A 6H-SiC MESFET was packaged and subjected to continuous electrical operation while residing in a 500 degree C oven in oxidizing air atmosphere for over 2400 hours. The transistor gain, saturation current (IDSS), and on-resistance (RDS) changed by less than 20% from initial values throughout the duration of the biased 500 degree C test. Another high-temperature packaged 6H-SiC MESFET was employed to form a simple one-stage high-temperature low-frequency voltage amplifier. This single-stage common-source amplifier demonstrated stable continuous electrical operation (negligible changes to gain and operating biases) for over 600 hours while residing in a 500 degree C air ambient oven. In both cases, increased leakage from annealing of the Schottky gate-to-channel diode was the dominant transistor degradation mechanism that limited the duration of 500 degree C electrical operation.

Neudeck, Philip G.↗
DOE OSTI2023

High-Performance Accelerator Modeling: Toward Improving Controls and Diagnostics for High-Brightness Beams in Experiment

One of the main charges of beam physics is to improve the overall beam quality delivered to experiments and applications, which is generally quantified by the beam brightness. High brightness beams have numerous applications, including x-ray free electron lasers and ultrafast electron diffraction (UED). To these ends, this thesis details efforts to use high-performance models --- high-fidelity models that execute quickly --- for controls and diagnostics. To demonstrate the generality of these techniques, this thesis focuses on three of the most successful photoinjector designs currently used worldwide: the UCLA/SLAC/BNL-type high-gradient S-band gun, the continuous-wave high-repetition rate VHF APEX gun and the L-band DESY-PITZ-type gun. Work will be shown from Pegasus (UCLA), HiRES (LBNL) and FAST (FNAL).Specifically, data-driven models for online virtual diagnostics are presented, in this case, in the context of UED at HiRES, leading to a temporal resolution improvement. Methods for improving the fidelity of physics-based models are discussed, with examples at HiRES, Pegasus and FAST. Markov-chain Monte Carlo analysis is applied to match simulations, in the context of photocathode studies at HiRES and Pegasus. Lastly, the augmentation of online model-based predictions with model-independent optimization is explored in a fluctuating environment at Pegasus and HiRES. A central theme of this dissertation is working with parameter fluctuations when modeling an accelerator beamline. Long-term drifts and shot-to-shot jitter exist in every accelerator to a varying degree and therefore play an important role in every chapter of this dissertation. This thesis attempts to address the issues associated with these fluctuations when trying to develop faithful model representations of the system.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Quantum advantage for differential equation analysis

Quantum algorithms for differential equation solving, data processing, and machine learning potentially offer an exponential speedup over all known classical algorithms. However, there also exist obstacles to obtaining this potential speedup in useful problem instances. The essential obstacle for quantum differential equation solving is that outputting useful information may require difficult postprocessing, and the essential obstacle for quantum data processing and machine learning is that inputting the data is a difficult task just by itself. In this study, we demonstrate that, when combined, these difficulties solve one another. We show how the output of quantum differential equation solving can serve as the input for quantum data processing and machine learning, allowing dynamical analysis in terms of principal components, power spectra, and wavelet decompositions. To illustrate this, we consider continuous-time Markov processes on epidemiological and social networks. These quantum algorithms provide an exponential advantage over existing classical Monte Carlo methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2004

Status of NASA Satellite, Field Observations, and Numerical Modeling Addressing the Impact of Urbanization on Short and Long Term Precipitation Variability

Howard (1833a) made the first documented observation of a temperature difference between an urban area and its rural environment. Manley (1958) termed this contrast the "urban heat island (UHI)". The UHI has now become a widely acknowledged, observed, and researched phenomenon because of its broad implications. It is estimated that by the year 2025, 60% of the world's population will live in cities (UNFP, 1999). In the United States, the current urban growth rate is approximately 12.5%, with 80% currently living in urban areas. As cities continue to grow, urban sprawl creates unique problems related to land use, transportation, agriculture, housing, pollution, and development for policymakers. Urban expansion and its associated urban heat islands also have measurable impacts on weather and climate processes.

Shepherd, J. Marshall↗
DOE OSTI2020

Fast and Accurate Core Analysis by the Full-Immersion Pressure-Pulse Decay: Part 2—Practice and Demonstration

Summary This study continues the work of presenting a novel approach for making petrophysical assessments of tight core samples. This method, the full-immersion pressure-pulse decay (Hannon 2019), involves applying a rapid increase in pressure in a chamber surrounding the entire outer surface area of a cylindrical sample, shutting the system in, and monitoring the pressure decay in the chamber as it reaches a new equilibrium. A precursor article covered the numerical simulator designed to model flow through the sample, demonstrating its performance and accuracy in addition to providing a first-order comparison between the speed and shape of the pressure-decay responses of the full-immersion method with those of other similar transient methods. This study covers the parameter-estimation procedure and experimental verification through a proof-of-concept laboratory investigation. The investigations provided here demonstrate that under appropriate, achievable experimental conditions, the pressure data can be analyzed in such a way that returns an estimate of the porosity and apparent permeabilities both parallel and perpendicular to bedding from a single test performed on a single cylindrical sample. After determining these experimental conditions (the uniqueness window), this report outlines a data-inversion strategy to estimate the petrophysical properties (porosity, horizontal permeability, and vertical permeability) from each test. This strategy is put to the test through comparisons with measurements performed by a commercial core laboratory. A common set of samples recovered from an outcrop of a tight-gas sandstone formation were investigated using the full-immersion method, and their results are compared with those from conventional steady-state measurement procedures performed by the commercial laboratory. Comparisons between petrophysical characterizations of these samples, which had permeabilities between 25 nd and 2.3 μd, demonstrated close agreement in most cases. However, whereas steady-state measurements performed at the professional laboratory required 4 to 5 hours of testing time per measurement of a single permeability, similar assessments using the full-immersion technique, requiring approximately 5 to 10 minutes to complete, returned estimates of the horizontal and vertical permeability simultaneously. Additional analyses are provided to determine principal reasons of discrepancies in instances where agreement was not as strong. Based on lessons learned from these experiences, the report closes with suggestions on areas of improvement in the experimental approach. Once complete, these developments should propel this technology to fill a critical need to determine petrophysical properties (porosity and permeability) of tight rocks in a time-efficient manner and in a way that does not compromise their accuracy.

Energy & Fuels↗
DOE OSTI2021

Vessel Fabrication Overview

LANL utilizes spherical vessels for numerous basic science, materials, and system performance experiments. These pressure vessels are designed and fabricated to meet the intent of ASME B&PV Section VIII, Division 3 and Code Case 2564. LANL has successfully contracted vessel construction in the past, utilizing HSLA-100 steel. As current vessel inventories approach design life, we are once again embarking upon a series of fabrication contracts to produce replacement components and vessel assemblies to support continued experiment cadence at LANL. LANL’s intent is to ramp up a continuous vessel supply chain for the DOE Complex.

42 ENGINEERING↗
NASA NTRS1973

Application of ERTS-1-data to the protection and management of New Jersey's coastal environment

The author has identified the following significant results. New Jersey's planned, regionalized network of sewage disposal facilities has been plotted on an ERTS-1 mosaic and circulation parameters for each of the planned outfall locations have been analyzed using the ERTS-1 imagery and comparative aircraft photography. Work is continuing on the circulation and dispersion of barge-dumped wastes in the New York Bight area. One of the largest remote sensing experiments ever attempted in this country was completed on April 7, 1973 during the ERTS-1 overpass. The test area included the northern portion of New Jersey and the Raritan Bay - New York Harbor area. Three NASA aircraft, two helicopters, nine surface vessels, 40 ground team personnel, and numerous oceanographic, radiometric, and meteorological equipment were deployed in an effort to characterize the surface and near-surface circulation dynamics in this 600 square mile area, during an entire tidal cycle. The analyses of these data in concert with all previous ERTS-1 overpasses will provide information that can lead to a better and more rational use of the nearshore marine environment. The data will be utilized to plan future outfall locations, regulating offshore disposal of wastes, etc.

Yunghans, R. S.↗
NASA NTRS1985

Finite-amplitude steady waves in plane viscous shear flows

Computations of two-dimensional solutions of the Navier-Stokes equations are carried out for finite-amplitude waves on steady unidirectional flow. Several cases are considered. The numerical method employs pseudospectral techniques in the streamwise direction and finite differences on a stretched grid in the transverse direction, with matching to asymptotic solutions when unbounded. Earlier results for Poiseuille flow in a channel are re-obtained, except that attention is drawn to the dependence of the minimum Reynolds number on the physical constraint of constant flux or constant pressure gradient. Attempts to calculate waves in Couette flow by continuation in the velocity of a channel wall fail. The asymptotic suction boundary layer is shown to possess finite-amplitude waves at Reynolds numbers orders of magnitude less than the critical Reynolds number for linear instability. Waves in the Blasius boundary layer and unsteady Rayleigh profile are calculated by employing the artifice of adding a body force to cancel the spatial or temporal growth. The results are verified by comparison with perturbation analysis in the vicinity of the linear-instability critical Reynolds numbers.

Milinazzo, F. A.↗
NASA NTRS1989

On the secondary instability in plane Poiseuille flow

Numerical experiments were performed to clarify apparent differences between experimental observations and a theoretical prediction of the secondary instability in plane Poiseuille flow. It is shown that subharmonic breakdown is unlikely in natural transition as a result of the initial growth of what we call the 'minus' modes and consequent forcing of Orr-Sommerfeld modes present in the background noise. Subharmonic breakdown was achieved only when these minus modes were continuously suppressed.

Kim, John↗
NASA NTRS1990

Vibration and wave propagation characteristics of multisegmented elastic beams

Closed form analytical solutions are derived for the vibration and wave propagation of multisegmented elastic beams. Each segment is modeled as a Timoshenko beam with possible inclusion of material viscosity, elastic foundation and axial forces. Solutions are obtained by using transfer matrix methods. According to these methods formal solutions are first constructed which relate the deflection, slope, moment and shear force of one end of the individual segment to those of the other. By satisfying appropriate continuity conditions at segment junctions, a global 4x4 matrix results which relates the deflection, slope, moment and shear force of one end of the beam to those of the other. If any boundary conditions are subsequently invoked on the ends of the beam one gets the appropriate characteristic equation for the natural frequencies. Furthermore, by invoking appropriate periodicity conditions the dispersion relation for a periodic system is obtained. A variety of numerical examples are included.

Nayfeh, Adnan H.↗
NASA NTRS1991

Programmable digital modem

The design of the Programmable Digital Modem (PDM) is outlined. The PDM will be capable of operating with numerous modulation techniques including: 2-, 4-, 8- and 16-ary phase shift keying (PSK), minimum shift keying (MSK), and 16-ary quadrature amplitude modulation (QAM), with spectral occupancy from 1.2x to 2x the data symbol rate. It will also be programmable for transmission rates ranging from 2.34 to 300 Mbit/s, where the maximum symbol rate is 75 Msymbol/s. Furthermore, these parameters will be executable in independent burst, dependent burst, or continuous mode. In dependent burst mode the carrier and clock oscillator sources are common from burst to burst. To achieve as broad a set of requirements as these, it is clear that the essential signal processing must be digital. In addition, to avoid hardware changes when the operational parameters are changed, a fixed interface to an analog intermediate frequency (IF) is necessary for transmission; and, common system level architectures are necessary for the modulator and demodulator. Lastly, to minimize size and power, as much of the design as possible will be implemented with application specific integrated circuit (ASIC) chips.

Poklemba, John J.↗
NASA NTRS1991

Analysis of metal-matrix composite structures. I - Micromechanics constitutive theory. II - Laminate analyses

The micromechanical constitutive theory is used to examine the nonlinear behavior of continuous-fiber-reinforced metal-matrix composite structures. Effective lamina constitutive relations based on the Abouli micromechanics theory are presented. The inelastic matrix behavior is modeled by the unified viscoplasticity theory of Bodner and Partom. The laminate constitutive relations are incorporated into a first-order deformation plate theory. The resulting boundary value problem is solved by utilizing the finite element method. Attention is also given to computational aspects of the numerical solution, including the temporal integration of the inelastic strains and the spatial integration of bending moments. Numerical results the nonlinear response of metal matrix composites subjected to extensional and bending loads are presented.

Arenburg, R. T.↗
NASA NTRS1993

Head-on collision of drops: A numerical investigation

The head-on collision of equal sized drops is studied by full numerical simulations. The Navier-Stokes equations are solved for fluid motion both inside and outside the drops using a front tracking/finite difference technique. The drops are accelerated toward each other by a body force that is turned off before the drops collide. When the drops collide, the fluid between them is pushed outward leaving a thin later bounded by the drop surface. This layer gets progressively thinner as the drops continue to deform and in several of the calculations this double layer is artificially removed once it is thin enough, thus modeling rupture. If no rupture takes place, the drops always rebound, but if the film is ruptured the drops may coalesce permanently or coalesce temporarily and then split again.

Nobari, M. R.↗
NASA NTRS2000

Forward Monte Carlo Computations of Polarized Microwave Radiation

Microwave radiative transfer computations continue to acquire greater importance as the emphasis in remote sensing shifts towards the understanding of microphysical properties of clouds and with these to better understand the non linear relation between rainfall rates and satellite-observed radiance. A first step toward realistic radiative simulations has been the introduction of techniques capable of treating 3-dimensional geometry being generated by ever more sophisticated cloud resolving models. To date, a series of numerical codes have been developed to treat spherical and randomly oriented axisymmetric particles. Backward and backward-forward Monte Carlo methods are, indeed, efficient in this field. These methods, however, cannot deal properly with oriented particles, which seem to play an important role in polarization signatures over stratiform precipitation. Moreover, beyond the polarization channel, the next generation of fully polarimetric radiometers challenges us to better understand the behavior of the last two Stokes parameters as well. In order to solve the vector radiative transfer equation, one-dimensional numerical models have been developed, These codes, unfortunately, consider the atmosphere as horizontally homogeneous with horizontally infinite plane parallel layers. The next development step for microwave radiative transfer codes must be fully polarized 3-D methods. Recently a 3-D polarized radiative transfer model based on the discrete ordinate method was presented. A forward MC code was developed that treats oriented nonspherical hydrometeors, but only for plane-parallel situations.

Battaglia, A.↗