Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pathway engineering”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Downregulation of Squalene Synthase Broadly Impacts Isoprenoid Biosynthesis in Guayule

Production of natural rubber by Parthenium argentaum (guayule) requires increased yield for economic sustainability. An RNAi gene silencing strategy was used to engineer isoprenoid biosynthesis by downregulation of squalene synthase (SQS), such that the pool of farnesyl diphosphate (FPP) substrate might instead be available to initiate natural rubber synthesis. Downregulation of SQS resulted in significantly reduced squalene and slightly increased rubber, but not in the same tissues nor to the same extent, partially due to an apparent negative feedback regulatory mechanism that downregulated mevalonate pathway isoprenoid production, presumably associated with excess geranyl pyrophosphate levels. A detailed metabolomics analysis of isoprenoid production in guayule revealed significant differences in metabolism in different tissues, including in active mevalonate and methylerythritol phosphate pathways in stem tissue, where rubber and squalene accumulate. New insights and strategies for engineering isoprenoid production in guayule were identified.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The bright frontiers of microbial metabolic optogenetics

In recent years, light-responsive systems from the field of optogenetics have been applied to several areas of metabolic engineering with remarkable success. By taking advantage of light's high tunability, reversibility, and orthogonality to host endogenous processes, optogenetic systems have enabled unprecedented dynamical controls of microbial fermentations for chemical production, metabolic flux analysis, and population compositions in co-cultures. In this article, we share our opinions on the current state of this new field of metabolic optogenetics. Furthermore, we make the case that it will continue to impact metabolic engineering in increasingly new directions, with the potential to challenge existing paradigms for metabolic pathway and strain optimization as well as bioreactor operation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Microstructural Transitions during Powder Metallurgical Processing of Solute Stabilized Nanostructured Tungsten Alloys

Exploiting grain boundary engineering in the design of alloys for extreme environments provides a promising pathway for enhancing performance relative to coarse-grained counterparts. Due to its attractive properties as a plasma facing material for fusion devices, tungsten presents an opportunity to exploit this approach in addressing the significant materials challenges imposed by the fusion environment. Here, we employ a ternary alloy design approach for stabilizing W against recrystallization and grain growth while simultaneously enhancing its manufacturability through powder metallurgical processing. Mechanical alloying and grain refinement in W-10 at.% Ti-(10,20) at.% Cr alloys are accomplished through high-energy ball milling with transitions in the microstructure mapped as a function of milling time. We demonstrate the multi-modal nature of the resulting nanocrystalline grain structure and its stability up to 1300 °C with the coarser grain size population correlated to transitions in crystallographic texture that result from the preferred slip systems in BCC W. Field-assisted sintering is employed to consolidate the alloy powders into bulk samples, which, due to the deliberately designed compositional features, are shown to retain ultrafine grain structures despite the presence of minor carbides formed during sintering due to carbon impurities in the ball-milled powders.

36 MATERIALS SCIENCE↗
DOE OSTI2024

DNA parts and gene constructs for plant biodesign

Plant biodesign requires the knowledge of DNA parts (e.g., genes, promoters, terminators), along with their combinations (as gene constructs) linked to engineered traits. DNA parts with validated or predicted functions in plants have been deposited in various online databases. However, these existing databases focus on basic biological functions of individual DNA parts, leaving a gap between basic knowledge and bioengineering applications. To fill this knowledge gap, we have created a user-friendly, open-ended database as a knowledge graph linking DNA parts to gene constructs to traits. This database contains experimentally validated DNA parts and gene constructs documented in peer-reviewed publications. The DNA parts include 1) molecular components with biological functions, such as genes involved in various biological processes (e.g., metabolic and signal transduction pathways) and 2) molecular components with technical functions, such as gene expression, genome engineering and sequence splicing. The gene constructs deposited in this database include both single-gene and multi-gene constructs. This database allows users to submit DNA parts and gene construct compositions linked to engineered traits described in peer-reviewed publications, providing a public digital repository for sharing the biodesign information among the researchers in the fields of plant biotechnology and plant synthetic biology.

