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 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Utilization of satellite data and regional scale numerical models in short range weather forecasting

Overwhelming evidence was developed in a number of studies of satellite data impact on numerical weather prediction that it is unrealistic to expect satellite temperature soundings to improve detailed regional numerical weather prediction. It is likely that satellite data over the United States would substantially impact mesoscale dynamical predictions if the effort were made to develop a composite moisture analysis system. The horizontal variability of moisture, most clearly depicited in images from satellite water vapor channels, would not be determined from conventional rawinsondes even if that network were increased by a doubling of both the number of sites and the time frequency.

Kreitzberg, C. W.↗

NASA NTRS1990

Three-dimensional numerical modelling of transport in horizontal MOCVD reactors

A systematic three-dimensional study for the MOCVD of gallium arsenide is presented. Numerical solutions, obtained for the system of equations governing the steady state heat, momentum, and mass transfer in horizontal reactors with large and small aspect ratio, are compared with experimental growth rate distributions. The effects on layer thickness uniformity of lateral thermal boundary conditions and susceptor tilt are explored. The validity of a parabolic flow approximation of the equations is also investigated by comparison with the complete three-dimensional treatment.

Ouazzani, Jalil↗

DOE OSTI2023

Comparison of Designs of Hydrogen Isotope Separation Columns by Numerical Modeling

Mixtures of gas-phase hydrogen isotopologues (diatomic combinations of protium, deuterium, and tritium) can be separated using columns containing a solid such as palladium that reversibly absorbs hydrogen. A temperature-swing process can transport hydrogen into or out of a column by inducing temperature-dependent absorption or desorption reactions. Here, we consider two designs: a thermal cycling absorption process, which moves hydrogen back and forth between two columns, and a simulated moving bed (SMB), where columns are in a circular arrangement. We present a numerical mass and heat transport model of absorption columns for hydrogen isotope separation. It includes a detailed treatment of the absorption–desorption reaction for palladium. By comparing the isotope concentrations within the columns as a function of position and time, we observe that SMB can lead to sharper separations for a given number of thermal cycles by avoiding the remixing of isotopes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

An immersed interface method for microstructure-scale electrochemical battery models: numerical formulation and performance portable implementation

We present the numerical formulation, verification, and performance portable implementation of an immersed interface method for microstructure scale electrochemical modeling of batteries. The innovation in this approach is the resolution of chemical species and electrostatic potential discontinuities at active interfaces without the use of interface conforming unstructured grids. A unified formulation on Cartesian grids for all domains (electrodes and electrolyte) is used with interfacial flux conditions applied using volume fraction or “color” function gradients. We have developed one dimensional and two dimensional test cases with analytic solutions for electrochemical modeling using which we verified the consistency and accuracy of our scheme. Our solver is also validated against solutions from a macroscale model and an unstructured multi-subdomain solver for a full lithium ion cell. We then demonstrated the utility of our solver on an image-based complex battery electrode microstructure at high charging rate. Our technique also exhibits good scalability on distributed memory architectures using central processing units (CPU), with problem sizes up to 1.8 billion degrees of freedom and with 5400 ranks. Initial performance studies of our open-source performance portable solver showed about 70 times speed up using a graphics processing unit (GPU) compared to single compute core for a problem with 4 million cells.

25 ENERGY STORAGE↗

NASA NTRS1988

Numerical Modeling Of Two-Phase, Reactive Flows

Capabilities of several computational schemes combined and enhanced. Advanced Rocket Injector-Combustor Code is computer program for numerical simulation of two-phase flows of chemically reacting or nonreacting components. Applied to sprays of fuel in injectors and combustors, flows in heat exchangers, evaporators, and mass-diffusion systems; and two-dimensional or axisymmetric flows containing liquid or solid suspensions. Used in research for development and verification of various physical models.

Liang, Pak-Yan↗

NASA NTRS1983

Numerical modeling of turbulent flow in a channel

Two-dimensional incompressible turbulent flow in a channel with a backward-facing step was studied numerically by Chorin's Random Vortex Method (RVM), an algorithm capable of tracing the action of elementary turbulent eddies and their cumulative effects without imposing any restrictions upon their motions. The step occurs in one side of a channel with otherwise flat, parallel walls; its height equals 1/3, 1/4 or 1/5 the width of the channel downstream. The main objective was to investigate the behavior of the large-scale turbulent eddies in a flow and the flow characteristics in the separated shear layer, the reattached zone, and the rebuilding boundary layer after reattachment. The unsteady vorticity field and the distribution of time-averaged turbulent statistics were obtained. The effects of expansion step height and initial boundary layer state were also studied. Comparisons were made with the available experimental results. The agreement is satisfactory in the velocity profiles and in the reattachment length, and fairly good in the turbulence profiles. Also a mechanism of the development of the reattaching turbulent flow was suggested by the numerical results.

