Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural simulation”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Structural uncertainty assessment for fire-engulfed objects in crosswind: Establishing credibility for a multiphysics wall-modeled large-eddy simulation paradigm

A structural uncertainty validation study for a large-scale, fire-engulfed, elevated object subjected to crosswind is presented to establish the credibility of a high-fidelity, low-Mach, turbulent reacting flow wall-modeled large-eddy simulation (WMLES) approach that includes multiphysics coupling to participating media radiation and conjugate heat transfer. To establish that WMLES can accurately predict surface quantities including drag and pressure coefficient in the low-Mach crosswind regime, a foundational elevated isothermal cylinder validation case is presented at a similar gap-to-diameter ratio of 0.25, spanning the subcritical to supercritical drag regime (Re 𝐷 = 1.1 × 10 5 and 4.3 × 10 5 , respectively). Here, this study exercised both static and dynamic coefficient LES (Smagorinsky and 𝑘 sgs ) with both local and exchange-based velocity sampling. Results showcase that the drag crisis (or the sudden drop in drag coefficient at increased Re 𝐷 ) is well captured when using an exchange-based dynamic coefficient WMLES methodology, while noting lack of mesh convergence and overall drag and pressure coefficient predictively when using a static coefficient, local velocity sampling WMLES. For the 𝒪⁡(10) m JP-8 liquid pool fire crosswind validation study presented, two experimental crosswind configurations (2 m/s and 9.5 m/s) are showcased for a fire-engulfed mock fuselage roughly 4 m in diameter. Using the best model-form practices identified in the isothermal study, dynamic coefficient 𝑘 sgs exchange-based WMLES fire validation findings demonstrate accurate peak irradiation and skin temperature predictions as a function of crosswind magnitude. Excessive yaw in the low-crosswind fuselage configuration, consistent with experimental findings, captured a significant predicted asymmetry in flame attachment and heat flux toward the downwind cylindrical cap—indicative of axial vortex structures transporting the flame along the upper and lower fuselage leeward surface. All fire mesh resolution simulations captured the experimental finding that as crosswind increased, predicted flame shape and peak irradiation magnitude onto the fuselage transitioned from a windward to a leeward cylinder location due to the migration of the upper- to lower-shear fuel/air mixing layer thereby demonstrating the novelty, significance, and credibility of this high-fidelity WMLES reacting flow framework.