plant biodesign synthetic biology gene constructs ↗
DOE OSTI2026

Final Technical Report for DE-SC0022206

This project developed foundational genetic, genomic, and epigenetic tools for anaerobic fungi (Neocallimastigomycota), a group of microorganisms with exceptional natural abilities to deconstruct lignocellulosic biomass. Efficient biomass deconstruction remains a major barrier to economical production of renewable fuels, chemicals, and materials from agricultural and forestry residues. The project sought to enable mechanistic studies and future engineering of anaerobic fungi by improving genomic resources, establishing methods for gene expression, and investigating epigenetic regulation of biomass-degrading pathways. Major accomplishments included generation of the first chromosome-scale genome assemblies for multiple anaerobic fungal species, providing publicly available genomic resources that support both engineering and fundamental biological research. The project established the first reproducible system for heterologous gene expression in anaerobic fungi and identified genomic features and mobile genetic elements that may support future development of stable transformation technologies. In parallel, the project demonstrated direct conversion of untreated lignocellulosic biomass into fuels and specialty chemicals through a fungal-yeast bioprocess and identified anaerobic fungal enzymes with utility for metabolic engineering. The research also revealed that epigenetic regulation plays an important role in controlling fungal gene expression and enzyme production, identifying potential strategies for enhancing biomass degradation. Collectively, this work established anaerobic fungi as a tractable emerging platform for bioenergy and biomanufacturing research, generated valuable public resources, trained the next generation of researchers, and advanced DOE-BER goals related to predictive biology, sustainable bioprocessing, and the circular bioeconomy.

