Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “building 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 703 records · Page 39

Coupled Induction Machine and HVAC Models for Simulating HVAC Performance Considering Grid Dynamics in Buildings

This paper presents the development of novel models that integrate induction machines with HVAC equipment, such as pumps, heat pumps, and chillers, to analyze the impact of electrical parameters on the operational performance of thermo-fluid systems. The proposed model employs a coupling technique that captures the dynamic interactions between induction machines and HVAC systems. By integrating electrical, thermal, and mechanical dynamics, the models provide a comprehensive framework for simulating real-world scenarios, including interactions with the electrical grid. This achievement was made possible through the development of a Computationally Efficient and Accurate Induction Machine (CEAIM) model. Implemented using the equation-based Modelica language, the CEAIM model has been validated against experimental results, manufacturer data sheets, and various operating conditions. Its performance has been compared with existing induction machine models in the Modelica Standard Library (MSL), demonstrating superior accuracy and computational efficiency. The CEAIM model predicts torque, speed, and power consumption with a coefficient of determination (R 2 ) ranging from 0.98 to 1 and a coefficient of variation of root mean square error (CVRMSE) between 0.27% and 6.67%. Additionally, CEAIM scales more efficiently than conventional MSL models, with a slower computational growth rate in large-scale simulations. After thorough validation of the CEAIM model, it was coupled with HVAC equipment as this approach provides a detailed multi-dimensional view of capturing electrical transients and mechanical performance. To support this, a case study was conducted to showcase its capabilities.

24 POWER TRANSMISSION AND DISTRIBUTION↗

AutoBEMGen

Generating energy models from scratch is a costly and time consuming task. A detailed energy model can have thousands of input parameters that are necessary for a complete simulation of a building. A solution to this problem is the AutoBEMGen framework. This Python-based framework requires a user supplied data file, currently in comma separated variable format, that defines a subset of building properties and parameters to quickly and easily generate a building energy model. This building properties include building type, 2-D footprint, number of floors, window-to-wall ratio, and more. AutoBEMGen then uses OpenStudio and EnergyPlus, two of DOE's flagship energy modeling tools, to create and simulate an OpenStudio model and EnergyPlus input file. Additionally, the AutoBEMGen framework is scalable to any number of user supplied buildings in the data field by leveraging multi-threading built into the Python programming language. By leveraging a small subset of building inputs required from end users, AutoBEMGen is able to generate complete building enervy models in a fraction of the time compared to generating them from scratch.

Adams, Mark B↗

Evaluate the impact of sensor accuracy on model performance in data-driven building fault detection and diagnostics using Monte Carlo simulation

The performance of data-driven fault detection and diagnostics (FDD) is heavily dependent on sensors. However, sensor inaccuracy and sensor faults are pervasive in building operation: inaccurate and missing sensor readings deteriorate FDD performance; sensor inaccuracy will also affect the selection of sensor for data-driven FDD in the model training process, which is another key factor of data-driven FDD performance. Sensor accuracy and sensor selection individually are well-studied research topics in this field, but the impact of sensor accuracy on sensor selection and its further impact on FDD performance has not been evaluated and quantified. In this paper, we developed a novel analysis methodology that comprehensively evaluates sensor fault on sensor selection and FDD accuracy. Monte Carlo simulation is applied to deal with multiple stochastic sensor inaccuracy and provide probabilistic analysis results of the impact of sensor inaccuracy on sensor selection and FDD accuracy. This methodology focuses on the net impact of fault states across a full sensor set. The developed methodology can be used for the early-stage sensor design and operation-stage sensor maintenance. Furthermore, a case study is conducted to demonstrate the analysis methodology using a commercial building model crated to Flexible Research Platform located at Oak Ridge National Laboratory, USA.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Large-Scale Simulation of Regional Demand Flexibility Implementation and Customer Economic Impact

