Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical model”

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 811 records · Page 45

Hydraulic characterization of a fault zone from fracture distribution

A quantitative assessment of how faults control the migration of geofluids is critical in many areas of geosciences. Here we integrated geological fieldwork, quantitative analysis of the fractures distribution and numerical modeling to build a geometrical representation of a fault zone and to characterize its hydraulic properties. Our target is a fault located in the Majella Mountain (Italy). We collected 21 scan lines across the fault profile in order to characterize its architecture. The numerical modeling of the fracture network of the damage zones and their hydraulic parameters was performed using both commercial (Move ® ) and open source software (dfnWorks and PFLOTRAN). Move ® was used to build a representative model of the fault zone using fracture spacing as a proxy, and to model the hydraulic parameters of the different fault domains. dfnWorks and PFLOTRAN were employed to infer the hydraulic parameters of the damage zones of the fault and then upscale these properties to an equivalent continuum domain, suitable for fluid flow simulations through the whole fault zone. Our findings show how even in a relatively small area it is possible to describe changes in terms of hydraulic properties of a fault zone and to build models capable to represent these variations.

58 GEOSCIENCES↗

Implementation of a Model Predictive Control Strategy to Regulate Temperature Inside Plug-Flow Solar Reactor With Countercurrent Flow

Abstract Solar-driven thermochemical energy storage systems are proven to be promising energy carriers (solar fuels) to utilize solar energy by using reactive solid-state pellets. However, the production of solar fuel requires a quasi-steady-state process temperature, which represents the main challenge due to the transient nature of solar power. In this work, an adaptive model predictive controller (MPC) is presented to regulate the temperature inside a tubular solar reactor to produce solid-state solar fuel for long-term thermal storage systems. The solar reactor system consists of a vertical tube heated circumferentially over a segment of its length by concentrated solar power, and the reactive pellets (MgMn2O4) are fed from the top end and flow downwards through the heated tube. A countercurrent flowing gas supplied from the lower end interacts with flowing pellets to reduce it thermochemically at a temperature range of 1000—1500 °C. A low-order physical model was developed to simulate the dynamics of the solar reactor including the reaction kinetics, and the proposed model was validated numerically by using a 7-kW electric furnace. The numerical model then was utilized to design the MPC controller, where the control system consists of an MPC code linked to an adaptive system identification code that updates system parameters online to ensure system robustness against external disturbances (sudden change in the flow inside the reactor), model mismatches, and uncertainty. The MPC controller parameters are tuned to enhance the system performance with minimum steady-state error and overshoot. The controller is tested to track different temperature ranges between 500 °C and 1400 °C with different particles/gas mass flowrates and ramping temperature profiles. Results show that the MPC controller successfully regulated the reactor temperature within ± 1 °C of its setpoint and maintained robust performance with minimum input effort when subjected to sudden changes in the amount of flowing media and the presence of chemical reaction.

Engineering↗

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗

Evaluation of Tritium Behavior in Forest Vegetation for the Purpose of Determining Appropriate Non-Zero Deposition Velocities for Tritium Oxide

Modeling and field measurements have been conducted between 2019 and 2021 to assess the practicality of using a non-zero deposition velocity for safety basis estimates of tritium oxide fate and transport modeling. A model was developed which used a complex deposition algorithm designed to assess how tritium oxide would mix based on the turbulent motions and wind speed effects that the forest canopy has on the atmosphere. The model was driven based on measurements of wind and turbulence taken from the Aiken AmeriFlux Tower which measures these properties at five levels located within and just above the forest. The model was then validated against a series of field experiments which were designed to test the model predictions and estimate the deposition velocity occurring over the forest environment at the Savannah River Site. The field releases used deuterium oxide as a surrogate for tritium oxide and was released as a fine mist which rapidly evaporated, creating a gaseous tracer in the atmosphere. Using air samplers, the elevations in deuterium concentration in the air relative to background measurements was assessed and then modeled. Generally, the numerical model tended to underpredict the amount of deuterium being mixed from above the canopy to the forest floor, indicating that the predictions it provides are still conservative relative to what was measured during the field experiments. Across a suite of modeling runs, the 95 th and 99 th percentile deposition velocities were estimated to be 1.2 and 0.7 cm s -1 , respectively. Estimated deposition velocities in the 2021 field experiments, which specifically assessed a release above the forest canopy and its mixing to the surface, predicted deposition velocities ranging from 1.75 to 6.61 cm s -1 . While these field releases do not cover all possible meteorological conditions, it seems appropriate to use a non-zero deposition velocity when performing safety-basis modeling of tritium oxide. The recommendation presented in this report is to use 1.0 cm s -1 . This is between the 95 th and 99 th percentile value estimated from the modeling study, suggesting it should be appropriate for the majority of release scenarios given the model’s apparent conservatism relative to field measurements.

