Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Closed coronal structures. III - Comparison of static models with X-ray, EUV, and radio observations

Numerical models of static coronal loops in energy balance are compared with high spatial resolution observations of extreme ultraviolet lines, broad-band X-ray emission, and interferometric observations at 2.8 cm of a solar active region. Difficulties of using scaling laws to test static models of coronal loops are reviewed. The theoretical model used for the comparison is summarized; the detailed X-ray, EUV, and microwave observations of the selected active region are presented; and the comparison of the model with the observations is performed. It is shown that simple static models with conductive flux vanishing at the loop base reproduce satisfactorily the observed properties in the upper portion of loop structures from compact, high-pressure loops in the core of the region to more extended, fainter loops and to large-scale loops interconnecting different active regions. Effects of changing loop parameters are investigated, and it is argued, that in contrast to the present approach, scaling laws cannot be used to discriminate between different static energy balance models. Some discrepancy is found between model predictions and observations for the lower sections of loop structures. Possible causes of the discrepancy are discussed.

Pallavicini, R.↗
NASA NTRS1982

A frost formation model and its validation under various experimental conditions

A numerical model that was used to calculate the frost properties for all regimes of frost growth is described. In the first regime of frost growth, the initial frost density and thickness was modeled from the theories of crystal growth. The 'frost point' temperature was modeled as a linear interpolation between the dew point temperature and the fog point temperature, based upon the nucleating capability of the particular condensing surfaces. For a second regime of frost growth, the diffusion model was adopted with the following enhancements: the generalized correlation of the water frost thermal conductivity was applied to practically all water frost layers being careful to ensure that the calculated heat and mass transfer coefficients agreed with experimental measurements of the same coefficients.

Dietenberger, M. A.↗
NASA NTRS1984

An in-depth review of photovoltaic system performance models

The features, strong points and shortcomings of 10 numerical models commonly applied to assessing photovoltaic performance are discussed. The models range in capabilities from first-order approximations to full circuit level descriptions. Account is taken, at times, of the cell and module characteristics, the orientation and geometry, array-level factors, the power-conditioning equipment, the overall plant performance, O and M effects, and site-specific factors. Areas of improvement and/or necessary extensions are identified for several of the models. Although the simplicity of a model was found not necessarily to affect the accuracy of the data generated, the use of any one model was dependent on the application.

Smith, J. H.↗
NASA NTRS1987

Analytical modeling of helium II in forced flow conditions

A numerical model was developed to examine the flow behavior of He II under forced flow conditions. The model basically solves the governing equations in the two-fluid model of He II, and can be used to assess the pump requirement for transferring He II in spaceborne systems. Predictions of liquid temperature and pressure drop agree well with empirical data obtained from small scale tests. At high transfer rates, the model predicts that the forced convection dominates the Gorter-Mellink conduction in heat transfer, and the classical turbulent flow correlation is sufficient for calculating the pressure drop across the transfer line.

Ng, Y. S.↗
NASA NTRS2018

WFIRST Coronagraph Flight Performance Modeling

As it has for the past few years, numerical modeling is being used to predict the on-orbit, high-contrast imaging performance of the WFIRST coronagraph, which was recently defined to be solely a technology demonstrator. A consequence has been a realignment of modeling priorities and revised applications of modeling uncertainty factors and margins, which apply to multiple factors such as pointing and wavefront jitter, thermally-induced deformations, polarization, and aberration sensitivities. At the same time, the models have increased in fidelity as additional parameters have been added, such as time-dependent pupil shear and mid-spatial-frequency deformations of the primary and secondary mirrors, detector effects, and reaction-wheel-speed-dependent pointing and wavefront jitter.

Trauger, John↗
DOE OSTI2021

A Computational Information Criterion for Particle-Tracking with Sparse or Noisy Data