The Distribution System Operator with Transactive (DSO+T) study conducted a large-scale simulation of over 60,000 customers in a region the size of Texas to demonstrate the effective coordination of distributed energy resources (DERs) in commercial and residential buildings. The integrated simulation included both the bulk (wholesale generation and transmission) and distribution systems. The DERs (including batteries, electric vehicles, air conditioning, and water heaters) participated in a transactive energy retail market that was integrated into an existing competitive wholesale market. The engineering and economic performance of the resulting demand flexibility was evaluated over annual simulations for both moderate and high renewable generation scenarios. A detailed parametric cost model was developed to enable detailed economic analysis of key stakeholders. In addition, fixed and dynamic customer tariffs were designed and applied to the customer population. This allowed the impact on annual customer bills to be analyzed for various building types (residential versus commercial; single- versus multi-family). This paper presents results showing the relative flexibility of batteries, electric vehicles, and building loads throughout the year and under different renewable scenarios. This feeds a detailed breakdown of the impact this flexibility has on the operating costs of the grid and the resulting net economic benefit. Finally, the study showed that practically all customer classes (including non-participating customers) save money under the proposed demand flexibility scheme. The study found overall net annual economic savings of $3.3-5.0B for a region the size of Texas equating to average customer bill savings of 10-16%.

Reeve, Hayden M.↗

Distributed model predictive control for coordinated, grid-interactive buildings

Continued focus on reducing carbon emissions and improving energy efficiency requires buildings to become grid-interactive and not just behave as static consumers. A distributed model predictive control (DMPC) algorithm known as Limited-Communication (LC) DMPC is modified to enable grid-interactive buildings. A grid-aggregator subsystem is added that allows for a bulk grid power reference signal to be followed while the individual building subsystems also achieve their local comfort objectives. The LC-DMPC algorithm is applied for the first time to systems with multiple buildings. Adequate power tracking is shown for different simulation scenarios involving heterogeneous buildings, and next steps are discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Understanding Uncertainty in Microstructure Evolution and Constitutive Properties in Additive Process Modeling

Coupled process–microstructure–property modeling, and understanding the sources of uncertainty and their propagation toward error in part property prediction, are key steps toward full utilization of additive manufacturing (AM) for predictable quality part development. The OpenFOAM model for process conditions, the ExaCA model for as-solidified grain structure, and the ExaConstit model for constitutive mechanical properties are used as part of the ExaAM modeling framework to examine a few of the various sources of uncertainty in the modeling workflow. In addition to “random” uncertainty (due to random number generation in the orientations and locations of grains present), the heterogeneous nucleation density N 0 and the mean substrate grain spacing S 0 are varied to examine their impact of grain area development as a function of build height in the simulated microstructure. While mean grain area after 1 mm of build is found to be sensitive to N 0 and S 0 , particularly at small N 0 and large S 0 (despite some convergence toward similar values), the resulting grain shapes and overall textures develop in a reasonably similar manner. As a result of these similar textures, ExaConstit simulation using ExaCA representative volume elements (RVEs) from various permutations of N 0 , S 0 , and location within the build resulted in similar yield stress, stress–strain curve shape, and stress triaxiality distributions. It is concluded that for this particular material and scan pattern, 15 layers is sufficient for ExaCA texture and ExaConstit predicted properties to become relatively independent of additional layer simulation, provided that reasonable estimates for N 0 and S 0 are used. However, additional layers of ExaCA will need to be run to obtain mean grain areas independent of build height and baseplate structure.

36 MATERIALS SCIENCE↗

Artificial Intelligence Thermostat to Detect Faults

Residential air conditioners and heat pumps often experience faults due to inadequate maintenance, which can severely reduce efficiency or even cause system failure. Common issues include dirty or clogged air filters and refrigerant leaks. These problems degrade performance and increase energy use and operating costs. This study presents a smart thermostat with embedded artificial intelligence to detect such faults and alert homeowners when maintenance is needed. The thermostat uses low-cost measurements—including return-air temperature, relative humidity, supply-air temperature, outdoor-air temperature, and condenser subcooling—to identify abnormal operations. Because different faults produce distinct response patterns, tailored algorithms are developed to recognize characteristic fault signatures. The investigation is built on a detailed co-simulation platform that couples EnergyPlus with the DOE/ORNL Heat Pump Design Model (HPDM). EnergyPlus represents the building’s dynamic environment, while HPDM is a high-fidelity, hardware-based model that can simulate fault-free performance as well as a wide range of faults, including gradual degradation such as minor refrigerant leakage. This platform provides a virtual training and testing environment that helps distinguish fault-induced behavior from normal operation and supports development of robust diagnostic algorithms. Using this framework, a Dynamic Bayesian Network was developed to identify two common faults—gradual refrigerant charge loss and indoor airflow blockage—and the AI-embedded thermostat was verified through annual building simulations.

