Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure 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 541 records · Page 30

Modeling deep convection in the Greenland Sea

The development of deep convective events in the high-latitude ocean is studied using a three-dimensional, coupled ice-ocean model. Oceanic mixing is described according to the level 2.5 turbulence closure scheme in which convection occurs in a continuous way, i.e., convective adjustment is not invoked. The model is forced by strong winds and surface cooling. Strong upwelling at the multilyear ice edge and consequent entrainment of warm Atlantic waters into the mixed layer is produced by winds parallel to the ice edge. Concomitant cooling drives deep convection and produces chimneylike structures. Inclusion of a barotropic mean flow over topography to the model provides important preconditioning and selects the location of deep convection. The most efficient preconditioning occurs at locations where the flow ascends a slope. In a stratified environment similar to the Greenland Sea with a 12 m/s wind the model simulations show that localized deep convection takes place after about 10 days to depths of 1000 m.

Hakkinen, S.↗

Integral boundary-layer models for turbulent separated flows

Two-dimensional, integral boundary-layer methods of the entrainment and dissipation type, with closure relationships which account for turbulent separated flows through the velocity and turbulence models, are developed. Kinematic shape factor relations used in the two methods are evaluated by comparison with recently obtained experimental data for separated flows. Prescribed-displacement-thickness computations are made with both methods and compared to experiment for a series of separated flows, including low speed and transonic shock-induced separation. Effects of upstream history on the turbulent shear stress development are shown to be important for the accurate prediction of the separated flows considered and the experimental trends are reasonably modeled by a rate equation derived from the turbulent kinetic energy equation.

Thomas, J. L.↗

Delamination Modeling of Composites for Improved Crash Analysis

Finite element crash modeling of composite structures is limited by the inability of current commercial crash codes to accurately model delamination growth. Efforts are made to implement and assess delamination modeling techniques using a current finite element crash code, MSC/DYTRAN. Three methods are evaluated, including a straightforward method based on monitoring forces in elements or constraints representing an interface; a cohesive fracture model proposed in the literature; and the virtual crack closure technique commonly used in fracture mechanics. Results are compared with dynamic double cantilever beam test data from the literature. Examples show that it is possible to accurately model delamination propagation in this case. However, the computational demands required for accurate solution are great and reliable property data may not be available to support general crash modeling efforts. Additional examples are modeled including an impact-loaded beam, damage initiation in laminated crushing specimens, and a scaled aircraft subfloor structures in which composite sandwich structures are used as energy-absorbing elements. These examples illustrate some of the difficulties in modeling delamination as part of a finite element crash analysis.

Fleming, David C.↗

Delamination Modeling of Composites for Improved Crash Analysis

Finite element crash modeling of composite structures is limited by the inability of current commercial crash codes to accurately model delamination growth. Efforts are made to implement and assess delamination modeling techniques using a current finite element crash code, MSC/DYTRAN. Three methods are evaluated, including a straightforward method based on monitoring forces in elements or constraints representing an interface; a cohesive fracture model proposed in the literature; and the virtual crack closure technique commonly used in fracture mechanics. Results are compared with dynamic double cantilever beam test data from the literature. Examples show that it is possible to accurately model delamination propagation in this case. However, the computational demands required for accurate solution are great and reliable property data may not be available to support general crash modeling efforts. Additional examples are modeled including an impact-loaded beam, damage initiation in laminated crushing specimens, and a scaled aircraft subfloor structures in which composite sandwich structures are used as energy-absorbing elements. These examples illustrate some of the difficulties in modeling delamination as part of a finite element crash analysis.

Fleming, David C.↗

Aerosol–Ice Formation Closure: A Southern Great Plains Field Campaign

Prediction of ice formation in clouds presents one of the grand challenges in the atmospheric sciences. Immersion freezing initiated by ice-nucleating particles (INPs) is the dominant pathway of primary ice crystal formation in mixed-phase clouds, where supercooled water droplets and ice crystals coexist, with important implications for the hydrological cycle and climate. However, derivation of INP number concentrations from an ambient aerosol population in cloud-resolving and climate models remains highly uncertain. We conducted an aerosol-ice formation closure pilot study using a field-observational approach to evaluate the predictive capability of immersion freezing INPs. The closure study relies on co-located measurements of the ambient size-resolved and single-particle composition and INP number concentrations. The acquired particle data serve as input in several immersion freezing parameterizations, that are employed in cloud-resolving and climate models, for prediction of INP number concentrations. We discuss in detail one closure case study in which a front passed through the measurement site, resulting in a change of ambient particle and INP populations. We achieved closure in some circumstances within uncertainties, but we emphasize the need for freezing parameterization of potentially missing INP types and evaluation of the choice of parameterization to be employed. Overall, this closure pilot study aims to assess the level of parameter details and measurement strategies needed to achieve aerosol-ice formation closure. The closure approach is designed to accurately guide immersion freezing schemes in models, and ultimately identify the leading causes for climate model bias in INP predictions.

ice formation↗

Closure measurements of naturally initiating small cracks

