Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic models”

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 667 records · Page 37

Design, dynamic modeling and wave basin verification of a Hybrid Wave–Current Energy Converter

Massive and high-density Marine and Hydrokinetic (MHK) energy is contained in the ocean, including waves, tidal streams, and ocean currents. Traditional MHK energy converters harvest energy from only a single MHK energy source, which does not fully exploit the energy potential that co-exists in multiple forms of MHK energy in the ocean. This paper presents the design and dynamics of a Hybrid Wave–Current Energy Converter (HWCEC) that can simultaneously convert both wave and current energy to electricity with a single Power Take-off (PTO) through the engagement and disengagement of the three one-way clutches. The critical design parameters are analyzed through modeling and simulation, including the transmission ratios, the electrical impedance, and the turbine-heave plate distance. Water basin tests in a wave–current tank were conducted, which shows the prototyped HWCEC can improve the electric power output by 38%–71% for regular waves and 79% for irregular waves, while providing a 70% reduction of the Peak to Average Ratio (PAR) of power compared to the baseline Wave Energy Converter (WEC). Here we further analyze the occurring percentages of the instant power of both WEC and HWCEC in irregular wave tests, shedding new insight on how HWCEC reduces the PAR.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A reactive molecular dynamics model for uranium/hydrogen containing systems

Uranium-based materials are valuable assets in the energy, medical, and military industries. However, understanding their sensitivity to hydrogen embrittlement is particularly challenging due to the toxicity of uranium and the computationally expensive nature of quantum-based methods generally required to study such processes. In this regard, we have developed a Chebyshev Interaction Model for Efficient Simulation (ChIMES) that can be employed to compute energies and forces of U and UH3 bulk structures with vacancies and hydrogen interstitials with accuracy similar to that of Density Functional Theory (DFT) while yielding linear scaling and orders of magnitude improvement in computational efficiency. Here, we show that the bulk structural parameters, uranium and hydrogen vacancy formation energies, and diffusion barriers predicted by the ChIMES potential are in strong agreement with the reference DFT data. We then use ChIMES to conduct molecular dynamics simulations of the temperature-dependent diffusion of a hydrogen interstitial and determine the corresponding diffusion activation energy. Our model has particular significance in studies of actinides and other high-Z materials, where there is a strong need for computationally efficient methods to bridge length and time scales between experiments and quantum theory.

36 MATERIALS SCIENCE↗

Dynamic modeling of spacecraft in a collisionless plasma

A new computational model is described which can simulate the charging of complex geometrical objects in three dimensions. Two sample calculations are presented. In the first problem, the capacitance to infinity of a complex object similar to a satellite with solar array paddles is calculated. The second problem concerns the dynamical charging of a conducting cube partially covered with a thin dielectric film. In this calculation, the photoemission results in differential charging of the object.

Katz, I.↗

Understanding methane/carbon dioxide partitioning in clay nano- and meso-pores with constant reservoir composition molecular dynamics modeling

The interactions among fluid species such as H 2 O, CO 2 , and CH 4 confined in nano- and meso-pores in shales and other rocks is of central concern to understanding the chemical behavior and transport properties of these species in the earth's subsurface and is of special concern to geological C-sequestration and enhanced production of oil and natural gas. The behavior of CO 2 , and CH 4 is less well understood than that of H 2 O. This paper presents the results of a computational modeling study of the partitioning of CO 2 and CH 4 between bulk fluid and nano- and meso-pores bounded by the common clay mineral montmorillonite. The calculations were done at 323 K and a total fluid pressure of 124 bars using a novel approach (constant reservoir composition molecular dynamics, CRC-MD) that uses bias forces to maintain a constant composition in the fluid external to the pore. This purely MD approach overcomes the difficulties in making stochastic particle insertion–deletion moves in dense fluids encountered in grand canonical Monte Carlo and related hybrid approaches. The results show that both the basal siloxane surfaces and protonated broken edge surfaces of montmorillonite both prefer CO 2 relative to CH 4 suggesting that methods of enhanced oil and gas production using CO 2 will readily displace CH 4 from such pores. This preference for CO 2 is due to its preferred interaction with the surfaces and extends to approximately 20 Å from them.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic Modeling using Output-Error Parameter Estimation based on Frequency Responses Estimated with Multisine Inputs

A method is developed for estimating model parameters, such as nondimensional stability and control derivatives, by fitting transfer function or state-space models to empirical frequency response data using the output-error approach. The frequency response data were computed using Fourier transforms of measured input and output data. The control surfaces were excited with periodic multisine inputs which facilitated time-ecient estimation of multiple-input multiple-output frequency responses. The method was applied to lateral data from a nonlinear flight dynamics simulation of the F-16 aircraft, and to longitudinal data from multiple repeated flight test maneuvers of the NASA T-2 subscale aircraft. Results using simulation data showed the frequency response method compared well to other standard methods for parameter estimation. In addition to including all the available inputs, outputs, and harmonic frequencies in the estimation, relatively small subsets of the measured data could also be used to focus on identifying specific parts of the model. Results from flight test data showed that parameter estimates and uncertainties determined from repeated maneuvers were accurate and in statistical agreement with each other.

