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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Development of Gridded Innovations and Observations Supplement to MERRA-2

Atmospheric reanalysis have become an important source of data for weather and climate research, owing to the continuity of the data, but especially because of the multitude of observational data included (radiosondes, commercial aircraft, retrieved data products and radiances). However, the presence of assimilated observations can vary based on numerous factors, and so it is difficult or impossible for a researcher to say with any degree of certainty how many and what type of observations contributed to the reanalysis data they are using at any give point in time or space. For example, quality control, transmission interruptions, and station outages can occasionally affect data availability. While orbital paths can be known, drift in certain instruments and the large number of available instruments makes it challenging to know which satellite is observing any region at any point in the diurnal cycle. Furthermore, there is information from the statistics generated by the data assimilation that can help understand the model and the quality of the reanalysis. Typically, the assimilated observations and their innovations are in observation-space data formats and have not been made easily available to reanalysis users.A test data set has been developed to make the MERRA-2 assimilated observations available for rapid and general use, by simplifying the data format. The observations are binned to a grid similar as MERRA-2 and saved as netCDF. This data collection includes the mean and number of observations in the bin as well as its variance. The data will also include the innovations from the data assimilation, the forecast departure and the analysis increment, as well as bias correction (for satellite radiances). We refer to this proof-of-concept data as the MERRA-2 Gridded Innovations and Observations (GIO). In this paper, we present the data format and its strengths and limitations with some initial testing and validation of the methodology.

SSMI↗
NASA NTRS2018

Development of Gridded Innovations and Observations Supplement to MERRA-2

Atmospheric reanalysis have become an important source of data for weather and climate research, owing to the continuity of the data, but especially because of the multitude of observational data included (radiosondes, commercial aircraft, retrieved data products and radiances). However, the presence of assimilated observations can vary based on numerous factors, and so it is difficult or impossible for a researcher to say with any degree of certainty how many and what type of observations contributed to the reanalysis data they are using at any give point in time or space, or their contribution to the eventual analyzed fields. For example, quality control, transmission interruptions, and station outages can occasionally affect data availability. While orbital paths are known, drift in certain instruments, cloud clearing, and the large number of available instruments makes it challenging to know which satellite is observing any region at any point in the diurnal cycle. Furthermore, there is information from the statistics generated by the data assimilation that can help understand the model and the quality of the reanalysis. Typically, the assimilated observations and their innovations are saved in observation-space data formats and are not easily made available to reanalysis users.Here, we present an early version of a data set has been developed to make the MERRA-2 assimilated observations available for rapid and general use, by simplifying the data format. The observations are binned to a grid similar as MERRA-2 and saved as netCDF. This data collection includes the mean and number of observations in the bin as well as its variance. The data will also include the innovations from the data assimilation, the forecast departure and the analysis increment, as well as bias correction (for satellite radiances). In this paper, we present the data format (called MERRA-2 Gridded Innovations and Observations or GIO) and its strengths and limitations with some initial testing and validation of the methodology.

Bosilovich, Michael G.↗
NASA NTRS1983

Higher modes of the Orr-Sommerfeld problem for boundary layer flows

The discrete spectrum of the Orr-Sommerfeld problem of hydrodynamic stability for boundary layer flows in semi-infinite regions is examined. Related questions concerning the continuous spectrum are also addressed. Emphasis is placed on the stability problem for the Blasius boundary layer profile. A general theoretical result is given which proves that the discrete spectrum of the Orr-Sommerfeld problem for boundary layer profiles (U(y), 0,0) has only a finite number of discrete modes when U(y) has derivatives of all orders. Details are given of a highly accurate numerical technique based on collocation with splines for the calculation of stability characteristics. The technique includes replacement of 'outer' boundary conditions by asymptotic forms based on the proper large parameter in the stability problem. Implementation of the asymptotic boundary conditions is such that there is no need to make apriori distinctions between subcases of the discrete spectrum or between the discrete and continuous spectrums. Typical calculations for the usual Blasius problem are presented.

Lakin, W. D.↗
DOE OSTI2023

A Refinement-by-Superposition -Method for (curl)- and (div)-Conforming Discretizations