Traditional probabilistic methods for the simulation of advection-diffusion equations (ADEs) often overlook the entropic contribution of the discretization, e.g., the number of particles, within associated numerical methods. Many times, the gain in accuracy of a highly discretized numerical model is outweighed by its associated computational costs or the noise within the data. Herein, we address the question of how many particles are needed in a simulation to best approximate and estimate parameters in one-dimensional advective-diffusive transport. To do so, we use the well-known Akaike Information Criterion (AIC) and a recently-developed correction called the Computational Information Criterion (COMIC) to guide the model selection process. Random-walk and mass-transfer particle tracking methods are employed to solve the model equations at various levels of discretization. Numerical results demonstrate that the COMIC provides an optimal number of particles that can describe a more efficient model in terms of parameter estimation and model prediction compared to the model selected by the AIC even when the data is sparse or noisy, the sampling volume is not uniform throughout the physical domain, or the error distribution of the data is non-IID Gaussian.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Characterizing Spatial and Temporal Variability of California’s Coastal Upwelling Using a Satellite-Derived CUTI Index

Ocean upwelling off California is a wind-driven coastal process with complex interactions between spatial wind patterns and sea surface temperatures. Large-scale spatial and temporal upwelling characterizations are traditionally based on numerical models, while quantifying upwelling with direct observations is not straightforward. To complement model-based upwelling characterizations, here we use satellite estimates of winds, sea surface temperatures, and ocean currents along the California coast to calculate the coastal upwelling transport index (CUTI), originally calculated from a numerical ocean model. The CUTI framework accounts for both Ekman and geostrophic offshore-directed ocean transport, the latter of which is often neglected in satellite-based upwelling assessments. The approach yields the CUTI over a 30-yr period at a daily temporal and 0.25° spatial resolution to analyze upwelling patterns, with focus on potential Morro Bay and Humboldt wind energy areas. Morro Bay (CUTI median of 0.63 m 2 s −1 ) exhibits stronger, more consistent upwelling throughout the year, whereas Humboldt (CUTI median of 0.45 m 2 s −1 ) exhibits a pronounced seasonal cycle with strongest upwelling in the summer and occasional downwelling in winter. The cluster analysis of 1993–2022 satellite maps of CUTI identifies four characteristic spatial upwelling patterns along the coast, corresponding to distinct spatial distributions of sea surface temperature anomalies, Ekman transport, and geostrophic ocean water transport. This dataset provides an observation-based foundation for characterizing natural coastal upwelling off the California coast, showing good agreement with the original model CUTI while being based on satellite data and offering increased spatial resolution.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Meshfree simulation and experimental validation of extreme thermomechanical conditions in friction stir extrusion

Abstract Friction stir extrusion (FSE) is a novel solid-phase processing technique that consolidates and extrudes metal powders, flakes, chips, or billets into high-performance parts by plastic deformation, which has the potential to save substantial processing time and energy. Currently, most studies on FSE are experimental and only a few numerical models have been developed to explain and predict the complex physics of the process. In this work, a meshfree simulation framework based on smoothed particle hydrodynamics (SPH) was developed for FSE. Unlike traditional grid-based methods, SPH is a Lagrangian particle-based method that can handle severe material deformations, capture moving interfaces and surfaces, and monitor the field variable histories explicitly without complicated tracking schemes. These aspects of SPH make it attractive for the FSE process, where in situ evolution of field variables is difficult to observe experimentally. To this end, a 3-D, fully thermomechanically coupled SPH model was developed to simulate the FSE of aluminum wires. The developed model was thoroughly validated by comparing the numerically predicted material flow, strain, temperature history, and extrusion force with experimental results for a certain set of process parameters. The validated SPH model can serve as an effective tool to predict and better understand the extreme thermomechanical conditions during the FSE process.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Benchmark study of the DTU OWC chamber with both two-way and one-way absorption

This paper reports on a benchmark study based on small-scale (1:50) measurements of a single, oscillating water column chamber mounted sideways in a long flume. The geometry of the OWC chamber is extracted from a barge-like, attenuator-type floating concept “KNSwing” with 40 chambers targeted for deployment in the Danish part of the North Sea. In addition to traditional two-way energy extraction we also consider one-way energy extraction with passive venting and compare chamber response, pressures and total absorbed energy between the two methods. A blind study was established for the numerical modeling, with participants applying several implementations of weakly nonlinear potential flow theory and commercial Navier–Stokes solvers (CFD). Both compressible and incompressible models were used for the air phase. Potential flow calculations predict more energy absorption near the chamber resonance for one-way absorption than for two-way absorption, but the opposite is found from the experimental measurements. This outcome is mainly attributed to energy losses in the experimental passive valve system, but this conclusion must be confirmed by better experimental measurements. Modeling the one-way valve in CFD proved to be very challenging and only one team was able to provide results which were generally closer to the experiments. The study illustrates the challenges associated with both numerical and experimental analysis of OWC chambers. Air compressibility effects were not found to be important at this scale, even with the large volume of additional air used for the one-way case.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Infrasonic directivity of monopole, dipole and bipole ground-surface reflected sources

