Search NASASearch

SEARCH · Search NASA

Results for “simulation software”

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.

510 records · Page 7

Natural Fiber Simulations and Analysis for Composites

As the U.S. government and the aviation industry strive to achieve net-zero carbon emissions by 2050, the Sustainable Manufacturing of Aircraft (SUMAC) project, supported by the Convergent Aeronautics Solutions (CAS) program of the National Aeronautics and Space Administration (NASA), aims to develop sustainably derived thermoplastic composites and manufacturing technologies for future aviation. A multidisciplinary team of experts drives this initiative through a collective effort. The team’s key focus areas are sustainably derived thermoplastic resins, natural fiber reinforcement, composite manufacturing, structural health monitoring, computational materials modeling, and systems analysis. This presentation will give an overview of the SUMAC project. Fiber-scale microscopy analysis characterized fiber cross sections finding that hemp and flax fibers are very deformable and non-circular (circularities of 0.2 were observed), and diameters range from submicron to 20 μm. Fiber shapes range from flattened, with aspect ratios Rlong/Rshort up to 5, to round, Rlong/Rshort =1. And Fiber core perimeters, which likely drives deformability, were calculated Pcore=7 μm. Explicit-fiber yarn simulations, parameterized from single-fiber measurements, resolved and identified slack as a cause for yarn strain-hardening tensile behavior. The bonded network model of PLA, PHA and PA11 resins captures tensile behavior, and can model an infused weave. Computational models of the weave unit cells matched the pattern (including plain, twill and triaxial braid), yarn ply (1to 2), and yarn twist from optical microscopy.

thermoplastic

Natural Fiber Simulations and Analysis for Composites

As the U.S. government and the aviation industry strive to achieve net-zero carbon emissions by 2050, the Sustainable Manufacturing of Aircraft (SUMAC) project, supported by the Convergent Aeronautics Solutions (CAS) program of the National Aeronautics and Space Administration (NASA), aims to develop sustainably derived thermoplastic composites and manufacturing technologies for future aviation. A multidisciplinary team of experts drives this initiative through a collective effort. The team’s key focus areas are sustainably derived thermoplastic resins, natural fiber reinforcement, composite manufacturing, structural health monitoring, computational materials modeling, and systems analysis. This presentation will give an overview of the SUMAC project. Fiber-scale microscopy analysis characterized fiber cross sections finding that hemp and flax fibers are very deformable and non-circular (circularities of 0.2 were observed), and diameters range from submicron to 20 μm. Fiber shapes range from flattened, with aspect ratios Rlong/Rshort up to 5, to round, Rlong/Rshort =1. And Fiber core perimeters, which likely drives deformability, were calculated Pcore=7 μm. Explicit-fiber yarn simulations, parameterized from single-fiber measurements, resolved and identified slack as a cause for yarn strain-hardening tensile behavior. The bonded network model of PLA, PHA and PA11 resins captures tensile behavior, and can model an infused weave. Computational models of the weave unit cells matched the pattern (including plain, twill and triaxial braid), yarn ply (1to 2), and yarn twist from optical microscopy.

thermoplastic

Simulation to a Newborn Supernova Remnant from a Low-mass Iron Core Star

Supernova remnant observations show a high degree of asymmetry, mixing, and inhomogeneity. These asymmetries are seeded during the early seconds of the explosion and are further enhanced and modified as the shock and ejecta move through the stellar progenitor and into the circumstellar medium. We present simulations of a 9.6 M⊙ zero-metallicity progenitor initialized after shock revival and evolved for several years when the ejecta is in the circumstellar medium. A suite of 1D and 2D simulations examines the effects of neutron-star wind and radioactive decay heating. In 1D, decay heating forms a low-density bubble that suppresses the reverse shock. While in 2D, the heating is localized to metal-rich pockets, inflating them and compressing the surrounding material into dense shells. In 3D, the neutron-star wind and decay heating modify the plume morphology, producing more large-scale structures. The extended plume morphology leads to an asymmetrical shock breakout. After breakout, the leading plumes cannot keep up with the shock front, resulting in deceleration and fragmentation by the reverse shock while retaining the large-scale asymmetry. The projected ejecta morphology and velocities are strongly viewing angle dependent. The relatively uniform metal-rich distribution does not resemble the strongly inhomogeneous ejecta structure of Cas A. The 160-isotope decay network shows that 24.4% of the radioactive heating comes from decay chains other than the canonical 56Ni chain. The low explosion energy, low 56Ni yield, and Ni/Fe ratio greater than unity suggest an observational signature similar to an electron capture supernova.