54 ENVIRONMENTAL SCIENCES↗

Analyzing Mechanical integrity of Microwave resonant Cavity flowmeter

The microwave resonant cavity flow sensor has been proposed to address the difficulties of measuring high-temperature molten salt flow in reactor vessels. The flow sensor consists of a metallic cylindrical cavity coupled to a microwave waveguide. The cylindrical cavity has one of the flat surfaces thin enough to undergo microscopic deflection due to dynamic pressure of the exerted by the flowing fluid. Deflection of the membrane causes a shift in the microwave resonant frequency. Fluid velocity can be estimated from the calibration curve. The objective of this report is to investigate mechanical integrity of the flowmeter’s membrane. The membrane is a thin stainless-steel circular plate with a clamped edge, subjected to uniform load. By calculating the stress on the plate due to deflection and, comparing the stress to the material ultimate tensile strength and yield strength, it can be estimated if the plate will fail. If the maximum stress is greater than the yield strength, the plate will undergo plastic deformation, and the sensor will be disabled. A stress greater than the ultimate tensile strength will cause the plate to fracture. The stress on the plate was calculated with an analytic closed form solution model, and with COMSOL Structural Mechanics Module which does not involve any approximations. Stresses in both the analytic model and COMSOL numerical model were calculated for stainless steel 316 thin circular plates under room temperature conditions, where the fluid is water. The diameter of the circular plate is approximately one inch, and thickness is 10mil. Fluid flow velocity is in the range from 0.5m/s to 1m/s. The results from the two models are similar. Both the analytic model and COMSOL model showed that maximum stresses on the plate, which are at the radial boundary of the plate, are three orders of magnitude smaller than the yield strength and ultimate tensile strength. This indicates that the sensor is at a low risk of mechanical failure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Microwave brightness temperature of a windblown sea

A mathematical model is developed for the apparent temperature of the sea at all microwave frequencies. The model is a numerical model in which both the clear water structure and white water are accounted for as a function of wind speed. The model produces results similar to Stogryn's model at 19.35 GHz for wind speeds less than 8 m/sec; it can use radiosonde data to calculate atmospheric effects and can incorporate an empirically determined antenna gain pattern. The corresponding computer program is of modular design and the logic of the main program is capable of treating a horizontally inhomogeneous surface or atmosphere. It is shown that a variation of microwave brightness temperature with zenith angle is necessary to produce the wind sensitivity of the horizontally polarized brightness temperature; the variation of sky temperature with frequency is sufficient to produce a frequency dependent wind sensitivity.

Hall, F. G.↗

Martian thermal boundary layers: Subhourly variations induced by radiative-conductive heat transfer within the dust-laden atmosphere-ground system

Radiative-conductive heat transfer has been investigated for the ground-atmosphere system of the planet Mars. The basic goal was the quantitative determination of time dependent vertical distributions of temperature and static stability for Southern-Hemispheric summer season and middle and polar latitudes, for both dust-free and dust-laden atmospheric conditions. The numerical algorithm which models at high spatial and temporal resolution the thermal energy transports in the dual ground-atmosphere system, is based on solution of the applicable heating rate equation, including radiative and molecular-conductive heat transport terms. The two subsystems are coupled by an internal thermal boundary condition applied at the ground-atmosphere interface level. Initial data and input parameters are based on Mariner 4, 6, 7, and 9 measurements and the JPL Mars Scientific Model. Numerical experiments were run for dust-free and dust-laden conditions in the midlatitudes, as well as ice-free and ice-covered polar regions. Representative results and their interpretation are presented. Finally, the theoretical framework of the generalized problem with nonconservative Mie scattering and explicit thermal-convective heat transfer is formulated, and applicable solution algorithms are outlined.

