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 613 records · Page 34

Numerical studies of motion and decay of vortex filaments

A computational code is developed for the integro-differential equations governing the motion of the centerlines of vortex filaments submerged in a background potential flow. These equations, which are derived from the method of matched asymptotic analysis, include the effect of decaying large-magnitude circumferential and axial velocity components in the vortical cores. Numerical examples are presented to assess the effect of large axial velocity and of nonsimilar initial profiles in vortical cores. The initial configurations of the filaments are chosen so as to fulfill the basic assumption of asymptotic analysis, which is the effective vortical core size is much smaller than all other length scales in the flowfield, e.g., the radius of curvature and interfilament distance. The computations are continued until the basic assumption is no longer valid, that is, when the merging or intersection of filaments have begun. Various types of local or global merging or intersection of filaments are classified and demonstrated by numerical examples.

Liu, C. H.↗

Steady inviscid three-dimensional flows

The present analysis combines some of the theoretical concepts suggested by Hawthorne (1955) with a numerical integration procedure suggested by Martin (1978). The resulting algorithm is for inviscid subsonic flows. Thus, it is restricted to high Reynolds number flows. Chang and Adamczyk (1983) have provided a detailed derivation of the present algorithm along with a discussion of its stability bounds. The present paper represents a summary of this work. The integration of the continuity equation is considered along with an evaluation of the entropy, total temperature, and vorticity field. Attention is given to the shear-flow algorithm construction, and an application to a shear flow in a turning channel. A description of numerical results is also provided. The discussed algorithm represents a new procedure for solving inviscid subsonic three-dimensional rotational flows.

Adamczyk, J. J.↗

Drop Breakup in Fixed Bed Flows as Model Stochastic Flow Fields

We examine drop breakup in a class of stochastic flow fields as a model for the flow through fixed fiber beds and to elucidate the general mechanisms whereby drops breakup in disordered, Lagrangian unsteady flows. Our study consists of two parallel streams of investigation. First, large scale numerical simulations of drop breakup in a class of anisotropic Gaussian fields will be presented. These fields are generated spectrally and have been shown in a previous publication to be exact representations of the flow in a dilute disordered bed of fibers if close interactions between the fibers and the drops are dynamically unimportant. In these simulations the drop shape is represented by second and third order small deformation theories which have been shown to be excellent for the prediction of drop breakup in steady strong flows. We show via these simulations that the mechanisms of drop breakup in these flows are quite different than in steady flows. The predominant mechanism of breakup appears to be very short lived twist breakups. Moreover, the occurrence of breakup events is poorly predicted by either the strength of the local flow in which the drop finds itself at breakup, or the degree of deformation that the drop achieves prior to breakup. It is suggested that a correlation function of both is necessary to be predictive of breakup events. In the second part of our research experiments are presented where the drop deformation and breakup in PDMS/polyisobutylene emulsions is considered. We consider very dilute emulsions such that coalescence is unimportant. The flows considered are simple shear and the flow through fixed fiber beds. Turbidity, small angle light scattering, dichroism and microscopy are used to interrogate the drop deformation process in both flows. It is demonstrated that breakup at very low capillary numbers occurs in both flows but larger drop deformation occurs in the fixed bed flow. Moreover, it is witnessed that breakup in the bed occurs continuously during flow and apparently with uniform probability through the bed length. The drop deformations witnessed in our experiments are larger than those predicted by the numerical simulations, and future plans to investigate these differences are discussed.

Shaqfeh, Eric S. G.↗

Primordial random motions and angular momenta of galaxies and galaxy clusters.

We study the decay of primordial random motions of galaxies and galaxy clusters in an expanding universe by solving a kinetic equation for the relaxation of differential energy spectra N(E, t). Systematic dissipative energy losses are included, involving gravitational drag by, and accretion of, intergalactic matter, as well as the effect of collisions with other systems. Formal and numerical solutions are described for two distinct modes of galaxy formation in a turbulent medium, corresponding to formation at a distinct epoch and to continuous formation of galaxies. We show that any primordial random motions of galaxies at the present epoch can amount to at most a few km/sec, and that collisions at early epochs can lead to the acquisition of significant amounts of primordial angular momentum.

Silk, J.↗

Transonic flow calculations

The development of relaxation methods to calculate transonic flows is discussed. Rather accurate predictions can be made for a number of flows of interest using the transonic potential flow equation, which may be derived from the Euler equations for inviscid compressible flow by introducing the assumption that the flow is irrotational. Once the choice of a mathematical model has been settled, the numerical procedure for actually computing a solution contains two main elements: the construction of a discrete approximation which converges to the solution of the continuous problem in the limit as the mesh width is reduced to zero, and the solution of the resulting set of nonlinear difference equations by a convergent iterative scheme. The choice of an appropriate coordinate system and its influence on the accuracy of the discrete approximation is discussed. Several applications of the general method are described.