Neopane, Sudarshan [University of Tennessee (UT)]

Unveiling the lithium-ion transport mechanism in Li{sub 2}ZrCl{sub 6} Solid-State Electrolyte {ital via} deep learning-accelerated molecular dynamics simulations.

Lithium zirconium chlorides (LZCs) present a promising class of cost-effective solid electrolytes for next-generation all-solid-state batteries. The unique crystal structure of LZCs plays a crucial role in facilitating lithium-ion mobility, which further affects the electrochemical performance. To understand the underlying mechanism governing ion transport, we employed deep learning-accelerated molecular dynamics simulation on Li2ZrCl6 (trigonal alpha- and monoclinic beta-LZC), focusing specifically on the zirconium coordination environment. Our results reveal that disordered alpha-LZC exhibits the highest ionic conductivity, while beta-LZC demonstrates significantly lower conductivity, closely aligning with experimental findings. The study confirms that across all phases, lithium migration proceeds via the site-to-site hopping mechanism, where variations in site residence times critically impact the overall ionic conductivity. In alpha-LZCs, lithium ions prefer to anisotropically diffuse across interlayers as the result of a lower energy barrier, driven primarily by collective diffusion. In contrast, lithium ions in beta-LZC primarily isotropically diffuse within the intralayer, hindered by higher energy barriers and determined by individual diffusion. The variation in ZrCl6 2- octahedral unit softening, induced by the specific layered arrangement of zirconium atoms, emerges as a critical determinant of the energy barriers across the LZC phases. These atomic-scale insights into the transport processes provide valuable guidance for the rational design and optimization of LZCs-based electrolytes, accelerating their practical application in advanced energy storage technologies.

Guo, Hanzeng

Spacecraft Fire Safety Predictions using Verified Saffire Model

A model developed using Fire Dynamics Simulator (FDS) that aimed to determine the effect of a fire in a spacecraft was validated by data collected during the Saffire campaign. The model used inlet and outlet temperatures and CO 2 concentrations of the Saffire payload where fire spread was taking place to determine the amount of heat and combustion products that made it into Northrop Grumman’s Cygnus vehicle. The model was then validated using six remote sensors in various places in the vehicle, as well as a far field device (FFD) in the open zenith section that was representative of average vehicle values. The current work focuses on using the model to predict fire safety scenarios. One simulation aimed to determine the fate of HCl, which sticks to surfaces. The model prediction showed that the HCl was removed from the atmosphere rapidly. This compared well against the FFD data in the Saffire VI campaign event where a bottle of 5% HCl was released into the vehicle. An additional simulation where the Environmental Control and Life Support System (ECLSS) was shut off once the FFD reached 5 ppm of HCl showed that HCl stayed in the atmosphere considerably longer. Continuing the simulation with the ECLSS activated and after temperatures returned to their initial conditions, resulted in a rapid removal of HCL similar to that observed in the original HCl release scenario model. Finally, a simulation that used the heat release rate from a lithium-ion battery test to determine the effect it would have on a spacecraft was performed. This simulation used a heat addition rate that is considerably higher than what was determined from the burning of solid fuels in the Saffire campaign and hence produced a non-trivial temperature increase in more locations within the vehicle.

Fire Safety

Simulated Lunar Gravity Testing of VIPER Loop Heat Pipe

NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) mission consists of a rover designed to explore the lunar south pole. One of the main challenges faced by the rover during the lunar polar exploration is the adverse thermal environment. Temperatures can fluctuate more than a 100°C between day and night; potentially dropping to −246°C in permanent shadow regions. To maintain the rover components within temperature limits, VIPER’s Thermal Management System (TMS) relies heavily on Loop Heat Pipes (LHPs). To assist the design of the thermal management system, an Engineering Design Unit (EDU) LHP has been tested in several opportunities under thermal vacuum (TVAC) environment. Of particular interest was the LHP performance in lunar gravity. To that end, the EDU LHP was tested at the nominal orientation within the rover assembly, vertical, and inclined such that the gravitational component acting on the LHP evaporator and condenser was 1/6g (lunar gravity is 1/6 of Earth gravity) when compared to the vertical configuration. This paper will examine the LHP conductance as one the key parameters to assess the effect of the gravity on the LHP performance.

Loop Heat Pipe

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics

A comparison of wrought and powder metallurgical FeCrAl claddings under simulated LWR accident transients