Here, we present refinement-by-superposition (RBS) hp-refinement infrastructure for computational electromagnetics (CEMs), which permits exponential rates of convergence. In contrast to dominant approaches to hp-refinement for continuous Galerkin methods, which rely on explicit constraint equations, the multilevel strategy presented drastically reduces the implementation complexity. Through the RBS methodology, enforcement of continuity occurs by construction, enabling arbitrary levels of refinement with ease, and without the practical (but not theoretical) limitations of constrained-node refinement. We outline the construction of the RBS hp-method for refinement with H (curl)- and H (div)-conforming finite cells. Numerical simulations for the 2-D finite element method (FEM) solution of the Maxwell eigenvalue problem demonstrate the effectiveness of RBS hp-refinement. As an additional goal of this work, we aim to promote the use of mixed-order (low- and high-order) elements in practical CEM applications.

42 ENGINEERING↗
DOE OSTI2026

Best practices in software development for robust and reproducible geoscientific models based on insights from the Global Carbon Budget's dynamic vegetation models

Computational models play an increasingly vital role in scientific research by enabling the numerical simulation of complex processes. Such models are also fundamental in geosciences. For instance, they offer critical insights into the impacts of global change on the Earth system today and in the future. Beyond their value as research tools, models are also software products and should therefore adhere to certain established software engineering standards. However, scientists are rarely trained as software developers, which can lead to potential deficiencies in software quality like unreadable, inefficient, or erroneous code. The complexity of models, coupled with their integration into broader workflows, also often makes it challenging to reproduce results, evaluate processes, and build upon them. In this paper, we review the state and current practices of the development processes of the state-of-the-art land surface models used by the Global Carbon Budget. We combine the experience of modelers from the respective research groups with the expertise of software engineers from tech companies to outline key principles and tools for improving software quality in research. We explore four main areas: (1) model testing and validation, (2) scientific, technical, and user documentation, (3) version control, continuous integration, and code review, and (4) the portability and reproducibility of workflows. Our review reveals that while modeling communities are incorporating many best practices, significant room for improvement remains in areas such as automated testing, automated documentation, and reproducibility. Therefore, we here identify and promote essential software engineering practices, including numerous examples of practices from within the community that can serve as guidelines for other models and could help streamline processes across the entire community. We conclude with an open-source example implementation of these principles, demonstrating portable and reproducible data flows, a continuous integration setup, and web-based visualizations. This example may serve as a practical resource for model developers, users, and all scientists engaged in scientific programming.

