Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible computation”

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

Transport and Dynamics in Toroidal Fusion Systems

The study entitled, "Transport and Dynamics in Toroidal Fusion Systems," (TDTFS) applied analytical theory and numerical computation to investigate topics of importance to confining plasma, the fourth state of matter, with magnetic fields. A central focus of the work is how non-thermal components of the ion particle distribution affect the "sawtooth" collective oscillation in the core of the tokamak magnetic configuration. Previous experimental and analytical research had shown and described how the oscillation frequency decreases and amplitude increases, leading to "monster" or "giant" sawteeth, when the non-thermal component is increased by injecting particle beams or by exciting ions with imposed electromagnetic waves. The TDTFS study applied numerical computation to self-consistently simulate the interaction between macroscopic collective plasma dynamics and the non-thermal particles. The modeling used the NIMROD code [Sovinec, Glasser, Gianakon, et al., J. Comput. Phys. 195, 355 (2004)] with the energetic component represented by simulation particles [Kim, Parker, Sovinec, and the NIMROD Team, Comput. Phys. Commun. 164, 448 (2004)]. The computations found decreasing growth rates for the instability that drives the oscillations, but they were ultimately limited from achieving experimentally relevant parameters due to computational practicalities. Nonetheless, this effort provided valuable lessons for integrated simulation of macroscopic plasma dynamics. It also motivated an investigation of the applicability of fluid-based modeling to the ion temperature gradient instability, leading to the journal publication [Schnack, Cheng, Barnes, and Parker, Phys. Plasmas 20, 062106 (2013)]. Apart from the tokamak-specific topics, the TDTFS study also addressed topics in the basic physics of magnetized plasma and in the dynamics of the reversed-field pinch (RFP) configuration. The basic physics work contributed to a study of two-fluid effects on interchange dynamics, where "two-fluid" refers to modeling independent dynamics of electron and ion species without full kinetic effects. In collaboration with scientist Ping Zhu, who received separate support, it was found that the rule-of-thumb criteria on stabilizing interchange has caveats that depend on the plasma density and temperature profiles. This work was published in [Zhu, Schnack, Ebrahimi, et al., Phys. Rev. Lett. 101, 085005 (2008)]. An investigation of general nonlinear relaxation with fluid models was partially supported by the TDTFS study and led to the publication [Khalzov, Ebrahimi, Schnack, and Mirnov, Phys. Plasmas 19, 012111 (2012)]. Work specific to the RFP included an investigation of interchange at large plasma pressure and support for applications [for example, Scheffel, Schnack, and Mirza, Nucl. Fusion 53, 113007 (2013)] of the DEBS code [Schnack, Barnes, Mikic, Harned, and Caramana, J. Comput. Phys. 70, 330 (1987)]. Finally, the principal investigator over most of the award period, Dalton Schnack, supervised a numerical study of modeling magnetic island suppression [Jenkins, Kruger, Hegna, Schnack, and Sovinec, Phys. Plasmas 17, 12502 (2010)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Computer–aided process intensification: Challenges, trends and opportunities

Chemicals based products and the processes that manufacture them play important roles in sustaining our society. Chemical, biochemical, petrochemical and related processes are characterized by either conversion of specific raw materials to desired chemicals or extraction of desired chemicals from specific raw materials. Well–known examples in the first category are the conversion of methanol and acetic acid to form methyl acetate and water, and, carbon dioxide hydrogenation to methanol via the reverse water gas shift reaction. In the second category, the use of de–methanizer, de–ethanizer, de–propanizer, and so on, separation steps in gas processing plants are well–known. In the conversion–based processes, as the conversions of raw materials are usually not 100% complete, energy intensive downstream separation techniques need to be employed to separate the desired products and recycle the unreacted raw materials. Note that the conversion process itself may also be energy intensive. In the extraction–based processes where product chemicals do not usually represent 100% of the raw material, energy intensive separation techniques are commonly employed to achieve the desired separations. This means that for both categories of processes, energy in different forms need to be supplied and the processes during operation have the potential to produce waste that may lead to negative environmental impacts and cause supply chain issues. Furthermore, planet earth is facing challenges with respect to energy, water, food, and environment, and significant improvements in process efficiencies of new and/or existing processes are needed.

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

Parafoveal Target Detectability Reversal Predicted by Local Luminance and Contrast Gain Control

This project is part of a program to develop image discrimination models for the prediction of the detectability of objects in a range of backgrounds. We wanted to see if the models could predict parafoveal object detection as well as they predict detection in foveal vision. We also wanted to make our simplified models more general by local computation of luminance and contrast gain control. A signal image (0.78 x 0.17 deg) was made by subtracting a simulated airport runway scene background image (2.7 deg square) from the same scene containing an obstructing aircraft. Signal visibility contrast thresholds were measured in a fully crossed factorial design with three factors: eccentricity (0 deg or 4 deg), background (uniform or runway scene background), and fixed-pattern white noise contrast (0%, 5%, or 10%). Three experienced observers responded to three repetitions of 60 2IFC trials in each condition and thresholds were estimated by maximum likelihood probit analysis. In the fovea the average detection contrast threshold was 4 dB lower for the runway background than for the uniform background, but in the parafovea, the average threshold was 6 dB higher for the runway background than for the uniform background. This interaction was similar across the different noise levels and for all three observers. A likely reason for the runway background giving a lower threshold in the fovea is the low luminance near the signal in that scene. In our model, the local luminance computation is controlled by a spatial spread parameter. When this parameter and a corresponding parameter for the spatial spread of contrast gain were increased for the parafoveal predictions, the model predicts the interaction of background with eccentricity.

Ahumada, Albert J., Jr.↗
DOE OSTI2023

Inferring the Energetics of CO 2 –Aniline Adduct Formation from Vibrational Spectroscopy

Control of atmospheric CO 2 is an important contemporary scientific and engineering challenge. Towards this goal, the reaction of CO 2 with amines to form carbamate bonds is an established method for CO 2 capture. However, controllable reversal of this reaction remains difficult and requires tuning the energetics of the carbamate bond. Through IR spectroscopy, we show that a characteristic frequency observed upon carbamate formation varies as a function of the substituent’s Hammett parameter for a family of para- substituted anilines. We present computational evidence that the vibrational frequency of the adducted CO 2 serves as a predictor of the energy of formation of the carbamate. Electron donating groups typically enhance the driving force of carbamate formation by transferring more charge to the adducted CO 2 and thus increasing the occupancy of the anti-bonding orbital in the carbon-oxygen bonds. Increased occupancy of the anti-bonding orbital within adducted CO 2 indicates a weaker bond, leading to a red shift in the characteristic carbamate frequency. As a result, our work serves the large field of CO 2 capture research where spectroscopic observables, such as IR frequencies, are more easily obtainable and can stand in as a descriptor of driving forces.

Amines↗
DOE OSTI2015

Similarity solution to three dimensional boundary layer flow of second grade nanofluid past a stretching surface with thermal radiation and heat source/sink

Development of human society greatly depends upon solar energy. Heat, electricity and water from nature can be obtained through solar power. Sustainable energy generation at present is a critical issue in human society development. Solar energy is regarded one of the best sources of renewable energy. Hence the purpose of present study is to construct a model for radiative effects in three-dimensional of nanofluid. Flow of second grade fluid by an exponentially stretching surface is considered. Thermophoresis and Brownian motion effects are taken into account in presence of heat source/sink and chemical reaction. Results are derived for the dimensionless velocities, temperature and concentration. Graphs are plotted to examine the impacts of physical parameters on the temperature and concentration. Numerical computations are presented to examine the values of skin-friction coefficients, Nusselt and Sherwood numbers. It is observed that the values of skin-friction coefficients are more for larger values of second grade parameter. Moreover the radiative effects on the temperature and concentration are quite reverse.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Specific Energy Analysis for the Manufacturing of Light-Weight Automobile Body

This manuscript analyses the specific energy requirements to produce two hybridised vehicle structures on current automotive assembly lines. It provides quantitative and specific electric (kWhr/vehicle) and fossil-fuel (MMBtu/vehicle) energy predictions of the body-in-white (BIW) manufacturing processes including body panel forming, welding, painting, and final assembly, in addition to the facility heating ventilation and air conditioning or HVAC consumption. Two hybridised BIW design criteria are analysed; the first is based on minimising the added cost per unit weight saved (+$/-kg), while the second is based on maximising the percentage of weight saved (% -kg). Also a new criterion; that is the added specific-energy per unit weight-saved or (+kWhr/vehicle/-kg) is proposed. The light-weight structures comprise body parts made of aluminium and steel. To compute energy consumptions, the study utilises the energy performance index (EPI) stochastic model in addition to complementing it with a panel-forming energy model. Cost analysis is also provided for all the major vehicle body-panels using an in-house Excel macro. The study reverse engineers a passenger vehicle BIW using a full-size coordinate measuring machine (CMM). The presented work further highlights the monetary value of the associated emissions when increasing the specific energy.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2013

Detector Simulation and WIMP Search Analysis for the Cryogenic Dark Matter Search Experiment

Astrophysical and cosmological measurements on the scales of galaxies, galaxy clusters, and the universe indicate that 85% of the matter in the universe is composed of dark matter, made up of non-baryonic particles that interact with cross-sections on the weak scale or lower. Hypothetical Weakly Interacting Massive Particles, or WIMPs, represent a potential solution to the dark matter problem, and naturally arise in certain Standard Model extensions. The Cryogenic Dark Matter Search (CDMS) collaboration aims to detect the scattering of WIMP particles from nuclei in terrestrial detectors. Germanium and silicon particle detectors are deployed in the Soudan Underground Laboratory in Minnesota. These detectors are instrumented with phonon and ionization sensors, which allows for discrimination against electromagnetic backgrounds, which strike the detector at rates orders of magnitude higher than the expected WIMP signal. This dissertation presents the development of numerical models of the physics of the CDMS detectors, implemented in a computational package collectively known as the CDMS Detector Monte Carlo (DMC). After substantial validation of the models against data, the DMC is used to investigate potential backgrounds to the next iteration of the CDMS experiment, known as SuperCDMS. Finally, an investigation of using the DMC in a reverse Monte Carlo analysis of WIMP search data is presented.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Stormtime transport of ring current and radiation belt ions

This is an investigation of stormtime particle transport that leads to formation of the ring current. Our method is to trace the guiding-center motion of representative ions (having selected first adiabatic invariants mu) in response to model substorm-associated impulses in the convection electric field. We compare our simulation results qualitatively with existing analytically tractable idealizations of particle transport (direct convective access and radial diffusion) in order to assess the limits of validity of these approximations. For mu approximately less than 10 MeV/G (E approximately less than 10 keV at L equivalent to 3) the ion drift period on the final (ring-current) drift shell of interest (L equivalent to 3) exceeds the duration of the main phase of our model storm, and we find that the transport of ions to this drift shell is appropriately idealized as direct convective access, typically from open drift paths. Ion transport to a final closed drift path from an open (plasma-sheet) drift trajectory is possible for those portions of that drift path that lie outside the mean stormtime separatrix between closed and open drift trajectories, For mu approximately 10-25 MeV/G (110 keV approximately less than E approximately less than 280 keV at L equivalent to 3) the drift period at L equivalent to 3 is comparable to the postulated 3-hr duration of the storm, and the mode of transport is transitional between direct convective access and transport that resembles radial diffusion. (This particle population is transitional between the ring current and radiation belt). For mu approximately greater than 25 MeV/G (radiation-belt ions having E approximately greater than 280 keV at L equivalent to 3) the ion drift period is considerably shorter than the main phase of a typical storm, and ions gain access to the ring-current region essentially via radial diffusion. By computing the mean and mean-square cumulative changes in 1/L among (in this case) 12 representative ions equally spaced in drift time around the steady-state drift shell of interest (L equivalent to 3), we have estimated (from both our forward and our time-reversed simulations) the time-integrated radial-diffusion coefficients D(sup sim)(sub LL) for particles having selected values of mu approximately greater than 15 MeV/G. The results agree surprisingly well with the predictions (D(sup ql)(sub LL)) of quasilinear radial diffusion theory, despite the rather brief duration (approximately 3 hrs) of our model storm and despite the extreme variability (with frequency) of the spectral-density function that characterizes the applied electric field during our model storm. As expected, the values of D(sup sim)(sub LL) deduced (respectively) from our forward and time-reversed simulations agree even better with each other and with D(sup sim)(sub LL) when the impulse amplitudes which characterize the individual substorms of our model storm are systematically reduced.

Chen, Margaret W.↗
NASA NTRS1999

A Numerical Analysis of Heat Transfer and Effectiveness on Film Cooled Turbine Blade Tip Models

A computational study has been performed to predict the distribution of convective heat transfer coefficient on a simulated blade tip with cooling holes. The purpose of the examination was to assess the ability of a three-dimensional Reynolds-averaged Navier-Stokes solver to predict the rate of tip heat transfer and the distribution of cooling effectiveness. To this end, the simulation of tip clearance flow with blowing of Kim and Metzger was used. The agreement of the computed effectiveness with the data was quite good. The agreement with the heat transfer coefficient was not as good but improved away from the cooling holes. Numerical flow visualization showed that the uniformity of wetting of the surface by the film cooling jet is helped by the reverse flow due to edge separation of the main flow.

Ameri, A. A.↗
DOE OSTI2018

All-Solid-State Synaptic Transistor with Ultralow Conductance for Neuromorphic Computing

Electronic synaptic devices are important building blocks for neuromorphic computational systems that can go beyond the constraints of von Neumann architecture. Although two-terminal memristive devices are demonstrated to be possible candidates, they suffer from several shortcomings related to the filament formation mechanism including nonlinear switching, write noise, and high device conductance, all of which limit the accuracy and energy efficiency. Electrochemical three-terminal transistors, in which the channel conductance can be tuned without filament formation provide an alternative platform for synaptic electronics. In this work, an all-solid-state electrochemical transistor made with Li ion–based solid dielectric and 2D α-phase molybdenum oxide (α-MoO 3 ) nanosheets as the channel is demonstrated. These devices achieve nonvolatile conductance modulation in an ultralow conductance regime (<75 nS) by reversible intercalation of Li ions into the α-MoO 3 lattice. Based on this operating mechanism, the essential functionalities of synapses, such as short- and long-term synaptic plasticity and bidirectional near-linear analog weight update are demonstrated. Simulations using the handwritten digit data sets demonstrate high recognition accuracy (94.1%) of the synaptic transistor arrays. These results provide an insight into the application of 2D oxides for large-scale, energy-efficient neuromorphic computing networks.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

PTM‐Psi : A python package to facilitate the computational investigation of p ost‐ t ranslational m odification on p rotein s tructures and their i mpacts on dynamics and functions

Abstract Post‐translational modification (PTM) of a protein occurs after it has been synthesized from its genetic template, and involves chemical modifications of the protein's specific amino acid residues. Despite of the central role played by PTM in regulating molecular interactions, particularly those driven by reversible redox reactions, it remains challenging to interpret PTMs in terms of protein dynamics and function because there are numerous combinatorially enormous means for modifying amino acids in response to changes in the protein environment. In this study, we provide a workflow that allows users to interpret how perturbations caused by PTMs affect a protein's properties, dynamics, and interactions with its binding partners based on inferred or experimentally determined protein structure. This Python‐based workflow, called PTM‐Psi , integrates several established open‐source software packages, thereby enabling the user to infer protein structure from sequence, develop force fields for non‐standard amino acids using quantum mechanics, calculate free energy perturbations through molecular dynamics simulations, and score the bound complexes via docking algorithms. Using the S ‐nitrosylation of several cysteines on the GAP2 protein as an example, we demonstrated the utility of PTM‐Psi for interpreting sequence–structure–function relationships derived from thiol redox proteomics data. We demonstrate that the S ‐nitrosylated cysteine that is exposed to the solvent indirectly affects the catalytic reaction of another buried cysteine over a distance in GAP2 protein through the movement of the two ligands. Our workflow tracks the PTMs on residues that are responsive to changes in the redox environment and lays the foundation for the automation of molecular and systems biology modeling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Design considerations and experimental observations for the TAMU air-cooled reactor cavity cooling system for the VHTR

The Reactor Cavity Cooling System (RCCS) is a promising passive decay heat removal system for the Very High Temperature Reactor (VHTR) to ensure reliability of the transfer of the core residual and decay heat to the environment under all off-normal circumstances. A small scale experimental test facility was constructed at Texas A and M University (TAMU) to study pertinent multifaceted thermal hydraulic phenomena in the air-cooled reactor cavity cooling system (RCCS) design based on the General Atomics (GA) concept for the Modular High Temperature Gas-Cooled Reactor (MHTGR). The TAMU Air-Cooled Experimental Test Facility is ⅛ scale from the proposed GA-MHTGR design. Groundwork for experimental investigations focusing into the complex turbulence mixing flow behavior inside the upper plenum is currently underway. The following paper illustrates some of the chief design considerations used in construction of the experimental test facility, complete with an outline of the planned instrumentation and data acquisition methods. Computational Fluid Dynamics (CFD) simulations were carried out to furnish some insights on the overall behavior of the air flow in the system. CFD simulations assisted the placement of the flow measurement sensors location. Preliminary experimental observations of experiments at 120oC inlet temperature suggested the presence of flow reversal for cases involving single active riser at both 5 m/s and 2.25 m/s, respectively and four active risers at 2.25 m/s. Flow reversal may lead to thermal stratification inside the upper plenum by means of steady state temperature measurements. A Particle Image Velocimetry (PIV) experiment was carried out to furnish some insight on flow patterns and directions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

The chemical imprint of silicate dust on the most metal-poor stars

We investigate the impact of dust-induced gas fragmentation on the formation of the first low-mass, metal-poor stars (<1 M {sub ☉}) in the early universe. Previous work has shown the existence of a critical dust-to-gas ratio, below which dust thermal cooling cannot cause gas fragmentation. Assuming that the first dust is silicon-based, we compute critical dust-to-gas ratios and associated critical silicon abundances ([Si/H]{sub crit}). At the density and temperature associated with protostellar disks, we find that a standard Milky Way grain size distribution gives [Si/H]{sub crit} = –4.5 ± 0.1, while smaller grain sizes created in a supernova reverse shock give [Si/H]{sub crit} = –5.3 ± 0.1. Other environments are not dense enough to be influenced by dust cooling. We test the silicate dust cooling theory by comparing to silicon abundances observed in the most iron-poor stars ([Fe/H] < -4.0). Several stars have silicon abundances low enough to rule out dust-induced gas fragmentation with a standard grain size distribution. Moreover, two of these stars have such low silicon abundances that even dust with a shocked grain size distribution cannot explain their formation. Adding small amounts of carbon dust does not significantly change these conclusions. Additionally, we find that these stars exhibit either high carbon with low silicon abundances or the reverse. A silicate dust scenario thus suggests that the earliest low-mass star formation in the most metal-poor regime may have proceeded through two distinct cooling pathways: fine-structure line cooling and dust cooling. This naturally explains both the carbon-rich and carbon-normal stars at extremely low [Fe/H].

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Anelastic sensitivity kernels with parsimonious storage for adjoint tomography and full waveform inversion

We introduce a technique to compute exact anelastic sensitivity kernels in the time domain using parsimonious disk storage. The method is based on a reordering of the time loop of time-domain forward/adjoint wave propagation solvers combined with the use of a memory buffer. It avoids instabilities that occur when time-reversing dissipative wave propagation simulations. The total number of required time steps is unchanged compared to usual acoustic or elastic approaches. The cost is reduced by a factor of 4/3 compared to the case in which anelasticity is partially accounted for by accommodating the effects of physical dispersion. We validate our technique by performing a test in which we compare the Kα sensitivity kernel to the exact kernel obtained by saving the entire forward calculation. This benchmark confirms that our approach is also exact. We illustrate the importance of including full attenuation in the calculation of sensitivity kernels by showing significant differences with physical-dispersion-only kernels.

Geochemistry & Geophysics↗
DOE OSTI2017

The quinic acid derivative KZ-41 prevents glucose-induced caspase-3 activation in retinal endothelial cells through an IGF-1 receptor dependent mechanism

Retinal microaneurysms, an early disease manifestation of diabetic retinopathy, are associated with retinal endothelial cell (REC) death and macular edema. We previously demonstrated that a quinic acid (QA) analog, KZ-41, promoted REC survival by blunting stress-induced p38 MAPK activation. Herein, we sought to expand our understanding of the prosurvival signal transduction pathways actuated by KZ-41. Using human RECs exposed to high glucose (25 mM, 72 hours), we demonstrated that KZ-41 blocks caspase-3 activation by triggering phosphorylation of the PI3K regulatory subunit (p85; Tyr458) and its downstream target Akt (Ser473). Akt signal transduction was accompanied by autophosphorylation of the receptor tyrosine kinase, insulin growth factor-1 receptor (IGF-1R). IGF-1R knockdown using either the tyrosine kinase inhibitor AG1024 or silencing RNA abolished KZ-41’s pro-survival effect. Under high glucose stress, caspase-3 activation correlated with elevated ERK1/2 phosphorylation and decreased insulin receptor substrate-1 (IRS-1) levels. KZ-41 decreased ERK1/2 phosphorylation and reversed the glucose-dependent reduction in IRS-1. To gain insight into the mechanistic basis for IGF-1R activation by KZ-41, we used molecular modeling and docking simulations to explore a possible protein:ligand interaction between the IGF-1R kinase domain and KZ-41. Computational investigations suggest two possible KZ-41 binding sites within the kinase domain: a region with high homology to the insulin receptor contains one potential allosteric binding site, and another potential site on the other side of the kinase domain, near the hinge domain. These data, together with previous proof-of-concept efficacy studies demonstrating KZ-41 mitigates pathologic retinal neovascularization in the murine oxygen-induced retinopathy model, suggests that QA derivatives may offer therapeutic benefit in ischemic retinopathies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Microseismic Focal Mechanisms and Implications for Changes in Stress during the 2014 Newberry EGS Stimulation

We adapt the relative polarity method from Shelly et al. (2016) to compute focal mechanisms for microearthquakes associated with the 2014 hydroshearing stimulation at the Newberry volcano geothermal site. We focus the analysis on events relocated by Aguiar and Myers (2018), who report that six event clusters predominantly comprise the 2014 sequence. Data quality allows focal mechanism analysis for four of the six event clusters. We use Hardebeck and Shearer (2002, 2003; hereafter HASH) to compute focal mechanisms based on first–motion polarities and S/P amplitude ratios. We manually determine P– and S–wave polarities for a well–recorded reference event in each cluster, then use waveform cross correlation to determine whether recordings of other events in the cluster are the same or reversed polarity at each network station. Here, most waveform polarities are consistent with the affiliated reference event, indicating similar focal mechanisms within each cluster. The deeper clusters are east–west–striking normal faults, whereas the shallower clusters, close to the top of the open–hole section of the borehole, are strike slip with east–west motion. Regional studies and prestimulation borehole breakouts find the maximum stress direction is vertical and maximum horizontal stress is approximately north–south. Fault geometry and focal mechanisms of microseismicity during the stimulation suggest that increased pressure from fluid injection predominantly caused changes in horizontal stress, consistent with predictions from numerical studies of stress change caused by fluid injection. At shallow depths, where previous studies suggest the difference between vertical and horizontal stress is lowest, injection appears to have rotated the direction of maximum stress from vertical to horizontal, resulting in strike–slip motion. At greater depth, vertical stress continued to be the dominant direction during the stimulation, but fault orientation indicates either reactivation of pre–existing fractures or rotation of the direction of maximum horizontal stress from approximately north–south to east–west.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Liquid Crystal Orientation and Shape Optimization for the Active Response of Liquid Crystal Elastomers

Liquid crystal elastomers (LCEs) are responsive materials that can undergo large reversible deformations upon exposure to external stimuli, such as electrical and thermal fields. Controlling the alignment of their liquid crystals mesogens to achieve desired shape changes unlocks a new design paradigm that is unavailable when using traditional materials. While experimental measurements can provide valuable insights into their behavior, computational analysis is essential to exploit their full potential. Accurate simulation is not, however, the end goal; rather, it is the means to achieve their optimal design. Such design optimization problems are best solved with algorithms that require gradients, i.e., sensitivities, of the cost and constraint functions with respect to the design parameters, to efficiently traverse the design space. In this work, a nonlinear LCE model and adjoint sensitivity analysis are implemented in a scalable and flexible finite element-based open source framework and integrated into a gradient-based design optimization tool. To display the versatility of the computational framework, LCE design problems that optimize both the material, i.e., liquid crystal orientation, and structural shape to reach a target actuated shapes or maximize energy absorption are solved. Multiple parameterizations, customized to address fabrication limitations, are investigated in both 2D and 3D. The case studies are followed by a discussion on the simulation and design optimization hurdles, as well as potential avenues for improving the robustness of similar computational frameworks for applications of interest.

42 ENGINEERING↗