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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Verification of numerical models for seismic fluid-structure interaction analysis of internal components in liquid-filled advanced reactors

Earthquake shaking of a liquid-filled advanced reactor induces fluid-structure interaction (FSI) between the reactor vessel (tank), its internal components, and the contained liquid. Seismic design, qualification, and risk assessment of such reactors must consider fluid-structure responses, their geometries and support conditions, and three-directional seismic inputs, all of which require numerical simulations. Numerical models used for analysis of safety-related nuclear equipment must be verified and validated. Here, this paper verifies numerical models of submerged components using analytical solutions. Seismic FSI analysis of the numerical models is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Prior analytical solutions are reworked in this paper, and calculation errors are identified and corrected. The solutions address frequencies of two concentric cylindrical pipes filled with liquid: the inner pipe is a submerged component, and the outer pipe contains the liquid. A unitless frequency coefficient is defined and presented for pipes of different materials, filled with different liquids, and with a range of dimensions for application to advanced reactors. The numerical models are verified here by comparing the lateral frequencies of submerged components with those calculated using the corrected analytical solutions. Recommendations for verification of numerical models of internal components in advanced reactors for seismic FSI analysis are provided. Although the unitless frequency coefficients and verification procedures are developed for application to advanced nuclear reactors, they are broadly applicable to FSI analysis of submerged components in liquid-filled vessels such as storage tanks, boilers, and steam generators.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1991

GEOSIM: A numerical model for geophysical fluid flow simulation

A numerical model which simulates geophysical fluid flow in a wide range of problems is described in detail, and comparisons of some of the model's results are made with previous experimental and numerical studies. The model is based upon the Boussinesq Navier-Stokes equations in spherical coordinates, which can be reduced to a cylindrical system when latitudinal walls are used near the pole and the ratio of latitudinal length to the radius of the sphere is small. The equations are approximated by finite differences in the meridional plane and spectral decomposition in the azimuthal direction. The user can specify a variety of boundary and initial conditions, and there are five different spectral truncation options. The results of five validation cases are presented: (1) the transition between axisymmetric flow and baroclinic wave flow in the side heated annulus; (2) the steady baroclinic wave of the side heated annulus; (3) the wave amplitude vacillation of the side heated annulus; (4) transition to baroclinic wave flow in a bottom heated annulus; and (5) the Spacelab Geophysical Fluid Flow Cell (spherical) experiment.

Butler, Karen A.↗

DOE OSTI2022

Combustion Characteristics of a Tier II Gasoline Certification Fuel and its Surrogate with Iso-butanol: experiments and detailed numerical modeling

This presentation will describe results of an experimental and numerical study of combustion of a Tier II gasoline certification fuel and its surrogate blended with iso-butanol as a potential biofuel additive. The experiments were carried out under conditions where gas transport originated entirely by fuel evaporation with external convection (both forced and natural) were minimized. The resulting near spherically symmetric droplet flames were amenable to detailed numerical modeling. The droplet initial diameters were in the narrow range of 0.55 mm to 0.62 mm while and burning is performed under standard atmosphere. The combustion characteristics such as evolution of droplet, flame and soot shell diameters for indolene/iso-butanol and surrogate/iso-butanol mixtures are compared. The presence of iso-butanol is experimentally found to reduce the formation of soot as the iso-butanol concentration was increased. Simulations are also reported using a detailed kinetic chemistry comprised of 324 species and 17608 reactions using the CRECK POLIMI mechanism. The detailed numerical model included unsteady gas and liquid transport, soot chemistry and radiative transport. Droplet burn rates were well predicted while flame diameters showed more variability depending on how the flame was defined in the numerical simulation (peak gas temperature and OH concentration). Simulations also show the influence of iso-butanol concentration on flame temperature, greenhouse gas emissions and particulate emissions.

Guo, Songtao↗

NASA NTRS1997

Numerical Modeling of HgCdTe Solidification: Effects of Phase Diagram, Double-Diffusion Convection and Microgravity Level

Melt convection, along with species diffusion and segregation on the solidification interface are the primary factors responsible for species redistribution during HgCdTe crystal growth from the melt. As no direct information about convection velocity is available, numerical modeling is a logical approach to estimate convection. Furthermore influence of microgravity level, double-diffusion and material properties should be taken into account. In the present study, HgCdTe is considered as a binary alloy with melting temperature available from a phase diagram. The numerical model of convection and solidification of binary alloy is based on the general equations of heat and mass transfer in two-dimensional region. Mathematical modeling of binary alloy solidification is still a challenging numericial problem. A Rigorous mathematical approach to this problem is available only when convection is not considered at all. The proposed numerical model was developed using the finite element code FIDAP. In the present study, the numerical model is used to consider thermal, solutal convection and a double diffusion source of mass transport.

