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,459 records · Page 81

The 3D modeling of high numerical aperture imaging in thin films

A modelling technique is described which is used to explore three dimensional (3D) image irradiance distributions formed by high numerical aperture (NA is greater than 0.5) lenses in homogeneous, linear films. This work uses a 3D modelling approach that is based on a plane-wave decomposition in the exit pupil. Each plane wave component is weighted by factors due to polarization, aberration, and input amplitude and phase terms. This is combined with a modified thin-film matrix technique to derive the total field amplitude at each point in a film by a coherent vector sum over all plane waves. Then the total irradiance is calculated. The model is used to show how asymmetries present in the polarized image change with the influence of a thin film through varying degrees of focus.

Flagello, D. G.↗

Convection phenomena in low-gravity processing - The GTE GaAs space experiment

A numerical model for the proposed GTE GaAs space experiment is presented. The model includes the actual ampule design and the experimentally obtained furnace thermal profiles. Results are presented based on both the traditional steady-state models and on fully time dependent simulations. It is demonstrated that the commonly presumed axisymmetric flows in numerical studies are physically accurate only when the local gravitational acceleration vector is parallel to the ampules's axis. Other gravitational orientations result in nontrivial three-dimensional flows which must be resolved using full three-dimensional models. Implications of the resulting flow and thermal fields on the growing crystal at conditions available during space processing are also discussed.

Arnold, W. A.↗

Cascade Models in Simulation of Extended Heavy Targets Irradiated by Accelerated Proton and Deuteron Beams

The paper presents a survey of the main numerical models used for simulation of interaction of accelerated particle beams with target nuclei. These models form the core of the software for simulation of various experiments and experimental facilities both for scientific and applied purposes. The beam and target parameters considered in detail in this study (protons and deuterons with energies from 0.66 to 4 AGeV and bulk U targets) cover the range of interest in development of new concepts of nuclear power production aided by accelerated particle beams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Simulation of the Impact of New Aircraft- and Satellite-based Ocean Surface Wind Measurements on Estimates of Hurricane Intensity

The Hurricane Imaging Radiometer (HIRAD) is a new airborne microwave remote sensor currently under development to enhance real-time hurricane ocean surface wind observations. HIRAD builds on the capabilities of the Stepped Frequency Microwave Radiometer (SFMR), which now operates on NOAA P-3, G-4, and AFRC C-130 aircraft. Unlike the SFMR, which measures wind speed and rain rate along the ground track directly beneath the aircraft, HIRAD will provide images of the surface wind and rain field over a wide swath (approximately 3 times the aircraft altitude). To demonstrate potential improvement in the measurement of peak hurricane winds, we present a set of Observing System Simulation Experiments (OSSEs) in which measurements from the new instrument as well as those from existing platforms (air, surface, and space-based) are simulated from the output of a high-resolution (approximately 1.7 km) numerical model. Simulated retrieval errors due to both instrument noise as well as model function accuracy are considered over the expected range of incidence angles, wind speeds and rain rates. Based on numerous simulated flight patterns and data source combinations, statistics are developed to describe relationships between the observed and true (from the model s perspective) peak wind speed. These results have implications for improving the estimation of hurricane intensity (as defined by the peak sustained wind anywhere in the storm), which may often go un-observed due to sampling limitations.

Uhlhorn, Eric↗

Forecasts of time averages with a numerical weather prediction model

Forecasts of time averages of 1-10 days in duration by an operational numerical weather prediction model are documented for the global 500 mb height field in spectral space. Error growth in very idealized models is described in order to anticipate various features of these forecasts and in order to anticipate what the results might be if forecasts longer than 10 days were carried out by present day numerical weather prediction models. The data set for this study is described, and the equilibrium spectra and error spectra are documented; then, the total error is documented. It is shown how forecasts can immediately be improved by removing the systematic error, by using statistical filters, and by ignoring forecasts beyond about a week. Temporal variations in the error field are also documented.

Roads, J. O.↗