Infrasound (acoustic waves below 20 Hz) can be used to detect, locate and quantify activity in the atmosphere such as volcanic eruptions and anthropogenic explosions. Attempts to quantify volcanic eruption parameters such as exit velocity, plume height and mass flow rate using infrasound data depend strongly on assumptions of the acoustic source type. Infrasonic sources may produce omnidirectional or directional wavefields, while propagation effects, such as interaction with topography, can induce further wavefield directivity that is measured by field instrumentation. Limited sampling of these wavefields can hinder our ability to infer the underlying source, and thus our understanding of the eruption characteristics. Equivalent sources are often used to represent acoustic source mechanisms and resultant wavefields. In this study, we review equivalent acoustic sources as they pertain to infrasonic scale and wavelengths commonly encountered in very local (⁠<5 km range) geophysical field deployments. We highlight the equivalent infrasonic bipole source that can be induced by ground-reflection of an elevated monopole; we are not aware of any prior infrasound studies that use the bipole source concept. We use analytical and numerical methods to explore source directivity of monopole, dipole and bipole ground-reflected sources at infrasonic frequencies as well as the additional directivity complications introduced by interactions with topography. We illustrate that for typical volcano-infrasound wavelengths, increasing height above the ground as well as increasing source frequency leads to increased wavefield directivity. Numerical modelling using a simple omnidirectional monopole source embedded in topography further illustrates that both horizontal and vertical infrasound directionality can be induced by topography at the distance scales appropriate for local volcano infrasound monitoring. Information summarized in this analytical and numerical exploration of infrasound directivity may be used to help guide future volcano-infrasound field deployments intended to estimate source parameters or quantify wavefield directivity. Analytic solutions for simple whole-space or half-space atmospheres provide useful formulations for planning or initially analysing geophysical field-scale experimental data; however, especially at very local distances from the source (⁠<5 km), 3-D simulations are necessary to account for complex topography commonly encountered in volcano-infrasound applications.

Infrasound↗
NASA NTRS2008

Towards Improved Forecasts of Atmospheric and Oceanic Circulations over the Complex Terrain of the Eastern Mediterranean

Forecasting atmospheric and oceanic circulations accurately over the Eastern Mediterranean has proved to be an exceptional challenge. The existence of fine-scale topographic variability (land/sea coverage) and seasonal dynamics variations can create strong spatial gradients in temperature, wind and other state variables, which numerical models may have difficulty capturing. The Hellenic Center for Marine Research (HCMR) is one of the main operational centers for wave forecasting in the eastern Mediterranean. Currently, HCMR's operational numerical weather/ocean prediction model is based on the coupled Eta/Princeton Ocean Model (POM). Since 1999, HCMR has also operated the POSEIDON floating buoys as a means of state-of-the-art, real-time observations of several oceanic and surface atmospheric variables. This study attempts a first assessment at improving both atmospheric and oceanic prediction by initializing a regional Numerical Weather Prediction (NWP) model with high-resolution sea surface temperatures (SST) from remotely sensed platforms in order to capture the small-scale characteristics.

Chronis, Themis↗
DOE OSTI2025

A modeling study of ocean thermal energy conversion resource and potential environmental effects around Kailua-Kona, Hawaii

Ocean Thermal Energy Conversion (OTEC) offers a promising renewable energy solution through a heat exchange process using the temperature difference between warm surface seawater and cold deep seawater. Because accurate resource characterization is critical for the optimal design and implementation of OTEC systems, a high-resolution numerical model is employed to better characterize the OTEC resource at Kona, Hawaii. Our model provides detailed spatial and temporal variability of the thermal gradient, which is essential for assessing the viability and efficiency of OTEC systems. The model results reveal distinct patterns and dynamics not captured by existing observations or models (e.g., lower-resolution information). These findings highlight the importance of using high-resolution models for accurate predictions of thermal gradient variability, ultimately supporting more efficient and sustainable OTEC deployment. Additionally, the study investigates the impacts of mixed water discharge from OTEC plants that can cause shock to organisms living in the surface water and potentially destabilize the water column. Understanding these effects is vital for minimizing any potential negative environmental consequences and ensuring the long-term viability of OTEC operations. Further, our model improves OTEC resource characterization, which can lead to optimal design and deployment of OTEC systems. The analysis of OTEC water discharge impacts can accelerate the development of OTEC technologies, overcoming permitting/consenting challenges. These findings contribute to the broader adoption of high-resolution modeling in ocean energy resource characterization, particularly for OTEC applications.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS2016