Dai, Y. W.↗

NASA NTRS2006

Numerical Modeling of Turbulence Effects within an Evaporating Droplet in Atomizing Sprays

A new approach to account for finite thermal conductivity and turbulence effects within atomizing liquid sprays is presented in this paper. The model is an extension of the T-blob and T-TAB atomization/spray model of Trinh and Chen (2005). This finite conductivity model is based on the two-temperature film theory, where the turbulence characteristics of the droplet are used to estimate the effective thermal diffhsivity within the droplet phase. Both one-way and two-way coupled calculations were performed to investigate the performance of this model. The current evaporation model is incorporated into the T-blob atomization model of Trinh and Chen (2005) and implemented in an existing CFD Eulerian-Lagrangian two-way coupling numerical scheme. Validation studies were carried out by comparing with available evaporating atomization spray experimental data in terms of jet penetration, temperature field, and droplet SMD distribution within the spray. Validation results indicate the superiority of the finite-conductivity model in low speed parallel flow evaporating spray.

Balasubramanyam, M. S.↗

NASA NTRS1992

Numerical Modeling of a Vortex Stabilized Arcjet

Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As the level of swirl and viscosity in the flowfield increased the mass flow rate and thrust decreased. The technique was used to predict the flow through a typical arcjet thruster geometry. Results indicate swirl and viscosity play an important role in the complex geometry of an arcjet.

Pawlas, Gary E.↗

NASA NTRS1994

Erosional dynamics, flexural isostasy, and long-lived escarpments: A numerical modeling study

Erosional escarpments common features of high-elevation rifted continets. Fission track data suffest that these escarpments form by base level lowering and/or marginal uplift during rifting, followed by lateral retreat of an erosion front across tens to hundreds of kioometers. Previous modeling studies have shown that this characteristic pattern of denudation can have a profound impact upon marginal isostatic uplift and the evolution of offshore sedimentary basins. Yet at present there is only a rudimentary understanding of the geomorphic mechanisms capable of driving such prolonged escarpment retreat. In this study we present a nonlinear, two-dimensional landscape evolution model tha tis used to asses the necessary and sufficient conditions for long-term retreat of a rift-generated escarpment. The model represents topography as a grid of cells, with drainage networkds evolving as water flows across the grid in the direction of steepest descent. The model accounts for sediment production by weathering, fluvial sediment transport, bedrock channel erosion, and hillslope sediment transport by diffusive mechanisms and by mass failure. Numerical experiments presented explore the effects of different combinations of erosion processes and of dynamic coupling between denudation and flexural isostatic uplift. Model results suggest that the necessary and sufficient conditions for long-term escarpment retreat are (1) incising bedrock channels in which the erosion rate increases with increasing drainage area, so that the channels steepen and propagate headward; (2) a low rate of sediment production relative to sediment transport efficiency, which promotes relief-generating processes over diffusive ones; (3) high continental elevation, which allows greater freedom for fluvial dissection; and (4) any process, including flexural isostatic uplift, that helps to maintain a drainage divide near an escarpment crest. Flexural isostatic uplift also facilitates escarpment, thereby increasing channel gradients and accelerating erosion which in turn generates additional isostatic uplift. Of all the above conditions, high continental elevation is common ot most rift margin escarpments and may ultimately be the most important factor.

Tucker, Gregory E.↗

NASA NTRS2001

Structural and Acoustic Numerical Modeling of a Curved Composite Honeycomb Panel

The finite and boundary element modeling of the curved section of a composite honeycomb aircraft fuselage sidewall was validated for both structural response and acoustic radiation. The curved panel was modeled in the pre-processor MSC/PATRAN. Geometry models of the curved panel were constructed based on the physical dimensions of the test article. Material properties were obtained from the panel manufacturer. Finite element models were developed to predict the modal parameters for free and supported panel boundary conditions up to a frequency of 600 Hz. Free boundary conditions were simulated by providing soft foam support under the four comers of the panel or by suspending the panel from elastic bands. Supported boundary conditions were obtained by clamping the panel between plastic tubing seated in grooves along the perimeter of a stiff and heavy frame. The frame was installed in the transmission loss window of the Structural Acoustic Loads and Transmission (SALT) facility at NASA Langley Research Center. The structural response of the curved panel due to point force excitation was predicted using MSC/NASTRAN and the radiated sound was computed with COMET/Acoustics. The predictions were compared with the results from experimental modal surveys and forced response tests on the fuselage panel. The finite element models were refined and updated to provide optimum comparison with the measured modal data. Excellent agreement was obtained between the numerical and experimental modal data for the free as well as for the supported boundary conditions. Frequency response functions (FRF) were computed relating the input force excitation at one panel location to the surface acceleration response at five panel locations. Frequency response functions were measured at the same locations on the test specimen and were compared with the calculated FRF values. Good agreement was obtained for the real and imaginary parts of the transfer functions when modal participation was allowed up to 3000 Hz. The validated finite element model was used to predict the surface velocities due to the point force excitation. Good agreement was obtained between the spatial characteristics of the predicted and measured surface velocities. The measured velocity data were input into the acoustic boundary element code to compute the sound radiated by the panel. The predicted sound pressure levels in the far-field of the panel agreed well with the sound pressure levels measured at the same location.

