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,063 records · Page 59

TPSAS-NF1676L-19616-DND

The present study has three goals: (1) perform an experiment where a simple laminate damage process can be characterized in high detail; (2) evaluate the performance of existing commercially available laminate damage simulation tools by modeling the experiment; (3) observe and understand the underlying physics of damage in a composite honeycomb sandwich structure subjected to low-velocity impact. A quasi-static indentation experiment has been devised to provide detailed information about a simple mixed-mode damage growth process. The test specimens consist of an aluminum honeycomb core with a cross-ply laminate facesheet supported on a stiff uniform surface. When the sample is subjected to an indentation load, the honeycomb core provides support to the facesheet resulting in a gradual and stable damage growth process in the skin. This enables real time observation as a matrix crack forms, propagates through a ply, and then causes a delamination. Finite element analyses were conducted in ABAQUS/Explicit(TradeMark) 6.13 that used continuum and cohesive modeling techniques to simulate facesheet damage and a geometric and material nonlinear model to simulate core crushing. The high fidelity of the experimental data allows a detailed investigation and discussion of the accuracy of each numerical modeling approach.

Mark McElroy↗

A model to determine open or closed cellular convection

The horizontal asymmetry of cellular convection is demonstrated to be detectable by consideration of the vertical asymmetry of the driving force. A two-dimensional numerical model for convection in an internally heated and cooled fluid is presented, based on equations for the conservation of temperature and vorticity. Attention is focused on the steady-state finite-amplitude solutions for fixed Rayleigh number, Prandtl number, and aspect ratio. Temperature in the model corresponds to potential temperature under dry conditions and the equivalent potential temperature under saturated conditions. Slowly varying convection driven by asymmetric boundary fluxes is expressed in terms of steady-state convection driven by an asymmetric internal heat source. It is shown that heating near the ground produces open convection patterns where most of the fluid is descending, while cooling near the top of the flow leads to closed cellular patterns with a preponderance of ascending fluid. Extension of the model to three dimensions is indicated.

Helfand, H. M.↗

Thermal emission from the dust coma of Comet Bowell and a model for the grains

An analysis of the physical characteristics of Comet Bowell is presented based on 10 and 20 microns emissions observed with the NASA 3 m IR telescope during 1981-82. Numerical models are defined for quantifying the number of dust particles at a given radius and at a given angular distance from the nucleus, the contribution of icy grains to the 20 microns flux and the average albedo, and the dust distribution in the coma. Large, dirty ice grains are judged to be responsible for production of OH at a rate of about 10 to the 29th mol/sec during the spring of 1981. The total surface area of the dirty ice grains is projected to be 10 to the 15th sq cm. The coma in 1982 would therefore have consisted of nonvolatile grains, i.e., old dirty ice grains. Further observations of other, similar comets are needed to identify the composition of the emitted grains and the heliocentric distance at which high water production is initiated.

Hanner, M. S.↗

Modeling flow at the nozzle of a solid rocket motor

The mechanical behavior of a rocket motor internal flow field results in a system of nonlinear partial differential equations which can be solved numerically. The accuracy and the convergence of the solution of the system of equations depends largely on how precisely the sharp gradients can be resolved. An adaptive grid generation scheme is incorporated into the computer algorithm to enhance the capability of numerical modeling. With this scheme, the grid is refined as the solution evolves. This scheme significantly improves the methodology of solving flow problems in rocket nozzle by putting the refinement part of grid generation into the computer algorithm.

Chow, Alan S.↗

Modeling of shallow and inefficient convection in the outer layers of the Sun using realistic physics

In an attempt to understand the properties of convective energy transport in the solar convective zone, a numerical model has been constructed for turbulent flows in a compressible, radiation-coupled, nonmagnetic, gravitationally stratified medium using a realistic equation of state and realistic opacities. The time-dependent, three-dimensional hydrodynamic equations are solved with minimal simplifications. The statistical information obtained from the present simulation provides an improved undserstanding of solar photospheric convection. The characteristics of solar convection in shallow regions is parameterized and compared with the results of Chan & Sofia's (1989) simulations of deep and efficient convection. We assess the importance of the zones of partial ionization in the simulation and confirm that the radiative energy transfer is negliglble throughout the region except in the uppermost scale heights of the convection zone, a region of very high superadiabaticity. When the effects of partial ionization are included, the dynamics of flows are altered significantly. However, we confirm the Chan & Sofia result that kinetic energy flux is nonnegligible and can have a negative value in the convection zone.

Kim, Yong-Cheol↗