Jameson, A.↗

Study of Advanced Propulsion Systems for Small Transport Aircraft Technology (STAT) Program

Definitions of takeoff gross weight, performance, and direct operating cost for both a 30 and 50 passenger airplane were established. The results indicate that a potential direct operating cost benefit, resulting from advanced technologies, of approximately 20 percent would be achieved for the 1990 engines. Of the numerous design features that were evaluated, only maintenance-related items contributed to a significant decrease in direct operating cost. Recommendations are made to continue research and technology programs for advanced component and engine development.

Baerst, C. F.↗

Large velocity shears and associated electrostatic waves and turbulence in the auroral F region

Broadband electrostatic waves at 10-1000 Hz have been observed with very large shears in the plasma flow velocity transverse to the ambient magnetic field in the auroral F region. The shears were detected through their perpendicular electric field signatures, which changed by as much as 200 mV/m over distances of only a few hundred meters. Transverse shears can be uniquely related to field-aligned currents through the current continuity equation, and the resulting field-aligned drift exceeds the threshold for excitation of current-driven electrostatic ion-cyclotron waves. A numerical simulation of this instability has been used to generate electric-field spectra in the rocket frame of reference, and these spectra are similar to the spectra generated form the actual rocket data.

Earle, G. D.↗

A new family of stable elements for the Stokes problem based on a mixed Galerkin/least-squares finite element formulation

Adding to the classical Hellinger-Reissner formulation, a residual form of the equilibrium equation, a new Galerkin/least-squares finite element method is derived. It fits within the framework of a mixed finite element method and is stable for rather general combinations of stress and velocity interpolations, including equal-order discontinuous stress and continuous velocity interpolations which are unstable within the Galerkin approach. Error estimates are presented based on a generalization of the Babuska-Brezzi theory. Numerical results (not presented herein) have confirmed these estimates as well as the good accuracy and stability of the method.

Franca, Leopoldo P.↗

Calculations of particle-laden flows by second-moment closures

A differential second-moment closure is applied to calculate nonswirling turbulent particle-laden jet flow. The closure is effected by extending constant-density single-phase turbulence models to the continuous (carrier) phase of two-phase flow. The dispersed (discrete) phase is treated by stochastic-Lagrangian technique. To reduce the effects of numerical (false) diffusion on the predicted results, a higher order differencing scheme, namely the flux-spline, is employed. In addition, an experimental study is conducted to provide data of good quality, especially near the inlet, for model assessment. The experiment consists of a confined coaxial jet flow of air with glass beads. The beads' diameter is 105 microns and is injected from the inner jet with a mass loading ratio of 0.9. The computed results are in good agreement with experimental data, especially for the carrier phase. However, the predicted turbulence fluctuations associated with the discrete phase are lower than the data at the downstream region of the flow field.

Nikjooy, M.↗

A High Order Discontinuous Galerkin Method for 2D Incompressible Flows

In this paper we introduce a high order discontinuous Galerkin method for two dimensional incompressible flow in vorticity streamfunction formulation. The momentum equation is treated explicitly, utilizing the efficiency of the discontinuous Galerkin method The streamfunction is obtained by a standard Poisson solver using continuous finite elements. There is a natural matching between these two finite element spaces, since the normal component of the velocity field is continuous across element boundaries. This allows for a correct upwinding gluing in the discontinuous Galerkin framework, while still maintaining total energy conservation with no numerical dissipation and total enstrophy stability The method is suitable for inviscid or high Reynolds number flows. Optimal error estimates are proven and verified by numerical experiments.

Liu, Jia-Guo↗

Mars Exploration 2003 to 2013 - An Integrated Perspective: Time Sequencing the Missions

The science goals for the Mars exploration program, together with the HEDS precursor environmental and technology needs, serve as a solid starting point for re-planning the program in an orderly way. Most recently, the community has recognized the significance of subsurface sampling as a key component in "following the water". Accessing samples from hundreds and even thousands of meters beneath the surface is a challenge that will call for technology development and for one or more demonstration missions. Recent mission failures and concerns about the complexity of the previously planned MSR missions indicate that, before we are ready to undertake sample return and deep sampling, the Mars exploration program needs to include: 1) technology development missions; and 2) basic landing site assessment missions. These precursor missions should demonstrate the capability for reliable & accurate soft landing and in situ propellant production. The precursor missions will need to carry out close-up site observations, ground-penetrating radar mapping from orbit and conduct seismic surveys. Clearly the programs should be planned as a single, continuous exploration effort. A prudent minimum list of missions, including surface rovers with ranges of more than 10 km, can be derived from the numerous goals and requirements; they can be sequenced in an orderly way to ensure that time is available to feed forward the results of the precursor missions. One such sequence of missions is proposed for the decade beginning in 2003.