Using 3-D Numerical Weather Data in Piloted Simulations

This report describes the process of acquiring and using 3-D numerical model weather data sets in NASA Langley's Research Flight Deck (RFD). A set of software tools implement the process and can be used for other purposes as well. Given time and location information of a weather phenomenon of interest, the user can download associated numerical weather model data. These data are created by the National Oceanic and Atmospheric Administration (NOAA) High Resolution Rapid Refresh (HRRR) model, and are then processed using a set of Mathworks' Matlab(TradeMark) scripts to create the usable 3-D weather data sets. Each data set includes radar re ectivity, water vapor, component winds, temperature, supercooled liquid water, turbulence, pressure, altitude, land elevation, relative humidity, and water phases. An open-source data processing program, wgrib2, is available from NOAA online, and is used along with Matlab scripts. These scripts are described with sucient detail to make future modi cations. These software tools have been used to generate 3-D weather data for various RFD experiments.

Daniels, Taumi S.↗
NASA NTRS1993

An adaptive multigrid model for hurricane track prediction

This paper describes a simple numerical model for hurricane track prediction which uses a multigrid method to adapt the model resolution as the vortex moves. The model is based on the modified barotropic vorticity equation, discretized in space by conservative finite differences and in time by a Runge-Kutta scheme. A multigrid method is used to solve an elliptic problem for the streamfunction at each time step. Nonuniform resolution is obtained by superimposing uniform grids of different spatial extent; these grids move with the vortex as it moves. Preliminary numerical results indicate that the local mesh refinement allows accurate prediction of the hurricane track with substantially less computer time than required on a single uniform grid.

Fulton, Scott R.↗
DOE OSTI2024

Performance on HPC Platforms Is Possible Without C++

Computing at large scales has become extremely challenging due to increasing heterogeneity in both hardware and software. More and more scientific workflows must tackle a range of scales and use machine learning and AI intertwined with more traditional numerical modeling methods, placing more demands on computational platforms. These constraints indicate a need to fundamentally rethink the way computational science is done and the tools that are needed to enable these complex workflows. The current set of C++-based solutions may not suffice, and relying exclusively upon C++ may not be the best option, especially because several newer languages and boutique solutions offer more robust design features to tackle the challenges of heterogeneity. In June 2023, we held a mini symposium that explored the use of newer languages and heterogeneity solutions that are not tied to C++ and that offer options beyond template metaprogramming and Parallel. For for performance and portability. In conclusion, we describe some of the presentations and discussion from the mini symposium in this article.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Development, Monitoring, and Control of Fracture Thermal Energy Storage in Crystalline Rock Formations (DEMO-FTES) [Abstract]

