Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EQUATIONS - STRESSES PLATES”

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

A one-equation turbulence transport model for high Reynolds number wall-bounded flows

A one-equation turbulence model that avoids the need for an algebraic length scale is derived from a simplified form of the standard k-epsilon model equations. After calibration based on well established properties of the flow over a flat plate, predictions of several other flows are compared with experiment. The preliminary results presented indicate that the model has predictive and numerical properties of sufficient interest to merit further investigation and refinement. The one-equation model is also analyzed numerically and robust solution methods are presented.

Baldwin, Barrett S.↗

A one-equation turbulence transport model for high Reynolds number wall-bounded flows

A one-equation turbulence model that avoids the need for an algebraic length scale is derived from a simplified form of the standard-k-epsilon model equations. After calibration based on well established properties of the flow over a flat plate, predictions of several other flows are compared with experiment. The preliminary results presented indicate that the model has predictive and numerical properties of sufficient interest to merit further investigation and refinement. The one-equation model is also analyzed numerically and robust solution methods are presented.

Baldwin, Barrett S.↗

Simulations Of Three-Dimensional Flow With Separation

Report discusses numerical solutions of Navier-Stokes equations for viscous supersonic flow over fin protruding from flat plate. Two types of configurations considered: one with wedge-shaped fin (sharp edge facing upstream), other with blunt (rectangular parallelepiped) fin. Only steady flows considered. Computations performed to study phenomena involved in separations of three-dimensional flows.

Hung, Ching-Mao↗

Evaluation of two-equation turbulence models for predicting transitional flows

An evaluation of the capabilities and limitations of four low-Reynolds number two-equation turbulence models for predicting bypass transition on a flat plate has been conducted. The ability of these models to reproduce the effect of the free-stream turbulence on transition was tested. The strengths and deficiencies of these models have been identified systematically. It is found that the k-epsilon models are capable of predicting the qualitative aspects of transition. However, the predictions are found to be sensitive to the initial conditions. Also, the start and end of transition were found to depend on the damping functions used in the low-Reynolds number versions of the k-epsilon models.

Abid, Ridha↗

Method of successive approximations for the solution of certain problems in aerodynamics

A method of successive approximations for the solution of problems in the fields of diffusion, boundary-layer flow, and heat-transfer is illustrated by solving problems in each of these fields. In most of the examples, the approximate solutions are compared with known accurate solutions and the agreement is shown to be good.

FLOW, TURBULENT - DIFFUSION↗

Validation of Heat Transfer and Film Cooling Capabilities of the 3-D RANS Code TURBO

The capabilities of the 3-D unsteady RANS code TURBO have been extended to include heat transfer and film cooling applications. The results of simulations performed with the modified code are compared to experiment and to theory, where applicable. Wilcox s k-turbulence model has been implemented to close the RANS equations. Two simulations are conducted: (1) flow over a flat plate and (2) flow over an adiabatic flat plate cooled by one hole inclined at 35 to the free stream. For (1) agreement with theory is found to be excellent for heat transfer, represented by local Nusselt number, and quite good for momentum, as represented by the local skin friction coefficient. This report compares the local skin friction coefficients and Nusselt numbers on a flat plate obtained using Wilcox's k-model with the theory of Blasius. The study looks at laminar and turbulent flows over an adiabatic flat plate and over an isothermal flat plate for two different wall temperatures. It is shown that TURBO is able to accurately predict heat transfer on a flat plate. For (2) TURBO shows good qualitative agreement with film cooling experiments performed on a flat plate with one cooling hole. Quantitatively, film effectiveness is under predicted downstream of the hole.

Shyam, Vikram↗

Learning robust parameter inference and density reconstruction in flyer plate impact experiments

Estimating physical parameters or material properties from experimental observations is a common objective in many areas of physics and material science. In many experiments, especially in shock physics, radiography is the primary means of observing the system of interest. However, radiography does not provide direct access to key state variables, such as density, which prevents the application of traditional parameter estimation approaches. Here we focus on flyer plate impact experiments on porous materials, and resolving the underlying parameterized equation of state (EoS) and crush porosity model parameters given radiographic observation(s). We use machine learning as a tool to demonstrate with high confidence that using only high impact velocity data does not provide sufficient information to accurately infer both EoS and crush model parameters, even with fully resolved density fields or a dynamic sequence of images. We thus propose an observable data set consisting of low and high impact velocity experiments/simulations that capture different regimes of compaction and shock propagation, and proceed to introduce a generative machine learning approach which produces a posterior distribution of physical parameters directly from radiographs. We demonstrate the effectiveness of the approach in estimating parameters from simulated flyer plate impact experiments, and show that the obtained estimates of EoS and crush model parameters can then be used in hydrodynamic simulations to obtain accurate and physically admissible density reconstructions. Finally, we examine the robustness of the approach to model mismatches, and find that the learned approach can provide useful parameter estimates in the presence of out-of-distribution radiographic noise and previously unseen physics, thereby promoting a potential breakthrough in estimating material properties from experimental radiographic images.

