Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Implementation”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Implementation of INDO/SCI with COSMO Implicit Solvation and Benchmarking for Solvatochromic Shifts

Accurate and rapid quantum mechanical prediction of solvatochromic shifts, particularly in systems where charge transfer plays a significant role, is important for many aspects of molecular and material design. Although the semiempirical INDO/SCI approach is computationally efficient and performs well for charge-transfer states, the availability of implicit solvent approaches has been limited. As such, we implement the COSMO solvent model with a perturbative state-specific correction to the excited-state energies with the INDO/SCI method. We show that for a series of prototypical π-conjugated molecules, our newly implemented INDO/SCI/COSMO model yields more accurate absorption energies and comparably accurate solvatochromic shifts to those computed using TD-ωB97XD and CIS with COSMO solvation at a substantially lower computational cost.

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

Sensitivity of Carbonaceous Aerosol Properties to the Implementation of a Dynamic Aging Parameterization in the Regional Climate Model RegCM

Freshly emitted soot is hydrophobic, but condensation of secondary aerosols and coagulation with other particles modify its hygroscopic optical properties. This conversion is referred to as “aerosol aging.” Many climate models represent this aging process with a fixed aging time scale, whereas in reality, it is a dynamic process that depends on environmental conditions. Here, we implement a dynamic aging parameterization scheme in the regional climate model RegCM4 in place of the fixed aging timescale of 1.15 days (~27.6 h) and examine its impact on the aerosol life cycle over the Indian subcontinent. The conversion from hydrophobic to hydrophilic aerosol is usually lower than 27.6 h over the entire landmass and lower than 10 h over the polluted Indo-Gangetic Basin (IGB), with seasonal variability. Due to the implementation of the dynamic aging scheme, the column burden and surface mass concentration of carbonaceous aerosols increase during the drier season (December–February) when washout is negligible. The burden is reduced during the wet season (June–September) due to a more efficient washout except over the IGB, where a reduction in precipitation as a result of radiative feedbacks increases the aerosol concentrations. Over the polluted IGB, surface dimming increases due to the dynamic aging scheme, with the top of the atmosphere forcing remaining mostly unchanged. As a result, atmospheric heating increases by at least 1.2 W/m 2 . Our results suggest that climate models should incorporate dynamic aging for a more realistic representation of aerosol simulations, especially in highly polluted regions.

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

Full Implementation of Matrix Approach to Biogeochemistry Module of CLM5

Abstract Earth system models (ESMs) have been rapidly developed in recent decades to advance our understanding of climate change‐carbon cycle feedback. However, those models are massive in coding, require expensive computational resources, and have difficulty in diagnosing their performance. It is highly desirable to develop ESMs with modularity and effective diagnostics. Toward these goals, we implemented a matrix approach to the Community Land Model version 5 (CLM5) to represent carbon and nitrogen cycles. Specifically, we reorganized 18 balance equations each for carbon and nitrogen cycles among the 18 vegetation pools in the original CLM5 into two matrix equations. Similarly, 140 balance equations each for carbon and nitrogen cycles among the 140 soil pools were reorganized into two additional matrix equations. The vegetation carbon and nitrogen matrix equations are connected to soil matrix equations via litterfall. The matrix equations fully reproduce simulations of carbon and nitrogen dynamics by the original model. The computational cost for forwarding simulation of the CLM5 matrix model was 26% more expensive than the original model, largely due to calculation of additional diagnostic variables, but the spin‐up computational cost was significantly saved. We showed a case study on modeled soil carbon storage under two forcing data sets to illustrate the diagnostic capability that the matrix approach uniquely offers to understand simulation results of global carbon and nitrogen dynamics. The successful implementation of the matrix approach to CLM5, one of the most complex land models, demonstrates that most, if not all, the biogeochemical models can be reorganized into the matrix form to gain high modularity, effective diagnostics, and accelerated spin‐up.

58 GEOSCIENCES↗
DOE OSTI2022

Cell-free prototyping enables implementation of optimized reverse β-oxidation pathways in heterotrophic and autotrophic bacteria