Grauer, Jared A.↗

Dynamic Modeling of Sequential Service Restoration in Islanded Single Master Microgrids

This paper presents a methodology for generating restoration sequences in distribution systems and microgrids with emphasis on deriving optimal switching times for remote control switches (RCS). Microgrid architecture considered are assumed to be operating in islanded single master mode. Control methods for handling the unbalance inherent in distribution systems using single-phase controlled PQ inverters were developed and incorporated in the proposed method. The problem of considering simultaneously both energy optimization and dynamic stability of the system, which is lacking in existing methods, was addressed by studying the transient behavior of microgrids based on extended EMTP simulations and then approximate and incorporate these behaviors into the optimization formulation for restoration. The approach is based on multi-objective optimization of final restored load and total restoration time. The methodology was studied on a modified IEEE 123 node test feeder. An EMTP simulation in PSCAD of the formed microgrids was used to validate the accuracy of the dynamic and steady-state characteristics of the resulting systems. By considering the transient behavior of DGs during restoration, a better restoration sequence that does not trigger protection relays could be realized.

42 ENGINEERING↗

Molecular Dynamics Modeling of the Structure and Na + -Ion Transport in Na 2 S + SiS 2 Glassy Electrolytes

Solid-state sodium batteries, a relatively safe and potentially cost-effective energy-storage technology, have attracted increasing scientific attention recently for application in stationary grid-scale energy storage. Identifying solid electrolytes with high electrochemical stability and high Na+-ion conductivity at room temperature is critically important to enable high energy densities with enhanced rate capabilities. We evaluated sodium sulfide silicon sulfide, xNa 2 S + (1 - x)SiS 2 , glasses as potential glassy solid electrolytes (GSEs) using molecular dynamics (MD) simulations. We employed ab initio MD to determine ion conduction mechanisms, to calculate energy barriers for ion hops, and to correlate these to the local short-range structure of 0.50Na 2 S + 0.50SiS 2 glass. To simulate much larger systems for accurately calculating the ionic conductivity, we parameterized empirical Buckingham-type potential and performed classical MD simulations. After validating these calculations by comparing the structure obtained from MD to that from X-ray scattering data, we calculated the ionic conductivity of these glasses for the range of 0.33 ≤ x ≤ 0.67 compositions. The calculated ionic conductivities at room temperature were in the range of similar to 10 -5 S/cm for the x = 0.50 composition and increased significantly with sodium sulfide (x) content. Furthermore, these calculations provide theoretical insights into the role of Na 2 S content on the ionic conductivity of GSEs aiding in the selection of specific compositions to enhance the ionic conductivity.

25 ENERGY STORAGE↗

In situ X-ray Scattering and Dynamical Modeling of Pd Catalyst Nanoparticles Formed in Flames

It has previously been demonstrated that organopalladium precursors can break down under combustion conditions, forming nanoparticles that catalyze ignition. Here, we use in situ small-angle X-ray scattering (SAXS) to probe the formation and growth of palladium nanoparticles in an ethylene fl ame doped with 28 ppm (mol) of Pd(acetate) 2 . The particles appear to nucleate in the fl ame front and are observed by SAXS to grow in size and mass in the high-temperature region of the fl ame ( ~ 2200 K) with median diameters that evolve from 1.5 to 3.0 nm. Transmission electron microscopy of particles collected on a grid located outside the fl ame shows that the particles are metallic palladium with sizes comparable to those determined by SAXS. Molecular dynamics simulation of particles of selected sizes indicates that at the fl ame temperature the particles are molten and the average mass density of the particle material is notably smaller than that of bulk, liquid Pd at the melting point. Both experimental and computational results point to homogeneous nucleation and particle - particle coalescence as mechanisms for particle formation and growth. Aerosol dynamics simulation reproduces the time evolution of the particle size distribution and suggests that a substantial fraction of the particles must be electrically charged during their growth process.

36 MATERIALS SCIENCE↗

X-57 Ground Dynamics Modeling and Analysis

This paper presents simulation work to assess the ground handling behavior of the X-57 airplane, prior to its first flight under electric propulsion. A tire deformation model was fit to manufacturer ground testing data to develop a detailed simulation model of previously unmodeled rolling drag forces on the tires. A special modeling effort corrected for a difference between the runway surface in the manufacturer data (grass) and in intended X-57 surface operations (asphalt/concrete). Analysis based on variation in the test conditions of the manufacturer test data showed that uncertainty in the new estimated rolling drag model should lead to no more than 1.5-percent uncertainty in the takeoff acceleration. Another batch simulation series swept across ground operating conditions and airplane configurations, characterizing turning radii and roll angles experienced for different open-loop steering commands in these various situations. This simulation allowed the identification of ground handling characteristics even beyond the regime of normal operations. As part of this task, a new, modified variant of the Fibonacci search algorithm identified the correct steering inputs for maintaining a straight path as a function of the operating condition. Finally, this paper describes a fixed-base piloted simulator assessment. This assessment supplements the previous simulated ground dynamics data with pilot-in-the-loop data and pilot feedback and confirms the overall acceptability of the ground handling qualities of the X-57 airplane prior to first flight.

