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,603 records · Page 89

Glass formulation and composition optimization with property models: A review

Abstract Glass is a versatile material with a remarkable history and many practical applications. It plays a critical role in our everyday lives, the advancement of science, and the development of many technologies. The Edisonian type trial‐and‐error method was commonly used for conventional design of glass compositions, which was time‐consuming and costly. With the urgent need to develop new glass compositions for technology applications rapidly, it has become necessary to develop precise property models with predictive powers using large databases and efficient formulation approaches. This paper reviews the design of glass compositions using these analytical and numerical models of composition–structure–property relations of glasses, some based on large databases and machine learning approaches. Aspects of data collection, model fitting, feature extraction, model evaluation, and uncertainty quantification will be covered. Furthermore, advances in the glass optimization framework and available tools are summarized with examples. The outlook and perspective for further glass property model development and formulation approaches are discussed.

Lu, Xiaonan↗

Innovating the next generation of commercial smart building software

Nearly 30% of commercial building energy use is wasted due to equipment faults and HVAC controls problems. The result is increased emissions, compromised comfort and productivity, and less reliable coordination of building power needs with a clean grid. The energy impact alone represents $17 billion in potential savings. Today’s smart building software provides a robust solution to address these operational deficiencies. Energy management and information systems (EMIS) are saving up to 9% on average, with two-year paybacks. They are being incorporated into energy management processes, commissioning services, and utility programs. As effective as they are, two barriers prevent even deeper benefits; limited personnel to fix problems once they are identified, and the expense and time to manually implement changes in control systems. In partnership with the research community, the EMIS industry is developing new capabilities to overcome these barriers. Moving beyond siloed products for either fault detection and diagnostics, or optimal control, these new capabilities empower users to not only automatically identify faults, but also to push corrective action, and control improvements to their buildings. In this paper, several areas for enhancements are documented: ‘one-time’ correction of faults such as setpoints, schedules, and economizer lockouts; short-term active testing for automated proportional integral derivative (PID) loop tuning and functional testing; and continuous supervisory control for demand flexibility and year-round efficiency. Results are presented from a pair of partner implementations out of a dozen providers integrating these enhancements into their products, including field tests from across the country, and insights into operator acceptance and integration into operations and maintenance practices.

Casillas, Armando↗

LEO high voltage solar array arcing response model

A series of mathematical models were developed that describe the electrical behavior of a large solar cell array floating electrically in the low Earth orbit (LEO) space plasma and struck by an arc at a point of negative bias. There are now three models in this series: ARCII, which is a fully analytical, linearized model; ARCIII, which is an extension of ARCIII that includes solar cell inductance as well as load reactance; Nonlinear ARC, which is a numerical model able to treat effects such as non-linearized, i.e., logarithmic solar cell I/V characteristics, conductance switching as a solar cell crosses plasma ground on a voltage excursion and non-ohmic plasma leakage current collection.

Metz, Roger N.↗

Flare energetics

In this investigation of flare energetics, researchers sought to establish a comprehensive and self-consistent picture of the sources and transport of energy within a flare. To achieve this goal, they chose five flares in 1980 that were well observed with instruments on the Solar Maximum Mission, and with other space-borne and ground-based instruments. The events were chosen to represent various types of flares. Details of the observations available for them and the corresponding physical parameters derived from these data are presented. The flares were studied from two perspectives, the impulsive and gradual phases, and then the results were compared to obtain the overall picture of the energics of these flares. The role that modeling can play in estimating the total energy of a flare when the observationally determined parameters are used as the input to a numerical model is discussed. Finally, a critique of the current understanding of flare energetics and the methods used to determine various energetics terms is outlined, and possible future directions of research in this area are suggested.

Wu, S. T.↗

Experimental and analytic studies of the triggered lightning environment of the F106B

The triggered lightning environment of the F106B aircraft is investigated. Scale modeling of the F106B with a metallized model was done to measure electric field enhancement factors on the aircraft and on canonically shaped conducting objects. These are then compared to numerically determined quantities. Detailed numerical modeling is done of the development of the triggered lightning channel. This is done using nonlinear air chemistry models to model a variety of physical phenomena which occur in a triggered lightning event. The effect of a triggered lightning strike on internal wires in the F106B is investigated using finite difference models and transmission line models to calculate the electromagnetic coupling of lightning currents through seams and joints of the aircraft to internal cables. Time domain waveforms are computed and compared to measured waveforms. The effect of thunderstorm particles on the initial triggering of a lightning strike is investigated. The electric field levels needed to cause air breakdown in the presence and absence of thunderstorm particles are calculated. This is done as a function of the size, shape, and density of the particles.

Rudolph, Terence↗

Weight and cost forecasting for advanced manned space vehicles