Analytic Couple Modeling Introducing Device Design Factor, Fin Factor, Thermal Diffusivity Factor, and Inductance Factor

A set of convenient thermoelectric device solutions have been derived in order to capture a number of factors which are previously only resolved with numerical techniques. The concise conversion efficiency equations derived from governing equations provide intuitive and straight-forward design guidelines. These guidelines allow for better device design without requiring detailed numerical modeling. The analytical modeling accounts for factors such as i) variable temperature boundary conditions, ii) lateral heat transfer, iii) temperature variable material properties, and iv) transient operation. New dimensionless parameters, similar to the figure of merit, are introduced including the device design factor, fin factor, thermal diffusivity factor, and inductance factor. These new device factors allow for the straight-forward description of phenomenon generally only captured with numerical work otherwise. As an example a device design factor of 0.38, which accounts for thermal resistance of the hot and cold shoes, can be used to calculate a conversion efficiency of 2.28 while the ideal conversion efficiency based on figure of merit alone would be 6.15. Likewise an ideal couple with efficiency of 6.15 will be reduced to 5.33 when lateral heat is accounted for with a fin factor of 1.0.

Seebeck effect↗

Axisymmetric Two-Dimensional Modeling of No Vent Filling of a Cryogenic Tank using Generalized Fluid System Simulation Program

This paper presents preliminary analysis using a finite volume based computational method for network flow analysis for the simulation of chill and fill of a cryogenic tank. An axisymmetric two-dimensional model of a spherical tank was developed using a nodal approach. No Vent Filling (NVF) was accomplished using a Thermodynamic Vent System assisted injector where the injector was cooled by the cold gas generated by Joule-Thompson expansion of the liquid propellant. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around spray droplets, and c) condensation of vapor at the interface of the cooled injector and ullage. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data and one-dimensional (1D) model predictions. It was demonstrated that the axisymmetric 2D model provided a more realistic prediction than the 1D model.

Cryogenics↗

Experiment 5: Science and Technology of Surface Controlled Oscillations: Report on USML-2 Results

Minuscule amounts (e.g., 1 part in 10,000) of a surface-active material in a liquid can drastically affect the surface behavior of the liquid, influencing how the material flows and mixes with other liquid and solid materials. In many respects, the science of surfactants has been empirical, with trial and error dominating over the ability to predict how surfactant type and concentration influence surface behavior. A program for the modeling of surfactant behavior has been established at Yale. This program combines experimental work performed both on the ground and in space, and theoretical and numerical modeling. By levitating a drop of liquid in air, away from solid container surfaces, and by manipulating the drop with acoustic radiation forces, we have been able to establish idealized conditions for surface behavior studies. The primary experiments involve the study of the free oscillations of initially deformed drops. In STS-73, the USML-2 mission of the Space Shuttle, we performed the following measurements: 1) the oscillation of a spherical drop in its quadrupole mode; 2) the oscillation of a drop about a deformed (oblate) shape; 3) the slow static squeezing of the drop from spherical to nearly flat; and 4) the superoscillations of drops when the radiation forces maintaining the drop in a flattened state are suddenly reduced. Analytic and numerical studies have enabled us to understand the physics of these oscillations and to extract material properties such as the dynamic surface tension and the surface viscosities (shear and dilatational). The relation to ground-based studies is essential, because the knowledge and understanding gleaned from our space studies enable us to interpret ground-based data.

Apfel, Robert E.↗

An evaluation of soundings, analyses and model forecasts derived from TIROS-N and NOAA-6 satellite data

TIROS-N and NOAA-6 temperature soundings over North America during three days in January 1980, and synoptic analyses and numerical-model forecasts derived from them, are compared with conventional data and analyses from NMC's limited-area fine-mesh model (LFM). The collocated sounding comparison revealed significant errors, especially near the surface and the tropopause. Satellite-derived thermal gradients were found to be weak, and thickness-analysis difference fields to propagate eastward, suggesting that sounding errors are correlated with synoptic patterns. The same pattern of anomalies is seen in the model forecasts. More detailed determinations of the correlation detected here could be used to optimize the assimilation of satellite soundings to conventional data.

Koehler, T. L.↗

New meteoroid model predictions for directional impacts on LDEF

An extensive body of data, from meteors, zodiacal light, spacecraft-borne impact detectors (Helios, Pioneer, Galileo, Ulysses), and other sources, forms the basis of a new numerical model for the distributions of interplanetary meteoroids. For each of the five populations in this model it is possible to evaluate meteoroid concentration and flux for oriented surfaces or detectors having arbitrary position and velocity in interplanetary space. For a spacecraft in geocentric orbit the effects of gravitational focussing and shielding by the Earth have been newly derived with full attention to the directionality of the particles, both on approach (i.e., relative to a massless Earth) and at the target. This modeling approach was exercised to provide an estimate of meteoroid fluence for each of several oriented surfaces on LDEF.

