Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation analytics”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

The evolution of cocoons surrounding light, extragalactic jets

If the mass density of supersonic, collimated material is less than that of the surrounding medium, a so-called light jet will be enveloped by a cocoon of overpressured shocked gas. Hydrodynamical simulations are used to understand the evolution of the cocoon. The cocoon's evolution is also compared to a simple analytic theory. To reconcile the theory with the simulations, the growth of the jet head must be taken into account. The overpressured cocoon stage exists for a relatively short astronomical time, after which only the region of the cocoon near the jet head remains overpressured. The spatial distribution of the optical emission often observed in distant extragalactic jet systems can be explained with this improved understanding of cocoon evolution.

Cioffi, Denis F.↗
DOE OSTI2025

Sedimentation and shear-induced dynamics of spheroids in fluids with spatial viscosity variations

A generalized reciprocal theorem is used to relate the force and torque induced on a particle in an inertia-less fluid with small variation in viscosity to integrals involving Stokes flow fields and the spatial dependence of viscosity. These resistivity expressions are analytically evaluated using spheroidal harmonics and then used to obtain the mobility of the spheroid during sedimentation, and in linear flows, of a fluid with linear viscosity stratification. The coupling between the rotational and translational motion induced by stratification rotates the spheroid’s centerline, creating a variety of rotational and translational dynamics dependent upon the particle’s aspect ratio, κ, and the component of the stratification unit vector in the gravity direction, d g . Spheroids with 0.55 ⪅ κ ⪅ 2.0 exhibit the largest variety of settling behaviors. Interestingly, this range covers most microplastics and typical microorganisms. One of the modes include a stable orientation dependent only on κ and d g , but independent of initial orientation, thus allowing for the potential control of settling angles and sedimentation rates. In a simple shear flow, cross-streamline migration occurs due to the stratification-induced force generated on the particle. Similarly, a particle no longer stays at the stagnation point of a uniaxial extensional flow. While fully analytical results are obtained for spheroids, numerical simulations provide a source of validation. These simulations also provide additional insights into the stratification-induced force- and torque-producing mechanisms through the stratification-induced stress, which is not accessed in the reciprocal theorem-based analytical calculations.

Geophysical and Geological Flows: Stratified flows↗
DOE OSTI2023

tinyIFD: A High-Throughput Binding Pose Refinement Workflow Through Induced-Fit Ligand Docking

A critical step in structure-based drug discovery is predicting whether and how a candidate molecule binds to a model of a therapeutic target. However, substantial protein side chain movements prevent current screening methods, such as docking, from accurately predicting the ligand conformations and require expensive refinements to produce viable candidates. Here, we present the development of a high-throughput and flexible ligand pose refinement workflow, called “tinyIFD”. The main features of the workflow include the use of specialized high-throughput, small-system MD simulation code mdgx.cuda and an actively learning model zoo approach. We show the application of this workflow on a large test set of diverse protein targets, achieving 66% and 76% success rates for finding a crystal-like pose within the top-2 and top-5 poses, respectively. We also applied this workflow to the SARS-CoV-2 main protease (M pro ) inhibitors, where we demonstrate the benefit of the active learning aspect in this workflow.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

High-performance multiblock PEMs containing a highly acidic fluorinated-hydrophilic domain for water electrolysis

The present paper describes the design and evaluation of novel hydrophilic–hydrophobic poly(arylene ether sulfone) (PAES) multiblock copolymers for their synergistic effects upon transport properties and their potential use in proton exchange membrane water electrolysis. The multiblock copolymers are prepared via a coupling reaction between (i) a hydrophilic segment consisting of a disulfonated quinone fluorinated biphenyl group that contains fluorine moieties next to the sulfonated groups to increase the acidity, and (ii) hydrophobic segments composed of non-sulfonated biphenyl sulfone to provide dimensional stability. Two different lengths (molecular weights; 5k and 10k, where k represents 10 3 g mol -1 ) of hydrophobic segments are used to investigate the effects of the membrane properties compared with those of Nafion® and PAES random copolymer (i.e., BPSH40). Atomic force microscopy images of the BPSH40 and multiblock membranes are shown to agree closely with a mesoscale simulation, thus confirming the importance of the morphological effect upon the transport properties. Furthermore, the multiblock copolymer with a higher proportion of hydrophilic segments (10 k–5k) was shown to provide enhanced performance (3.41 A cm -2 at 1.9 V) compared to the multiblock copolymer with equal proportions of hydrophilic and hydrophobic segments (10 k–10 k) due to the greater continuity of nano-sized ionic channels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Multiscale water drop contact angles at selected silica surfaces