Shen, Bo [ORNL] (ORCID:0000000336600393)↗

Building wakes in disturbed layers

The wind profiles in the wakes of buildings in a disturbed shear layer are investigated for conditions experienced during aircraft take-off and landing. Experiments were conducted in a pipe flow apparatus in which an equilibrium flow is disturbed by pipe roughness, simulating obstacles upwind of a building, with the disturbed flow made to flow past a second obstacle, simulating a building. Velocity and turbulence profiles of disturbed and equilibrium flows upstream and downstream of the obstacles are presented. It is concluded that the integrated momentum flow of a boundary layer behind an obstacle can be greater than the upstream momentum flow, if the exponent of the upstream velocity profile is sufficiently high and the boundary layer has been disturbed by upstream obstacles. The rise of the turbulence level behind the obstacle is also dependent on the profile exponent in a disturbed layer.

Logan, E.↗

Model America - data and models of every U.S. building

The 5-year goal of the 'Model America' concept was to generate a model of every building in the United States. This data repository delivers on that goal. Oak Ridge National Laboratory (ORNL) has developed the Automatic Building Energy Modeling (AutoBEM) software suite to process multiple types of data, extract building-specific descriptors, generate building energy models, and simulate them on High Performance Computing (HPC) resources. For more information, see AutoBEM-related publications (bit.ly/AutoBEM). There were 125,714,640 buildings detected in the United States and this dataset contains 122,930,327 (97.8%) buildings which resulted in a successful simulation. Future, annual updates have been proposed that may include additional buildings, data improvements, or other algorithmic enhancements. This dataset of 122.9 million buildings includes: Models (state_county.zip) - OpenStudio (v3.1.0) and EnergyPlus (v9.4) building energy models. Please note that the download requires the free Globus Connect Personal (https://www.globus.org/globus-connect-personal); Each model has approximately 3,000 building input descriptors that can be extracted. Please see the EnergyPlus(v9.4) 2,784-page Input/Output Reference Guide (https://energyplus.net/sites/all/modules/custom/nrel_custom/pdfs/pdfs_v9.4.0/InputOutputReference.pdf) for everything that can be retrieved or simulated from these models. These models were derived from the following metadata, which is not included in this dataset: 1. ID - unique building ID 2. County - county name 3. State - state name 4. CZ - ASHRAE Climate Zone designation 5. Clim_Zone - text label of climate zone 6. est_year - estimated year of construction 7. est_commercial - estimated building type (0=residential, 1=commercial) 8. Centroid - building center location in latitude/longitude (from Footprint2D) 9. Footprint2D - building polygon of 2D footprint (lat1/lon1_lat2/lon2_...) 10. Height - building height (meters) 11. Area2D - footprint area (ft2) 12. BuildingType - DOE prototype building designation (IECC=residential) as implemented by OpenStudio-standards 13. WWR_surfaces - percent of each facade (pair of points from Footprint2D) covered by fenestration/windows (average 14.5% for residential, 40% for commercial buildings) 14. NumFloors - number of floors (above-grade) 15. Area - estimate of total conditioned floor area (ft2) 16. Standard - building vintage. These models are made free and openly available in hopes of stimulating any simulation-informed use case. Data is provided as-is with no warranties, express or implied, regarding fitness for a particular purpose. We wish to thank our sponsors which include Oak Ridge National Laboratory (ORNL) Laboratory Directed Research and Development (LDRD), U.S. Dept. of Energy's (DOE) Building Technologies Office (BTO), Office of Electricity (OE), Biological and Environmental Research (BER), and National Nuclear Security Administration (NNSA). This research used resources of the Argonne Leadership Computing Facility, which is a DOE Office of Science User Facility supported under Contract DE-AC02-06CH11357. Please cite as: New, Joshua R., Adams, Mark, Bass, Brett, Berres, Anne, and Clinton, Nicholas (2021). 'Model America - data and models of every U.S. building. [Data set].' Constellation, doi.ccs.ornl.gov/ui/doi/339, April 14, 2021

24 POWER TRANSMISSION AND DISTRIBUTION↗

Effect of Solvents on Lignin–Surface Interactions via Molecular Dynamics Simulations

Lignin, an essential building block of lignocellulosic biomass, is a potential abundant source of aromatic monomers for the polymer and chemical industry. Reductive catalytic fractionation (RCF) is one promising process that can produce high yields of phenolic monomers and oligomers from lignin under different catalytic conditions. An important choice in optimizing RCF is the selection of solvent; however, detailed insights into solvent effects on lignin behaviors and interactions remain limited. Here, in this work, we perform all-atom molecular dynamics simulations to study the solvation of lignin, solvent-mediated conformational changes, and the interaction of solvated lignin oligomers with model surfaces. We focus on the behavior of an oligomeric lignin model compound in methanol, ethanol, a binary mixture of ethanol and water, and water at both the RCF reaction temperature (473 K) and room temperature. Analysis of structural features of lignin suggests that these three organic solvent systems favorably solvate lignin, resulting in a more extended conformation suitable for catalytic conversion to valuable chemicals. We further introduce model palladium (Pd) and carbon (C) surfaces to understand how solvent choice impacts adsorption onto a representative catalytic surface and support, and to quantify the competition among the reactant and solvent molecules for the surface. Unbiased simulations suggest strong adsorption of lignin on both Pd and C surfaces at 473 K, with notable solvent-mediated differences in adsorption energies. Additionally, our findings indicate that lignin adsorption is promoted by the entropy change resulting from the displacement of solvent molecules from the surface. This study provides a molecular perspective of adsorption of lignin onto varying surfaces, which is a step towards understanding and optimizing the catalytic conversion of lignin into valuable chemicals.

adsorption↗

Energy Evaluation of a New Construction Pilot Community: Fresno, California

A new construction pilot community was constructed by builder-partner Wathen-Castanos Hybrid Homes based on a single occupied test house that was designed to achieve greater than 30% energy savings with respect to the Building America House Simulation Protocols developed by NREL. Builders face several key problems when implementing a whole-house systems integrated measures package from a single test house into multiple houses. This report addresses factors of implementation and scalability at the community scale and proposes methodologies by which community-scale energy evaluations can be performed based on results at the occupied test house level. Research focused on the builder and trade implementation of a measures package and the actual utility usage in the houses at the community scale of production. Five occupants participated in this research by providing utility bills and information on occupancy and miscellaneous gas and electric appliance use for their houses. IBACOS used these utility data and background information to analyze the actual energy performance of the houses. The actual utility bill readings were compared to projected energy consumption using BEopt with actual weather and thermostat set points for normalization.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A Study of the Role of Clouds in the Relationship Between Land Use/Land Cover and the Climate and Air Quality of the Atlanta Area

The goal of Project ATLANTA is to derive a better scientific understanding of how land cover changes associated with urbanization affect local and regional climate and air quality. Clouds play a significant role in this relationship. Using GOES images, we found that in a 63-day period (5 July-5 September 1996) there were zero days which were clear for the entire daylight period. Days which are cloud-free in the morning become partly cloudy with small cumulus clouds in the afternoon in response to solar heating. This result casts doubt on the applicability of California-style air quality models which run in perpetual clear skies. Days which are clear in the morning have higher ozone than those which are cloudy in the morning. Using the RAMS model, we found that urbanization increases the skin surface temperature by about 1.0-1.5 C on average under cloudy conditions, with an extreme of +3.5 C. Clouds cool the surface due to their shading effect by 1.5-2.0 C on average, with an extreme of 5.0 C. RAMS simulates well the building stage of the cumulus cloud field, but does poorly in the decaying phase. Next year's work: doing a detailed cloud climatology and developing improved RAMS cloud simulations.

Kidder, Stanley Q.↗

Zinc Corrosion in Borated, Buffered Post-LOCA Solutions

This current research provides a broad assessment of recently evaluated zinc corrosion in simulated non-buffered and buffered post-LOCA solution as a function of pH, temperature, velocity, boron concentration and surface area. Corrosion rates obtained can be used for further zinc borate formation analysis or to plan further zinc corrosion studies. Experiments are conducted utilizing a three-electrode electrochemical system in a glass cell. Aqueous solutions used in testing simulate reactor containment building (RCB) pool water under various LOCA scenarios that are representative of US PWR RCB pool water: boric acid only, boric acid with trisodium phosphate, boric acid with sodium hydroxide, and boric acid with sodium tetraborate. The specimens used in this work are machined from 99.6% pure zinc rods that are used as a surrogate for the high purity zinc layer of galvanized steel. Prior to the experiment, the specimens are first polished and then left exposed to the atmosphere for approximately one day at room temperature. Prior to each test, the zinc specimen are sequentially rinsed with deionized water and acetone, and then dried by evaporation in air at room temperature. Prepared samples are then mounted in the glass cell test apparatus and rotated on a spindle at a desired velocity for a specified duration. Instantaneous corrosion rate as a function of time and total zinc mass corroded, which may either be released into solution or incorporated into surface scale, is determined by electrochemical impedance spectroscopy (EIS) analysis. Corrosion rates after 24 hours of immersion are determined from a full potential dynamic scan and from weight loss measurements of descaled post-test samples. The solution zinc concentration resulting from zinc corrosion is determined by inductively coupled plasma optical emission spectroscopy (ICP-OES). Preliminary results show that zinc corrosion with a zinc surface area-to-liquid volume ratio of 0.0003 m2/L in a borated, buffered solution (55 deg. C, pH 7.2) exposed to a velocity of 0.5 ft/s is between 4.9 X to 2.3 X less than that obtained by the Kyrk et al study, which examines a lower temperature, lower pH with a 5.5 X larger surface area-to-volume ratio. This initial finding suggests that the corrosion of zinc, which supplies a zinc source term to form zinc borates, may not be insignificant under conditions of interest (borated, buffered post-LOCA solutions). Further experiments are currently underway that evaluate the influence of pH, velocity, temperature, buffer type, boron concentration and surface area-to-volume ratio on zinc corrosion. Preliminary results suggest two of the three buffer types examined as conditions of interest to US PWRs likely produce zinc corrosion within the range observed by Kyrk et al. Current preliminary results indicate that buffer type, temperature, pH and velocity each have an impact on the zinc corrosion rate. Quantification of zinc mass lost to solution and incorporated into scale is currently underway, as is the characterization of velocity and boric acid effects on zinc corrosion.

36 MATERIALS SCIENCE↗

Characteristics and selection of near-fault simulated earthquake ground motions for nonlinear analysis of buildings

Earthquake-induced ground shaking near rupturing faults is highly sensitive to the rupture characteristics, seismic wave propagation patterns and site conditions, and field recordings of near-fault shaking are relatively sparse. These challenges complicate the assessment of the seismic performance of near-fault structures. A common approach to representing near-fault ground motion in engineering analysis is to explicitly consider and select records with strong directivity pulses (pulse records). We use three-dimensional high-resolution physics-based earthquake simulations to test this approach in the context of scenario-based ground motion record selection, and to study the important characteristics of near-fault ground shaking. We highlight the deficiencies associated with classifying near-fault simulated records as “pulse” or “non-pulse,” based on the presence of a single dominating pulse in the velocity time history. We show that this approach is inadequate for characterizing near-fault shaking on soft soils which can be dominated by both forward rupture directivity and basin amplification effects. We conduct ground motion selection experiments for the analysis of near-fault structures with and without explicit classification of the pulse features in the records, and evaluate the bias in the predicted structural demands. We find that the maximum interstory drift demands on building structures imposed by unscaled site-specific simulated ground motion records selected based on relevant spectral shape features are not sensitive to the classification of records as pulse/non-pulse. Therefore, with regard to predicting the maximum interstory drifts in near-fault buildings, we do not find justification for the binary pulse classification of near-fault records.

42 ENGINEERING↗

Hydronic Heating Coil Versus Propane Furnace, Rehoboth Beach, Delaware (Fact Sheet)

Insight Homes constructed two houses in Rehoboth Beach, Delaware, with identical floor plans and thermal envelopes but different heating and domestic hot water (DHW) systems. Each house is 1,715-ft 2 with a single story, three bedrooms, two bathrooms, and the heating, ventilation, and air conditioning (HVAC) systems and ductwork located in conditioned crawlspaces. The standard house, which the builder offers as its standard production house, uses an air source heat pump (ASHP) with supplemental propane furnace heating. The Building America test house uses the same ASHP unit with supplemental heat provided by the DHW heater (a combined DHW and hydronic heating system, where the hydronic heating element is in the air handler). Both houses were occupied during the test period. Results indicate that efficiency of the two heating systems was not significantly different. Three issues dominate these results; lower system design performance resulting from the indoor refrigerant coil selected for the standard house, an incorrectly functioning defrost cycle in the standard house, and the low resolution of the natural gas monitoring equipment. The thermal comfort of both houses fell outside the ASHRAE Standard 55 heating range but was within the ACCA room-to-room temperature range when compared to the thermostat temperature. The monitored DHW draw schedules were input into EnergyPlus to evaluate the efficiency of the tankless hot water heater model using the two monitored profiles and the Building America House Simulation Protocols. The results indicate that the simulation is not significantly impacted by the draw profiles.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Comparison of effectiveness and energy use of airborne pathogen mitigation measures to meet clean air targets in a prototypical office building

Organizations such as ASHRAE and the Centers for Disease Control and Prevention (CDC) have proposed guidelines for controlling infectious aerosols in buildings, which can be met through measures such as modified operation of the heating, ventilation, and air-conditioning (HVAC) system or incorporating air-cleaning technologies. However, more research is needed to understand the trade-offs between health, energy, and comfort aspects when designing measures for these guidelines. To address this gap, this paper presents an analysis using new models for air-cleaning technologies, including in-duct and in-room germicidal ultraviolet (GUV) systems and portable air cleaners (PACs). These models are incorporated into an existing prototypical office building model and six measures are designed to meet ASHRAE Standard 241 and CDC clean air targets: MERV 13 HVAC filtration, maximum outdoor air supplied to the building, PACs, and in-duct, upper-room, and whole-room GUV. The measures are simulated for an office building in a cool and humid climate compared against a baseline simulation using MERV 8 filtration. The results show that all measures, except for the maximum outdoor air case, can meet the ASHRAE 241 standard without significant impacts on energy or comfort. The HVAC system measures were not able to meet the CDC target with the default system sizing and lead to significant energy increases, while the in-room measures were able to meet the CDC target with small impacts on energy consumption. Finally, this paper consolidates the simulation findings and provides practical guidance for building operators to meet clean air targets while limiting energy and comfort impacts.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Computational Simulation of an Additively Manufactured Marine Component

Early work in finite-element analysis of additive manufacturing (AM) components was not able to include distortion and residual stress of the entire build, yielding uncertainty in the computational model predictions. To address this, current versions of recently developed computational tools are being used to understand their role in validation at the Naval Surface Warfare Center, Carderock Division, for laser powder-bed fusion (L-PBF) AM. A 316L stainless steel L-PBF build containing multiple parts, which exhibited the same local part failure over multiple AM fabrication runs, was selected for computational analysis. Single-part simulations analyzed optimization of computational speed and prediction accuracy. Full-build simulations exposed the smallest re-coater tolerances early in part fabrication as a result of part shrinkage and subsequent cupping. This cupping was predicted to change the powder spreading and result in the non-consolidation of the support structures, which paralleled the failures seen in the physical builds. Part redesign of the conventionally produced part to align with AM fabrication requirements was recommended to alleviate the build failures. This project is part of a larger program focused on validation of the computational tools to determine best practices for insertion of AM simulation into the greater Naval enterprise.

36 MATERIALS SCIENCE↗

End-Use Analysis of ASHRAE Standard 90.1-2022

These data dashboards present the data and findings from PNNL's detailed analysis of ASHRAE Standard 90.1 on US commercial new construction buildings. The dashboards support simulation-based assessment of U.S. commercial buildings compliant with ASHRAE Standard 90.1-2004, 2007, 2009, 20210, 2013, 2016, 2019, and 2022.

Nambiar, Chitra↗