Abstract Carbon-negative synthesis of biochemical products has the potential to mitigate global CO 2 emissions. An attractive route to do this is the reverse β-oxidation (r-BOX) pathway coupled to the Wood-Ljungdahl pathway. Here, we optimize and implement r-BOX for the synthesis of C4-C6 acids and alcohols. With a high-throughput in vitro prototyping workflow, we screen 762 unique pathway combinations using cell-free extracts tailored for r-BOX to identify enzyme sets for enhanced product selectivity. Implementation of these pathways into Escherichia coli generates designer strains for the selective production of butanoic acid (4.9 ± 0.1 gL −1 ), as well as hexanoic acid (3.06 ± 0.03 gL −1 ) and 1-hexanol (1.0 ± 0.1 gL −1 ) at the best performance reported to date in this bacterium. We also generate Clostridium autoethanogenum strains able to produce 1-hexanol from syngas, achieving a titer of 0.26 gL −1 in a 1.5 L continuous fermentation. Our strategy enables optimization of r-BOX derived products for biomanufacturing and industrial biotechnology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Performant implementation of the atomic cluster expansion (PACE) and application to copper and silicon

The atomic cluster expansion is a general polynomial expansion of the atomic energy in multi-atom basis functions. Here we implement the atomic cluster expansion in the performant C++ code that is suitable for use in large-scale atomistic simulations. We briefly review the atomic cluster expansion and give detailed expressions for energies and forces as well as efficient algorithms for their evaluation. We demonstrate that the atomic cluster expansion as implemented in shifts a previously established Pareto front for machine learning interatomic potentials toward faster and more accurate calculations. Moreover, general purpose parameterizations are presented for copper and silicon and evaluated in detail. We show that the Cu and Si potentials significantly improve on the best available potentials for highly accurate large-scale atomistic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structural differences between human and mouse neurons and their implementation in generative AIs

Mouse and human brains have different functions that depend on their neuronal networks. We analyzed nanometer-scale three-dimensional structures of brain tissues of the mouse medial prefrontal cortex and compared them with structures of the human anterior cingulate cortex. The obtained results indicated that mouse neuronal somata are smaller and neurites are thinner than those of human neurons. We implemented these characteristics of mouse neurons in convolutional layers of a generative adversarial network (GAN) and a denoising diffusion implicit model (DDIM), which were then subjected to image generation tasks using photo datasets of cat faces, cheese, human faces, birds, and automobiles. The mouse-mimetic GAN outperformed a standard GAN in the image generation task using the cat faces and cheese photo datasets, but underperformed for human faces and birds. The mouse-mimetic DDIM gave similar results, suggesting that the nature of the datasets affected the results. Analyses of the five datasets indicated differences in their image entropy, which should influence the number of parameters required for image generation. The preferences of the mouse-mimetic AIs coincided with the impressions commonly associated with mice. The relationship between the neuronal network and brain function should be investigated by implementing other biological findings in artificial neural networks.

generative AI↗
DOE OSTI2021

Computational general relativistic force-free electrodynamics: I. Multi-coordinate implementation and testing

General relativistic force-free electrodynamics is one possible plasma-limit employed to analyze energetic outflows in which strong magnetic fields are dominant over all inertial phenomena. The amazing images of black hole (BH) shadows from the Galactic Center and the M87 galaxy provide a first direct glimpse into the physics of accretion flows in the most extreme environments of the universe. The efficient extraction of energy in the form of collimated outflows or jets from a rotating BH is directly linked to the topology of the surrounding magnetic field. We aim at providing a tool to numerically model the dynamics of such fields in magnetospheres around compact objects, such as BHs and neutron stars. To do so, we probe their role in the formation of high energy phenomena such as magnetar flares and the highly variable teraelectronvolt emission of some active galactic nuclei. In this work, we present numerical strategies capable of modeling fully dynamical force-free magnetospheres of compact astrophysical objects. Here, we provide implementation details and extensive testing of our implementation of general relativistic force-free electrodynamics in Cartesian and spherical coordinates using the infrastructure of the E INSTEIN T OOLKIT . The employed hyperbolic/parabolic cleaning of numerical errors with full general relativistic compatibility allows for fast advection of numerical errors in dynamical spacetimes. Such fast advection of divergence errors significantly improves the stability of the general relativistic force-free electrodynamics modeling of BH magnetospheres.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Development and Implementation of an Optimization Model for Hydropower and Total Dissolved Gas in the Mid-Columbia River System