In this paper, multiscale advancing contact angles for glycerol/water drops at silica surfaces are reported for millidrops, submicron-drops, and nanodrops. Selected silica surfaces were muscovite, silicon, and talc. The contact angles for millidrops (1–2 mm) were determined by the traditional sessile drop technique. For submicron-drops (0.1–1.0 μm), a hollow tip Atomic Force Microscope (AFM) procedure was used. The contact angles for nanodrops (~7 nm) were examined from Molecular Dynamics (MD) simulation. The results were compared to evaluate the effect of drop size on the contact angle. In the case of the hydrophobic talc surface, the 75° advancing contact angle did not vary significantly with drop size. For the hydrophilic muscovite surface, the water drop wet the surface and an advancing contact angle of about 10° was found for the millidrops and submicron-drops. However, for the MD simulated nanodrops, attachment and spreading of the ~7 nm drop created a 2D film of molecular dimensions, the contact angle of which was difficult to define and varied from 0° to 17°. Perhaps of equal interest from the MD simulation results was that the spreading of the glycerol/water nanodrop at the muscovite surface resulted in crystallographic directional transport of water molecules to the extremities of the 2D film. Such separation and segregation left the center of the film with an increased concentration of glycerol. Based on these results, the line tension, which has been found in other investigations to account for contact angle decrease with a decrease in drop size, does not seem to be a significant factor in this study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Reflection and transmission of electromagnetic pulses at a planar dielectric interface -- theory and quantum lattice simulations

There is considerable interest in the application of quantum information science to advance computations in plasma physics. A particular point of curiosity is whether it is possible to take advantage of quantum computers to speed up numerical simulations relative to conventional computers. Many of the topics in fusion plasma physics are classical in nature. In order to implement them on quantum computers it will require couching a classical problem in the language of quantum mechanics. Electromagnetic waves are routinely used in fusion experiments to heat a plasma or to generate currents in the plasma. The propagation of electromagnetic waves is described by Maxwell equations with an appropriate description of the plasma as a dielectric medium. Before advancing to the tensor dielectric of a magnetized plasma, this paper considers electromagnetic wave propagation in a one-dimensional inhomogeneous scalar dielectric. The classic theory of scattering of plane electromagnetic waves at a planar interface, separating two different dielectric media, leads to Fresnel equations for reflection and transmission coefficients. In contrast to plane waves, this paper is on the reflection and transmission of a spatially confined electromagnetic pulse. Following an analytical formulation for the scattering of a Gaussian pulse, it is deduced that the maximum transmission coefficient for a pulse is $\sqrt{n_2/n_1}$ times that for a plane wave; the incident and transmitted pulses propagating in dielectric media with refractive indices $n_1$ and $n_2$, respectively. The analytical theory is complemented by numerical simulations using a quantum lattice algorithm for Maxwell equations. The algorithm, based on the Riemann-Silberstein-Weber representation of the electromagnetic fields and expressed in term of qubits, is an interleaved sequence of entangling operators at each lattice site and unitary streaming operators which transmit information from one site to an adjacent lattice site. Besides substantiating results from the theory for Gaussian pulses, numerical simulations show their validity for non-Gaussian pulses. Apart from their time-asymptotic forms, the simulations display an interplay between the incident, reflected, and transmitted pulses in the vicinity of the transition region between two dielectric media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

Ion Dynamics Model for Collisionless Radio Frequency Sheaths

Full scale reactor model based on fluid equations is widely used to analyze high density plasma reactors. It is well known that the submillimeter scale sheath in front of a biased electrode supporting the wafer is difficult to resolve in numerical simulations, and the common practice is to use results for electric field from some form of analytical sheath model as boundary conditions for full scale reactor simulation. There are several sheath models in the literature ranging from Child's law to a recent unified sheath model [P. A. Miller and M. E. Riley, J. Appl. Phys. 82, 3689 (1997)l. In the present work, the cold ion fluid equations in the radio frequency sheath are solved numerically to show that the spatiotemporal variation of ion flux inside the sheath, commonly ignored in analytical models, is important in determining the electric field and ion energy at the electrode. Consequently, a semianalytical model that includes the spatiotemporal variation of ion flux is developed for use as boundary condition in reactor simulations. This semianalytical model is shown to yield results for sheath properties in close agreement with numerical solutions.