The initiation and growth of cracks in smooth 2024-T3 aluminum specimens have been investigated using acetate replicas. Crack opening displacements were measured as a function of load at several positions along the crack as it grew across the thickness of the specimen and became a through-the-thickness crack of a few millimeters in length. Tests run for R-ratios of 0.5, 0.0 and -1.0 showed that closure loads were well predicted by Newman's (1984) model. It is noted that small cracks grow slightly faster than would be predicted by long-crack data, and it is suggested that effects other than plasticity-induced closure must be taken into account.

Sharpe, W. N., Jr.↗

A unified view of convective transports by stratocumulus clouds, shallow cumulus clouds, and deep convection

A bulk planetary boundary layer (PBL) model was developed with a simple internal vertical structure and a simple second-order closure, designed for use as a PBL parameterization in a large-scale model. The model allows the mean fields to vary with height within the PBL, and so must address the vertical profiles of the turbulent fluxes, going beyond the usual mixed-layer assumption that the fluxes of conservative variables are linear with height. This is accomplished using the same convective mass flux approach that has also been used in cumulus parameterizations. The purpose is to show that such a mass flux model can include, in a single framework, the compensating subsidence concept, downgradient mixing, and well-mixed layers.

Randall, David A.↗

A physically based fractional cloudiness parameterization

A bulk PBL model with a simple internal vertical structure and a simple second-order closure has been developed for use as a PBL parameterization in a large-scale model. A method to determine the value of sigma, the fractional area covered by rising motion, is given. Large-eddy simulations of Moeng (1984, 1986) are used to validate various aspects of the parameterization.

Randall, David A.↗

Effects of stores on longitudinal aerodynamic characteristics of a fighter at supersonic speeds

Experimental investigations of single and twin stores representative of advanced, elliptical cross section missile concepts were made at Mach numbers from 1.60 to 2.16 to substantiate theoretically predicted results. The stores were mounted on the fuselage of a model representing a fighter configuration. Store base closure effects in the carriage condition were also obtained through tests with and without base closure fairings.

Dollyhigh, S. M.↗

Numerical analysis of turbulent wakes of axial flow turbomachinery rotor blades

The wakes of rotor blades are turbulent, three-dimensional, and are subjected to curvature and rotation effects. The objective of this study is to predict the development of such wakes and compare the predictions with the existing data. A finite difference procedure is employed in the numerical analysis of the wake utilizing the continuity, momentum, and turbulence closure equations in the rotating curvilinear and non-orthogonal coordinate system. The turbulence closure is affected by the two equations model and the simplified Reynolds stress model. The predictions from both the models agree well with the mean velocity and the turbulence data.

Hah, C.↗

Navier-Stokes solutions for nonaxisymmetric nozzle flows

The compressible Navier-Stokes equations for two-dimensional and three-dimensional flow are solved for nonaxisymmetric nozzle configurations. Both internal flow alone and combined internal and external flow cases are considered. An implicit approximate factorization algorithm is employed to solve the governing equations which are cast in terms of a body-fitted coordinate system. Turbulence closure is accomplished with a composite algebraic model. Two-dimensional results for a wedge plug nozzle flow in the form of surface pressures compare favorably with experimental data. Three-dimensional calculations for two nonaxisymmetric nozzles are also presented and compared with data.

Swanson, R. C.↗

Computation of aircraft component flow fields at transonic Mach numbers using a three-dimensional Navier-Stokes algorithm

A computer analysis was developed for calculating steady (or unsteady) three-dimensional aircraft component flow fields. This algorithm, called ENS3D, can compute the flow field for the following configurations: diffuser duct/thrust nozzle, isolated wing, isolated fuselage, wing/fuselage with or without integrated inlet and exhaust, nacelle/inlet, nacelle (fuselage) afterbody/exhaust jet, complete transport engine installation, and multicomponent configurations using zonal grid generation technique. Solutions can be obtained for subsonic, transonic, or hypersonic freestream speeds. The algorithm can solve either the Euler equations for inviscid flow, the thin shear layer Navier-Stokes equations for viscous flow, or the full Navier-Stokes equations for viscous flow. The flow field solution is determined on a body-fitted computational grid. A fully-implicit alternating direction implicit method is employed for the solution of the finite difference equations. For viscous computations, either a two layer eddy-viscosity turbulence model or the k-epsilon two equation transport model can be used to achieve mathematical closure.

Shrewsbury, George D.↗

Critical comparison of second-order closures with direct numerical simulations of homogeneous turbulence