Managing energy, water, and environmental priorities and constraints within a cascade hydropower system is a challenging multiobjective optimization effort that requires advanced modeling and forecasting tools. Within the mid-Columbia River system, there is currently a lack of specific solutions for predicting how coordinated operational decisions can mitigate the impacts of total dissolved gas (TDG) supersaturation while satisfying multiple additional policy and hydropower generation objectives. In this study, a reduced-order TDG uptake equation is developed that predicts tailrace TDG at seven hydropower facilities on the mid-Columbia River. The equation is incorporated into a general multiobjective river, reservoir, and hydropower optimization tool as a prioritized operating goal within a broader set of system-level objectives and constraints. A test case is presented to assess the response of TDG and hydropower generation when TDG supersaturation is optimized to remain under state water-quality standards. Satisfaction of TDG as an operating goal is highly dependent on whether constraints that limit TDG uptake are implemented at a higher priority than generation requests. According to the model, an opportunity exists to reduce TDG supersaturation and meet hydropower generation requirements by shifting spillway flows to different time periods. In conclusion, a coordinated effort between all project owners is required to implement systemwide optimized solutions that satisfy the operating policies of all stakeholders.

13 HYDRO ENERGY↗
DOE OSTI2024

Simulation Evaluation of a Large-Scale Implementation of Virtual-Phase Link-Based Model Predictive Control

Traffic congestion is a serious problem in the US, and traffic signal control is one of the effective solutions to congestion. Previous research on model predictive control (MPC)-based traffic signal control showed substantial benefits over conventional methods. This study focused on implementing MPC over a large-scale network with complex intersections and the impact of cycle length, network size, and imperfect state estimation on performances. This study implemented a virtual phase link (VPL)-based model predictive control method which used the number of vehicles in each VPL as input state variables and was suitable for National Electrical Manufacturing Association (NEMA) ring-barrier control. To test the impact of network size, the performance of distributed MPC (36 intersections in the network are divided into five subnetworks) was compared with that of MPC over the full network for a set of cycle lengths. To test the impact of imperfect state estimation, we synthetically infused estimation error and developed two scenarios, MPC-error and MPC-error narrow, which had higher and lower estimation errors, respectively. The performance of these MPC methods was compared with that of the existing time-of-day (TOD) method and an offline method that used Webster's method for split and MULTIBAND for cycle length and offset optimization. Trajectory and linkwise signal performance measures were collected from the simulation to evaluate performance. The distributed MPC method with perfect state estimation had the lowest delay and highest energy efficiency of all the methods. The performance of MPC decreased as the prediction inaccuracy increased. MPC-error had 7% and 11% more delay than MPC-error narrow in the morning and evening peaks, respectively. Overall, simulation results suggest that even with imperfect state estimation, MPC methods will outperform offline methods significantly.

large-scale simulation↗
DOE OSTI2010

Relative unitary implementability of perturbed quantum field dynamics on de Sitter space

We study the quantum dynamics of a Klein-Gordon field on de Sitter space based on the Euclidean vacuum. We prove time evolution is not unitarily implementable. We also consider a Klein-Gordon field perturbed by a local potential V. In this case we prove that the deviation from the V=0 dynamics is unitarily implementable.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2010

A Parallel Implementation of a Smoothed Particle Hydrodynamics Method on Graphics Hardware Using the Compute Unified Device Architecture

The smoothed particle hydrodynamics (SPH), which is a class of meshfree particle methods (MPMs), has a wide range of applications from micro-scale to macro-scale as well as from discrete systems to continuum systems. Graphics hardware, originally designed for computer graphics, now provide unprecedented computational power for scientific computation. Particle system needs a huge amount of computations in physical simulation. In this paper, an efficient parallel implementation of a SPH method on graphics hardware using the Compute Unified Device Architecture is developed for fluid simulation. Comparing to the corresponding CPU implementation, our experimental results show that the new approach allows significant speedups of fluid simulation through handling huge amount of computations in parallel on graphics hardware.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2010

Image processing with cellular nonlinear networks implemented on field-programmable gate arrays for real-time applications in nuclear fusion