Iron-chromium-aluminum (FeCrAl) alloys are potential accident tolerant fuel (ATF) cladding candidates for light-water reactors but are difficult to fabricate as thin-walled tubes via conventional cast-and-wrought routes. Powder metallurgy (PM) offers a manufacturing alternative with improved compositional control, but its transient accident performance has not been directly benchmarked against wrought variants. This study evaluates the burst behavior of commercially developed PM-processed FeCrAl alloys, PM-C26M (Fe-12Cr-6Al-2Mo) and the high precipitate density FA-SMT (Fe-22Cr-5Al-3Mo), under simulated light-water reactor accident transient conditions. Burst testing was conducted using the Severe Accident Test Station with heating rates of 5 °C/s and 50 °C/s and internal pressures ranging from 25 MPa to 100 MPa. PM-C26M reproduced wrought C26M burst behavior within 7–37 °C across the stress range, indicating that PM processing does not compromise transient strength. FA-SMT exhibited markedly higher burst temperatures and reduced heating-rate sensitivity, consistent with its engineered precipitate strengthening. FA-SMT rupture exhibited axial "unzipping" rather than the lateral tearing characteristic of PM- and wrought C26M. Post-test EBSD and fractography indicate that this behavior is strongly correlated with strain-gated intergranular void nucleation associated with the dense precipitate architecture of FA-SMT, a response absent in the comparatively clean PM-C26M matrix and consistent with rupture morphologies reported for oxide-dispersion strengthened (ODS) FeCrAl of similar base-matrix chemistry to PM-C26M. These findings highlight the potential of powder metallurgy as a viable fabrication route for ATF claddings from an accident performance standpoint.

Bell, Sam [ORNL] (ORCID:0000000251905657)

Sintering effects on NpO 2 grain size and morphology: The role of precursor Np phase

Neptunium dioxide (NpO 2 ) is a key phase in nuclear material processing as a target material for the production of plutonium-238 ( 238 Pu) and has been historically synthesized via the calcination of a Np oxalate precursor. Alternative synthesis methods for NpO 2 are now more prevalent, necessitating their study and comparison with the more common oxalate route. The purpose of this work was to investigate the microstructural properties of NpO 2 synthesized from a nitrate-based Np precursor phase via the assessment of NpO 2 particle size and morphology as a function of calcination temperature and time. Scanning electron microscopy (SEM) was used to probe the primary grain size and morphology of NpO 2 after calcination at temperatures ranging from 700 to 1100 °C and hold times ranging from 1 to 10 h. Post-image analysis using ImageJ software enabled the quantification of mean particle diameter. This analysis indicated that particle diameter increases with both increasing calcination temperature and hold time. Primary particles were shown to be clumped in irregular patterns into the overall rough, blocky aggregates, but this macroscopic morphology was not affected by calcination time or temperature. Although trends in primary grain size of NpO 2 were consistent with available literature from other Np precursor phases, the macroscopic morphology of the NpO 2 aggregates was quite different than reported for other precursors. Through comparison with historical literature on Np oxalate, this work emphasizes the importance of Np precursor on the physical properties of NpO 2 .

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA

Cryogenic Spray Quenching of A Simulated Propellant Storage Tank Wall With Heat Transfer Enhancement By A Thin-Film Coating and Flow Pulsing in Microgravity

Human space exploration to the Moon, Mars, and possibly asteroids is NASA’s biggest challenge for the new millennium. One of the critical elements to this mission is the effective, sufficient, and reliable supply of cryogenic propellant fluids. Future lower-earth-orbiting (LEO) propellant fuel depots and human-carrying orbital transfer spacecraft flying to the moon and Mars will have to utilize the high thrust and high efficiency of liquid cryogenic chemical propulsion or nuclear thermal propulsion. Efficient in-space tank-to-tank propellant transfer (propellant fuel depot to orbital transfer spacecraft) of cryogenic propellants is an enabling technology for the planned Crewed Mars Surface Mission. The transfer of cryogenic propellants in space, however, has yet to be accomplished, solely due to the unavailability of cryogenic quenching heat transfer data during chilldown (quenching) and filling of the propellant receiver tank in reduced gravity and microgravity as liquid propellant cannot be stored in a required liquid state until the tank is quenched down to the liquid temperature. Therefore, highly energy efficient thermal-fluid management breakthrough concepts to conserve and minimize the cryogen consumption during propellant transfer have become the focus of research and engineering development, especially for the deep-space mission to Mars. In this paper, we introduce such concepts and demonstrate their feasibility for cryogenic storage tank chilldown in parabolic flights under a simulated space microgravity condition. In order to maximize the storage tank chilldown efficiency for the least amount of cryogen consumption, the technology adopted included cryogenic spray cooling, Teflon thin-film coating of the simulated tank surface, and spray flow pulsing. The completed flight experiments successfully demonstrated that spray cooling is the most efficient cooling method for the tank chilldown in microgravity. In microgravity, Teflon coating alone can improve the efficiency up to 72% and the efficiency can be improved up to 59% by flow pulsing alone. However, Teflon coating together with flow pulsing was found to substantially enhance the chilldown efficiency in microgravity for up to 113%.