97 MATHEMATICS AND COMPUTING↗

Finite Element and Plate Theory Modeling of Acoustic Emission Waveforms

A comparison was made between two approaches to predict acoustic emission waveforms in thin plates. A normal mode solution method for Mindlin plate theory was used to predict the response of the flexural plate mode to a point source, step-function load, applied on the plate surface. The second approach used a dynamic finite element method to model the problem using equations of motion based on exact linear elasticity. Calculations were made using properties for both isotropic (aluminum) and anisotropic (unidirectional graphite/epoxy composite) materials. For simulations of anisotropic plates, propagation along multiple directions was evaluated. In general, agreement between the two theoretical approaches was good. Discrepancies in the waveforms at longer times were caused by differences in reflections from the lateral plate boundaries. These differences resulted from the fact that the two methods used different boundary conditions. At shorter times in the signals, before reflections, the slight discrepancies in the waveforms were attributed to limitations of Mindlin plate theory, which is an approximate plate theory. The advantages of the finite element method are that it used the exact linear elasticity solutions, and that it can be used to model real source conditions and complicated, finite specimen geometries as well as thick plates. These advantages come at a cost of increased computational difficulty, requiring lengthy calculations on workstations or supercomputers. The Mindlin plate theory solutions, meanwhile, can be quickly generated on personal computers. Specimens with finite geometry can also be modeled. However, only limited simple geometries such as circular or rectangular plates can easily be accommodated with the normal mode solution technique. Likewise, very limited source configurations can be modeled and plate theory is applicable only to thin plates.

Prosser, W. H.↗

Theoretical modeling of a bottom-raised oscillating surge wave energy converter structural loadings and power performances

Here, this study presents theoretical formulations to evaluate the fundamental parameters and performance characteristics of a bottom-raised oscillating surge wave energy converter (OSWEC) device. Employing a flat plate assumption and potential flow formulation in elliptical coordinates, closed-form equations for the added mass, radiation damping, and excitation forces/torques in the relevant pitch-pitch and surge-pitch directions of motion are developed and used to calculate the system's response amplitude operator and the forces and moments acting on the foundation. The model is benchmarked against numerical simulations using WAMIT and WEC-Sim, showcasing excellent agreement. The sensitivity of plate thickness on the analytical hydrodynamic solutions is investigated over several thickness-to-width ratios ranging from 1:80 to 1:10. The results show that as the thickness of the benchmark OSWEC increases, the deviation of the analytical hydrodynamic coefficients from the numerical solutions grows from 3% to 25%. Differences in the excitation forces and torques, however, are contained within 12%. While the flat plate assumption is a limitation of the proposed analytical model, the error is within a reasonable margin for use in the design space exploration phase before a higher-fidelity (and thus more computationally expensive) model is employed. A parametric study demonstrates the ability of the analytical model to quickly sweep over a domain of OSWEC dimensions, illustrating the analytical model's utility in the early phases of design.

13 HYDRO ENERGY↗

Mechanical loading prediction through accelerometry data during walking and running

ABSTRACT Currently, there is no way to assess mechanical loading variables such as peak ground reaction forces (pGRF) and peak loading rate (pLR) in clinical settings. The purpose of this study was to develop accelerometry‐based equations to predict both pGRF and pLR during walking and running. One hundred and thirty one subjects (79 females; 76.9 ± 19.6 kg) walked and ran at different speeds (2–14 km·h −1 ) on a force plate–instrumented treadmill while wearing accelerometers at their ankle, lower back and hip. Regression equations were developed to predict pGRF and pLR from accelerometry data. Leave‐one‐out cross‐validation was used to calculate prediction accuracy and Bland–Altman plots. Our pGRF prediction equation was compared with a reference equation previously published. Body mass and peak acceleration were included for pGRF prediction and body mass and peak acceleration rate for pLR prediction. All pGRF equation coefficients of determination were above 0.96, and a good agreement between actual and predicted pGRF was observed, with a mean absolute percent error (MAPE) below 7.3%. Accuracy indices from our equations were better than previously developed equations. All pLR prediction equations presented a lower accuracy compared to those developed to predict pGRF. Walking and running pGRF can be predicted with high accuracy by accelerometry‐based equations, representing an easy way to determine mechanical loading in free‐living conditions. The pLR prediction equations yielded a somewhat lower prediction accuracy compared with the pGRF equations.

Veras, Lucas↗