A mass and cost estimating computerized methology for predicting advanced manned space vehicle weights and costs was developed. The user friendly methology designated MERCER (Mass Estimating Relationship/Cost Estimating Relationship) organizes the predictive process according to major vehicle subsystem levels. Design, development, test, evaluation, and flight hardware cost forecasting is treated by the study. This methodology consists of a complete set of mass estimating relationships (MERs) which serve as the control components for the model and cost estimating relationships (CERs) which use MER output as input. To develop this model, numerous MER and CER studies were surveyed and modified where required. Additionally, relationships were regressed from raw data to accommodate the methology. The models and formulations which estimated the cost of historical vehicles to within 20 percent of the actual cost were selected. The result of the research, along with components of the MERCER Program, are reported. On the basis of the analysis, the following conclusions were established: (1) The cost of a spacecraft is best estimated by summing the cost of individual subsystems; (2) No one cost equation can be used for forecasting the cost of all spacecraft; (3) Spacecraft cost is highly correlated with its mass; (4) No study surveyed contained sufficient formulations to autonomously forecast the cost and weight of the entire advanced manned vehicle spacecraft program; (5) No user friendly program was found that linked MERs with CERs to produce spacecraft cost; and (6) The group accumulation weight estimation method (summing the estimated weights of the various subsystems) proved to be a useful method for finding total weight and cost of a spacecraft.

Williams, Raymond↗

On the structure of gaseous confined laminar diffusion flames: Numerical investigation

The structure and characteristics of gaseous confined laminar diffusion flames are investigated by numerically solving the time-dependent two-dimensional axisymmetric conservation equations. The numerical model accounts for the important chemical and physical processes involved, including axial diffusion, viscous effects, radial convection, and finite-rate chemistry. The numerical results clearly show that the flame has a finite thickness and leakage of fuel vapor into the flame zone is possible. The effect of heat release is found to induce some radial flow. Predicted flame shape and dimensions are compared to the classical Burke-Schumann flame. The numerically calculated flame is observed to be about 15 percent taller and 5 percent narrower than that of the Burke-Schumann solution under the same conditions.

Mawid, M. A.↗

A Boundary Condition for Simulation of Flow Over Porous Surfaces

A new boundary condition is presented.for simulating the flow over passively porous surfaces. The model builds on the prior work of R.H. Bush to eliminate the need for constructing grid within an underlying plenum, thereby simplifying the numerical modeling of passively porous flow control systems and reducing computation cost. Code experts.for two structured-grid.flow solvers, TLNS3D and CFL3D. and one unstructured solver, USM3Dns, collaborated with an experimental porosity expert to develop the model and implement it into their respective codes. Results presented,for the three codes on a slender forebody with circumferential porosity and a wing with leading-edge porosity demonstrate a good agreement with experimental data and a remarkable ability to predict the aggregate aerodynamic effects of surface porosity with a simple boundary condition.

Frink, Neal T.↗

Multiyear Simulations of the Martian Water Cycle with the Ames General Circulation Model

Mars atmosphere is carbon dioxide dominated with non-negligible amounts of water vapor and suspended dust particles. The atmospheric dust plays an important role in the heating and cooling of the planet through absorption and emission of radiation. Small dust particles can potentially be carried to great altitudes and affect the temperatures there. Water vapor condensing onto the dust grains can affect the radiative properties of both, as well as their vertical extent. The condensation of water onto a dust grain will change the grain s fall speed and diminish the possibility of dust obtaining high altitudes. In this capacity, water becomes a controlling agent with regard to the vertical distribution of dust. Similarly, the atmosphere s water vapor holding capacity is affected by the amount of dust in the atmosphere. Dust is an excellent green house catalyst; it raises the temperature of the atmosphere, and thus, its water vapor holding capacity. There is, therefore, a potentially significant interplay between the Martian dust and water cycles. Previous research done using global, 3-D computer modeling to better understand the Martian atmosphere treat the dust and the water cycles as two separate and independent processes. The existing Ames numerical model will be employed to simulate the relationship between the Martian dust and water cycles by actually coupling the two cycles. Water will condense onto the dust, allowing the particle's radiative characteristics, fall speeds, and as a result, their vertical distribution to change. Data obtained from the Viking, Mars Pathfinder, and especially the Mars Global Surveyor missions will be used to determine the accuracy of the model results.

R M Haberle↗

Representing Functions in n Dimensions to Arbitrary Accuracy

A method of approximating a scalar function of n independent variables (where n is a positive integer) to arbitrary accuracy has been developed. This method is expected to be attractive for use in engineering computations in which it is necessary to link global models with local ones or in which it is necessary to interpolate noiseless tabular data that have been computed from analytic functions or numerical models in n-dimensional spaces of design parameters.

Scotti, Stephen J.↗

Addition of Electrostatic Forces to EDEM with Applications to Triboelectrically Charged Particles

Tribocharging of particles is common in many processes including fine powder handling and mixing, printer toner transport and dust extraction. In a lunar environment with its high vacuum and lack of water, electrostatic forces are an important factor to consider when designing and operating equipment. Dust mitigation and management is critical to safe and predictable performance of people and equipment. The extreme nature of lunar conditions makes it difficult and costly to carryout experiments on earth which are necessary to better understand how particles gather and transfer charge between each other and with equipment surfaces. DEM (Discrete Element Modeling) provides an excellent virtual laboratory for studying tribocharging of particles as well as for design of devices for dust mitigation and for other purposes related to handling and processing of lunar regolith. Theoretical and experimental work has been performed pursuant to incorporating screened Coulombic electrostatic forces into EDEM Tm, a commercial DEM software package. The DEM software is used to model the trajectories of large numbers of particles for industrial particulate handling and processing applications and can be coupled with other solvers and numerical models to calculate particle interaction with surrounding media and force fields. In this paper we will present overview of the theoretical calculations and experimental data and their comparison to the results of the DEM simulations. We will also discuss current plans to revise the DEM software with advanced electrodynamic and mechanical algorithms.