spray

Ice in Convective Storm Environments: Assessing Impacts of Microphysics on Simulated Reflectivities

Satellite radar retrievals of clouds and precipitation rely on assumptions about the microphysical properties of hydrometeors such as particle size distributions (PSDs) and massdimensional relationships, which have been shown to vary in distinct environments. A deeper understanding of linkages between mass-dimensional relationships, ice crystal shape, and radar reflectivities across cloud and precipitation regimes is necessary for developing new retrievals for upcoming satellite radar missions, including the INvestigation of Convective UpdraftS (INCUS), as well as future missions that will feature space-borne radars observing clouds, convection, and precipitation. This study investigates ice microphysical characteristics through aircraft in-situ and remote sensing observations from two convective events during the Mid-latitude Continental Convective Clouds Experiment (MC3E) in the Southern Great Plains, and tests implications for ice particle habit in forward modeled radar observations. During MC3E, the Ka/Ku- band High Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP) was onboard the NASA ER-2 aircraft, while in-situ measurements were taken onboard the University of North Dakota Citation. In-situ PSDs and particle shape information from optical array probe data are used to forward model radar reflectivities from particle scattering databases and are compared against HIWRAP reflectivities. Various particle combinations can be used to find agreement with HIWRAP observations, and the forward modeled reflectivities were sensitive to particle type. Varying the riming on the aggregate and graupel particles resulted in larger than 10 dB differences in reflectivities. The in-situ and radar observations imply the occurrence of aggregation, size sorting and lofting particles.

Julia A. Shates

WRF-Comfort: simulating microscale variability in outdoor heat stress at the city scale with a mesoscale model

Abstract. Urban overheating and its ongoing exacerbation due to global warming and urban development lead to increased exposure to urban heat and increased thermal discomfort and heat stress. To quantify thermal stress, specific indices have been proposed that depend on air temperature, mean radiant temperature (MRT), wind speed, and relative humidity. While temperature and humidity vary on scales of hundreds of meters, MRT and wind speed are strongly affected by individual buildings and trees and vary on the meter scale. Therefore, most numerical thermal comfort studies apply microscale models to limited spatial domains (commonly representing urban neighborhoods with building blocks) with resolutions on the order of 1 m and a few hours of simulation. This prevents the analysis of the impact of city-scale adaptation and/or mitigation strategies on thermal stress and comfort. To solve this problem, we develop a methodology to estimate thermal stress indicators and their subgrid variability in mesoscale models – here applied to the multilayer urban canopy parameterization BEP-BEM within the Weather Research and Forecasting (WRF) model. The new scheme (consisting of three main steps) can readily assess intra-neighborhood-scale heat stress distributions across whole cities and for timescales of minutes to years. The first key component of the approach is the estimation of MRT in several locations within streets for different street orientations. Second, mean wind speed and its subgrid variability are downscaled as a function of the local urban morphology based on relations derived from a set of microscale LES and RANS simulations across a wide range of realistic and idealized urban morphologies. Lastly, we compute the distributions of two thermal stress indices for each grid square, combining all the subgrid values of MRT, wind speed, air temperature, and absolute humidity. From these distributions, we quantify the high and low tails of the heat stress distribution in each grid square across the city, representing the thermal diversity experienced in street canyons. In this contribution, we present the core methodology as well as simulation results for Madrid (Spain), which illustrate strong differences between heat stress indices and common heat metrics like air or surface temperature both across the city and over the diurnal cycle.

Geology

Thermalization and criticality on an analogue–digital quantum simulator

Abstract Understanding how interacting particles approach thermal equilibrium is a major challenge of quantum simulators 1,2 . Unlocking the full potential of such systems towards this goal requires flexible initial state preparation, precise time evolution and extensive probes for final state characterization. Here we present a quantum simulator comprising 69 superconducting qubits that supports both universal quantum gates and high-fidelity analogue evolution, with performance beyond the reach of classical simulation in cross-entropy benchmarking experiments. This hybrid platform features more versatile measurement capabilities compared with analogue-only simulators, which we leverage here to reveal a coarsening-induced breakdown of Kibble–Zurek scaling predictions 3 in theXYmodel, as well as signatures of the classical Kosterlitz–Thouless phase transition 4 . Moreover, the digital gates enable precise energy control, allowing us to study the effects of the eigenstate thermalization hypothesis 5–7 in targeted parts of the eigenspectrum. We also demonstrate digital preparation of pairwise-entangled dimer states, and image the transport of energy and vorticity during subsequent thermalization in analogue evolution. These results establish the efficacy of superconducting analogue–digital quantum processors for preparing states across many-body spectra and unveiling their thermalization dynamics.