Solomon, Kevin [University of Delaware] (ORCID:000↗
DOE OSTI2020

Development of a clostridia-based cell-free system for prototyping genetic parts and metabolic pathways

Gas fermentation by autotrophic bacteria, such as clostridia, offers a sustainable path to numerous bioproducts from a range of local, highly abundant, waste and low-cost feedstocks, such as industrial flue gases or syngas generated from biomass or municipal waste. Unfortunately, designing and engineering clostridia remains laborious and slow. The ability to prototype individual genetic part function, gene expression patterns, and biosynthetic pathway performance in vitro before implementing designs in cells could help address these bottlenecks by speeding up design. Unfortunately, a high-yielding cell-free gene expression (CFE) system from clostridia has yet to be developed. Here in this paper, we report the development and optimization of a high-yielding (236 ± 24 μg/mL) batch CFE platform from the industrially relevant anaerobe, Clostridium autoethanogenum. A key feature of the platform is that both circular and linear DNA templates can be applied directly to the CFE reaction to program protein synthesis. We demonstrate the ability to prototype gene expression, and quantitatively map aerobic cell-free metabolism in lysates from this system. We anticipate that the C. autoethanogenum CFE platform will not only expand the protein synthesis toolkit for synthetic biology, but also serve as a platform in expediting the screening and prototyping of gene regulatory elements in non-model, industrially relevant microbes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Engineering-scale Modeling of High-Temperature Creep and Creep Crack Growth in Alloy 316H

This document demonstrates completion of the goals described in the technical narrative of the Department of Energy’s Industry Funding Opportunity Announcement (iFOA) project entitled “Modeling and Simulation Development Pathways to Accelerating KP-FHR Licensing,” which relates to the development and demonstration of capabilities for conducting engineering-scale simulations of Alloy 316H components under high-temperature conditions. This work encompassed two major aspects: integrating and testing constitutive models for the creep response of 316H at high temperatures, and developing tools for modeling creep crack growth (CCG) in 316H. Two classes of constitutive models were used in this effort: phenomenological models based on behavior observed at the engineering scale, and reduced-order models (ROMs) that represent the nonlinear response of mesoscale models that capture the sensitivity to material microstructure and processing. Likewise, the approaches employed for CCG modeling considered both simplified engineering approaches and detailed simulations of creep and damage ahead of the crack tip. These developments, which were performed by utilizing the Grizzly code as well as the open-source libraries it depends on, strengthen Grizzly’s ability to support licensing and safety analyses of high-temperature reactor components.

316H↗
DOE OSTI2022

Improved genome editing by an engineered CRISPR-Cas12a

CRISPR-Cas12a is an RNA-guided, programmable genome editing enzyme found within bacterial adaptive immune pathways. Unlike CRISPR-Cas9, Cas12a uses only a single catalytic site to both cleave target double-stranded DNA (dsDNA) (cis-activity) and indiscriminately degrade single-stranded DNA (ssDNA) (trans-activity). To investigate how the relative potency of cis- versus trans-DNase activity affects Cas12a-mediated genome editing, we first used structure-guided engineering to generate variants of Lachnospiraceae bacterium Cas12a that selectively disrupt trans-activity. The resulting engineered mutant with the biggest differential between cis - and trans -DNase activity in vitro showed minimal genome editing activity in human cells, motivating a second set of experiments using directed evolution to generate additional mutants with robust genome editing activity. Notably, these engineered and evolved mutants had enhanced ability to induce homology-directed repair (HDR) editing by 2–18-fold compared to wild-type Cas12a when using HDR donors containing mismatches with crRNA at the PAM-distal region. Finally, a site-specific reversion mutation produced improved Cas12a (iCas12a) variants with superior genome editing efficiency at genomic sites that are difficult to edit using wild-type Cas12a. This strategy establishes a pipeline for creating improved genome editing tools by combining structural insights with randomization and selection. The available structures of other CRISPR-Cas enzymes will enable this strategy to be applied to improve the efficacy of other genome-editing proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Engineering Relaxor Behavior in (BaTiO 3 ) n /(SrTiO 3 ) n Superlattices

Complex-oxide superlattices provide a pathway to numerous emergent phenomena because of the juxtaposition of disparate properties and the strong interfacial interactions in these unit-cell-precise structures. This is particularly true in superlattices of ferroelectric and dielectric materials, wherein new forms of ferroelectricity, exotic dipolar textures, and distinctive domain structures can be produced. Here, in this study, relaxor-like behavior, typically associated with the chemical inhomogeneity and complexity of solid solutions, is observed in (BaTiO 3 ) n /(SrTiO 3 ) n (n = 4–20 unit cells) superlattices. Dielectric studies and subsequent Vogel–Fulcher analysis show significant frequency dispersion of the dielectric maximum across a range of periodicities, with enhanced dielectric constant and more robust relaxor behavior for smaller period n. Bond-valence molecular-dynamics simulations predict the relaxor-like behavior observed experimentally, and interpretations of the polar patterns via 2D discrete-wavelet transforms in shorter-period superlattices suggest that the relaxor behavior arises from shape variations of the dipolar configurations, in contrast to frozen antipolar stripe domains in longer-period superlattices (n = 16). Moreover, the size and shape of the dipolar configurations are tuned by superlattice periodicity, thus providing a definitive design strategy to use superlattice layering to create relaxor-like behavior which may expand the ability to control desired properties in these complex systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bulk nanocrystalline Al alloys with hierarchical reinforcement structures via grain boundary segregation and complexion formation

Grain size engineering, particularly reducing grain size into the nanocrystalline regime, offers a promising pathway to further improve the strength-to-weight ratio of Al alloys. Unfortunately, the fabrication of nanocrystalline metals often requires non-equilibrium processing routes, which typically limit the specimen size and require large energy budgets. In this study, multiple dopant elements in ternary Al alloys are deliberately selected to enable segregation to the grain boundary region and promote the formation of amorphous complexions. Three different fully dense bulk nanocrystalline Al alloys (Al-Mg-Y, Al-Fe-Y, and Al-Ni-Y) with small grain sizes were successfully fabricated using a simple powder metallurgy approach, with full densification connected directly to the onset of amorphous complexion formation. All the compositions demonstrate densities above 99% with grain sizes <60 nm following consolidation via hot pressing at 585 °C. The very fine grain structure results in excellent mechanical properties, as evidenced by nanoindentation hardness values in the range of 2.2-2.8 GPa. Detailed microstructural characterization verifies the segregation of all dopant species to grain boundaries as well as the formation of amorphous complexions, which suggests their influential role in aiding effective consolidation and endowing thermal stability in the alloys. Moreover, nanorods with a core-shell structure are also observed at the grain boundaries, which likely contribute to the stabilization of the grain structure while also strengthening the materials. Lastly, intermetallic particles with sizes of hundreds of nanometers form in all systems. As a whole, the results presented here demonstrate a general alloy design strategy of segregation and boundary evolution pathway that enables the fabrication of multiple nanocrystalline Al alloys with hierarchical microstructures and improved performance.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An elastin-like polymer targeting vascular endothelial growth factor receptor-1 reduces survival in serum-starved endothelial cells

Peptides often exhibit biological activity that depends on the context in which they are displayed and delivered. Understanding and controlling these contextual effects on peptide function is critical for designing targeted and responsive peptide-based biomaterials and therapeutics. Genetically engineered protein polymers such as elastin-like polypeptides (ELPs) can incorporate bioactive peptide motifs and are attractive candidates for biomaterials used in tissue engineering and targeted drug delivery. They also present an opportunity for investigating and modulating cell signaling pathways by presenting a peptide ligand in various defined chemical and physical environments. Vascular endothelial growth factor receptor-1 (VEGFR1) signaling plays important and complex roles in cell survival and angiogenesis, but polymeric materials that interact with this signaling axis are scarce. In this study, a novel genetically engineered elastin-like polymer that targets VEGFR1 is characterized. This polymer, termed R1B-ELP, binds to human endothelial cells in a manner dependent on its VEGFR1-targeting motif and, based on cell proliferation and cytotoxicity assays, demonstrates activity consistent with disrupting pro-survival signaling necessary for endothelial cell function under conditions of environmental stress. Notably, these findings indicate that ELP fusion alters the functional behavior of the targeting peptide. Modulators of VEGFR1 signaling have potential applications in basic studies of angiogenesis as well as in therapeutic applications targeting vascular or inflammatory diseases.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Direct time-resolved detection and quantification of key reactive intermediates in diethyl ether oxidation at T = 450–600 K

High-pressure multiplexed photoionization mass spectrometry (MPIMS) with tunable vacuum ultraviolet (VUV) ionization radiation from the Lawrence Berkeley Labs Advanced Light Source is used to investigate the oxidation of diethyl ether (DEE). Kinetics and photoionization (PI) spectra are simultaneously measured for the species formed. Several stable products from DEE oxidation are identified and quantified using reference PI cross-sections. In addition, we directly detect and quantify three key chemical intermediates: peroxy (ROO˙), hydroperoxyalkyl peroxy (˙OOQOOH), and ketohydroperoxide (HOOP[double bond, length as m-dash]O, KHP). These intermediates undergo dissociative ionization (DI) into smaller fragments, making their identification by mass spectrometry challenging. Furthermore, with the aid of quantum chemical calculations, we identify the DI channels of these key chemical species and quantify their time-resolved concentrations from the overall carbon atom balance at T = 450 K and P = 7500 torr. This allows the determination of the absolute PI cross-sections of ROO˙, ˙OOQOOH, and KHP into each DI channel directly from experiment. The PI cross-sections in turn enable the quantification of ROO˙, ˙OOQOOH, and KHP from DEE oxidation over a range of experimental conditions that reveal the effects of pressure, O 2 concentration, and temperature on the competition among radical decomposition and second O 2 addition pathways.

42 ENGINEERING↗
DOE OSTI2020

Hydrodynamic and frictional modulation of deformations in switchable colloidal crystallites

Displacive transformations in colloidal crystals may offer a pathway for increasing the diversity of accessible configurations without the need to engineer particle shape or interaction complexity. To date, binary crystals composed of spherically symmetric particles at specific size ratios have been formed that exhibit floppiness and facile routes for transformation into more rigid structures that are otherwise not accessible by direct nucleation and growth. There is evidence that such transformations, at least at the micrometer scale, are kinetically influenced by concomitant solvent motion that effectively induces hydrodynamic correlations between particles. Here, we study quantitatively the impact of such interactions on the transformation of binary bcc-CsCl analog crystals into close-packed configurations. We first employ principal-component analysis to stratify the explorations of a bcc-CsCl crystallite into orthogonal directions according to displacement. We then compute diffusion coefficients along the different directions using several dynamical models and find that hydrodynamic correlations, depending on their range, can either enhance or dampen collective particle motions. These two distinct effects work synergistically to bias crystallite deformations toward a subset of the available outcomes.

42 ENGINEERING↗
DOE OSTI2020

Internal Radiation Dose Evaluation For An Unruptured Post Release Tristructural Istropic Fuel Particle for Advanced and Micro-reactor Applications

There are unique benefits from advanced/micro reactor designs and fuel types that offer safety features in the case of an accident that may reduce environmental consequences compared to conventional reactors and fuels. Several micro reactor designs, for both civilian and military applications, are proposed using tristructural-isotropic (TRISO) fuel. TRISO particles are a robust advanced nuclear fuel type that leads to the unique question of how unruptured, activated TRISO particles will interact with humans. Human exposure to nuclear fuel, or radionuclides, includes all pathways of exposure. However, TRISO particles are 0.9 mm in size, and that particle size restricts internal dose assessment to the ingestion pathway; eliminating the inhalation and dermal absorption pathways. Furthermore, activity of the TRISO particle was established by High Temperature Engineering Test Reactor simulations assuming the TRISO particle was one day post-release from the reactor. International Commission of Radiological Protection publications 100 and 133 were used to develop the internal dosimetry model for an unruptured TRISO particle for the International Commission of Radiological Protection reference male and female dose assessment. The rectosigmoid colon had the highest internal organ dose for both reference male (2.1 Sv) and reference female (2.3 Sv). The internal dose from ingestion of the scenario-defined TRISO particle was 0.25 Sv for the reference male and 0.29 Sv for the refence female, which exceeds the annual occupational effective dose limit of 0.05 Sv in the code of federal regulations, 10.CFR.20 Subpart C.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models

Mathematical modeling of plant metabolism allows the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In current research using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For 13 C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Report on the Integration of Experimental and Modeling Data for Initial Equivalence Study of Mechanical Performance in Irradiated LPBF 316 Stainless Steel

To ensure the rapid development, deployment, and use of advanced nuclear technologies, faster qualification approaches are needed. Typically, the primary pathway uses traditional data packages consistent with the American Society of Mechanical Engineers Boiler and Pressure Vessel Code, which does not consider the environmental effects the material will experience such as corrosion and radiation damage. Examining radiation effects requires a significant amount of space in US facilities at the Advanced Test Reactor and the High Flux Isotope Reactor (HFIR) and suffers from natural gradients in temperature and neutron flux profiles. Ion irradiation may enable rapid assessment of radiation-induced damage to a material and is proposed as part of an accelerated materials qualification framework through the Advanced Materials and Manufacturing Technologies program. To enhance the utility of ion irradiation as an examination tool, this report provides the initial assessment of engineering-relevant properties of microstructures produced from ion irradiation in the near-surface volume. Nanoindentation, Vickers hardness, and known tensile properties were brought together with simple mathematical models and experimental data for an initial equivalence study of the mechanical performance of irradiated laser powder bed fusion (LPBF) 316 stainless steels across length scales. Direct observation of the calculated ion irradiation yield stress and measured neutron irradiation yield stress at 2 dpa showed that both datasets exhibit the same trend with irradiation temperature and overlap within an acceptable band of stress values. Ion irradiations at 10 dpa serve as a prediction of properties to compare to postirradiation examination of HFIR-irradiated LPBF 316H further in the program. This work is a significant demonstration of the Licensing Approach with Ions and Neutrons, which uses ion irradiations to generate mechanical property information more rapidly than through neutron irradiations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Coupled Microbial-Conversion and Computational-Fluid-Dynamics (CFD) Models for Butanediol Production in Micro-Aerated Reactors

Microbial conversion of substrates to macromolecules has been widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries. These bioreactions are also being investigated in the production of low-value commodities such as biofuels [1]. Gas-liquid mass-transfer and transport-reaction coupling are important challenges when designing and scaling up these reactor systems. Experiments in wellmixed small-scale reactors have enabled characterization of microbial reactivity, while their coupling with macroscale transport remains relatively unexplored. In this work, we use a coupled metabolic-CFD model to study the action of a genetically engineered microbe Zymomonas mobilis [2] on sugars to produce 2,3-Butanediol (BDO). BDO is an important hydrocarbon intermediate that can be catalytically upgraded to several fuels and chemicals [3]. An important aspect to this particular microbial conversion is the need for micro-aerated environments as opposed to traditional aerobic fermentation. Slight variations in oxygen concentration can result in competing reaction pathways that disable BDO production. Hence, gas-liquid mass transfer and transport need to be optimized in large-scale reactors to maximize BDO production, for which CFD is a valuable tool. The aerobic-fermentation CFD model previously developed by the authors [5] for simulating bubble-column and airlift reactors at scale was used in this study. The Reynolds-averaged mass, momentum and energy transport equations for interpenetrating gas and liquid phase are solved in this model along with the transport and interphase mass transfer of oxygen. Our previous work used a phenomenological model for microbial oxygen uptake that neglected microbial growth and other reaction pathways. In this work, a detailed metabolic model enabled prediction of product formation and inhibition pathways. In order to manage computational cost, we used a subcycling technique [6] that takes advantage of the clear separation in transport (~ 200 sec) and reaction (~ 2-3 hours) timescales. The CFD model is first solved to steady state, after which the metabolic model is advanced at every cell in the computational domain using the local oxygen concentration. The CFD model is then run to achieve a new steady state that provides a new oxygen distribution for the metabolic model. This process, where reaction and fluid updates are interleaved together, is iterated until reactants are completely exhausted. This work will examine the performance of different reactor designs such as bubble column and airlift reactors at scale (250-500 m3). Oxygen mass-transfer coefficient and distribution are critically analyzed among reactors, and optimization studies pertaining to aeration is presented. Furthermore, it has been observed in experiments that high BDO production may be achieved by manipulating the aerobic environment over the course of reaction, such that oxygen concentration is high during the growth phase, and very low as sugar is depleted. This characteristic will be addressed by our simulations for which a timedependent scheduling strategy for aeration is presented that maximizes BDO production. [1] Humbird, D., Davis, R., and McMillan, J., Aeration costs in stirred-tank and bubble column bioreactors, Biochemical Engineering Journal, 127, 161—166, 2017 [2] Yang, S., Mohagheghi, A., Franden, M. A., Chou, Y.-C., Chen, X., Dowe, N., Himmel, M. E., and Zhang, M., Metabolic engineering of zymomonas mobilis for 2, 3-butanediol production from lignocellulosic biomass sugars. Biotechnology for biofuels, 9(1):189, 2016 [3] Kim, S. J., Sim, H. J., Kim, J. W., Lee, Y. G., Park, Y. C., and Seo, J. H., Enhanced production of 2,3-butanediol from xylose by combinatorial engineering of xylose metabolic pathway and cofactor regeneration in pyruvate decarboxylase-deficient Saccharomyces cerevisiae. Bioresource Technology, 245:1551–1557, 2017 [4] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998

29 ENERGY PLANNING, POLICY, AND ECONOMY↗