Hogue, Michael D.↗

Discrete Element Modeling of Triboelectrically Charged Particles

Tribocharging of particles is common in many processes including fine powder handling and mixing, printer toner transport and dust extraction. In a lunar environment with its high vacuum and lack of water, electrostatic forces are an important factor to consider when designing and operating equipment. Dust mitigation and management is critical to safe and predictable performance of people and equipment. The extreme nature of lunar conditions makes it difficult and costly to carry out experiments on earth which are necessary to better understand how particles gather and transfer charge between each other and with equipment surfaces. DEM (Discrete Element Modeling) provides an excellent virtual laboratory for studying tribocharging of particles as well as for design of devices for dust mitigation and for other purposes related to handling and processing of lunar regolith. Theoretical and experimental work has been performed pursuant to incorporating screened Coulombic electrostatic forces into EDEM, a commercial DEM software package. The DEM software is used to model the trajectories of large numbers of particles for industrial particulate handling and processing applications and can be coupled with other solvers and numerical models to calculate particle interaction with surrounding media and force fields. While simple Coulombic force between two particles is well understood, its operation in an ensemble of particles is more complex. When the tribocharging of particles and surfaces due to frictional contact is also considered, it is necessary to consider longer range of interaction of particles in response to electrostatic charging. The standard DEM algorithm accounts for particle mechanical properties and inertia as a function of particle shape and mass. If fluid drag is neglected, then particle dynamics are governed by contact between particles, between particles and equipment surfaces and gravity forces. Consideration of particle charge and any tribocharging and electric field effects requires calculation of the forces due to these effects.

Hogue, Michael D.↗

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↗

Benchmark of the LAMMPS code on CRESCO and SUMMIT HPC systems

In this work we present the benchmark results of a numerical model where the solid Deuterium grows under gas pressure. The two-phase molecular Deuterium systems (gas + solid) is modeled using classical molecular dynamics approach and the LAMMPS code. The performance between different modern supercomputers are analyzed and compared. In addition, the numerical simulations reach high accuracy to produce results that match with the data collected through laboratory experiments.

Lupo Pasini, Massimiliano↗

Spatiotemporal parallelization of an analytical heat conduction model for additive manufacturing via a hybrid OpenMP + MPI approach

The ability to do thermal simulations for entire additive manufacturing builds is a key computational problem facing the additive manufacturing community; however, complex numerical models considering multiple physical phenomena currently do not have the capacity for simulations at this scale. To this end, conduction only analytic models offer a viable approach due to the massive drop in computational expense. In this work, we extend an existing implementation which uses a governing equation which can be evaluated at any point in space and time. This implementation already utilizes OpenMP with a spatial decompositions scheme stemming from a melt pool tracking algorithm. Furthermore, we then combine this with a parallel in time (PinT) approach to make the problem highly parallelizable. The new scheme, which uses MPI for internode communication and OpenMP for intranode communication, is shown to scale very well across multiple computational nodes. This approach results in the ability to simulate the 3D solidification conditions for entire layers of additively manufactured parts in minutes making part scale thermal simulations more practical.

36 MATERIALS SCIENCE↗

Where, when, and why: Occurrence of fast-pairwise collective neutrino oscillation in three-dimensional core-collapse supernova models

Fast-pairwise collective neutrino oscillation represents a key uncertainty in the theory of core-collapse supernova (CCSN). Despite the potentially significant impact on CCSN dynamics, it is usually neglected in numerical models of CCSN because of the formidable technical difficulties of self-consistently incorporating this physics. In this paper, we investigate the prospects for the occurrence of fast flavor conversion by diagnosing electron neutrino lepton number (ELN) crossing in more than a dozen state-of-the-art three-dimensional CCSN models. ELN crossings provide a necessary condition for triggering flavor conversion. Although only zeroth and first angular moments are available from the simulations, our new method enables us to look into the angular distributions of neutrinos in momentum space and provide accurate insight into ELN crossings. Our analysis suggests that fast flavor conversion generally occurs in the postshock region of CCSNe, and that explosive models provide more favorable conditions for the flavor conversion than failed CCSNe. Here we also find that there are both common and progenitor-dependent characteristics. Classifying ELN crossings into two types, we analyze the generation mechanism of each case by scrutinizing the neutrino radiation field and matter interactions. We find key ingredients of CCSN dynamics driving the ELN crossings: proto-neutron star convection, asymmetric neutrino emission, neutrino absorptions, and scatterings. This study suggests that we need to accommodate fast flavor conversions in realistic CCSN models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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.↗