Pallmann, A. J.↗

Dynamics of the North Sea pole tide reconsidered

The anomalous tide in the North Sea, driven by the Chandler wobble, is here re-examined. A previously published solution in which these 'pole tide' observations were explained as a forced co-oscillation from the deep water tide is shown to be untenable, because it is very unlikely that the deep ocean tide is anything but equilibrium. Using the same physics as in the previous solution (viscous, linear, quasi-geostrophic dynamics), it is now suggested that the only plausible explanation is that there is a coincidental resonance between the long-wavelength branch topographic Rossby wave of the North Sea and the direct driving by the pole tide potential. This resonance is not a basin-resonance, but that of a free wave driven by a travelling forcing function. Simplified models reproduce the observations in qualitative fashion. If this explanation is the correct one, then the Arctic seas north of Asia and North America appear to be likely candidates as other regions where there may be a large pole tide response and, hence, a locus of wobble dissipation. Given the difficulties of observation, the long times required, and the general instability of the dissipation calculation with the simplified analytical model, probably a numerical model will be required for further progress.

Wunsch, Carl↗

Microphysical fundamentals governing cirrus cloud growth: Modeling studies

For application to Global Climate Models, large scale numerical models of cirrus cloud formation and maintenance need to be refined to more reliably simulate the effects and feedbacks of high level clouds. A key aspect is how ice crystal growth is initiated in cirrus, which has started a cloud microphysical controversy between camps either believing that heterogeneous or homogeneous drop freezing is predominantly responsible for cold cirrus ice crystal nucleation. In view of convincing evidence for the existence of highly supercooled cloud droplets in the middle and upper troposphere, however, it is concluded that active ice nuclei are rather scarce at cirrus cloud altitudes, and so a new understanding of cirrus cloud formation is needed. This understanding is sought through an examination of cirrus cloud growth models.

Sassen, Kenneth↗

The impact of mesoscale data on the simulation of a mesoscale convective weather system

The objectives are: (1) to demonstrate the sensitivity of a numerical simulation of a mesoscale convective weather system to the initial conditions; and (2) to provide further evidence of the need for a high resolution observing system that is compatible with numerical-model initial data requirements. To this end, a series of nine numerical-model sensitivity experiments were conducted in which one or more variables from one or more observations (soundings) were omitted and/or adjusted in the model initialization. Two types of surroundings were available for manipulation in the sensitivity experiments. Specifically, for the particular event being simulated (the 1977 Johnstown flash flood), an independent mesoanalysis was available from Bosart and Sanders (1981). In his analysis, Bosart produced a fine-resolution three-dimensional array of data from which soundings that helped to define mesoscale features could be extracted. These 'nonconventional' soundings were added to the conventional sounding data routinely available for model initialization.

Fritsch, J. M.↗

Summary of Cumulus Parameterization Workshop

A workshop on cumulus parameterization took place at the NASA Goddard Space Flight Center from December 3-5, 2001. The major objectives of this workshop were (1) to review the problem of representation of moist processes in large-scale models (mesoscale models, Numerical Weather Prediction models and Atmospheric General Circulation Models), (2) to review the state-of-the-art in cumulus parameterization schemes, and (3) to discuss the need for future research and applications. There were a total of 31 presentations and about 100 participants from the United States, Japan, the United Kingdom, France and South Korea. The specific presentations and discussions during the workshop are summarized in this paper.

Tao, Wei-Kuo↗

The Mars Dust Cycle: Investigating the Effects of Radiatively Active Water Ice Clouds on Surface Stresses and Dust Lifting Potential with the NASA Ames Mars General Circulation Model