Measurement of flow properties coupled to experimental and numerical analyses of dense, granular flows for solar thermal energy storage

Granular flows of sintered bauxite proppants were researched along an inclined plane for solar thermal energy storage applications. Granular flow properties needed to drive numeric granular models were measured for improved numerical model predictions for Carbobead CP 50/140, 40/100, and 30/60 particles. Particle shape and size distributions were determined by coupling optical microscopy to an in-house image processing algorithm. The impulse excitation technique was used to measure elastic and shear moduli, and compute Poisson’s ratio. The coefficient of static sliding friction was measured using the slip-stick method, and the static rolling friction was determined from measured shear on particles positioned between two hot-pressed plates. The coefficient of restitution was measured by dropping particles on a surface and determining the kinetic energy before and after impact with the surface using high resolution particle tracking velocimetry. Particle size did not significantly impact the coefficients of restitution and static rolling friction, however, particle shape distribution resulted in a large variation in measurements. An inclined flow experiment was performed to characterize granular flows of Carbobead CP 30/60 particles using particle image velocimetry. Numerical models of the experiment using discrete element method were generated with the measured mechanical properties as inputs for comparison with experimental results. A constant directional torque rolling friction model best predicted bulk granular flow behavior. Furthermore, good agreement between the model and experiment was achieved at ambient, steady state conditions, with average velocity differences <10%.

14 SOLAR ENERGY↗

An investigation to improve the Menhaden fishery prediction and detection model through the application of ERTS-A data

The author has identified the following significant results. Linear regression of secchi disc visibility against number of sets yielded significant results in a number of instances. The variability seen in the slope of the regression lines is due to the nonuniformity of sample size. The longer the period sampled, the larger the total number of attempts. Further, there is no reason to expect either the influence of transparency or of other variables to remain constant throughout the season. However, the fact that the data for the entire season, variable as it is, was significant at the 5% level, suggests its potential utility for predictive modeling. Thus, this regression equation will be considered representative and will be utilized for the first numerical model. Secchi disc visibility was also regressed against number of sets for the three day period September 27-September 29, 1972 to determine if surface truth data supported the intense relationship between ERTS-1 identified turbidity and fishing effort previously discussed. A very negative correlation was found. These relationship lend additional credence to the hypothesis that ERTS imagery, when utilized as a source of visibility (turbidity) data, may be useful as a predictive tool.

Maughan, P. M.↗

Measurement of spray combustion processes

A free jet configuration was chosen for measuring noncombusting spray fields and hydrocarbon-air spray flames in an effort to develop computational models of the dynamic interaction between droplets and the gas phase and to verify and refine numerical models of the entire spray combustion process. The development of a spray combustion facility is described including techniques for laser measurements in spray combustion environments and methods for data acquisition, processing, displaying, and interpretation.

Peters, C. E.↗

Speculations on the origin and evolution of the Utopia-Elysium lowlands of Mars

This paper proposes a qualitative model for the origin of the Utopia-Elysium northern lowlands on eastern Mars in terms of the long-term evolution of two large overlapping impact basins. The model, which is consistent with both the observed geologic constraints and more quantitative results obtained by numerical modeling of smaller (Orientale-size) impact basins, is shown to qualitatively account for the major topographic variation seen in the Utopia-Elysium region, including the overall 'lowness' of the area and localized depressions.

Frey, Herbert V.↗

Effect of thermal convection on the shape of a solid-liquid interface

The effect of thermal convection on the shape of solid-liquid interface was investigated in experiments conducted in a transparent Bridgman-type directional solidification furnace. The relationship was numerically modeled using a standard 2D finite-difference approach, with the solid-liquid deformable interface approximated by a blocking-off technique. The directional solidification furnace was used with pure succinonitrile (which is also transparent) contained in a long square ampoule made of borosilicate glass. With the furnace in the vertical configuration, a flat interface was observed, in agreement with the model. On the other hand, a highly distorted interface was obtained in the horizontal configuration; the numerical results showed a strong recirculating cell in front of the interface due to natural thermal convection. The results indicate that thermal convection is responsible for the interface distortion.