Bune, Andris V.↗

DOE OSTI2024

Exploring Failure of Adhesively-Bonded Joints with Different Void Sizes at the Same Void Volume Fraction through A Micro-scale Numerical Modeling

This paper studied the effects of void sizes on the failure behavior of adhesively-bonded materials under global shear via micro-scale computational modeling. The numerical results indicate that large void sizes with increased distances between them at the same void volume fraction can alleviate the reduced joint strength of a weaker adherend-adhesive interface caused by interfacial voids. However, this is not the case for an adherend-adhesive interface with the interfacial tensile strength being sufficiently higher than the adhesive strength, showing the negligible effect on the joint strength due to various void sizes and the amount of interfacial voids. This preliminary investigation provides insightful information for scaling up the void size in the macro-scale computational modeling of joints, and enhances the understanding of the micro-mechanical adhesion in adherend-adhesive interfaces with different levels of bonding.

Qiao, Yao↗

DOE OSTI2022

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow $Q^{avg}_{Ao}$, maximum aortic root flow $Q^{max}_{Ao}$, average LVAD flow $Q^{avg}_{VAD}$, and maximum LVAD flow $Q^{max}_{VAD}$; (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2000

Numerical Modeling of Drying Residual RP-1 in Rocket Engines

When a Rocket Engine shuts down under a fuel rich environment, a significant amount of unburned RP-1 is trapped In the engine. It is necessary to clean the residual RP-1 prior to subsequent firing to avoid any explosion due to detonation. The conventional method is to dry RP-1 with inert gas such as Nitrogen or Helium. It is difficult to estimate the drying time unless the engine is adequately equipped with instruments to measure the trace of RP-1 during the drying process. Such instrumentation in flight hardware is often impractical and costly. On the other hand numerical modeling of the drying process can provide a good insight for a satisfactory operation of the process. A numerical model can provide answer to questions such as a) how long it takes to dry, b) which fluid is a better dryer for RP-1, c) how to reduce drying time etc. The purpose of the present paper is to describe a numerical model of drying RP-1 trapped in a cavity with flowing nitrogen or helium. The numerical model assumes one dimensional flow of drying fluid in contact with liquid pool of RP-1. An evaporative mass transfer takes place across the contact surface.

Majumdar, Alok↗

NASA NTRS1981

A mesoscale sixth-order numerical modelling system

A numerical simulation system is currently under development for NASA which is intended to improve the modeling of subsynoptic and mesoscale adjustments associated with cyclogenesis, severe storm development and atmospheric transport processes. The model utilizes a standard hydrostatic sigma-p coordinate primitive equation set, with x,y-space differencing accurate to eighth order. A three-step dynamic initialization procedure is employed between the analysis of real-time data and grid interpolation. Results of an 18-hour simulation during which synoptic scale cyclogenesis, subsynoptic scale jet streak adjustments, mesoscale convergence zones and tornadic storms were observed have shown the present model to have the potential for simulating the fine-scale structure of features associated with cyclogenesis and intense squall-line development. The mesoscale model was also found to produce less truncation than the NWS LFM model, although a frictionless version of the mesoscale model somewhat overdeepens and overaccelerates features.

Kaplan, M. L.↗

DOE OSTI2023

A hybrid analytical and numerical model for cross-over and performance decay in a unit cell vanadium redox flow battery

Developing an accurate and efficient model for cross-over is critical for improving the long-term performance of redox flow batteries (RFBs). Here, this work presents a hybrid analytical and numerical model that combines a two-dimensional analytical solution to the active species, a one-dimensional analytical model for cross-over mechanisms, and a zero-dimensional numerical model for outlet concentrations of reactants. By comparing with experiment over 41 cycles (ca. 144 h), the model reported a mean voltage difference of 0.0089 V, mean time difference of 14 s per cycle (of 3.5 h), maximum relative difference for capacity and energy of 1.34% and 1.63%, respectively. The predicted mean concentrations for V 2+ , V 3+ , and VO$_{2}^{+}$ in membrane are 65%~77% of measured values from the literature. Upon validation, the model reproduced behaviors in electrolyte imbalance similar to those observed in experiments and numerical models, and revealed the control of cross-over, self-discharge, stoichiometry of side reactions, and Coulombic efficiency on electrolyte imbalance. The model also demonstrates excellent computational efficiency for simulating 41 cycles (around 37,000 points) within 3~4 s. The demonstrated efficiency and accuracy for predicting cross-over and its impacts on voltage, capacity&energy decay, membrane concentrations, and electrolyte imbalance makes it a reliable tool for optimizing RFBs’ long-term performance.

25 ENERGY STORAGE↗

DOE OSTI2024

Wave-powered water pump for upwelling in aquaculture: Numerical model and ocean test