Gregor, Konstantin [Technical Univ. of Munich (Ger↗
DOE OSTI2026

A variational mimetic finite difference method for elliptic interface problems on non-matching polytopal meshes with geometric interface inconsistencies

A new variational mimetic finite difference method for elliptic interface problems with perfect and imperfect thermal contacts on non-matching polytopal meshes with geometric interface inconsistencies is developed and analyzed theoretically and numerically. The method is defined on multiple non-matching submeshes with gaps and overlaps along their interfaces. The discrete equations are derived from a minimization problem for the augmented Dirichlet functional. For a perfect thermal contact, the functional uses a modified mimetic gradient with extended stencil which couples unknowns from both sides of an interface, as well as penalty terms to enforce weak continuity of temperature across the interface. The method leads to a symmetric positive definite matrix for any scaling of the penalty terms. For an imperfect thermal contact, the Dirichlet functional is supplemented with a quadratic jump term along the interface related to the interface thermal resistance. We prove that the method conserves the total heat flux across each interface. In conclusion, the obtained results are verified with numerical experiments showing convergence in the discrete L 2 and L ∞ norms.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Theoretical and experimental studies of baroclinic processes

The geophysical fluid flow simulation code, GEOSIM, is being used to study the phenomenon of vacillation in the baroclinic annulus. Having verified that the code predicts vacillation for the same points as the experiments, the work is aiming toward explaining the mechanics of vacillation and pointing out some of the sensitivities of the results to the numerical method. Researchers are finding that there is a structural change associated with amplitude vacillation, where the structural changes are in the vertical. The results disagree with the premise of Lindzen et al, that the vacillation is due to constructive and destructive interference of neutral modes with different phase speeds. The researchers are continuing to study the Spacelab 3 Geophysical Fluid Flow Cell (GFFC) results with horizontal temperature gradients and heating from below. GEOSIM has been used to compute a wide range of cases, and these are being compared with the observations. The computations and observations compare well, and the model is being used to extend the results beyond cases studied in the experiments and to study the mechanics and predictability of the flows. The study of fully nonlinear baroclinic instability using the GFFC apparatus is proceeding with the numerical code. While the first instability that occurs is of planetary scale, secondary instabilities consisting of small-scale, penetrative convection occurs where cold fluid flows over a warm surface. The simultaneous modeling of the planetary scale and the convective scale is possible because of the nonhydrostatic formulation of the model. Some of these results have been animated on the Stardent computer, which shows the explosive nature of the small-scale convection.

Miller, Tim↗
DOE OSTI2019

Large Area Airborne Contamination Monitoring

Current D&D operations at Hanford have demonstrated a flaw in the current state of the art capability of defining airborne contamination boundaries. Airborne particulate emissions of Pu-239 from CM2H operations on the Hanford Pu Finishing Plant were detected well beyond areas controlled for airborne plutonium, putting numerous workers at risk for radiological assimilations. Forty-two Hanford workers had positive bioassay results for plutonium uptake. Work areas determined to be free of airborne contamination as measured by state-of-the-art continuous alpha monitoring (CAM) systems were later discovered not to be. Airborne contamination was scattered to a such degree that CAM systems were surveying insufficient volumes of air to accurately predict where respiratory protection was required. A $\$ 19,000$ state of the art portable CAM system used at both Hanford and at SRS screens air for alpha activity at the rate of 1 CFM. This research project will explore the feasibility of using commercially available air purifying systems to screen significantly larger volumes of air for airborne actinide contamination. Commercial household HEPA based air purifying units screen large volumes of air, a simple $\$ 300$ unit being looked at for this study could screen a volume of air that would require over 300 of the $\$ 19,000$ CAMs. The goal of the project was to determine the feasibility of building a detection system capable of screening the air purifiers pre-filters and HEPA filters for any actinide bearing hot particles that have been captured. The simplest technologies available to Field-Radiological-Protection personnel are hand-held alpha field survey meters. With the size and the corrugated structure of these HEPA filters, hand screening is not practical. This project researched alternatives capable of automating this function. Two potential technologies were examined. One technology evaluated was the measurement of xray emissions from the actinides, taking advantage of the higher emission rates of the L x-rays compared to the gamma emissions of isotopes such as Pu-238 or Pu-239. High resolution to low resolution systems were evaluated. Branching ratio data on the actinides in the literature is incomplete, and a study was initiated to measure branching ratios of the various L emissions of a number of the plutonium isotopes. A Monte Carlo simulation was also run to calculate the effect the attenuating HEPA media would have on x-ray measurements on spectrometers centered in each quadrant above the surface of the HEPA. The second technology evaluated was measuring the light from a scintillating media applied to the surface of the HEPA filter using a ThorLabs Scientific-Grade Digital Camera Model 8051M. A light tight enclosure capable of holding the 2’x2’x1” HEPA filter was fabricated. Numerous scintillation media were evaluated for both the sensitivity and for resolution of measuring alpha activity. Figure 1 shows images of a 4”x4” section of HEPA filter with a 3 x 3 grid of Pu-238 contamination loaded inside the crevices of the HEPA.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Performance and wake flow characterization of a 1:8.7-scale reference USDOE MHKF1 hydrokinetic turbine to establish a verification and validation test database

As hydrokinetic turbine technologies continue to advance towards commercialization, public datasets on the performance characteristics for these devices and their flow field effects are invaluable to advance our understanding of these technologies and to validate analytical and numerical models. The Applied Research Laboratory at The Pennsylvania State University (ARL Penn State) collaborated with Sandia National Laboratories and the University of California at Davis to design, fabricate (at a 1:8.7 scale), and experimentally test a novel hydrokinetic turbine rotor design to provide an open platform and dataset for further study and development. The water tunnel test of this three-bladed, horizontal-axis rotor recorded power production, blade loading, the near-wake flow, cavitation effects, and noise generation. These state-of-the-art measurements demonstrate much of the complex physics associated with the flow through an unducted, horizontal-axis turbine, and they elucidate the performance characteristics and flow field effects at an unprecedented fidelity, accuracy and resolution. Measurements of powering coefficients (power, torque and thrust) as a function of tip-speed-ratio were performed. The dataset also includes unsteady measurements of driveshaft loading, blade strain, tower pressures, and radiated noise. Detailed flow mapping using laser Doppler velocimetry, and planar and stereo particle image velocimetry includes measurements of mean velocity and Reynolds stresses. Although the wake measurements are limited to less than half a diameter, they reveal the complex flow patterns in the near-wake structure of the rotor. The full database, available at the United States Department of Energy’s marine and hydrokinetic data repository, includes tunnel and model Computer Aided Design geometry files and inflow data sufficient for a “Model-the-Test” computational Verification and Validation study.

13 HYDRO ENERGY↗
DOE OSTI2025

Absorbing boundary conditions in material point method adopting perfectly matched layer theory

This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work integrates the Perfectly Matched Layer (PML) theory into the implicit MPM framework. The proposed approach introduces absorbing particles surrounding the computational domain that efficiently absorb outgoing waves and reduce reflections, allowing for accurate modeling of wave propagation and its further impact on geotechnical slope stability analysis. The study also includes several benchmark tests to validate the effectiveness of the proposed method, such as several types of impulse loading and symmetric and asymmetric base shaking. The conducted numerical tests also demonstrate the ability to handle large deformation problems, including the failure of elasto-plastic soils under gravity and dynamic excitations. The findings extend the capability of MPM in simulating continuous analysis of earthquake-induced landslides, from shaking to failure.

58 GEOSCIENCES↗
DOE OSTI2021

Wide Dust Gaps in Protoplanetary Disks Induced by Eccentric Planets: A Mass-eccentricity Degeneracy

The tidal perturbation of embedded protoplanets on their natal disks has been widely attributed to be the cause of gap-ring structures in submillimeter images of protoplanetary disks around T Tauri stars. Numerical simulations of this process have been used to propose scaling of characteristic dust-gap width/gap-ring distance with respect to planet mass. Applying such scaling to analyze observed gap samples yields a continuous mass distribution for a rich population of hypothetical planets in the range of several Earth to Jupiter masses. In contrast, the conventional core-accretion scenario of planet formation predicts a bimodal mass function due to (1) the onset of runaway gas accretion above ~20 Earth masses and (2) suppression of accretion induced by gap opening. Here, we examine the dust disk response to the tidal perturbation of eccentric planets as a possible resolution of this paradox. Based on simulated gas and dust distributions, we show the gap-ring separation of Neptune-mass planets with small eccentricities might become comparable to that induced by Saturn-mass planets on circular orbits. Furthermore, this degeneracy may obliterate the discrepancy between the theoretical bimodal mass distribution and the observed continuous gap width distribution. Despite damping due to planet–disk interaction, modest eccentricity may be sustained either in the outer regions of relatively thick disks or through resonant excitation among multiple super Earths. Moreover, the ring-like dust distribution induced by planets with small eccentricities is axisymmetric even in low viscosity environments, consistent with the paucity of vortices in Atacama Large Millimeter/submillimeter Array images.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

Aerodynamic influence coefficient method using singularity splines.

A new numerical formulation with computed results, is presented. This formulation combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of the loading function methods. The formulation employs a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfies all of the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise (termed 'spline'). Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral.

Mercer, J. E.↗
NASA NTRS1980

The radiation of sound by the instability waves of a compressible plane turbulent shear layer

The problem of acoustic radiation generated by instability waves of a compressible plane turbulent shear layer is solved. The solution provided is valid up to the acoustic far-field region. It represents a significant improvement over the solution obtained by classical hydrodynamic-stability theory which is essentially a local solution with the acoustic radiation suppressed. The basic instability-wave solution which is valid in the shear layer and the near-field region is constructed in terms of an asymptotic expansion using the method of multiple scales. This solution accounts for the effects of the slightly divergent mean flow. It is shown that the multiple-scales asymptotic expansion is not uniformly valid far from the shear layer. Continuation of this solution into the entire upper half-plane is described. The extended solution enables the near- and far-field pressure fluctuations associated with the instability wave to be determined. Numerical results show that the directivity pattern of acoustic radiation into the stationary medium peaks at 20 degrees to the axis of the shear layer in the downstream direction for supersonic flows. This agrees qualitatively with the observed noise-directivity patterns of supersonic jets.

Tam, C. K. W.↗
NASA NTRS1985

A Survey of Composite Grid Generation for General Three-dimensional Sections

The generation and use of composite grids for general three-dimensional physical boundary configurations is discussed, and the availability of several codes or procedures is noted. With the composite framework, the physical region is segmented into sub-regions, each bounded by six curved sides, and a grid is generated in each sub-region. These grids may be joined at the interfaces between the sub-regions with various degrees of continuity. This structure allows codes to be constructed to operate on rectangular blocks in computational space, so that existing solution procedures can be readily incorporated in the construction of codes for general configurations. Numerical grid generation is an integral part of the numerical solution of partial differential equations and is one of the pacing items in the development of codes for general configurations. The numerically generated grid frees the computational simulation from restriction to certain boundary shapes and allows general codes to be written in which the boundary shape is specified simply by input. The numerically generated grid allows all computation to be done on a fixed square grid in the computational space, which is always rectangular by construction.

Thompson, J. F.↗
NASA NTRS1985

NASA's space station planning - An update

The Manned Space Station (MSS) program is summarized, along with the program guidelines, activities, plans, management approach, international involvement and the present and future private sector role in space. The MSS, the first continually-manned facility in space, will be configured to grow modularly to meet future needs. Comprising both manned and unattached unmanned elements, the program is to advance the state of robotics and automation. Numerous uses of the MSS have been defined in support of and as the primary platform for astronomical observations, experimentation on prototype space-based industrial processes, and as a staging area for the launch and refurbishment of spacecraft. Current studies are targeted at defining the final design of the initial operating configuration, the hardware to be developed by the nations and organizations which are sharing in the work, and developing interest among potential commercial users of the MSS capabilities.

Hoban, F. T.↗
NASA NTRS1989

An integrating matrix formulation for buckling of rotating beams including the effects of concentrated masses

The integrating matrix technique of computational mechanics is extended to include the effects of concentrated masses. The stability of a flexible rotating beam with discrete masses is analyzed to determine the critical rotational speeds for buckling in the inplane and out-of-plane directions. In this problem, the beam is subjected to compressive centrifugal forces arising from steady rotation about an axis which does not pass through the clamped end of the beam. To determine the eigenvalues from which stability is assessed, the differential equations of motion are solved numerically by combining the extended integrating matrix method with an eigenanalysis. Stability boundaries for a discrete mass representation of a uniform beam are shown to asymptotically approach the stability boundaries for the corresponding continuous mass beam as the number of concentrated masses is increased. An error in the literature is also noted for the discrete mass problem concerning the behavior of the critical rotational speed for inplane buckling as the radius of rotation of the clamped end of the beam is reduced.

Lakin, William D.↗
NASA NTRS1994

Instability of streamwise vortices in plane channel flows

We present analysis and numerical experiments on the instability of streamwise vortices in 'minimal channel' flows and argue that this instability is a key feature in the observed intermittent cycle of formation, break-up, and re-formation of these structures. The base flow is a three-component, two-dimensional pair of counter-rotating rolls with axes aligned along the direction of the mean shear. While it is not a steady solution to the Navier-Stokes equations, we show numerically that this flow is unstable on a fast time scale to a secondary, three-dimensional Floquet mode. The growth of the secondary instability does not saturate in a new equilibrium, but continues until highly unstable local shear layers form and the entire flow breaks down into turbulence. Our analysis is motivated in part by the strong similarities between the intermittent turbulent cycle in minimal channel flows and one studied, both experimentally and in computations, in Couette-Taylor flow.

Coughlin, K.↗
NASA NTRS2000

Microgravity Research Results and Experiences from the NASA Mir Space Station Program

The Microgravity Research Program Office (MRPO) participated aggressively in Phase I of the International Space Station Program using the Russian Mir Space Station. The Mir Station offered an otherwise unavailable opportunity to explore the advantages and challenges to long duration microgravity space research. Payloads with both NASA and commercial backing were included as well as cooperative research with the Canadian Space Agency (CSA). From this experience, much was learned about dealing with long duration on orbit science utilization and developing new working relationships with our Russian partner to promote efficient planning, operations, and integration to solve complexities associated with a multiple partner program. Microgravity participation in the NASA Mir Program began with the first joint NASA Mir flight to the Mir Space Station. The earliest participation setup acceleration measurement capabilities that were used throughout the Program. Research, conducted by all Microgravity science disciplines, continued on each subsequent increment for the entire three-year duration of the Program. The Phase I Program included the Microgravity participation of over 30 Fluids, Combustion, Materials, and Biotechnology Sciences and numerous commercially sponsored research payloads. In addition to the research gained from Microgravity investigations, long duration operation of facility hardware was tested. Microgravity facilities operated on Mir included the Space Acceleration Measurement System (SAMS), the Microgravity Glovebox (MGBX), the Biotechnology System (BTS) and the Canadian Space Agency sponsored Microgravity Isolation Mount (MIM). The Russian OPTIZONE Furnace was also incorporated into our material science research. All of these efforts yielded significant and useful scientific research data. This paper focuses on the microgravity research conducted onboard the Mir space station. It includes the Program preparation and planning necessary to support this type of cross increment research experience; the payloads which were flown; and summaries of significant microgravity science findings. Most importantly this paper highlights the various disciplines of microgravity research conducted during the International Space Station, Phase 1 Program onboard the Mir Station. A capsulation of significant research and the applicability of our findings are provided. In addition, a brief discussion of how future microgravity science gathering capabilities, hardware development and payload operations techniques have enhanced our ability to conduct long duration microgravity research.

Schagheck, R. A.↗