Numerical simulation of crossing/turbulent boundary layer interaction at Mach 8.3 comparison of zero and two-equation turbulence models

A 3D hypersonic crossing shock wave/turbulent boundary layer interaction is examined numerically. The test geometry consists of a pair of opposing sharp fins of angle alpha = 15 deg mounted on a flat plate. The freestream Mach number is 8.28. Two theoretical models are evaluated. The full 3D Reynolds-averaged Navier-Stokes equations are solved using the Baldwin-Lomax algebraic turbulent eddy viscosity model and the Rodi turbulence model. Computed results for both cases show good agreement with experiment for flat plate surface pressure and for pitot pressure and yaw angle profiles in the flowfield. General agreement is obtained for surface flow direction. Fair to poor agreement is obtained for surface heat transfer, indicating a need for more accurate turbulence models. The overall flowfield structure is similar to that observed in previous crossing shock interaction studies.

Narayanswami, N.↗

Investigation of vortex breakdown on a delta wing using Euler and Navier-Stokes equations

A numerical investigation of leading-edge vortex breakdown on a delta wing at high angles of attack is presented. The analysis has been restricted to low speed flows on a flat plate wing with sharp leading edges. Both Euler and Navier-Stokes (assuming fully laminar and fully turbulent flows) equations have been used in this study and the results are compared against experimental data. Predictions of vortex breakdown progression with angle of attack with both Euler and Navier-Stokes equations are shown to be consistent with the experimental data. However, the Navier-Stokes predictions show significant improvements in breakdown location at angles of attack where the vortex breakdown approaches the wing apex. The location of the primary vortex and the level of vorticity in the pre-breakdown regions in these flowfield solutions are affected very little by the viscous effects, even though the Navier-Stokes solutions exhibit viscous phenomena such as secondary and tertiary vortices which are not present in the Euler solutions. In the post-breakdown regions, however, the levels of vorticity in the primary vortex have increased differences between the Euler and Navier-Stokes solutions at comparable locations. Navier-Stokes indicates the presence of a secondary vortex even after the primary vortex is burst. The predicted trajectories of the primary vortex are in very good agreement with the test data, the laminar solutions providing the overall best comparison. The Euler shows a small displacement of the primary vortex, relative to experiment, due to the lack of secondary vortices. The turbulent Navier-Stokes, in general, fall between the Euler and laminar solutions. These findings are based on solutions from meshes that are not usually considered fine enough for resolving vortical flows. To further understand the vortex breakdown phenomenon and the fine details of the vortical flow structure, further calculations with the finer meshes and improved turbulence models are necessary.

S. Agrawal↗

An initial investigation into methods of computing transonic aerodynamic sensitivity coefficients

The primary accomplishments of the project are as follows: (1) Using the transonic small perturbation equation as a flowfield model, the project demonstrated that the quasi-analytical method could be used to obtain aerodynamic sensitivity coefficients for airfoils at subsonic, transonic, and supersonic conditions for design variables such as Mach number, airfoil thickness, maximum camber, angle of attack, and location of maximum camber. It was established that the quasi-analytical approach was an accurate method for obtaining aerodynamic sensitivity derivatives for airfoils at transonic conditions and usually more efficient than the finite difference approach. (2) The usage of symbolic manipulation software to determine the appropriate expressions and computer coding associated with the quasi-analytical method for sensitivity derivatives was investigated. Using the three dimensional fully conservative full potential flowfield model, it was determined that symbolic manipulation along with a chain rule approach was extremely useful in developing a combined flowfield and quasi-analytical sensitivity derivative code capable of considering a large number of realistic design variables. (3) Using the three dimensional fully conservative full potential flowfield model, the quasi-analytical method was applied to swept wings (i.e. three dimensional) at transonic flow conditions. (4) The incremental iterative technique has been applied to the three dimensional transonic nonlinear small perturbation flowfield formulation, an equivalent plate deflection model, and the associated aerodynamic and structural discipline sensitivity equations; and coupled aeroelastic results for an aspect ratio three wing in transonic flow have been obtained.

Carlson, Leland A.↗

One-dimensional analysis of plane and radial thin film flows including solid-body rotation

The flow of a thin liquid film with a free surface along a horizontal plate which emanates from a pressurized vessel is examined by integrating the equations of motion across the thin liquid layer and discretizing the integrated equations using finite difference techniques. The effects of 0-g and solid-body rotation will be discussed. The two cases of interest are plane flow and radial flow. In plane flow, the liquid is considered to be flowing along a channel with no change in the width of the channel, whereas in radial flow the liquid spreads out radially over a disk, so that the area changes along the radius. It is desired to determine the height of the liquid film at any location along the plate of disk, so that the heat transfer from the plate or disk can be found. The possibility that the flow could encounter a hydraulic jump is accounted for.