In the past years cameras have become increasingly common tools in scientific applications. They are now quite systematically used in magnetic confinement fusion, to the point that infrared imaging is starting to be used systematically for real-time machine protection in major devices. However, in order to guarantee that the control system can always react rapidly in case of critical situations, the time required for the processing of the images must be as predictable as possible. The approach described in this paper combines the new computational paradigm of cellular nonlinear networks (CNNs) with field-programmable gate arrays and has been tested in an application for the detection of hot spots on the plasma facing components in JET. The developed system is able to perform real-time hot spot recognition, by processing the image stream captured by JET wide angle infrared camera, with the guarantee that computational time is constant and deterministic. The statistical results obtained from a quite extensive set of examples show that this solution approximates very well an ad hoc serial software algorithm, with no false or missed alarms and an almost perfect overlapping of alarm intervals. The computational time can be reduced to a millisecond time scale for 8 bit 496x560-sized images. Moreover, in our implementation, the computational time, besides being deterministic, is practically independent of the number of iterations performed by the CNN - unlike software CNN implementations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Conceptual Design, Implementation and Commissioning of Data Acquisition and Control System for Negative Ion Source at IPR

Negative ion Experimental facility has been setup at IPR. The facility consists of a RF based negative ion source (ROBIN)--procured under a license agreement with IPP Garching, as a replica of BATMAN, presently operating in IPP, 100 kW 1 MHz RF generators and a set of low and high voltage power supplies, vacuum system and diagnostics. 35 keV 10A H- beam is expected from this setup. Automated successful operation of the system requires an advanced, rugged, time proven and flexible control system. Further the data generated in the experimental phase needs to be acquired, monitored and analyzed to verify and judge the system performance. In the present test bed, this is done using a combination of PLC based control system and a PXI based data acquisition system. The control system consists of three different Siemens PLC systems viz. (1) S-7 400 PLC as a Master Control, (2) S-7 300 PLC for Vacuum system control and (3) C-7 PLC for RF generator control. Master control PLC directly controls all the subsystems except the Vacuum system and RF generator. The Vacuum system and RF generator have their own dedicated PLCs (S-7 300 and C-7 respectively). Further, these two PLC systems work as a slave for the Master control PLC system. Communication between PLC S-7 400, S-7 300 and central control room computer is done through Industrial Ethernet (IE). Control program and GUI are developed in Siemens Step-7 PLC programming software and Wincc SCADA software, respectively. There are approximately 150 analog and 200 digital control and monitoring signals required to perform complete closed loop control of the system. Since the source floats at high potential ({approx}35 kV); a combination of galvanic and fiber optic isolation has been implemented. PXI based Data Acquisition system (DAS) is a combination of PXI RT (Real time) system, front end signal conditioning electronics, host system and DAQ program. All the acquisition signals coming from various sub-systems are connected and acquired by the PXI RT system, through only fiber optics link for signal conditioning, electrical isolation and better noise immunity. Real time and Host application programs are developed in LabVIEW and the data shall be stored with a facility of online display of selected parameters. Mathematical calculations and report generation will take place at the end of each beam shot. The paper describes in detail about the design approach, implementation strategy, program development, commissioning and operational test result of ROBIN through a data acquisition and control system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2011

Results of the implementation on a mock-up of the full 3dB hybrid matching option of the ITER ICRH system

Each of the two ICRH antennas for ITER must couple 20MW to the plasma in the 40-55MHz band via an array of 24 radiating shorted straps fed by four generators. The matching system must provide automatic matching control on the mean load provided by the plasma and be resilient (parallel {Gamma}{sub G} parallel <0.2) to a wide range of fast antenna load excursions occurring in ELMy plasmas. Furthermore, good control of the current distribution in the strap array must be possible for the various heating and current drive scenarios. Two load resilient matching options have been considered for ITER: the 4 'Conjugate-T'(CT) and the 4 hybrids ones, the first being presently considered as a back-up option. Automatic control of these 2 options has been developed, and tested for optimization on a low-powered scaled mock-up. Successful implementation of the simultaneous feedback control of 11 actuators for the matching of the 4 CT and for the control of the toroidal phasing has already been achieved. The matching and the array current control of the 3dB hybrid option are provided by simultaneous feedback control of the decouplers and double stub tuners (in total 23 actuators) and this has also been successfully achieved for the full array. The paper discusses the circuit implementation and presents the obtained results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Implementation of moire-schlieren deflectometry on a small scale fast capillary plasma discharge