Briggs, G.↗

Impact of GPM Rainrate Data Assimilation on Simulation of Hurricane Harvey (2017)

Built upon Tropical Rainfall Measuring Mission (TRMM) legacy for next-generation global observation of rain and snow. The GPM was launched in February 2014 with Dual-frequency Precipitation Radar (DPR) and GPM Microwave Imager (GMI) onboard. The GPM has a broad global coverage approximately 70deg S -70deg N with a swath of 245/125-km for the Ka (35.5 GHz)/Ku (13.6 GHz) band radar, and 850-km for the 13-channel GMI. GPM also features better retrievals for heavy, moderate, and light rain and snowfall To develop methodology to assimilate GPM surface precipitation data with Grid-point Statistical Interpolation (GSI) data assimilation system and WRF ARW model To investigate the potential and the value of utilizing GPM observation into NWP for operational environment The GPM rain rate data has been successfully assimilated using the GSI rain data assimilation package. Impacts of rain rate data have been found in temperature and moisture fields of initial conditions. 2.Assimilation of either GPM IMERG or GPROF rain product produces significant improvement in precipitation amount and structure for Hurricane Harvey (2017) forecast. Since IMERG data is available half-hourly, further forecast improvement is expected with continuous assimilation of IMERG data

Precipitation↗

Leveraging Pre-storm Soil Moisture Estimates for Enhanced Land Surface Model Calibration in Ungauged Hydrologic Basins

Despite long-standing efforts, hydrologists still lack robust tools for calibrating land surface model (LSM) streamflow estimates within ungauged basins. Using surface soil moisture estimates from the Soil Moisture Active Passive Level 4 Soil Moisture (L4_SM) product, precipitation observations, and streamflow gauge measurements for 617 medium-scale (200-10,000 km2) basins in the contiguous United States, we measure the temporal (Spearman) rank correlation between antecedent (i.e., pre-storm) surface soil moisture (ASM) and the storm-scale runoff coefficient (RC; the fraction of storm-scale precipitation accumulation converted into streamflow). In humid and semi-humid basins, this rank correlation is shown to be sufficiently strong to allow for the substitution of storm-scale RC observations (available only in basins that are both lightly regulated and gauged) with high-quality ASM values (available quasi-globally from L4_SM) in streamflow calibration procedures. Using this principle, we define a new, basin-wise LSM streamflow calibration approach based on L4_SM alone and successfully apply it to identify LSM configurations that produce a high rank correlation with observed RC. However, since the approach cannot detect RC bias, it is less successful in identifying LSM configurations with low mean-absolute error. Plain Text Summary Accurately forecasting the fraction of rainfall that runs off into streams, as opposed to infiltrates into the soil, is critical for flash-flood prediction, water-resource monitoring, and tracking the transport of nutrients from agricultural fields into local waterways. Such forecasting is typically performed by hydrologic models that attempt to represent the physical processes responsible for surface runoff generation. However, to provide accurate streamflow forecasts, these models typically need to be calibrated against actual streamflow observations. This is problematic given the relatively poor, and declining, global availability of stream gauges. This paper presents a novel model calibration strategy that uses soil moisture from remote sensing and numerical modeling in place of streamflow observations during calibration. This transition has significant practical advantages because, unlike streamflow observations, the soil moisture data are continuously available across space. Our results demonstrate that this new approach can significantly improve hydrologic models within humid and semi-humid basins lacking sufficient ground-based instrumentation for traditional streamflow calibration.

hydrologic models↗

Transport Phenomena in Thin Rotating Liquid Films Including: Nucleate Boiling

In this grant, experimental, numerical and analytical studies of heat transfer in a thin liquid film flowing over a rotating disk have been conducted. Heat transfer coefficients were measured experimentally in a rotating disk heat transfer apparatus where the disk was heated from below with electrical resistance heaters. The heat transfer measurements were supplemented by experimental characterization of the liquid film thickness using a novel laser based technique. The heat transfer measurements show that the disk rotation plays an important role on enhancement of heat transfer primarily through the thinning of the liquid film. Experiments covered both momentum and rotation dominated regimes of the flow and heat transfer in this apparatus. Heat transfer measurements have been extended to include evaporation and nucleate boiling and these experiments are continuing in our laboratory. Empirical correlations have also been developed to provide useful information for design of compact high efficiency heat transfer devices. The experimental work has been supplemented by numerical and analytical analyses of the same problem. Both numerical and analytical results have been found to agree reasonably well with the experimental results on liquid film thickness and heat transfer Coefficients/Nusselt numbers. The numerical simulations include the free surface liquid film flow and heat transfer under disk rotation including the conjugate effects. The analytical analysis utilizes an integral boundary layer approach from which