The dust cycle is a critically important component of Mars' current climate system. Dust is present in the atmosphere of Mars year-round but the dust loading varies with season in a generally repeatable manner. Dust has a significant influence on the thermal structure of the atmosphere and thus greatly affects atmospheric circulation. The dust cycle is the most difficult of the three climate cycles (CO2, water, and dust) to model realistically with general circulation models. Until recently, numerical modeling investigations of the dust cycle have typically not included the effects of couplings to the water cycle through cloud formation. In the Martian atmosphere, dust particles likely provide the seed nuclei for heterogeneous nucleation of water ice clouds. As ice coats atmospheric dust grains, the newly formed cloud particles exhibit different physical and radiative characteristics. Thus, the coupling between the dust and water cycles likely affects the distributions of dust, water vapor and water ice, and thus atmospheric heating and cooling and the resulting circulations. We use the NASA Ames Mars GCM to investigate the effects of radiatively active water ice clouds on surface stress and the potential for dust lifting. The model includes a state-of-the-art water ice cloud microphysics package and a radiative transfer scheme that accounts for the radiative effects of CO2 gas, dust, and water ice clouds. We focus on simulations that are radiatively forced by a prescribed dust map, and we compare simulations that do and do not include radiatively active clouds. Preliminary results suggest that the magnitude and spatial patterns of surface stress (and thus dust lifting potential) are substantial influenced by the radiative effects of water ice clouds.

Kahre, Melinda A.↗

Tensile Creep Behavior of SiC(f)/SiC Ceramic Matrix Minicomposites

Single fiber-tow minicomposites represent the major load-bearing element of woven and laminate ceramic matrix composites (CMCs). To understand the effects of fiber type, fiber content, and matrix cracking on tensile creep in SiCf/SiC CMCs, single-tow SiCf/SiC minicomposites with different fiber types and contents were investigated. The minicomposites studied contained either Hi-Nicalon™ or Hi-Nicalon™ Type S SiC fibers with a boron nitride (BN) interphase and a chemical-vapor-infiltrated-silicon-carbide (CVI-SiC) matrix. Tensile creep was performed at 1200 °C in air. A bottom-up creep modeling approach was applied where creep parameters of the fibers and matrix were obtained separately at 1200 °C. Next, a theoretical model based on the rule of mixtures was derived to model the fiber and matrix creep-time-dependent stress redistribution. Fiber and matrix creep parameters, load transfer model results, and numerical modeling were used to construct a creep strain model to predict creep damage evolution of minicomposites with different fiber types and contents.

SiC/SiC Minicomposites Creep↗

Progressive collapse analysis of stiffened panel structure with consideration of actual welding distortion and residual stresses

To date, inclusion of welding distortion and residual stresses in determining the strength of marine structures is based on empirical formulations due to welding complexity. In this study, a more feasible numerical model was developed from the perspective of inherent strain theory. The accuracy of this simplified numerical model was confirmed by transient thermo-elastic–plastic finite element analysis using a novel dynamic mesh refining and iterative substructure method which is 50 times faster than conventional finite element method. Welding distortion by the two numerical approaches was validated with experimental measurement. The collapse analyses indicated that the inherent strain model predicted a compressive strength and load-displacement curve with high accuracy, as validated by the transient thermo-mechanical model. The global stiffness and the ultimate strength of the welded structure would be overestimated by about 22% if no welding effects were considered. The proposed computation procedure for determining the ultimate strength of welded structures is accurate, simple, and practical and enables improvement of ultimate strength through optimization of welding joint design to reduce the welding residual stress and distortion.

Huang, Hui↗

Retrieval of Microphysical Characteristics of Particles in Atmospheres of Distant Comets from Ground-Based Polarimetry

We summarize unique aperture data on the degree of linear polarization observed for distant comets C/2010 S1, C/2010 R1, C/2011 KP36, C/2012 J1, C/2013 V4, and C/2014 A4 with heliocentric distances exceeding 3 AU. Observations have been carried out at the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences (Nizhnij Arkhyz, Russia) during the period from 2011 to 2016. The measured negative polarization proves to be significantly larger in absolute value than what is typically observed for comets close to the Sun. We compare the new observational data with the results of numerical modeling performed with the T-matrix and superposition T-matrix methods. In our computer simulations, we assume the cometary coma to be an optically thin cloud containing particles in the form of spheroids, fractal aggregates composed of spherical monomers, and mixtures of spheroids and aggregate particles. We obtain a good semi-quantitative agreement between all polarimetric data for the observed distant comets and the results of numerical modeling for the following models of the cometary dust: (i) a mixture of submicrometer water-ice oblate spheroids with aggregates composed of submicrometer silicate monomers; and (ii) a mixture of submicrometer water-ice oblate spheroids and aggregates consisting of both silicate and organic monomers. The microphysical parameters of these models are presented and discussed.