Divine, Neil↗

A one-dimensional modeling study of carbonaceous haze effects on the springtime Arctic environment

The effect of carbonaceous Arctic haze on the surface energy balance was investigated using a specially developed one-dimensional three-component numerical model which accounts for the transfer of IR and solar radiation through the atmosphere, turbulent mixing within the boundary layer, and heat conduction through snow-capped sea ice. The results of calculations indicate that, when haze is present, the surface temperature increases, provided the relative humidity through the column does not change. The increase is due to an increase in the absorbed radiation at the surface.

Emery, Christopher A.↗

Experiments And Model Development For The Investigation Of Sooting And Radiation Effects In Microgravity Droplet Combustion

This study involves flight experiments (for droplets between 1.5 to 5 mm) and supportive ground-based experiments, with concurrent numerical model development and validation. The experiments involve two fuels: n-heptane, and ethanol. The diagnostic measurements include light extinction for soot volume fraction, two-wavelength pyrometry and thin-filament pyrometry for temperature, spectral detection for OH chemiluminescence, broadband radiometry for flame emission, and thermophoretic sampling with subsequent transmission electron microscopy for soot aerosol property calculations.

Yozgatligil, Ahmet↗

Bubble Formation on a Wall in Cross-Flowing Liquid and Surrounding Fluid Motion,With and Without Heating

Application of gas-liquid two-phase flow systems for space-based thermal management and for the HEDS program demands a precise control of bubble size distribution in liquid. The necessity of bulk liquid motion for controlling bubble size and frequency in the space environment has been suggested by recent studies on pool, forced convection boiling and bubble formation in flowing liquid. The present work, consisting of two parts, explores bubble generation at wall in a cross-flowing liquid, i.e., in a forced convection boiling configuration. A schematic is shown. The first part looks into the bubble formation process under isothermal conditions in a reduced gravity environment, by injecting gas through a hole in the wall of a flowing liquid channel. In the latter part with channel wall heating, flow and temperature fields near a single bubble are studied under normal (1-g) and micro-gravity (mu-g) conditions. The objective of the isothermal experiments is to experimentally investigate the effects of liquid cross-flow velocity, gas flow rate, and orifice diameter on bubble formation. Data were taken mainly under reduced gravity conditions but some data were taken in normal gravity for comparison. The reduced gravity experiment was conducted aboard the NASA DC-9 Reduced Gravity Aircraft. The results show that the process of bubble formation and detachment depends on gravity, the orifice diameter (D(sub N)), the gas flow rate (Q(sub g)), and the liquid cross-flow velocity (U(sub L)). The reduced gravity data are shown. The data are analyzed based on a force balance, and two different detachment mechanisms are identified. When the gas momentum is large, the bubble detaches from the injection orifice as the gas momentum overcomes the attaching effects of liquid drag and inertia. The surface tension force is much reduced because a large part of the bubble pinning edge at the orifice is lost as the bubble axis is tilted by the liquid flow. When the gas momentum is small, the force balance in the liquid flow direction is important, and the bubble detaches when the bubble axis inclination exceeds a certain angle. With wall heating, liquid motion around an air bubble in cross-flowing 2cs silicone oil is experimentally investigated in 1-g. A spectral element based steady 2D numerical model is also developed. The traces of particles from experimental flow visualization and the corresponding computed streamlines are shown. At the upstream side of the bubble facing the cross-flow, thermocapillary and forced convection create liquid motion away from the wall, up along the surface. At the downstream side, a competing interaction between the two creates a recirculation cell, causing the bulk liquid to stagnate on the surface and separate thereafter. The important dimensionless parameters are - Surface tension and local cross-flow Reynolds numbers R(sub sigma)(U(sub ref)a/v) and Re(sub loc)(U(sub L)a/v), respectively based on reference thermocapillary U(sub ref)(sigma(sub T)Delta.T/mu, Delta T=T(sub wall)-T(sub liquid)) and local cross-flow velocity U(sub L), Prandtl number P(sub r) and Grashoff number Gr(rho.g.beta.DeltaTa(sup 3)/mu.v). Variation of the stagnation point with R(sub sigma) and Re(sub loc) is shown. Also shown are good agreement between experimental and numerical results in 1-g. The computational model is extended to mu-g condition to investigate temperature and velocity on the bubble surface, stagnation and reattachment points of the recirculation cell and wall heat transfer. It is observed that wall heating significantly alters the flow field around the bubble and thus the forces acting on the bubble, which control its detachment. Thus a combination of heating and liquid cross-flow can be utilized to precisely control bubble formation in a mu-g environment.