Domino, Stefan Paul [Sandia National Laboratories ↗
DOE OSTI2021

SHarD: A beam dynamics simulation code for dielectric laser accelerators based on spatial harmonic field expansion

In order to demonstrate acceleration of electrons to relativistic scales by an on chip dielectric laser accelerator (DLA), a ponderomotive focusing scheme capable of capturing and transporting electrons through nanometer-scale apertures over extended interaction lengths has been proposed. Here we present a Matlab-based numerical code (SHarD) utilizing a spatial harmonic expansion of the fields within the dielectric structure to simulate the evolution of the beam phase space distribution in this scheme. The code can be used to optimize key-parameters for the accelerator performance such as the final energy, transverse spot size evolution and total number of electrons accelerated through currently fabricated structures. Eventually, the simulation model will be applied to inform the phase mask profile to be added to a pulse front tilt drive laser pulse using a liquid crystal mask in the experimental setup being assembled at UCLA Pegasus Laboratory.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Group structure selection with random forests

Choosing an appropriate group structure for multigroup transport is far from an exact science. For some applications, one blindly uses a group structure developed years ago by forgotten methods. Furthermore, one sometimes uses the same group structure for a variety of problems, even if the group structure was originally developed with a certain application in mind. In this work, we create optimized group structures with simulated annealing for critical assembly test problems and apply a random forest regressor with bagging to choose the best group structure based on parameters of the different test problems. The optimized group structures were generated using a simulated annealing optimizer for several simple, spherical, and unreflected problems. The optimization was performed to minimize a cost function that included fission rate, absorption rate, leakage, and k{sub eff}. A random forest regressor was then trained on a set of International Criticality Safety Benchmark Evaluation Project inputs and used to select one of these six group structures. The trained machine learning model chose the best group structure 65% of the time, and one of the three best 89% of the time. Furthermore, it decreased the L2 error over all test problems by a factor of 25 when compared to the standard Los Alamos 70-group structure. In other words, the model chose group structure that were far more appropriate for the test problems than the traditional LANL group structure. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Machine Learning assisted optimization and parameter space exploration dataset of spin ice Hamiltonian

This repository contains both simulated and experimental structure factor data for the data challenge involving the inverse scattering problem. The simulated data were generated during a machine-learning-assisted optimization routine described in ref[1]. The experimental structure factor was measured on a rare-earth oxide, Dy2Ti2O7 using diffuse neutron scattering from time-of-flight techniques on the CORELLI instrument at the Spallation Neutron Source, Oak Ridge National Laboratory. A Metropolis Monte Carlo code implemented to run in a High-performance computing setting was used to calculated simulated structure factors for the spin-ice Hamiltonian at 680 mK, which is the same temperature as for the experimental data. The total size of all the files in this repository is 5.12 GB. A detailed description of the files is given below. ExperimentalData_630mK.dat – A linearized version of 3-dimensional experimental data of size 61×81×21. This data was processed to remove an estimation of non-magnetic background, including nuclear scattering signal and instrumentation background. Parameters.dat – 6700 samples were evaluated over the 4-dimensional parameter space (J_1, J_2, J_3 and J_(3^' )). There is an additional parameter, D in the spin Hamiltonian to mimic the dipolar interaction between magnetic ions. However, this parameter, D was fixed to a value determined by prior work. This file contains five columns for the parameters J_1, J_2, J_3, J_(3^' ) and D respectively. 3D_Simulation_Data.dat – The simulated structure factor, S(Q) data are included in this file. Each raw contains a linearized array of 3D volumes of S(Q) calculated for the parameter set given in the corresponding row of the file Parameters.dat. The size of the volume data was matched to the experimental data. Qx(h,-h,0).dat, Qy(k,k,-2k).dat, Qz(l,l,l).dat – These files contain the h, k, and l values along with the reciprocal vectors [h,-h,0], [k,k,-2k] and [l,l,l] respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Structural Origin of Chiroptical Properties in Perovskite Nanocrystals with Chiral Organic Ligands

Abstract The authors investigate how chiral ligands attached to perovskite nanocrystal (PNC) surfaces structurally distort the perovskite lattice. Chiral electro‐optical properties of the resulting PNCs are demonstrated through the fabrication of a circularly polarized light (CPL) detector with a discrimination of up to 14% between left‐ and right‐handed CPL. Both experimental and electronic‐structure‐based simulations are combined to provide insights into the interactions (both structural and electronic) between chiral organic ligands and PNCs. The major finding is a centro‐asymmetric distortion of the surface lattice that penetrates up to five atomic unit cells deep into the PNCs, which is the likely cause of the chiral‐optical properties. Spin‐polarized transport through chiral‐PNCs results from the chiral‐induced spin selectivity effect and amplifies the discrimination between left and right‐handed CPL as is experimentally demonstrated in the detectors.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Automated Transformation of Lignin Topologies into Atomic Structures with LigninBuilder

Lignin is an abundant aromatic heteropolymer found in secondary plant cell walls and is a potential feedstock for conversion into bioderived fuels and chemicals. Lignin chemical diversity complicates traditional structural studies, and so, relatively little experimental evidence exists for how lignin structure exists in aqueous solution or how lignin polymers respond to changes in their chemical environment. Molecular modeling can address these concerns; however, prior computational structural lignin models typically did not capture lignin heterogeneity, as only a few polymers were considered. LigninBuilder creates a framework for building structural libraries for lignin from existing topological libraries, permitting significantly greater diversity of lignin structures to be sampled at atomic detail. As a demonstration of its capabilities, LigninBuilder was applied to three libraries of lignin from hardwood, softwood, and grass, and the resulting polymer structures were simulated in an aqueous environment. The lignins adopted compact globular structures, as would be expected for polymers in poor solvents. The differences between the libraries were largest when quantifying the packing of nonadjacent aromatic residues, with greater branching within the polymer resulting in poorer aromatic packing. Individual lignin polymers were also found to undergo rapid conformational changes, with the dwell time within a state growing as the square of the molecular weight. This first application of LigninBuilder demonstrates the potential for atomic-level insight into lignin interactions. LigninBuilder is distributed as a plugin to the visualization software VMD, lowering the barrier for modeling lignin structure in diverse environments.

09 BIOMASS FUELS↗
DOE OSTI2022

Quantum Simulation of Molecular Electronic States with a Transcorrelated Hamiltonian: Higher Accuracy with Fewer Qubits

Simulation of electronic structure is one of the most promising applications on noisy intermediate-scale quantum (NISQ) era devices. However, NISQ devices suffer from a number of challenges like limited qubit connectivity, short coherence times, and sizable gate error rates. Thus, desired quantum algorithms should require shallow circuit depths and low qubit counts to take advantage of these devices. Here, we attempt to reduce quantum resource requirements for molecular simulations on a quantum computer while maintaining the desired accuracy with the help of classical quantum chemical theories of canonical transformation and explicit correlation. In this work, compact ab initio Hamiltonians are generated classically, in the second quantized form, through an approximate similarity transformation of the Hamiltonian with (a) an explicitly correlated two-body unitary operator with generalized pair excitations that remove the Coulombic electron–electron singularities from the Hamiltonian and (b) a unitary one-body operator to efficiently capture the orbital relaxation effects required for accurate description of the excited states. The resulting transcorrelated Hamiltonians are able to describe both the ground and the excited states of molecular systems in a balanced manner. Using the variational quantum eigensolver (VQE) method based on the unitary coupled cluster with singles and doubles (UCCSD) ansatz and only a minimal basis set (ANO-RCC-MB), we demonstrate that the transcorrelated Hamiltonians can produce ground state energies comparable to the reference CCSD energies with the much larger cc-pVTZ basis set. This leads to a reduction in the number of required CNOT gates by more than 3 orders of magnitude for the chemical species studied in this work. Furthermore, using the quantum equation of motion (qEOM) formalism in conjunction with the transcorrelated Hamiltonian, we are able to reduce the deviations in the excitation energies from the reference EOM-CCSD/cc-pVTZ values by an order of magnitude. In conclusion, the transcorrelated Hamiltonians developed here are Hermitian and contain only one- and two-body interaction terms and thus can be easily combined with any quantum algorithm for accurate electronic structure simulations.

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

Hydroxyl radical yields in the heavy ion radiolysis of water

Scavenging experiments coupled with Monte Carlo track structure and kinetics simulations have been used to determining and understand the yields of hydroxyl radicals in the radiolysis of water with protons and carbon ions. Carbon dioxide production from aerated formic acid solutions with concentrations ranging from 10-3 to 1 mol./dm3 was used as a probe of the time dependence of the hydroxyl radical yield from about 5 ns to 1 ?s. Numerical extrapolation of the results at slightly different ion energies in combination with data from gamma radiolysis experiments allows a systematic analysis of both track average and track segment yields, with considerable track chemistry found to occur on the nanosecond to microsecond time scale. Track chemistry calculations utilizing Monte Carlo track structure simulations and nonhomogeneous independent reaction times diffusion-reaction modeling of product yields reproduce experimental observations. Subsequently, the track simulation methods were used to predict hydroxyl radical kinetics in deaerated and aerated water at neutral pH. Taken with literature data for helium ion beam studies, the results offer insight into the potential for indirect radiation induced damage to biological systems for light ions with very different track structures.

38 - RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCL↗
DOE OSTI2022

A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation

Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.

36 MATERIALS SCIENCE↗
DOE OSTI2019

The simulation on absorption properties of metamaterial/GaAs/electrode layer hybrid structure based Terahertz photoconductive detector

Terahertz (THz) photodetectors have attracted great attention from scientists worldwide for their application in security checking, biomedical treatment and astronomical observation of remote stars and distant galaxies. As a typical THz detector, extrinsic GaAs based photoconductive detector is facing critical technical bottlenecks in the epitaxial growth of sufficiently thick and high-quality GaAs absorption layer. In this work, a novel THz photoconductive detector based on metamaterial/GaAs/electrode layer hybrid structure was designed and simulated. By setting the periodic split ring resonator (SRR) structure as 88 μm pitch with 8 μm width, the absorption peaks exist at the wavelength of about 142 and 367 μm, which originate from the resonant cavity and the SRR dipole resonance effect, and the novel device shows a significant enhancement compared with the conventional GaAs photoconductive detector. Thus, the necessary thickness of GaAs absorption-layer is largely reduced, and the resonant absorption peak can be modulated by changing the thickness of absorption layer. This work provides a novel device structure which can solve the critical epitaxial growth bottleneck of GaAs photoconductive detector and used for the astronomical observation, security check, etc.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Cooperative Education

Los Alamos National Laboratory (LANL) is a multidisciplinary national laboratory that conducts research and development in national security, engineering, materials science, computational modeling, and advanced manufacturing. The laboratory develops innovative technologies to address complex scientific and engineering challenges. This project focuses on the development and evaluation of high-performance absorbing structures through computational design, simulation, and engineering analysis. Absorbing structures are used in applications where damage mitigation, structural protection, and material efficiency are critical performance requirements. The increasing demand for lightweight, high-strength, and highly efficient structural systems has created a need for improved design methodologies capable of maximizing absorption while minimizing weight and material usage. The project utilizes advanced engineering software, including 3D CAD software and FEA, to generate and optimize structural concepts. Computational simulations are performed to evaluate structural behavior under loading conditions, while mathematical analyses are conducted using Python-based tools as well as established analytical equations from material and structural mechanics. The project benefits LANL by supporting the development of advanced design methodologies and improving the understanding of material and structural performance. During the internship term, a significant portion of the design development, simulation, and data analysis activities will be completed. Success of the project depends on collaboration among engineering mentors and technical staff members. Work will be conducted at Los Alamos National Laboratory using laboratory computing resources and engineering software.

42 ENGINEERING↗
DOE OSTI2019

Pair distribution function analysis of neutron-irradiated silicon carbide

We have employed x-ray total scattering to investigate the structure of polycrystalline 3C-silicon carbide following neutron irradiation. Here, the structure as a function of irradiation temperature and dose was quantified by analyzing pair distribution functions. Although the SiC matrix retains its crystal structure after irradiation, a significant increase in the diffuse scattering component is observable indicating that neutron irradiation leads to changes in both the short- and medium-range order. These changes include both an irradiation dose- and temperature-dependent increase in the vacancy concentration leading to an increase in the Si and C atomic displacement parameters. A dose-dependent decrease in the size of defect free material is also quantified from the structural refinements due to an increase in the number of defects and defect clusters. Evidence of additional correlations in the short-range order (up to ~4 Å) from differential pair distribution function analysis indicate that combinations of atomistic defects including anti-site defects, vacancies and defect clusters are present after these irradiation conditions. Such structural information will be valuable for direct comparison of experimental and simulated atomic structures of irradiated silicon carbide.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Valence shell photoelectron angular distributions and vibrationally resolved spectra of imidazole: A combined experimental–theoretical study

Linearly polarized synchrotron radiation has been used to record polarization dependent valence shell photoelectron spectra of imidazole in the photon energy range 21–100 eV. These have allowed the photoelectron angular distributions, as characterized by the anisotropy parameter β, and the electronic state intensity branching ratios to be determined. Complementing these experimental data, theoretical photoionization partial cross sections and β-parameters have been calculated for the outer valence shell orbitals. The assignment of the structure appearing in the experimental photoelectron spectra has been guided by vertical ionization energies and spectral intensities calculated by various theoretical methods that incorporate electron correlation and orbital relaxation. Strong orbital relaxation effects have been found for the 15a', nitrogen lone-pair orbital. The calculations also predict that configuration mixing leads to the formation of several low-lying satellite states. The vibrational structure associated with ionization out of a particular orbital has been simulated within the Franck–Condon model using harmonic vibrational modes. The adiabatic approximation appears to be valid for the X 2A" state, with the β-parameter for this state being independent of the level of vibrational excitation. However, for all the other outer valence ionic states, a disparity occurs between the observed and the simulated vibrational structure, and the measured β-parameters are at variance with the behavior expected at the level of the Franck–Condon approximation. Furthermore, these inconsistencies suggest that the excited electronic states may be interacting vibronically such that the nuclear dynamics occur over coupled potential energy surfaces.

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

Predicting initial dissolution rates using structural features from molecular dynamics simulations

Predicting chemical durability of glass materials is important for various applications from daily life such as drink glass and kitchen ware to advanced technologies such as nuclear waste disposal and biomedicine. In this work, we explored prediction of initial dissolution rate through structural features from molecular dynamics (MD) simulations for a wide range of glass compositions (total 28) including borosilicate and aluminosilicate glasses, ZrO 2 -containing and V 2 O 5 -containing boroaluminosilicate glasses. The initial dissolution rates (r 0 ) measured experimentally at 90 °C with varying solution conditions were correlated with structural features (e.g., polyhedral linkages and non-bridging oxygen species) obtained from MD simulations, either from this study or from literature. Since hydrolysis of the glass network through breaking of the network former linkages (e.g., Si-O-Si, Si-O-Al, etc.) is a critical step of network glass dissolution, the statistics of these linkages obtained from MD were correlated to r0 through linear regression, where the coefficient of determination (R 2 ) and root mean square error are found to be 0.949 and 0.681, respectively. This model was compared and discussed with existing models developed by various approaches including machine learning, the kinetic rate equation, topological constraint theory, and other descriptors from MD simulations. The discussion provides insights on future model improvements to predict glass dissolution. In addition, the impact of V 2 O 5 on the glass dissolution was examined in detail, implicating that the impact is not the same across all glass compositions and test conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Structure and Vibrational Properties of Potassium-Promoted Tungsten Oxide Catalyst Monomeric Sites Supported on Alumina (K 2 O/WO 3 /Al 2 O 3 ) Characterized Using Periodic Density Functional Theory

Al 2 O 3 -supported tungsten oxide catalysts have been instrumental in many industrially relevant reactions and their reactivity is controlled by their molecular structure. In turn, their molecular structure has primarily been derived via Raman measurements with assignments made using model compounds of known local (molecular) coordination. Here, the structure and simulated Raman spectra of unpromoted and K + -promoted tungsten oxide catalyst monomeric sites supported on γ-Al 2 O 3 (110) (K 2 O/WO x /γ-Al 2 O 3 (110)) were studied using periodic DFT methods. Two different WO x -grafted monomers on the γ-Al 2 O 3 (110) surface were identified with a total energy difference of 0.17 eV between both structures. Importantly, both structures showed the presence of W═O and W–OH moieties, thus providing additional insights into experimental Raman data, which typically describe only W═O moieties. The grafted WO 3 species were stabilized when the present W–OH groups were oriented toward the alumina surface leading to the formation of H-bonds, calculated at 1039 cm –1 , for example in the vicinity of the W═O vibrations. The W═O bond length was altered in the presence of K + , as shown experimentally and theoretically in blue shifting of the Raman band corresponding to W═O. The W vibrations were well localized in the calculated spectra, and little shifts were observed upon the different WO x molecular geometry, explaining why a single Raman peak is mostly observed experimentally. The acidity of the lowest energy catalyst structures was investigated by simulated NH 3 adsorption vibrational frequency and binding energy calculations. Results suggested that NH 3 prefers to bind in a Lewis-like structure with no proton donation from either W–OH or Al–OH moiety. This challenged some of the literature observations where Brønsted acid sites have been suggested to exist on near-monolayer coverage WO x species on γ-Al 2 O 3 prepared by calcining at moderate (400 °C) to high temperatures (700 °C). Finally, this work provided new insights into the molecular structure of WO x /γ-Al 2 O 3 and K 2 O/WO x /γ-Al 2 O 3 catalysts not immediately available from experimental measurements alone.

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

DFT Accurate Interatomic Potential for Molten NaCl from Machine Learning

Molten alkali chloride salts are a critical component in concentrated solar power and nuclear applications. Despite their ubiquity, the extreme chemical reactivity of molten alkali chlorides at high temperatures has presented a significant challenge in characterizing atomic structures and dynamic properties experimentally. In this work, we we investigate molten NaCl by performing high temperature molecular dynamics simulations using a Gaussian Approximation Potential (GAP) trained on Density Functional Theory (DFT) datasets. Our GAP model, trained with a meager 1000 atomic configurations, arrives at near ab initio accuracy with a mean absolute error of 1.5 meV/atom thus enabling fast analysis of high temperature salt properties on large length (5000 ion pairs) and time (> 1ns) scales currently inaccessible to ab initio simulations. Calculated structure factors and diffusion constants from our GAP model simulations show excellent agreement with experiments. Our results indicate that GAP models are able to capture the many-body interactions required to accurately model ionic-systems.

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

Structural Basis for Binding of Allosteric Drug Leads in the Adenosine A 1 Receptor

Despite intense interest in designing positive allosteric modulators (PAMs) as selective drugs of the adenosine A 1 receptor (A 1 AR), structural binding modes of the receptor PAMs remain unknown. Using the first X-ray structure of the A 1 AR, we have performed all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) technique to determine binding modes of the A 1 AR allosteric drug leads. Two prototypical PAMs, PD81723 and VCP171, were selected. Each PAM was initially placed at least 20 Å away from the receptor. Extensive GaMD simulations using the AMBER and NAMD simulation packages at different acceleration levels captured spontaneous binding of PAMs to the A 1 AR. The simulations allowed us to identify low-energy binding modes of the PAMs at an allosteric site formed by the receptor extracellular loop 2 (ECL2), which are highly consistent with mutagenesis experimental data. Furthermore, the PAMs stabilized agonist binding in the receptor. In the absence of PAMs at the ECL2 allosteric site, the agonist sampled a significantly larger conformational space and even dissociated from the A 1 AR alone. In summary, the GaMD simulations elucidated structural binding modes of the PAMs and provided important insights into allostery in the A 1 AR, which will greatly facilitate the receptor structure-based drug design.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Modeling of Charge Amplification in LArTPCs

Single-phase liquid argon time projection chambers (LArTPCs) are an important part of Fermilab's neutrino experimental program. Current LArTPCs are able to detect interactions producing energy deposits above the MeV scale, and focus on neutrinos at GeV energies. The Liquid Argon Charge Amplification Device, or LArCADe, program aims to employ modified anode geometries in order to obtain electron avalanches that would help reduce detection thresholds, opening the door to new physics opportunities at the keV scale (e.g. CEvNS and dark-matter WIMP searches). To model and quantify this phenomenon, a simulation of electron transport in LAr was built, able to simulate excitations and ionizations in non-uniform fields. Using COMSOL simulations of the anode structures, the built simulation produced quantitative and realistic estimates of LAr charge amplification near different tip geometries. These results will inform the next generation of tips produced for LArCADe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