Wave-powered upwelling can increase the productivity and survivability of several aquaculture species. This enhancement is due to transporting cold, nutrient-rich ocean water, typically found lower in the water column, to the surface. Macroalgaes, like kelp, exhibit increased growth from these altered conditions. The University of New Hampshire’s (UNH) wave-powered water pump (wave pump) is a point absorber wave energy converter (WEC) that uses ocean waves to create relative motion between a spar buoy and a concentric float which drives an internal pump. A numerical model of the wave pump was developed using WEC-Sim to predict device performance in the ocean. Wave pump performance was evaluated during a five day ocean test near Appledore Island in Maine in March 2023, where volumetric flow rate, relative distance between spar and float, and wave conditions were measured. These data were then used for numerical model validation. The ocean deployment recorded the device’s performance in a variety of sea states, with average significant wave heights up to 0.7 m. The ocean test data were compared to the WEC-Sim numerical model of the device with favorable results. Average values of device stroke period, stroke height, and flow rate agreed between the ocean test and model data to within approximately 16 to 22%. Furthermore, the validated numerical model provides a valuable tool for improving the design and developing a commercial-scale, wave-powered water pump for use in aquaculture.

16 TIDAL AND WAVE POWER↗

DOE OSTI2021

Validation of numerical models for seismic fluid-structure-interaction analysis of nuclear, safety-related equipment

Seismic design, qualification, and risk assessment of nuclear safety-related vessels filled with liquid will have to consider the interaction between the vessel (tank), the contained liquid, and submerged components, if any. Seismic fluid-structure-interaction (FSI) analysis of nuclear vessels will rely on numerical models, which are required to be verified and validated. This study validates previously-verified numerical models using test data generated from earthquake-simulator experiments involving a cylindrical tank. Two solvers in LS-DYNA for FSI simulations are used for the numerical analysis: Arbitrary-Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD). Numerical and test results are compared for responses critical to seismic design of advanced reactor vessels: hydrodynamic pressures on the tank wall, reactions at the support, and wave heights of the contained liquid. Analysis is performed for one-, two-, and three-directional seismic inputs with a range of intensity, and rocking motions. The accuracy of the numerical results and different methods for outputting wave heights are discussed. Recommendations for validation of seismic FSI numerical models of advanced reactors are provided. As a result, the validation exercise presented in this study is broadly applicable to cylindrical tanks, regardless of industry sector.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

NASA NTRS1990

Numerical modeling of an arcjet thruster

A numerical model of an arcjet thruster was developed and the results of calculations using this model are compared to the experimental data from a 30 kW arcjet using ammonia as a propellant. The model contains equations for the conservation of mass, radial, axial, and azimuthal momentum, energy, and the radial component of the magnetic field and provides a solution for three components of velocity, temperature, and the axial and radial current as a function of position in an axisymmetric flowfield. The model predicts 12 percent more specific impulse than the experiment when the mass flow and power input are constrained to the experimental values. A major conclusion of this study is that improved procedures for calculating transport properties are necessary if the accuracy of the model is to be improved.

Rhodes, Robert P.↗

NASA NTRS1995

Dynamics of the Antarctic Circumpolar Current. Evidence for Topographic Effects from Altimeter Data and Numerical Model Output

Geosat altimeter data and numerical model output are used to examine the circulation and dynamics of the Antarctic Circumpolar Current (ACC). The mean sea surface height across the ACC has been reconstructed from height variability measured by the altimeter, without assuming prior knowledge of the geoid. The results indicate locations for the Subantarctic and Polar Fronts which are consistent with in situ observations and indicate that the fronts are substantially steered by bathymetry. Detailed examination of spatial and temporal variability indicates a spatial decorrelation scale of 85 km and a temporal e-folding scale of 34 days. Empirical Orthogonal Function analysis suggests that the scales of motion are relatively short, occuring on 1000 km length-scales rather than basin or global scales. The momentum balance of the ACC has been investigated using output from the high resolution primitive equation model in combination with altimeter data. In the Semtner-Chervin quarter-degree general circulation model topographic form stress is the dominant process balancing the surface wind forcing. In stream coordinates, the dominant effect transporting momentum across the ACC is bibarmonic friction. Potential vorticity is considered on Montgomery streamlines in the model output and along surface streamlines in model and altimeter data. (AN)

Gille, Sarah T.↗

NASA NTRS2007

Numerical Modeling of Propellant Boil-Off in a Cryogenic Storage Tank