Bose, Deepak↗
DOE OSTI2020

Predicting chemical reaction equilibria in molten carbonate fuel cells via molecular simulations

It has been recently suggested that hydroxide ions can be formed in the electrolyte of molten carbonate fuel cells when water vapor is present. The hydroxide can replace carbonate in transporting electrons across the electrolyte, thereby reducing the CO 2 separation efficiency of the fuel cell although still producing electricity. In this work, we obtain the equilibrium concentration of hydroxide in five molten alkali carbonate salts from molecular simulations. The results reveal that there can be a substantial amount of hydroxide in the electrolyte at low partial pressures of CO 2 . In addition, we find that the equilibrium concentration of molecular water dissolved in the electrolyte is over two orders of magnitude higher than that of CO 2 . Increasing the size and polarizability (or in other words reducing the “hardness”) of the cations present in the electrolyte can reduce the hydroxide fraction, but at the cost of lowering ionic conductivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

MARTINI-Compatible Coarse-Grained Model for the Mesoscale Simulation of Peptoids

Peptoids (poly-N-substituted glycines) are a class of synthetic polymers that are regioisomers of peptides (poly-C-substituted glycines), in which the point of side-chain connectivity is shifted from the backbone C to the N atom. Peptoids have found diverse applications as peptidomimetic drugs, protein mimetic polymers, surfactants, and catalysts. Computational modeling is valuable in the understanding and design of peptoid-based nanomaterials. In this work, we report the bottom-up parameterization of coarse-grained peptoid force fields based on the MARTINI peptide force field against all-atom peptoid simulation data. Our parameterization pipeline iteratively refits coarse-grained bonded interactions using iterative Boltzmann inversion and nonbonded interactions by matching the potential of mean force for chain extension. We assure good sampling of the amide bond cis/trans isomerizations in the all-atom simulation data using parallel bias metadynamics. We develop coarse-grained models for two representative peptoids—polysarcosine (poly(N-methyl glycine)) and poly(N-((4-bromophenyl)ethyl)glycine)—and show their structural and thermodynamic properties to be in excellent accord with all-atom calculations but up to 25-fold more efficient and compatible with MARTINI force fields. Here, this work establishes a new rigorously parameterized coarse-grained peptoid force field for the understanding and design of peptoid nanomaterials at length and time scales inaccessible to all-atom calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Locating fast-varying line disturbances with the frequency mismatch

In an attempt to provide an efficient method for line disturbance identification in complex networks of diffusively coupled agents, we recently proposed to leverage the frequency mismatch. The frequency mismatch filters out the intricate combination of interactions induced by the network structure and quantifies to what extent the trajectory of each agent is affected by the disturbance. In this previous work, we provided some analytical evidence of its efficiency when the perturbation is assumed to be slow. Here we claim that the frequency mismatch performs actually well for most disturbance regimes. This is shown through a series of simulations and is backed up by an analytical argument. Therefore, we argue that the frequency mismatch is an efficient and elegant tool for line disturbance location in complex networks of diffusively coupled agents.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Shock state distributions in porous tantalum and characterization with multipoint velocimetry

Heterogenous materials under shock compression can be expected to reach different shock states throughout the material according to local differences in microstructure and the history of wave propagation. Here, in this work, a compact, multiple-beam focusing optic assembly is used with high-speed velocimetry to interrogate the shock response of porous tantalum films prepared through thermal-spray deposition. The distribution of particle velocities across a shocked interface is compared to results obtained using a set of defocused interferometric beams that sampled the shock response over larger areas. The two methods produced velocity distributions along the shock plateau with the same mean, while a larger variance was measured with narrower beams. The finding was replicated using three-dimensional, mesoscopically resolved hydrodynamics simulations of solid tantalum with a pore structure mimicking statistical attributes of the material and accounting for radial divergence of the beams, with agreement across several impact velocities. Accounting for pore morphology in the simulations was found to be necessary for replicating the rise time of the shock plateau. The validated simulations were then used to show that while the average velocity along the shock plateau could be determined accurately with only a few interferometric beams, accurately determining the width of the velocity distribution, which here was approximately Gaussian, required a beam dimension much smaller than the spatial correlation lengthscale of the velocity field, here by a factor of ∼30×, with implications for the study of other porous materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2005