Loren J. Newton↗

Molecular dynamics modelling of EGCG clusters on ceramide bilayers

A novel method of atomistic modelling and characterization of both pure ceramide and mixed lipid bilayers is being developed, using only the General Amber ForceField. Lipid bilayers modelled as pure ceramides adopt hexagonal packing after equilibration, and the area per lipid and bilayer thickness are consistent with previously reported theoretical results. Mixed lipid bilayers are modelled as a combination of ceramides, cholesterol, and free fatty acids. This model is shown to be stable after equilibration. Green tea extract, also known as epigallocatechin-3-gallate, is introduced as a spherical cluster on the surface of the mixed lipid bilayer. It is demonstrated that the cluster is able to bind to the bilayers as a cluster without diffusing into the surrounding water.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computational Fluid Dynamics Modeling of Bubbling in a Viscous Fluid for Validation of Waste Glass Melter Modeling

At the Hanford site, radioactive waste stored in underground tanks is slated for vitrification for final disposal. A comprehensive knowledge of the glass batch melting process will be useful in optimizing the process, which could potentially reduce the cost and duration of this multi-billion dollar cleanup effort. We are developing a high-fidelity heat transfer model of a Joule-heated ceramic lined melter to improve the understanding of the complex, inter-related processes occurring with the melter. The glass conversion rates in the cold cap layer are dependent on promoting efficient heat transfer. In practice, heat transfer is augmented by inserting air bubblers into the molten glass. However, the computational simulations must be validated to provide confidence in the solutions. As part of a larger validation procedure, it is beneficial to split the physics of the melter into smaller systems to validate individually. The substitution of molten glass for a simulant liquid with similar density and viscosity at room temperature provides a way to study mixing through bubbling as an isolated effect without considering the heat transfer dynamics. The simulation results are compared to experimental data obtained by the Vitreous State Laboratory at the Catholic University of America using bubblers placed within a large acrylic tank that is similar in scale to a pilot glass waste melter. Comparisons are made for surface area of the rising air bubbles between experiments and CFD simulations for a variety of air flow rates and bubble injection depths. Also, computed bubble rise velocity is compared to a well-accepted expression for bubble terminal velocity.

bubbling↗

Computational Fluid Dynamics Modeling of Bubbling in a Viscous Fluid for Validation of Waste Glass Melter Modeling

At the Hanford site, radioactive waste stored in underground tanks is slated for vitrification for final disposal. A comprehensive knowledge of the glass batch melting process will be useful in optimizing the process, which could potentially reduce the cost and duration of this multi-billion dollar cleanup effort. We are developing a high-fidelity heat transfer model of a Joule-heated ceramic lined melter to improve the understanding of the complex, inter-related processes occurring with the melter. The glass conversion rates in the cold cap layer are dependent on promoting efficient heat transfer. In practice, heat transfer is augmented by inserting air bubblers into the molten glass. However, the computational simulations must be validated to provide confidence in the solutions. As part of a larger validation procedure, it is beneficial to split the physics of the melter into smaller systems to validate individually. The substitution of molten glass for a simulant liquid with similar density and viscosity at room temperature provides a way to study mixing through bubbling as an isolated effect without considering the heat transfer dynamics. The simulation results are compared to experimental data obtained by the Vitreous State Laboratory at the Catholic University of America using bubblers placed within a large acrylic tank that is similar in scale to a pilot glass waste melter. Comparisons are made for surface area of the rising air bubbles between experiments and CFD simulations for a variety of air flow rates and bubble injection depths. Also, computed bubble rise velocity is compared to a well-accepted expression for bubble terminal velocity (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A stochastic-dynamic model for global atmospheric mass field statistics

A model that yields the spatial correlation structure of atmospheric mass field forecast errors was developed. The model is governed by the potential vorticity equation forced by random noise. Expansion in spherical harmonics and correlation function was computed analytically using the expansion coefficients. The finite difference equivalent was solved using a fast Poisson solver and the correlation function was computed using stratified sampling of the individual realization of F(omega) and hence of phi(omega). A higher order equation for gamma was derived and solved directly in finite differences by two successive applications of the fast Poisson solver. The methods were compared for accuracy and efficiency and the third method was chosen as clearly superior. The results agree well with the latitude dependence of observed atmospheric correlation data. The value of the parameter c sub o which gives the best fit to the data is close to the value expected from dynamical considerations.

Ghil, M.↗