A numerical model to predict boil-off of stored propellant in large spherical cryogenic tanks has been developed. Accurate prediction of tank boil-off rates for different thermal insulation systems was the goal of this collaboration effort. The Generalized Fluid System Simulation Program, integrating flow analysis and conjugate heat transfer for solving complex fluid system problems, was used to create the model. Calculation of tank boil-off rate requires simultaneous simulation of heat transfer processes among liquid propellant, vapor ullage space, and tank structure. The reference tank for the boil-off model was the 850,000 gallon liquid hydrogen tank at Launch Complex 39B (LC- 39B) at Kennedy Space Center, which is under study for future infrastructure improvements to support the Constellation program. The methodology employed in the numerical model was validated using a sub-scale model and tank. Experimental test data from a 1/15th scale version of the LC-39B tank using both liquid hydrogen and liquid nitrogen were used to anchor the analytical predictions of the sub-scale model. Favorable correlations between sub-scale model and experimental test data have provided confidence in full-scale tank boil-off predictions. These methods are now being used in the preliminary design for other cases including future launch vehicles

Majumdar, A. K.↗

DOE OSTI2023

n-Propylbenzene Droplet Combustion: Experiments and Numerical Modeling

n-Propylbenzene (PB, C9H12, 432K boiling point) is important as a constituent in surrogates for some transportation fuels but little is understood about its droplet burning characteristics. We present an experimental and numerical study of a fuel droplet burning under conditions where the sole means of gas transport is fuel evaporation at the droplet surface thereby creating a one-dimensional gas transport that is amenable to detailed numerical modeling. The one-dimensional gas transport that results is amenable to detailed numerical modeling. Initial droplet diameters (Do) ranged from 0.3 mm to about 5 mm to better understand the mechanisms that control burning over such a wide range of Do. For Do < 1 mm droplets were burned while anchored to 14 m SiC fibers within a sealed container in free-fall during which time the droplets were ignited by spark discharge and the burning process recorded by video cameras. Larger and free-floating droplets were studied in the spaced-based platform of the orbiting International Space Station for Do > 1.0 mm which provided unlimited times for experimental observation. Measurements of droplet, flame and soot shell diameters were made and are discussed in the presentation. The experimental results showed formation of optically thick soot clouds developing during burning that required a robust image analysis algorithm to enable measurement of droplet diameter. The measured droplet diameters were well predicted by a numerical model that included soot formation, a PB kinetic mechanism comprising 576 species and 27369 reactions, unsteady gas and liquid transport, variable properties and radiative transport. Flame diameters showed more variability compared to the simulations, depending on how the flame was defined (peak gas temperature and OH concentration). Droplets burned slower as Do was increased and the PB sooting propensity increased with Do. Flame extinction was observed for Do > 3.3 mm indicating a strong effect of radiation as Do was increased compared to the smaller droplet diameters examined where extinction was not experimentally observed. For the largest droplets examined the possibility that the observed extinction is indicative of a transition to cool flame burning is still under investigation.

Guo, Songtao↗

NASA NTRS1976

Interactive applications of satellite observations and mesoscale numerical models

The impact of numerical weather prediction (NWP) and satellite meteorology on operational weather forecasting has become overwhelming in the past few years. The paper looks toward the merger of these technologies in making short range 6-18 h forecasts through the use of mesoscale NWP models. A short-range (2-18 h) mesoscale forecast system envisioned for the near future is described that includes four components: hydrodynamic numerical models, large-scale and mesoscale; satellites, polar orbiting for high latitudes and geostationary for low latitudes; mesoclimatology, derived in large part from satellite data; and special-purpose simple models and empirical relations. It is important that the components of the forecast system be developed in parallel rather than in series if the system is to be completed within five years. There is enough evidence to substantiate the revolution in the mesoscale weather prediction in the next five years.

Kreitzberg, C. W.↗

DOE OSTI2021

A process to verify numerical models for seismic fluid-structure interaction in advanced reactor vessels

Seismic design and qualification of a liquid-filled advanced nuclear reactor will have to account for fluid-structure interaction (FSI). Interaction between the tank, internal components, and contained liquid will rely on analysis of numerical models that must be verified and validated. Here this study demonstrates a verification process for models of a base-supported cylindrical tank by comparing numerical predictions and analytical solutions. The numerical models are consistent with the assumptions made to derive analytical solutions, namely, either a rigid or a linear elastic tank, ideal fluid, and small-amplitude, unidirectional, horizontal inputs. One software platform is used to illustrate the process. Seismic FSI analysis is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Reported responses are those used for design, including hydrodynamic pressures on the tank wall, shear forces and moments at the tank base, and wave heights of the contained liquid. The accuracy of the numerical results is discussed. The numerical models are verified for calculating the pressures on the tank wall and reactions at its base. Accurate simulation of wave action is challenging for both solvers. Recommendations for modeling, code development, and steps for verification are provided. Although focused on reactor vessels and one software platform, the verification process described herein is broadly applicable to liquid-filled vessels and other finite element codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