Approximately half of global energy consumption is used for heating and cooling. Because fossil fuels are used to meet most of this demand, heating and cooling of buildings produces a large portion of global greenhouse gas (GHG) emissions. The proposed project seeks to demonstrate the feasibility of a seasonal thermal energy storage concept called fracture thermal energy storage (FTES), which has the potential to dramatically lower energy requirements for heating and cooling and improve the resilience of building energy systems. FTES is a technique for building a highly efficient heat exchanger by creating a carefully designed set of fractures in the ground below a building. This heat exchanger allows thermal energy to be stored over seasonal timeframes—for example, hot thermal energy that is easily captured in the summer or energy from waste heat sources—to meet heating needs during colder months. FTES offers a solution with a relatively small footprint and lower cost compared to currently operating aquifer thermal energy storage (ATES) systems and could be used in the many areas without aquifers suitable for ATES. If commercialized, FTES systems could exceed the 2.5 TWh of energy storage per year from the more than 2,800 ATES systems in operation worldwide, which range from 0.1 MW to 30 MW. The savings in CO2 emissions per year are also expected to match or exceed individual ATES projects, with the largest, a 30 MW system used to provide heating and cooling to the University of Technology in Eindhoven in the Netherlands, estimated to save 13,000 tons of CO2 emissions per year. FTES utilizes mature drilling and fracturing technology and therefore has the potential to be rapidly commercialized once demonstrated. The ability to construct and establish flow through an FTES heat exchanger has been demonstrated by a previous project and the potential for efficient, large-capacity energy storage has been shown using numerical models. However, no experimental validation of these numerical estimates of thermal energy storage has been made. The first crucial need to advance FTES technology is to identify the sensitivity of key metrics such as thermal energy storage and production rates, capacities, and efficiencies to design parameters such as the number of fractures, depth/temperature of fractures, size of fractures, and circulation rates. The second crucial need is experimental testing of achievable thermal performance with optimized system design parameters. The proposed scope of work seeks to systematically address these two critical needs through a highly complementary international collaboration spanning theory, laboratory, and mesoscale field evaluation. The proposed work plan calls for using dimensional analysis and existing state-of-the-art numerical simulators to design carefully scaled laboratory and 10-meter-scale field tests of the thermal efficiency of FTES. The existing advanced laboratory and intermediate-scale field testbeds that will be used for this project will allow for detailed monitoring of the system performance during the test and of how the performance changes across time and length scales. These results will determine the feasibility of full- scale FTES systems. If the thermal performance is consistent with model predictions, the results will provide a strong economic justification for rapid commercialization of FTES technology in a wide range of geographical areas.

25 ENERGY STORAGE↗
NASA NTRS2003

A New Objective Technique for Verifying Mesoscale Numerical Weather Prediction Models

This report presents a new objective technique to verify predictions of the sea-breeze phenomenon over east-central Florida by the Regional Atmospheric Modeling System (RAMS) mesoscale numerical weather prediction (NWP) model. The Contour Error Map (CEM) technique identifies sea-breeze transition times in objectively-analyzed grids of observed and forecast wind, verifies the forecast sea-breeze transition times against the observed times, and computes the mean post-sea breeze wind direction and speed to compare the observed and forecast winds behind the sea-breeze front. The CEM technique is superior to traditional objective verification techniques and previously-used subjective verification methodologies because: It is automated, requiring little manual intervention, It accounts for both spatial and temporal scales and variations, It accurately identifies and verifies the sea-breeze transition times, and It provides verification contour maps and simple statistical parameters for easy interpretation. The CEM uses a parallel lowpass boxcar filter and a high-order bandpass filter to identify the sea-breeze transition times in the observed and model grid points. Once the transition times are identified, CEM fits a Gaussian histogram function to the actual histogram of transition time differences between the model and observations. The fitted parameters of the Gaussian function subsequently explain the timing bias and variance of the timing differences across the valid comparison domain. Once the transition times are all identified at each grid point, the CEM computes the mean wind direction and speed during the remainder of the day for all times and grid points after the sea-breeze transition time. The CEM technique performed quite well when compared to independent meteorological assessments of the sea-breeze transition times and results from a previously published subjective evaluation. The algorithm correctly identified a forecast or observed sea-breeze occurrence or absence 93% of the time during the two- month evaluation period from July and August 2000. Nearly all failures in CEM were the result of complex precipitation features (observed or forecast) that contaminated the wind field, resulting in a false identification of a sea-breeze transition. A qualitative comparison between the CEM timing errors and the subjectively determined observed and forecast transition times indicate that the algorithm performed very well overall. Most discrepancies between the CEM results and the subjective analysis were again caused by observed or forecast areas of precipitation that led to complex wind patterns. The CEM also failed on a day when the observed sea- breeze transition affected only a very small portion of the verification domain. Based on the results of CEM, the RAMS tended to predict the onset and movement of the sea-breeze transition too early and/or quickly. The domain-wide timing biases provided by CEM indicated an early bias on 30 out of 37 days when both an observed and forecast sea breeze occurred over the same portions of the analysis domain. These results are consistent with previous subjective verifications of the RAMS sea breeze predictions. A comparison of the mean post-sea breeze winds indicate that RAMS has a positive wind-speed bias for .all days, which is also consistent with the early bias in the sea-breeze transition time since the higher wind speeds resulted in a faster inland penetration of the sea breeze compared to reality.

Case, Jonathan L.↗