Faghri, Amir↗

Resolution requirements for velocity gradients in turbulence

Since high resolution numerical simulations of turbulent flows, or experiments at high Reynolds numbers, represent a substantial investment in resources, the estimation of the minimum resolution required for the study of a given property has been the subject of continued interest. Different properties require, in general, different resolutions, and the present paper is dedicated to the requirements for the measurement of the probability distribution functions of the velocity gradients and, in particular, of their low order moments. The deviation of these quantities from the values corresponding to a Gaussian distribution was one of the first indications of the presence of Reynolds number-dependent intermittency and has been the object of recent interest as numerical simulations have become able to explore the distribution of gradients in the low Reynolds number range, while new experiments have extended the range to increasingly high Reynolds numbers. We will use progressive filtering of the results of numerical simulations of isotropic turbulence as a model for the effect of a sensor of finite size. The numerical issues will be addressed first to ensure that the simulations are fully resolved from the point of view of the velocity gradients. This will also give us an estimate for the numerical resolution required for the different quantities.

Jimenez, Javier↗

Technical Report Series on Global Modeling and Data Assimilation, Volume 43: Initial Evaluation of the Climate - MERRA-2

The years since the introduction of MERRA have seen numerous advances in the GEOS-5 Data Assimilation System as well as a substantial decrease in the number of observations that can be assimilated into the MERRA system. To allow continued data processing into the future, and to take advantage of several important innovations that could improve system performance, a decision was made to produce MERRA-2, an updated retrospective analysis of the full modern satellite era. One of the many advances in MERRA-2 is a constraint on the global dry mass balance; this allows the global changes in water by the analysis increment to be near zero, thereby minimizing abrupt global interannual variations due to changes in the observing system. In addition, MERRA-2 includes the assimilation of interactive aerosols into the system, a feature of the Earth system absent from previous reanalyses. Also, in an effort to improve land surface hydrology, observations-corrected precipitation forcing is used instead of model-generated precipitation. Overall, MERRA-2 takes advantage of numerous updates to the global modeling and data assimilation system. In this document, we summarize an initial evaluation of the climate in MERRA-2, from the surface to the stratosphere and from the tropics to the poles. Strengths and weaknesses of the MERRA-2 climate are accordingly emphasized.

MERRA-2↗

Numerical Model Studies of the Martian Mesoscale Circulations

Studies concerning mesoscale topographical effects on Martian flows examined low-level jets in the near equatorial latitudes and the dynamical intensification of flow by steep terrain. Continuation of work from previous years included evaluating the dissipation of cold air mass outbreaks due to enhanced sensible heat flux, further sensitivity and scaling evaluations for generalization of the characteristics of Martian mesoscale circulation caused by horizontal sensible heat-flux gradients, and evaluations of the significance that non-uniform surface would have on enhancing the polar CO2 ice sublimation during the spring. The sensitivity of maximum and minimum atmospheric temperatures to changes in wind speed, surface albedo, and deep soil temperature was investigated.

Segal, M.↗

Unified approach for incompressible flows

A unified approach for solving incompressible flows has been investigated in this study. The numerical CTVD (Centered Total Variation Diminishing) scheme used in this study was successfully developed by Sanders and Li for compressible flows, especially for the high speed. The CTVD scheme possesses better mathematical properties to damp out the spurious oscillations while providing high-order accuracy for high speed flows. It leads us to believe that the CTVD scheme can equally well apply to solve incompressible flows. Because of the mathematical difference between the governing equations for incompressible and compressible flows, the scheme can not directly apply to the incompressible flows. However, if one can modify the continuity equation for incompressible flows by introducing pseudo-compressibility, the governing equations for incompressible flows would have the same mathematical characters as compressible flows. The application of the algorithm to incompressible flows thus becomes feasible. In this study, the governing equations for incompressible flows comprise continuity equation and momentum equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. Thus, the CTVD schemes can be implemented. In addition, the physical and numerical boundary conditions are properly implemented by the characteristic boundary conditions. Accordingly, a CFD code has been developed for this research and is currently under testing. Flow past a circular cylinder was chosen for numerical experiments to determine the accuracy and efficiency of the code. The code has shown some promising results.

Chang, Tyne-Hsien↗