Mennetrier, C.↗

Evolution of Nonlinear Internal Waves in China Seas

Synthetic Aperture Radar (SAR) images from ERS-I have been used to study the characteristics of internal waves of Taiwan in the East China Sea, and east of Hainan Island in the South China Sea. Rank-ordered packets of internal solitons propagating shoreward from the edge of the continental shelf were observed in the SAR images. Based on the assumption of a semidiurnal tidal origin, the wave speed can be estimated and is consistent with the internal wave theory. By using the SAR images and hydrographic data, internal waves of elevation have been identified in shallow water due to a thicker mixed layer as compared with the bottom layer on the continental shelf. The generation mechanism includes the influences of the tide and the Kuroshio intrusion across the continental shelf for the formations of elevation internal waves. The effects of water depth on the evolution of solitons and wave packets are modeled by nonlinear Kortweg-deVries (KdV) type equation and linked to satellite image observations. The numerical calculations of internal wave evolution on the continental shelf have been performed and compared with the SAR observations. For a case of depression waves in deep water, the solitons first disintegrate into dispersive wave trains and then evolve to a packet of elevation waves in the shallow water area after they pass through a turning point of approximately equal layer depths has been observed in the SAR image and simulated by numerical model.

Liu, Antony K.↗

The Mode Localization Phenomenon: Visual Test Results

The high degree of manufacturing precision required in space structures may make them moresusceptible to mode localization. In this demonstration, the mode localization phenomenon isinvestigated using an existing full scale antenna-like structure having the complexity, manufacturingprecision, and dynamic characteristics of space structures. The actual imperfections were measuredand incorporated in a numerical model of the structure. The impact of several parameters thatprovoke the phenomenon were first studied by numerical simulations, then experimentally. Theresults show excellent agreement, thereby confirming the importance of mode localization in thedesign of actual space structure.

Levine-West, M.↗

Aircraft Flaps Modeling in OpenMDAO

The goal of this project was to develop a model for a single subsystem in the aerodynamics discipline, in this case the flaps of an aircraft. Flaps are high-lift devices used by planes to allow for quicker takeoffs, and slower landings. A computer-based aircraft model of the flaps of an aircraft was developed using the Python based open-source framework OpenMDAO. OpenMDAO is used to develop multi-disciplinary aircraft models using gradient-based optimization; design optimization (MDO) is concerned with solving design problems involving numerical models of complex engineering systems. There were 4 components in the model; each has input values, output variables, and equations to calculate said outputs. The variables and equations are sourced from NASA Fortran code from the 1970s, in a project called the General Aviation Synthesis Program (GASP). These variables and equations which create the model are being converted to Python for ease of use. The flaps model developed will be integrated into a larger model of a conventional aircraft’s flight phases. All subsystems of the model will first be built using the parameters of a Boeing 737 MAX-8, to validate its functionality and accuracy. Then, the aircraft model will be used for hybrid-electric research; running optimizations to improve efficiency, minimize fuel burn, and advance hybrid-electric technology in the aerospace field.

Computer-based aircraft modeling↗

Normal modes of the world's oceans: A numerical investigation using Proudman functions

The numerical modeling of the normal modes of the global oceans is addressed. The results of such modeling could be expected to serve as a guide in the analysis of observations and measurements intended to detect these modes. The numerical computation of normal modes of the global oceans is a field in which several investigations have obtained results during the past 15 years. The results seem to be model-dependent to an unsatisfactory extent. Some modeling areas, such as higher resolution of the bathymetry, inclusion of self-attraction and loading, the role of the Arctic Ocean, and systematic testing by means of diagnostic models are addressed. The results show that the present state of the art is such that a final solution to the normal mode problem still lies in the future. The numerical experiments show where some of the difficulties are and give some insight as to how to proceed in the future.

Sanchez, Braulio V.↗

A simple model describing the nonlinear dynamics of the dusk/dawn asymmetry in the high-latitude thermospheric flow