Science & Technology - Other Topics

Interactions Between Climate Policy and Technology-influenced Travel Behavior: Mitigating Induced Demand from CACC

Advances in vehicle technology have influenced the development of automated vehicle systems, where vehicles that do not require human intervention are already deployed in the roadway networks. While these advances are proved to increase roadway safety and highway capacity, more research is needed to understand the long-term and regional-level impacts on mobility, land use, energy consumption, and emissions. This study proposes a multi-model approach to analyze the effect of vehicle automation and deep decarbonization policies over a period from 2020 to 2040 in Austin, Texas. We use the Global Change Analysis Model (GCAM) to develop internally the scenarios that are then passed to the SMART Mobility modeling workflow, a large-scale simulation framework combining the POLARIS activity-based travel demand model and mesoscopic traffic simulator with the Autonomie vehicle energy consumption model and the UrbanSim land use simulator. Results suggest that the introduction of vehicles with advanced automation could increase fuel consumption when no decarbonization policies are implemented. Also, advances in vehicle technology research and development could lead to a decline in energy use in the long-term. Energy pricing and vehicle electrification incentives could help reduce the impact of vehicle automation. Finally, our analysis indicates the relevance of introducing land use processes in longterm vehicle automation studies.

land use

The impact of aerosol mixing state on immersion freezing: insights from classical nucleation theory and particle-resolved simulations

Immersion freezing, initiated by ice-nucleating particles (INPs) in supercooled aqueous droplets, plays an important role in the formation of ice crystals within clouds. The efficiency of immersion freezing depends strongly on INP composition and, crucially, on the mixing state – how chemical species are distributed across the particle population. Here, we quantify the impact of aerosol mixing state on immersion freezing using a combined theoretical and particle-resolved modeling approach. We derive analytical expressions for the frozen fraction of internally and externally mixed INP populations based on classical nucleation theory, showing that the frozen fraction is sensitive to whether ice-active species are present in all particles or only in a subset of the population. We introduce a multi-species immersion freezing scheme into the particle-resolved model PartMC, using the water activity-based immersion freezing model (ABIFM) to compute freezing probabilities for mixed-composition particles. To improve computational efficiency, we implement a Binned Tau-Leaping algorithm and demonstrate an order-of-magnitude speedup with minimal accuracy loss. Simulations reproduce the analytical trends in limiting cases and extend the analysis to more general aerosol populations, where mixing state continues to exert a substantial control on frozen fraction. Sensitivity analyses across particle size, species type, and cooling condition reveal that the mixing state effect is most pronounced when small amounts of highly efficient INPs are mixed with less efficient materials. These findings underscore the need to represent aerosol mixing state explicitly in models of heterogeneous ice nucleation to reduce uncertainty in cloud-phase partitioning.

54 ENVIRONMENTAL SCIENCES

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE

Spike-Nosed Bodies and Forward Injected Jets in Supersonic Flow

The paper contains new numerical simulation and experimental test results of blunt body drag reduction using thin spikes mounted in front of a body and one- or two-phase jets injected against a supersonic flow. Numerical simulations utilizing the NASA CFL3D code were conducted at the Hampton University Fluid Mechanics and Acoustics Laboratory (FM&AL) and experimental tests were conducted using the facilities of the IM/MSU Aeromechanics and Gas Dynamics Laboratory. Previous results were presented at the 37 t h AIAA/ASME/SAE/ASEE Joint Propulsion Conference. Those results were based on some experimental and numerical simulation tests for supersonic flow around spike-nosed or shell-nosed bodies, and numerical simulations were conducted only for a single spike-nosed or shell-nosed body at zero attack angle, alpha=0 o . In this paper, experimental test results of gas, liquid and solid particle jet injection against a supersonic flow are presented. In addition, numerical simulation results for supersonic flow around a multiple spike-nosed body with non-zero attack angles and with a gas and solid particle forward jet injection are included. Aerodynamic coefficients: drag, C D , lift, C L , and longitudinal momentum, M z , obtained by numerical simulation and experimental tests are compared and show good agreement.

M Gilinsky