Bhunia, Avijit↗

A model of the wall boundary layer for ducted propellers

The objective of the present study is to include a representation of a wall boundary layer in an existing finite element model of the propeller in the wind tunnel environment. The major consideration is that the new formulation should introduce only modest alterations in the numerical model and should still be capable of producing economical predictions of the radiated acoustic field. This is accomplished by using a stepped approximation in which the velocity profile is piecewise constant in layers. In the limit of infinitesimally thin layers, the velocity profile of the stepped approximation coincides with that of the continuous profile. The approach described here could also be useful in modeling the boundary layer in other duct applications, particularly in the computation of the radiated acoustic field for sources contained in a duct.

Eversman, Walter↗

The role of entrainment by falling raindrops in microbursts

The numerical model of Krueger et al. (1986) for dry microburst simulations is used to study the role of entrainment by falling raindrops. Two series of numerical simulations were conducted: a control series, and a series with the raindrop fall speed set to zero so that the rain moved with the air instead of falling through it. The results show that entrainment due to falling raindrops helps microbursts with large raindrop mixing ratios to form in stable stratifications. It is found that entrainment appears to contribute to the small spatial and temporal scales that characterize microburst outflows.

Krueger, Steven K.↗

High Strain Rate Response of Adhesively Bonded Fiber-Reinforced Composite Joints – A Computational Study to Guide Experimental Design

Adhesively bonded carbon fiber-reinforced epoxy composite laminates are widely used in aerospace applications. During a high energy impact event, these laminates are often subjected to high strain rate loading. However, the influence of high strain rate loading on the response of these composite joints is not well understood. Computational finite element (FE) modeling and simulations are conducted to guide the design of high strain rate experiments. Two different experimental designs based on split Hopkinson bar were numerically modeled to simulate Mode I and Mode II types loading in the composite. In addition, the computational approach adopted in this study helps in understanding the high strain rate response of adhesively bonded composite joints subjected to nominally Mode I and Mode II loading. The modeling approach consists of a ply-level 3D FE model, a progressive damage constitutive model for the composite material behavior and a cohesive tie-break contact element for interlaminar delamination.

Ravindran, Suraj↗

Modeling the Gas Dynamics Environment in a Subscale Solid Rocket Test Motor

Subscale test motors are often used for the evaluation of solid rocket motor component materials such as internal insulation. These motors are useful for characterizing insulation performance behavior, screening insulation material candidates and obtaining material thermal and ablative property design data. One of the primary challenges associated with using subscale motors however, is the uncertainty involved when extrapolating the results to full-scale motor conditions. These uncertainties are related to differences in such phenomena as turbulent flow behavior and boundary layer development, propellant particle interactions with the wall, insulation off-gas mixing and thermochemical reactions with the bulk flow, radiation levels, material response to the local environment, and other anomalous flow conditions. In addition to the need for better understanding of physical mechanisms, there is also a need to better understand how to best simulate these phenomena using numerical modeling approaches such as computational fluid dynamics (CFD). To better understand and model interactions between major phenomena in a subscale test motor, a numerical study of the internal flow environment of a representative motor was performed. Simulation of the environment included not only gas dynamics, but two-phase flow modeling of entrained alumina particles like those found in an aluminized propellant, and offgassing from wall surfaces similar to an ablating insulation material. This work represents a starting point for establishing the internal environment of a subscale test motor using comprehensive modeling techniques, and lays the groundwork for improving the understanding of the applicability of subscale test data to full-scale motors. It was found that grid resolution, and inclusion of phenomena in addition to gas dynamics, such as two-phase and multi-component gas composition are all important factors that can effect the overall flow field predictions.

Eaton, Andrew M.↗

Topology of large-scale structure. IV - Topology in two dimensions

In a recent series of papers, an algorithm was developed for quantitatively measuring the topology of the large-scale structure of the universe and this algorithm was applied to numerical models and to three-dimensional observational data sets. In this paper, it is shown that topological information can be derived from a two-dimensional cross section of a density field, and analytic expressions are given for a Gaussian random field. The application of a two-dimensional numerical algorithm for measuring topology to cross sections of three-dimensional models is demonstrated.

Melott, Adrian L.↗