Numerical modeling↗

Solar Physics

The areas of emphasis are: (1) develop theoretical models of the transient release of magnetic energy in the solar atmosphere, e.g., in solar flares, eruptive prominences, coronal mass ejections, etc.; (2) investigate the role of the Sun's magnetic field in the structuring of solar corona by the development of three-dimensional numerical models that describe the field configuration at various heights in the solar atmosphere by extrapolating the field at the photospheric level; (3) develop numerical models to investigate the physical parameters obtained by the ULYSSES mission; (4) develop numerical and theoretical models to investigate solar activity effects on the solar wind characteristics for the establishment of the solar-interplanetary transmission line; and (5) develop new instruments to measure solar magnetic fields and other features in the photosphere, chromosphere transition region and corona. We focused our investigation on the fundamental physical processes in solar atmosphere which directly effect our Planet Earth. The overall goal is to establish the physical process for the Sun-Earth connections.

Wu, S. T.↗

Overview of the Field Phase of the NASA Tropical Cloud Systems and Processes (TCSP)Experiment

The Tropical Cloud Systems and Processes experiment is sponsored by the National Aeronautics and Space Administration (NASA) to investigate characteristics of tropical cyclone genesis, rapid intensification and rainfall using a three-pronged approach that emphasizes satellite information, suborbital observations and numerical model simulations. Research goals include demonstration and assessment of new technology, improvements to numerical model parameterizations, and advancements in data assimilation techniques. The field phase of the experiment was based in Costa Rica during July 2005. A fully instrumented NASA ER-2 high altitude airplane was deployed with Doppler radar, passive microwave instrumentation, lightning and electric field sensors and an airborne simulator of visible and infrared satellite sensors. Other assets brought to TCSP were a low flying uninhabited aerial vehicle, and a surface-based radiosonde network. In partnership with the Intensity Forecasting Experiment of the National Oceanic and Atmospheric Administration (NOAA) Hurricane Research Division, two NOAA P-3 aircraft instrumented with radar, passive microwave, microphysical, and dropsonde instrumentation were also deployed to Costa Rica. The field phase of TCSP was conducted in Costa Rica to take advantage of the geographically compact tropical cyclone genesis region of the Eastern Pacific Ocean near Central America. However, the unusual 2005 hurricane season provided numerous opportunities to sample tropical cyclone development and intensification in the Caribbean Sea and Gulf of Mexico as well. Development of Hurricane Dennis and Tropical Storm Gert were each investigated over several days in addition to Hurricane Emily as it was close to Saffir-Simpson Category 5 intensity. An overview of the characteristics of these storms along with the pregenesis environment of Tropical Storm Eugene in the Eastern Pacific will be presented.

Hood, Robbie E.↗

Piezoelectric Vibration Damping Study for Rotating Composite Fan Blades

Resonant vibrations of aircraft engine blades cause blade fatigue problems in engines, which can lead to thicker and aerodynamically lower performing blade designs, increasing engine weight, fuel burn, and maintenance costs. In order to mitigate undesirable blade vibration levels, active piezoelectric vibration control has been investigated, potentially enabling thinner blade designs for higher performing blades and minimizing blade fatigue problems. While the piezoelectric damping idea has been investigated by other researchers over the years, very little study has been done including rotational effects. The present study attempts to fill this void. The particular objectives of this study were: (a) to develop and analyze a multiphysics piezoelectric finite element composite blade model for harmonic forced vibration response analysis coupled with a tuned RLC circuit for rotating engine blade conditions, (b) to validate a numerical model with experimental test data, and (c) to achieve a cost-effective numerical modeling capability which enables simulation of rotating blades within the NASA Glenn Research Center (GRC) Dynamic Spin Rig Facility. A numerical and experimental study for rotating piezoelectric composite subscale fan blades was performed. It was also proved that the proposed numerical method is feasible and effective when applied to the rotating blade base excitation model. The experimental test and multiphysics finite element modeling technique described in this paper show that piezoelectric vibration damping can significantly reduce vibrations of aircraft engine composite fan blades.

Min, James B.↗