The Influence of Clocking Angle of the Projectile on the Simulated Impact Response of a Shuttle Leading Edge Wing Panel

An analytical study was conducted to determine the influence of clocking angle of a foam projectile impacting a space shuttle leading edge wing panel. Four simulations were performed using LS-DYNA. The leading edge panels are fabricated of multiple layers of reinforced carbon-carbon (RCC) material. The RCC material was represented using Mat 58, which is a material property that can be used for laminated composite fabrics. Simulations were performed of a rectangular-shaped foam block, weighing 0.23-lb., impacting RCC Panel 9 on the top surface. The material properties of the foam were input using Mat 83. The impact velocity was 1,000 ft/s along the Orbiter X-axis. In two models, the foam impacted on a corner, in one model the foam impacted the panel initially on the 2-in.-long edge, and in the last model the foam impacted the panel on the 7-in.- long edge. The simulation results are presented as contour plots of first principal infinitesimal strain and time history plots of contact force and internal and kinetic energy of the foam and RCC panel.

Jackson, Karen E.↗
NASA NTRS1978

Control of wind turbine generators connected to power systems

A unique simulation model based on a Mode-O wind turbine is developed for simulating both speed and power control. An analytical representation for a wind turbine that employs blade pitch angle feedback control is presented, and a mathematical model is formulated. For Mode-O serving as a practical case study, results of a computer simulation of the model as applied to the problems of synchronization and dynamic stability are provided. It is shown that the speed and output of a wind turbine can be satisfactorily controlled within reasonable limits by employing the existing blade pitch control system under specified conditions. For power control, an additional excitation control is required so that the terminal voltage, output power factor, and armature current can be held within narrow limits. As a result, the variation of torque angle is limited even if speed control is not implemented simultaneously with power control. Design features of the ERDA/NASA 100-kW Mode-O wind turbine are included.

Hwang, H. H.↗
DOE OSTI2020

Zeolite-templated carbons as effective sorbents to remove methylsiloxanes and derivatives: A computational screening

Though widely used in our daily lives, volatile methylsiloxanes and derivatives are emerging contaminants and becoming a high-priority environment and public health concern. Developing effective sorbent materials can remove siloxanes in a cost-effective manner. Herein, by means of Grand Canonical Monte Carlo (GCMC) simulations, we evaluated the potentials of the recently proposed 68 stable zeolite-templated carbons (ZTCs) for the removal of four linear methylsiloxanes and derivatives as well as two cyclic methylsiloxanes by the calculated average loading and average adsorption energy values. Four ZTCs, namely ISV, FAU1, FAU3, and H8326836, were identified with the top 50% adsorption performance toward all the six targeted contaminants, which outperform activated carbons. Further first principles computations revealed that steric hindrance, electrostatic interactions (further enhanced by charge transfer), and CH-π interactions account for the outstanding adsorption performance of these ZTCs. This work provides a quick procedure to computationally screen promising ZTCs for siloxane removal, and help guide future experimental and theoretical investigations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Molecular-based analysis of nanoparticle solvation: Classical density functional approach

Proper statistical mechanics understanding of nanoparticle solvation processes requires an accurate description of the molecular structure of the solvent. Achieving this goal with standard molecular dynamics (MD) simulation methods is challenging due to large length scales. An alternative approach to this problem can be formulated using classical density functional theory (cDFT), where a full configurational description of the positions of all the atoms is replaced by collective atomic site densities in the molecule. Using an example of the negatively charged silica-like system in an aqueous polar environment represented by a two-site water model, we demonstrate here that cDFT can reproduce MD data at a fraction of the computational cost. An important implication of this result is the ability to understand how the solvent molecular features may affect the system’s properties at the macroscopic scale. A concrete example highlighted in this work is the analysis of nanoparticle interactions with sizes of up to 100 nm in diameter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1991

Small body mass determination from a flyby

This paper investigates the accuracy with which the mass of a small body can be determined from a flyby. Due to the nearly linear nature of hyperbolic flybys of small bodies, analytic solutions can be derived. Four analytic solutions are developed, each of increasing complexity with different applications. These solutions are compared to a complex computer simulation for accuracy. When these methods are applied to nearly linear trajectories, the analytic solutions agree with the computer simulation results to within 4.0 percent except in one boundary case which is discussed in the text. Also described is the linear range for which these equations are applicable. These solutions are then applied to the proposed asteroid flyby of the CRAF mission to predict the mass estimate accuracy.