A simple nonlinear, axisymmetric, shallow-water numerical model has been used to study the asymmetry in the neutral flow between the dusk and dawn sides of the auroral oval. The results indicate that the Coriolis force and the curvature terms are nearly in balance on the evening side and require only a small pressure gradient to effect adjustment. The result is smaller neutral velocities near dawn and larger velocities near dusk than would be the case for a linearized treatment. A consequence is that more gravity wave energy is produced on the morning side than on the evening side.

Gundlach, J. P.↗

Radiation breakage of DNA: a model based on random-walk chromatin structure

Monte Carlo computer software, called DNAbreak, has recently been developed to analyze observed non-random clustering of DNA double strand breaks in chromatin after exposure to densely ionizing radiation. The software models coarse-grained configurations of chromatin and radiation tracks, small-scale details being suppressed in order to obtain statistical results for larger scales, up to the size of a whole chromosome. We here give an analytic counterpart of the numerical model, useful for benchmarks, for elucidating the numerical results, for analyzing the assumptions of a more general but less mechanistic "randomly-located-clusters" formalism, and, potentially, for speeding up the calculations. The equations characterize multi-track DNA fragment-size distributions in terms of one-track action; an important step in extrapolating high-dose laboratory results to the much lower doses of main interest in environmental or occupational risk estimation. The approach can utilize the experimental information on DNA fragment-size distributions to draw inferences about large-scale chromatin geometry during cell-cycle interphase.

Chromatin/radiation effects↗

Crack opening calculation in phase-field modeling of fluid-filled fracture: A robust and efficient strain-based method

The phase-field method has become popular for the numerical modeling of fluid-filled fractures, thanks to its ability to represent complex fracture geometry without algorithms. However, the algorithm-free representation of fracture geometry poses a significant challenge in calculating the crack opening (aperture) of phase-field fracture, which governs the fracture permeability and hence the overall hydromechanical behavior. Although several approaches have been devised to compute the crack opening of phase-field fracture, they require a sophisticated algorithm for post-processing the phase-field values or an additional parameter sensitive to the element size and alignment. Here, we develop a novel method for calculating the crack opening of fluid-filled phase-field fracture, which enables one to obtain the crack opening without additional algorithms or parameters. Here we transform the displacement-jump-based kinematics of a fracture into a continuous strain-based version, insert it into a force balance equation on the fracture, and apply the phase-field approximation. Through this procedure, we obtain a simple equation for the crack opening which can be calculated with quantities at individual material points. We verify the proposed method with analytical and numerical solutions obtained based on discrete representations of fractures, demonstrating its capability to calculate the crack opening regardless of the element size or alignment.

58 GEOSCIENCES↗

Development and Experimental Validation of a Path-Dependent Spin Forming Finite Element Model

Spin forming is an advanced manufacturing process widely used in the aerospace and defense sectors to produce lightweight, high-strength cylindrical components with tight dimensional tolerances. This study explores the applicability of the path-dependent Mechanical Threshold Stress (MTS) constitutive model by simulating the evolution of geometry, machining forces, and plastic deformation during the spin forming of a 10-mm thick 6061-O aluminum cylinder. While numerical modeling of spin forming has advanced substantially over the past decade, systematic verification and experimental validation of material models remain limited, particularly in predicting through-thickness process evolution. The MTS model, incorporating a Voce hardening rule, is employed for its ability to represent cyclic loading, rapidly varying temperature fields, and strain rates characteristic of spin forming. Numerical convergence analysis indicates discretization uncertainties between 0.3% and 9.2% for key quantities of interest. Experimental validation demonstrates that the MTS model, when implemented with a verified mesh, accurately reproduces both elastic and plastic behavior of 6061-O aluminum, predicting peak roller loads within 11–18% of measurements, geometric tolerances within 3%, and plastic strain distributions within 10% of experimental values. Collectively, these results establish a validated computational framework for predictive spin-forming simulations with quantified confidence, providing a foundation for extension to other alloys, geometries, and forming conditions.

Spin forming↗