Grosveld, Ferdinand W.↗

NASA NTRS1997

Science Support for Space-Based Droplet Combustion: Drop Tower Experiments and Detailed Numerical Modeling

This program supports the engineering design, data analysis, and data interpretation requirements for the study of initially single component, spherically symmetric, isolated droplet combustion studies. Experimental emphasis is on the study of simple alcohols (methanol, ethanol) and alkanes (n-heptane, n-decane) as fuels with time dependent measurements of drop size, flame-stand-off, liquid-phase composition, and finally, extinction. Experiments have included bench-scale studies at Princeton, studies in the 2.2 and 5.18 drop towers at NASA-LeRC, and both the Fiber Supported Droplet Combustion (FSDC-1, FSDC-2) and the free Droplet Combustion Experiment (DCE) studies aboard the shuttle. Test matrix and data interpretation are performed through spherically-symmetric, time-dependent numerical computations which embody detailed sub-models for physical and chemical processes. The computed burning rate, flame stand-off, and extinction diameter are compared with the respective measurements for each individual experiment. In particular, the data from FSDC-1 and subsequent space-based experiments provide the opportunity to compare all three types of data simultaneously with the computed parameters. Recent numerical efforts are extending the computational tools to consider time dependent, axisymmetric 2-dimensional reactive flow situations.

Marchese, Anthony J.↗

NASA NTRS1981

Numerical modeling of diffusive physical vapor transport in cylindrical ampoules

Diffusion limited physical vapor transport in closed cylindrical ampoules of aspect ratio (length/radius) between 0.5 and 10 was studied numerically. The transport of a crystal forming component A through an inert component B, which undergoes zero net transport, was considered. A case in which the partial pressure of B is comparable to that of A, as well as systems with very low partial pressures of B (impurity cases) were treated. Concentrations and velocity distributions were obtained from the coupled transport equations for momentum, mass and components. It was found that the interdiffusion of components in viscous interaction with the container walls leads to (1) recirculation of component B, even in the absence of gravity; (2) concentration gradients normal to the main transport direction in the vapor space; and (3) non-uniform interfacial concentration gradients and, hence, non-uniform crystal growth rates.

Greenwell, D. W.↗

NASA NTRS1988

Numerical modeling of multidimensional flow in seals and bearings used in rotating machinery

The rotordynamic behavior of turbomachinery is critically dependent on fluid dynamic rotor forces developed by various types of seals and bearings. The occurrence of self-excited vibrations often depends on the rotor speed and load. Misalignment and rotor wobbling motion associated with differential clearance were often attributed to stability problems. In general, the rotative character of the flowfield is a complex three dimensional system with secondary flow patterns that significantly alter the average fluid circumferential velocity. A multidimensional, nonorthogonal, body-fitted-grid fluid flow model is presented that describes the fluid dynamic forces and the secondary flow pattern development in seals and bearings. Several numerical experiments were carried out to demonstrate the characteristics of this complex flowfield. Analyses were performed by solving a conservation form of the three dimensional Navier-Stokes equations transformed to those for a rotating observer and using the general-purpose computer code PHOENICS with the assumptions that the rotor orbit is circular and that static eccentricity is zero. These assumptions have enabled a precise steady-state analysis to be used. Fluid injection from ports near the seal or bearing center increased fluid-film direct dynamic stiffness and, in some cases, significantly increased quadrature dynamic stiffness. Injection angle and velocity could be used for active rotordynamic control; for example, injection, when compared with no injection, increased direct dynamic stiffness, which is an important factor for hydrostatic bearings.

Hendricks, R. C.↗

NASA NTRS1996

Numerical modelling of the circulation of superrotating atmospheres: Venus and Titan

Atmospheric superrotation is known to exist on Venus where the cloud cover, at 70 km above the surface, rotates about 60 times faster than the solid planet. Superrotation is strongly suspected on Titan on the basis of numerical studies and observations. The numerical simulation of superrotation is described. It is shown how the general circulation models, originally developed for terrestrial meteorology and climatology, can be used to infer the circulation in such atmospheres. The emphasis is on the case of Titan for which a specific general circulation model was developed.

