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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

NASA NTRS2018

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↗

NASA NTRS2022

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↗

NASA NTRS2005

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↗

NASA NTRS1994

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↗

NASA NTRS2015

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↗

NASA NTRS1996

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.↗

NASA NTRS1995

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↗

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.↗

DOE OSTI2020

On differentiable local bounds preserving stabilization for Euler equations

This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A differentiable local bounds preserving method has been developed, which combines a Rusanov artificial diffusion operator and a differentiable shock detector. Nonlinear stabilization schemes are usually stiff and highly nonlinear. This issue is mitigated by the differentiability properties of the proposed method. Moreover, in order to further improve the nonlinear convergence, we also propose a continuation method for a subset of the stabilization parameters. The resulting method has been successfully applied to steady and transient problems with complex shock patterns. Numerical experiments show that it is able to provide sharp and well resolved shocks. Furthermore, the importance of the differentiability is assessed by comparing the new scheme with its non-differentiable counterpart. Numerical experiments suggest that, for up to moderate nonlinear tolerances, the method exhibits improved robustness and nonlinear convergence behavior for steady problems. Additionally, in the case of transient problem, we also observe a reduction in the computational cost.

42 ENGINEERING↗

DOE OSTI2026

Dynamically reconfigurable topological routing in nonlinear photonic systems

The propagation path of topologically protected states is bound to the interface between regions with different topology, and as such, the functionality of linear photonic devices leveraging these states is fixed during fabrication. Here, we propose a mechanism for dynamic control over a driven dissipative system’s local topology, yielding reconfigurable topological interfaces and thus tunable paths for protected routing. We illustrate our approach in non-resonantly pumped polariton lattices, where the nonlinear interaction between the polaritons and the exciton reservoir due to non-resonant pumping can yield picosecond-scale changes in the propagation paths of the chiral edge states. To analytically confirm the numerically observed topological dynamics, we generalize the spectral localizer framework to non-linear non-Hermitian Chern materials and apply this framework to a continuous model of the polariton system based on a driven-dissipative Gross-Pitaevskii equation. In doing so, we show that the local changes in the polariton lattice’s topology are captured by a local Chern marker. Looking forward, we anticipate such reconfigurable topological routing will enable the realization of novel classes of topological photonic devices.

Wong, Stephan [Sandia National Laboratories (SNL-N↗

DOE OSTI2023

National Transmission Planning Study - Long-Term Modeling Results

The National Transmission Planning Study (NTP Study) analyzed the transformation needed to ensure the U.S. transmission system continues to reliably meet electricity demand as the power sector evolves. The study linked several long- and short-term power systems models to test numerous interregional and regional transmission buildout scenarios. The results included in this dataset are from the long-term capacity expansion modeling, completed with NREL's ReEDS model. This includes nearly 100 future transmission scenarios with a wide range of economic, reliability, and resilience conditions. The dataset includes all the results visualized in the Tableau Public site, as linked below.

Array↗

DOE OSTI2025

Model-based analysis of solute transport and potential carbon mineralization in the active layer of a hillslope underlain by permafrost with seasonal variability and climate change

Permafrost carbon, stored in frozen organic matter across vast Arctic and sub-Arctic regions, represents a substantial and increasingly vulnerable carbon reservoir. As global temperatures rise, the accelerated thawing of permafrost releases greenhouse gases, exacerbating climate change. However, freshly thawed permafrost carbon may also experience lateral transport by groundwater flow to surface water recipients such as rivers and lakes, increasing the terrestrial-to-aquatic transfer of permafrost carbon. Mobilization and subsurface transport mechanisms are poorly understood and not accounted for in global climate models, leading to high uncertainties in the predictions of the permafrost carbon feedback. Here, we focus on a hillslope in Endalen Valley, Svalbard, as a representative example of a high-Arctic hillslope underlain by continuous permafrost. We analyze solute transport in the form of a non-reactive tracer representing dissolved organic carbon (DOC) using a physics-based numerical model with the objective to study governing cryotic and hydrodynamic transport mechanisms relevant for warming permafrost regions. We first analyze transport times for DOC pools at different locations within the active layer under present-day climatic conditions and proceed to study susceptibility for deeper ancient carbon release in the upper permafrost due to thaw under different warming scenarios. Results suggest that DOC in the active layer near the permafrost table experiences rapid lateral transport upon thaw due to saturated conditions and lateral flow, while DOC close to the ground surface experiences slower transport due to flow in unsaturated soil. Deeper permafrost carbon release exhibits vastly different transport behaviors depending on warming and thaw rate. Gradual warming leads to small fractions of DOC being mobilized every year, while the majority moves vertically through percolation and cryosuction. Abrupt thaw resulting from a single very warm year leads to faster lateral transport times, similar to active layer DOC released in saturated conditions. Lastly, we analyze the potential susceptibility of DOC to mineralization to CO 2 prior to export due to soil moisture and temperature conditions. We find that high liquid saturation during transport coincides with very low mineralization rates and potentially inhibits mineralization into CO 2 before export. Overall, the results highlight the importance of subsurface hydrologic and thermal conditions for the retention and lateral export of permafrost carbon by subsurface flow.

Hamm, Alexandra [Stockholm Univ. (Sweden)] (ORCID:↗

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.↗