Recently, several models have been proposed for closing the second-moment equations, in which the velocity-pressure gradient tensor and the dissipation rate tensor are two of the most important terms. In the literature, these correlation tensors are usually decomposed into a so-called rapid term and a return-to-isotropy term. Models of these terms have been used in global flow calculations together with other modeled terms. However, their individual behaviors in different flows have not been fully examined because they are unmeasurable in the laboratory. Recently, the development of direct numerical simulation (DNS) of turbulence has given us the possibility to do this kind of study. With direct numerical simulation, we may use the solution to calculate exactly the values of these correlation terms and then directly compare them with the values from their modeled formulations. In this paper, we make direct comparisons of five representative rapid models and eight return-to-isotropy models using the DNS data of 45 homogeneous flows, which were done by Rogers et al. (1986) and Lee and Reynolds (1985). The purpose of these direct comparisons is to explore the performance of these models in different flows and identify the ones that give the best performance. The paper also describes the modeling procedure, model constraints, and the various evaluated models. The detailed results of the direct comparisons are discussed, and a few concluding remarks on turbulence models are given.

Shih, Tsan-Hsing↗

Observations of Circumstellar Disks with Infrared Interferometry

Star formation is arguably the area of astrophysics in which infrared interferometry has had the biggest impact. The optically thick portion of T Tauri and Herbig Ae/Be disks DO NOT extend to a few stellar radii of the stellar surface. Emission is coming from near the dust sublimation radius, but not all from a single radius. The Herbig Ae stars can be either flared or self-shadowed but very massive (early Be) stars are geometrically thin. The Herbig Ae stars can be either flared or self-shadowed but very massive (early Be) stars are geometrically thin. Observational prospects are rapidly improving: a) Higher spectral resolution will allow observations of the gas: jets, winds, accretion. b) Closure phase and imaging will help eliminate model uncertainties/dependencies.

young steller objects↗

Hybrid Reynolds-Averaged/Large Eddy Simulation of the Flow in a Model SCRamjet Cavity Flameholder

Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. Experimental data available for this configuration include velocity statistics obtained from particle image velocimetry. Several turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged/large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This e ort was undertaken to not only assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community, but to also begin to understand how this capability can best be used to augment standard Reynolds-averaged simulations. The numerical errors were quantified for the steady-state simulations, and at least qualitatively assessed for the scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results displayed a high degree of variability when comparing the flameholder fuel distributions obtained from each turbulence model. This prompted the consideration of applying the higher-fidelity scale-resolving simulations as a surrogate "truth" model to calibrate the Reynolds-averaged closures in a non-reacting setting prior to their use for the combusting simulations. In general, the Reynolds-averaged velocity profile predictions at the lowest fueling level matched the particle imaging measurements almost as well as was observed for the non-reacting condition. However, the velocity field predictions proved to be more sensitive to the flameholder fueling rate than was indicated in the measurements.

Baurle, R. A.↗

Galactic cosmic-ray mediation of a spherical solar wind flow. 1: The steady state cold gas hydrodynamical approximation

Realistic models of the outer heliosphere should consider that the interstellar cosmic-ray pressure becomes comparable to pressures in the solar wind at distances more than 100 AU from the Sun. The cosmic-ray pressure dynamically affects solar wind flow through deceleration. This effect, which occurs over a scale length of the order of the effective diffusion length at large radial distances, has important implications for cosmic-ray modulation and acceleration. As a first step toward solution of this nonlinear problem, a steady state numerical model was developed for a relatively cold spherical solar wind flow which encounters the confining isotropic pressure of the surrounding Galactic medium. This pressure is assumed to be dominated by energetic particles (Galactic cosmic rays). The system of equations, which are solved self-consistently, includes the relevant hydrodynamical equations for the solar wind flow and the spherical cosmic-ray transport equation. To avoid the closure parameter problem of the two-fluid model, the latter equation is solved for the energy-dependent cosmic-ray distribution function.

Le Roux, J. A.↗

Nonlinear fracture of concrete and ceramics

The nonlinear fracture process zones in an impacted unnotched concrete bend specimen, a prenotched ceramic bend specimen, and an unnotched ceramic/ceramic composite bend specimen were estimated through hybrid experimental numerical analysis. Aggregate bridging in concrete, particulate bridging in ceramics, and fiber bridging in ceramic/ceramic composite are modeled by Barenblatt-type cohesive zones which are incorporated into the finite-element models of the bend specimens. Both generation and propagation analyses are used to estimate the distribution of crack closure stresses in the nonlinear fracture process zones. The finite-element models are then used to simulate fracture tests consisting of rapid crack propagation in an impacted concrete bend specimen, and stable crack growth and strain softening in a ceramic and ceramic/ceramic composite bend specimens.

Kobayashi, Albert S.↗

On the divergence of the Hall current in the westward traveling surge

Recent observations suggest that the Hall current tends to be highly divergence-free in the westward traveling surge. It is shown in a model calculation that the blockage of the Hall current from closure along field lines leads to a significant reduction of the divergence of the Hall current, but is not divergence-free. The westward electrojet in this models is found to rotate counterclock wise around the leading edge of the surge and merge into the eastward electrojet. The difference in strength between the two electrojets at the head of the surge goes into the upward field-aligned current. It is also shown that the conductivity enhanced by the discrete auroral precipitation does not contribute to the divergence of the Hall current, provided that the electric fields are coplanar along a field line. The above results strongly suggest that the Hall current tends to be highly divergence-free in regions of bright auroral forms, especially in the westward traveling surge.

Kan, J. R.↗