Hourding, F.↗

NASA NTRS1992

Numerical models for the collapse and fragmentation of centrally condensed molecular cloud cores

The gravitational collapse and fragmentation of centrally condensed molecular cloud cores are investigated using a new hydrodynamical code. The numerical scheme is second-order accurate and uses explicit finite difference methods to advance the fluid variables on a 3D Cartesian grid. Two initial power-law density profiles, rho varies as r exp -1 and rho varies as r exp -2, are considered, as well as two initial density perturbations in the azimuthal coordinate theta, rho-i - rho(1 + a cos 2theta) where a = 0.1 and 0.5. Fragmentation is found to be possible in these centrally condensed cores if the initial conditions also include differential rotation. Models which collapse with initial uniform rotation do not produce fragments. If molecular cloud cores are indeed centrally condensed, as suggested by observations of star-forming regions and by studies of ambipolar diffusion, then differential rotation may be a mechanism for producing binary protostars during gravitational collapse.

Myhill, Elizabeth A.↗

DOE OSTI2025

A Numerical Modeling Framework for Flocculation and Cohesive Sediment Transport in the Wave Bottom Boundary Layer

Flocculation, a critical process in coastal and estuarine systems, plays a significant role in sediment transport, nutrient cycling, and ecological health. This study develops a cohesive sediment transport modeling framework tailored to the wave bottom boundary layer under dilute and equilibrium conditions, explicitly incorporating flocculation effects via a Population Balance Equation (PBE). Using Direct Numerical Simulation, six baseline cases, each with a distinct sediment concentration profile resulting from a constant settling velocity and critical erosion shear stress, are generated to drive the PBE flocculation model for given floc yield strength and stickiness. Results reveal that flocculation significantly influences sediment concentration profiles promoting three distinct stages, well‐mixed, transition to lutocline, and well‐developed lutocline. At low concentrations with well‐mixed profiles, cohesive floc properties are less significant, and turbulence is a main flocculation driver. In contrast, as concentration increases, cohesive floc properties become crucial, facilitating lutocline formation. Here, the analysis also highlights limitations of depth‐averaged settling velocity as a parameterization. It is suitable for well‐mixed and transitional profiles but fails in well‐developed lutoclines, where empirical formulations that explicitly incorporate turbulent shear rate and sediment concentration better capture variability. This study underscores the necessity of incorporating flocculation effects into sediment transport models to enhance predictions of sediment dynamics in wave bottom boundary layers.

Penaloza‐Giraldo, Jorge A. [Oak Ridge National Lab↗

NASA NTRS1990

Two-dimensional atmospheric transport and chemistry model - Numerical experiments with a new advection algorithm

Extensive testing of the advective scheme proposed by Prather (1986) has been carried out in support of the California Institute of Technology-Jet Propulsion Laboratory two-dimensional model of the middle atmosphere. The original scheme is generalized to include higher-order moments. In addition, it is shown how well the scheme works in the presence of chemistry as well as eddy diffusion. Six types of numerical experiments including simple clock motion and pure advection in two dimensions have been investigated in detail. By comparison with analytic solutions, it is shown that the new algorithm can faithfully preserve concentration profiles, has essentially no numerical diffusion, and is superior to a typical fourth-order finite difference scheme.

Shia, Run-Lie↗

NASA NTRS2001

Numerical Modeling of Liquid-Vapor Phase Change

We implemented a two- and three-dimensional finite difference/front tracking technique to solve liquid-vapor phase change problems. The mathematical and the numerical features of the method were explained in great detail in our previous reports, Briefly, we used a single formula representation which incorporated jump conditions into the governing equations. The interfacial terms were distributed as singular terms using delta functions so that the governing equations would be the same as conventional conservation equations away from the interface and in the vicinity of the interface they would provide correct jump conditions. We used a fixed staggered grid to discretize these equations and an unstructured grid to explicitly track the front. While in two dimensions the front was simply a connection of small line segments, in three dimensions it was represented by a connection of small triangular elements. The equations were written in conservative forms and during the course of computations we used regriding to control the size of the elements of the unstructured grid. Moreover, we implemented a coalescence in two dimensions which allowed the merging of different fronts or two segments of the same front when they were sufficiently close. We used our code to study thermocapillary migration of bubbles, burst of bubbles at a free surface, buoyancy-driven interactions of bubbles, evaporation of drops, rapid evaporation of an interface, planar solidification of an undercooled melt, dendritic solidification, and a host of other problems cited in the reference.

Esmaeeli, Asghar↗