Thomas, S.↗

Numerical simulation of boundary layers. Part 1: Weak formulation and numerical method

A numerical method designed to solve the time-dependent, three-dimensional, incompressible Navier-Stokes equations in boundary layers is presented. The fluid domain is the half-space over a flat plate, and periodic conditions are applied in the horizontal directions. The discretization is spectral. The basis functions are divergence-free and a weak formulation of the momentum equation is used, which eliminates the pressure term. An exponential mapping and Jacobi polynomials are used in the semi-infinite direction, with the irrotational component receiving special treatment. Issues related to the accuracy, stability and efficiency of the method are discussed. Very fast convergence is demonstrated on some model problems with smooth solutions. The method has also been shown to accurately resolve the fine scales of transitional and turbulent boundary layers.

Spalart, P. R.↗

High-precision shock equation of state measurements for metallic fluid carbon between 15 and 20 Mbar

Diamond is an efficient ablator material to convert the energy of high-power giant lasers into ablation pressure with applications for High-Energy-Density (HED) science, planetary science, and Inertial Confinement Fusion (ICF) research at the National Ignition Facility (NIF). Unfortunately, current theoretical equation of state models cannot reproduce all the observed experimental data in the multi-megabar regime particularly relevant for HED and ICF research. New experimental data on the behavior of carbon at extreme pressures and temperatures are, therefore, essential to improve our predictive capability to design and analyze dynamic compression experiments for HED or ICF research and build improved equation of state models in the future. In this paper, we report high-precision laser-driven shock compression measurements on diamond single crystals at the Omega Laser Facility. Using ultrafast Doppler optical Velocimetry Interferometer System for Any Reflector (VISAR) to track the leading shock front and a quartz plate as an in situ reference, we obtain relative pressure-density shock equation-of-state measurements between 15 and 20 Mbar with an impedance-matching procedure. We also report shock-and-release measurements in a spherical geometry at the NIF. The new data provide tight constraints on the compressibility of warm dense carbon along the Hugoniot of full density diamond, allowing us to discriminate between existing theoretical equation-of-state models. We find that both LLNL LEOS 9061 and LANL Sesame 7835 models capture well the shock compressibility in the explored range. LANL Sesame 7835 also reproduces well the observed shock-and-release behavior of diamond near 10–20 Mbar.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Applications of the theory of optimal control of distributed-parameter systems to structural optimization

An extension of classical methods of optimal control theory for systems described by ordinary differential equations to distributed-parameter systems described by partial differential equations is presented. An application is given involving the minimum-mass design of a simply-supported shear plate with a fixed fundamental frequency of vibration. An optimal plate thickness distribution in analytical form is found. The case of a minimum-mass design of an elastic sandwich plate whose fundamental frequency of free vibration is fixed. Under the most general conditions, the optimization problem reduces to the solution of two simultaneous partial differential equations involving the optimal thickness distribution and the modal displacement. One equation is the uniform energy distribution expression which was found by Ashley and McIntosh for the optimal design of one-dimensional structures with frequency constraints, and by Prager and Taylor for various design criteria in one and two dimensions. The second equation requires dynamic equilibrium at the preassigned vibration frequency.

Armand, J. P.↗

Multi-metallic conduction cooled superconducting radio-frequency cavity with high thermal stability

Superconducting radio-frequency cavities are commonly used in modern particle accelerators for applied and fundamental research. Such cavities are typically made of high-purity, bulk Nb and with cooling by a liquid helium bath at a temperature of ~2 K. The size, cost and complexity of operating a particle accelerator with a liquid helium refrigerator make the current cavity technology not favorable for use in industrial-type accelerators. We have developed a multi-metallic 1.495 GHz elliptical cavity conductively cooled by a cryocooler. The cavity has a ~2 µm thick layer of Nb 3 Sn on the inner surface, exposed to the rf field, deposited on a ~3 mm thick bulk Nb shell and a bulk Cu shell, of thickness $\geqslant\! 5$ mm deposited on the outer surface by electroplating. A bolt-on Cu plate 1.27 cm thick was used to thermally connect the cavity equator to the second stage of a Gifford-McMahon cryocooler with a nominal capacity of 2 W at 4.2 K. The cavity was tested initially in liquid helium at 4.3 K and reached a peak surface magnetic field of ~36 mT with a quality factor of 2 × 10 9 . The cavity cooled by the cryocooler achieved a peak surface magnetic field of ~29 mT, equivalent to an accelerating gradient of 6.5 MV m –1 . The conduction-cooled cavity could be operated in continuous-wave with as high as 5 W dissipation in the cavity for 1 h without any thermal breakdown, because of the Cu outer layer with high thermal conductivity. This result represents a paradigm shift in the technology of superconducting accelerator cavities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