Sybert, Catherine B.↗
DOE OSTI2020

Electrolytically Assisted Surface Decontamination (EASD{sup TM}) for POCO Operations - 20282

NNL in collaboration with C-Tech Innovation Ltd and Sellafield Ltd has been exploring innovative technologies to enable a significant reduction of radiological hazards within facilities during the Post Operational Clean Out (POCO) phase of a nuclear plant's life cycle. Reducing the hazard by effectively decontaminating plants in-situ during POCO delivers huge cost reductions for future decommissioning operations. These cost reduction benefits are realized by reducing the number and complexity of remote operations as well as lowering the long-term waste disposal costs. Whilst chemical decontamination can achieve the desired level of decontamination, applying aggressive chemical reagents is hazardous, potentially difficult to control and requires there to be complimentary effluent treatment and waste routes. This work has aimed to develop flexible and controllable decontamination processes which could be operated without the additional complexity and issues associated with chemical decontamination. The processes needed to be relatively fast and effective to minimize the time operators would spend in an active area. In addition, there was a driver to produce decontamination methodologies which generated a secondary waste compatible with current routes and which takes advantage of current waste capacity. The strategy therefore was to come up with a solution that works with, rather than against, the science and engineering behind process plants at Sellafield. Electrolytically Assisted Surface Decontamination (EASD{sup TM}) is an innovative electrochemical decontamination process (developed by NNL and C-Tech Innovation Ltd) that can remove activity from contaminated metal in very short time periods. Application of the patented electric waveform to the surface when treating contaminated metal with nitric acid has been shown to significantly enhance the decontamination performance. When compared to proposed baseline washout procedures at Sellafield, the only change to the process is the applied electrical waveform which is temporary and controllable. The current causes dissolution of the metal surface being decontaminated which leads to activity transferring from the plant item into the nitric acid effluent stream. This innovation has the potential to transform the POCO process and allow a nitric acid based washout approach to deliver POCO quickly and cost effectively within the existing infrastructure. Several devices have now been designed that incorporate this EASD{sup TM} technology and enable the decontamination process to be applied to a range items commonly identified as being contaminated during nuclear decommissioning programmes such as pipework, tanks and hotspots of walls/floors. Inactive and active laboratory-scale testing has been completed using both simulated and 'real-life' contaminated (low level waste) metallic items retrieved from different nuclear sites. Active trials demonstrated that contaminated items could be treated to free release levels within minutes. Engineering scale tests are currently being performed, with guidance from Sellafield's system engineers and plant managers, to provide the necessary re-assurance the designed devices could be successfully deployed in an on-plant scenario. The final stage of this development work is to carry out an active demonstration of the decontamination technology on Thorp at Sellafield. This paper aims to highlight progress made to date with the EASD{sup TM} technology and, more specifically, the development of an in-situ decontamination device for deployment within radioactive pipework. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Experimental and numerical studies of seismic fluid-structure interaction in a base-supported cylindrical vessel

Seismic design and qualification of advanced reactors will rely heavily on the use of verified and validated numerical models capable of capturing the interaction of the vessel, its contained fluid, and the internal equipment: fluid-structure interaction (FSI) analysis. Analytical solutions can be used for preliminary sizing and design of such vessels but their application is limited to simple geometries and boundary conditions, and small amplitude, translational (and rotational) inputs. To validate numerical models for seismic FSI analysis in finite element codes, a comprehensive set of experiments was performed on a liquid-filled cylindrical vessel, using a 6 degree-of-freedom earthquake simulator. Results in terms of sloshing frequency, damping ratio in sloshing modes, and hydrodynamic responses (wave height, hydrodynamic pressure, base shear, and base moment) for multi-directional earthquake simulator inputs are reported and compared with analytical solutions for liquid-filled vessels. The impact of seismic (base) isolation on hydrodynamic responses was studied using earthquake simulator inputs generated using a virtual isolation system. Data from the experiments are used to validate a numerical model of the fluid-filled vessel using the Arbitrary Lagrangian Eulerian (ALE) solver in the commercial finite element program LS-DYNA. Validation studies are presented for multi-directional seismic inputs, including rocking motions. Lagrangian modeling approaches using an elastic material formulation for the fluid are also investigated and their limitations and possible applications are identified. Here, the results are broadly applicable to the seismic response of base supported, liquid-filled vessels.

58 GEOSCIENCES↗