We present the results of an implementation of a refractive diagnostic to study fast dynamics in capillary discharges. It consists of a moire-schlieren deflectometry technique that provides a quantitative analysis of the refractive index gradients. The technique is composed of an angular deflection mapping system (moire deflectometry) and a spatial Fourier filter (schlieren). Temporal resolution of 12 ps, 50 {mu}m of spatial resolution and minimum detectable gradient of ({nabla}n{sub e}){sub min}=6x10{sup 18}cm{sup -4} were obtained. With these parameters, a large aspect ratio capillary discharge of 15 ns half period current was studied; the diagnostic was implemented axially along the alumina tube of 1.6 mm inner diameter and 21 mm length. The detectable electron density for these conditions was 1x10{sup 17}cm{sup -3}. From the interpretation of the fringe displacement, we are able to measure the velocity of the radial compression wave and the compression ratio due to the Lorentz force. On axis, electron densities of the order of 5x10{sup 17}cm{sup -3} were obtained at the time of maximum soft x-ray emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Parallel implementation of three-dimensional molecular dynamic simulation for laser-cluster interaction

The objective of this article is to report the parallel implementation of the 3D molecular dynamic simulation code for laser-cluster interactions. The benchmarking of the code has been done by comparing the simulation results with some of the experiments reported in the literature. Scaling laws for the computational time is established by varying the number of processor cores and number of macroparticles used. The capabilities of the code are highlighted by implementing various diagnostic tools. To study the dynamics of the laser-cluster interactions, the executable version of the code is available from the author.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Molecular properties of excited electronic state: Formalism, implementation, and applications of analytical second energy derivatives within the framework of the time-dependent density functional theory/molecular mechanics

This work extends our previous works [J. Liu and W. Z. Liang, J. Chem. Phys. 135, 014113 (2011); J. Liu and W. Z. Liang, J. Chem. Phys. 135, 184111 (2011)] on analytical excited-state energy Hessian within the framework of time-dependent density functional theory (TDDFT) to couple with molecular mechanics (MM). The formalism, implementation, and applications of analytical first and second energy derivatives of TDDFT/MM excited state with respect to the nuclear and electric perturbations are presented. Their performances are demonstrated by the calculations of adiabatic excitation energies, and excited-state geometries, harmonic vibrational frequencies, and infrared intensities for a number of benchmark systems. The consistent results with the full quantum mechanical method and other hybrid theoretical methods indicate the reliability of the current numerical implementation of developed algorithms. The computational accuracy and efficiency of the current analytical approach are also checked and the computational efficient strategies are suggested to speed up the calculations of complex systems with many MM degrees of freedom. Finally, we apply the current analytical approach in TDDFT/MM to a realistic system, a red fluorescent protein chromophore together with part of its nearby protein matrix. The calculated results indicate that the rearrangement of the hydrogen bond interactions between the chromophore and the protein matrix is responsible for the large Stokes shift.

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

Efficient implementation of the three-dimensional reference interaction site model method in the fragment molecular orbital method

The three-dimensional reference interaction site model (3D-RISM) method was efficiently implemented in the fragment molecular orbital (FMO) method. The method is referred to as the FMO/3D-RISM method, and allows us to treat electronic structure of the whole of a macromolecule, such as a protein, as well as the solvent distribution around a solute macromolecule. The formalism of the FMO/3D-RISM method, for the computationally available form and variational expressions, are proposed in detail. A major concern leading to the implementation of the method was decreasing the computational costs involved in calculating the electrostatic potential, because the electrostatic potential is calculated on numerous grid points in three-dimensional real space in the 3D-RISM method. In this article, we propose a procedure for decreasing the computational costs involved in calculating the electrostatic potential in the FMO method framework. The strategy involved in this procedure is to evaluate the electrostatic potential and the solvated Fock matrix in different manners, depending on the distance between the solute and the solvent. The electrostatic potential is evaluated directly in the vicinity of the solute molecule by integrating the molecular orbitals of monomer fragments of the solute molecule, whereas the electrostatic potential is described as the sum of multipole interactions when an analog of the fast multipole method is used. The efficiency of our method was demonstrated by applying it to a water trimer system and three biomolecular systems. The FMO/3D-RISM calculation can be performed within a reasonable computational time, retaining the accuracy of some physical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