Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “operator splitting”

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 739 records · Page 41

Advanced Design for the WIQ Magnet With Steering Corrector Function

The Facility for Rare Isotopes Beams (FRIB) delivers heavy-ion primary beams at energies of up to 300 MeV/u at 10 kW of beam power to generate rare isotope beams for experiments and will eventually operate at beam power of 400 kW. The preseprator of the Advanced Rare Isotope Separator (ARIS) is equipped with six warm-iron quadrupole (WIQ) singlets and two dipoles integrated right after the production target. They have a compact structure and operate in a high radiation vacuum environment within a hot cell having remote handling capabilities for installation and maintenance. Due to asymmetry with respect to the quadrupole poles, nested sextupole excitations in WIQs induce vertical dipoles that offset the centroid trajectory; Magnet misalignments also result in trajectory offsets. Such offsets degrade separator performance but can be minimized by changing the current distribution on sextupole and octupole coils. In this work, we show how modifications to the WIQ coil design can allow superimposed dipole fields to be included to the octupole and sextupole windings, as well as addition of dipole components by splitting coil currents over groups with separator power supplies. Adjusting the group currents can cancel the sextupole-induced vertical dipole component which can be as high as 0.012 Tm. Octupole coil changes may superimpose a horizontal dipole integrated strength as high as 0.0332 Tm. Unwanted higher harmonics induced as a side effect of the new design are kept to a minimum such that separator performance is preserved as much as possible.

Accelerator magnets↗

Optimization Studies of Radiation Shielding for PIP-II Project at Fermilab

The Proton Improvement Plan-II (PIP-II) at Fermilab represents a significant advancement in the quest to answer some of the most profound questions about our universe using the world's most intense high-energy neutrino beam. The project requires the construction of a new addition to the Fermilab accelerator complex – an 800-MeV high-intensity superconducting linear accelerator. Ensuring the safety and regulatory compliance of this ambitious project is paramount, necessitating thorough dose rate assessments under both normal operational and accidental scenarios to align with the Fermilab Radiological Control Manual (FRCM) standards. Our approach included a shielding optimization used for the simulations with the Monte Carlo code MARS [1,2,3] to incorporate new magnet and collimator designs, essential for reflecting the current state of PIP-II infrastructure. The implementation of high-resolution detector planes, despite their computational demands, enabled us to gather detailed radiation field data crucial for optimizing shielding configurations. To overcome the significant computational demands, we developed a branching code that drastically reduced simulation runtimes while maintaining statistical integrity. This was achieved through geometry splitting and the application of Russian Roulette techniques, tailored to prioritize regions of interest based on predefined importances and weight limits.

43 PARTICLE ACCELERATORS↗

Optimization Studies of Radiation Shielding for PIP-II Project at Fermilab

The Proton Improvement Plan-II (PIP-II) at Fermilab represents a significant advancement in the quest to answer some of the most profound questions about our universe using the world's most intense high-energy neutrino beam. The project requires the construction of a new addition to the Fermilab accelerator complex – an 800-MeV high-intensity superconducting linear accelerator. Ensuring the safety and regulatory compliance of this ambitious project is paramount, necessitating thorough dose rate assessments under both normal operational and accidental scenarios to align with the Fermilab Radiological Control Manual (FRCM) standards. Our approach included a shielding optimization used for the simulations with the Monte Carlo code MARS [1,2,3] to incorporate new magnet and collimator designs, essential for reflecting the current state of PIP-II infrastructure. The implementation of high-resolution detector planes, despite their computational demands, enabled us to gather detailed radiation field data crucial for optimizing shielding configurations. To overcome the significant computational demands, we developed a branching code that drastically reduced simulation runtimes while maintaining statistical integrity. This was achieved through geometry splitting and the application of Russian Roulette techniques, tailored to prioritize regions of interest based on predefined importances and weight limits.

Makovec, Alajos↗

The Slope Imaging Multi-Polarization Photon-Counting Lidar: Development and Performance Results

The Slope Imaging Multi-polarization Photon-counting Lidar is an airborne instrument developed to demonstrate laser altimetry measurement methods that will enable more efficient observations of topography and surface properties from space. The instrument was developed through the NASA Earth Science Technology Office Instrument Incubator Program with a focus on cryosphere remote sensing. The SIMPL transmitter is an 11 KHz, 1064 nm, plane-polarized micropulse laser transmitter that is frequency doubled to 532 nm and split into four push-broom beams. The receiver employs single-photon, polarimetric ranging at 532 and 1064 nm using Single Photon Counting Modules in order to achieve simultaneous sampling of surface elevation, slope, roughness and depolarizing scattering properties, the latter used to differentiate surface types. Data acquired over ice-covered Lake Erie in February, 2009 are documenting SIMPL s measurement performance and capabilities, demonstrating differentiation of open water and several ice cover types. ICESat-2 will employ several of the technologies advanced by SIMPL, including micropulse, single photon ranging in a multi-beam, push-broom configuration operating at 532 nm.

Dabney, Phillip↗

Probing the flavour structure of dimension-6 EFT operators in multilepton final states in proton-proton collisions at $\sqrt{s}=13$ TeV

An analysis of the flavour structure of dimension-6 effective field theory (EFT) operators in multilepton final states is presented, focusing on the interactions of quarks with Z bosons. For the first time, the flavour structure of these operators is disentangled by simultaneously probing the interactions with different quark generations. The analysis targets the associated production of a top quark pair and a Z boson, as well as diboson processes in final states with at least three leptons, which can be electrons or muons. The data were recorded by the CMS experiment in the years 2016–2018 in proton-proton collisions at a centre-of-mass energy of 13 TeV and correspond to an integrated luminosity of 138 fb −1 . Consistency with the standard model of particle physics is observed and limits are set on the selected Wilson coefficients, split into couplings to light- and heavy-quark generations.

Flavour Physics↗

Lattice Mie resonances and emissivity enhancement in mid-infrared iron pyrite metasurfaces

High-refractive-index antennas with characteristic dimensions comparable to wavelength have a remarkable ability to support pronounces electric and magnetic dipole resonances. Furthermore, periodic arrangements of such resonant antennas result in narrow and strong lattice resonances facilitated by the lattice. We design iron pyrite antennas operating in the mid-infrared spectral range due to the material’s low-energy bandgap and high refractive index. We utilize Kirchhoff’s law, stating that emissivity and absorptance are equal to each other in equilibrium, and we apply it to improve the thermal properties of the iron pyrite metasurface. Through the excitation of collective resonances and manipulation of the antenna lattice’s period, we demonstrate our capacity to control emissivity peaks. These peaks stem from the resonant excitation of electric and magnetic dipoles within proximity to the Rayleigh anomalies. In the lattice of truncated-cone antennas, we observe Rabi splitting of electric and magnetic dipole lattice resonances originating from the antennas’ broken symmetry. We demonstrate that the truncated-cone antenna lattices support strong out-of-plane magnetic dipole lattice resonances at oblique incidence. We show that the truncated-cone antennas, as opposed to disks or cones, facilitate a particularly strong resonance and bound state in the continuum at the normal incidence. Our work demonstrates the effective manipulation of emissivity peaks in iron pyrite metasurfaces through controlled lattice resonances and antenna design, offering promising avenues for mid-infrared spectral engineering.

Islam, Md Sakibul↗

Investigation of Dynamic Force/Vibration Transmission Characteristics of Four-Square Type Gear Durability Test Machines

In this study, design requirements for a dynamically viable, four-square type gear test machine are investigated. Variations of four-square type gear test machines have been in use for durability and dynamics testing of both parallel- and cross-axis gear set. The basic layout of these machines is illustrated. The test rig is formed by two gear pairs, of the same reduction ratio, a test gear pair and a reaction gear pair, connected to each other through shafts of certain torsional flexibility to form an efficient, closed-loop system. A desired level of constant torque is input to the circuit through mechanical (a split coupling with a torque arm) or hydraulic (a hydraulic actuator) means. The system is then driven at any desired speed by a small DC motor. The main task in hand is the isolation of the test gear pair from the reaction gear pair under dynamic conditions. Any disturbances originated at the reaction gear mesh might potentially travel to the test gearbox, altering the dynamic loading conditions of the test gear mesh, and hence, influencing the outcome of the durability or dynamics test. Therefore, a proper design of connecting structures becomes a major priority. Also, equally important is the issue of how close the operating speed of the machine is to the resonant frequencies of the gear meshes. This study focuses on a detailed analysis of the current NASA Glenn Research Center gear pitting test machine for evaluation of its resonance and vibration isolation characteristics. A number of these machines as the one illustrated has been used over last 30 years to establish an extensive database regarding the influence of the gear materials, processes surface treatments and lubricants on gear durability. This study is intended to guide an optimum design of next generation test machines for the most desirable dynamic characteristics.

Kahraman, Ahmet↗

Modulation of the coordination environment enhances the electrocatalytic efficiency of Mo single atoms toward water splitting

Here, enhancing the catalytic efficiency through engineering active site environments is expected to work pronouncedly for single atom catalysts (SACs) because of intense atomic scale interactions involved between SAs and their coordination environments. Taking Mo SACs for catalyzation of the hydrogen evolution reaction (HER) as an example, three SACs of different coordination environments, namely Mo-O 2 N 2 , Mo-O 2 N 1 C 1 , and Mo-O 2 C 2 , were successfully created for demonstration. The HER performances are in an increasing order of Mo-O 2 N 2 , Mo-O 2 N 1 C 1 , and Mo-O 2 C 2 , exhibiting η 10 of 98, 71, and 61 mV, η 500 of 340, 248, and 200 mV, Tafel slopes of 95.8, 39.6, and 33.8 mV dec -1 , and current density decays of 9, 6, and 6% after a 50 hour operation at an initial current density of 100 mA cm -2 , respectively. Substituting C with N in the coordination environment results in inferior HER catalytic efficiency and stability. Density functional theory calculations reveal that replacing carbon with nitrogen for coordination with the Mo SA on a carbon substrate of a higher N-doping level shifts the d-band center of Mo more negatively from the Fermi level, thereby increasing the hydrogen adsorption energy and thus decelerating the hydrogen desorption kinetics, giving consequent inferior HER activities.

25 ENERGY STORAGE↗

Broadband integrated polarization splitter and rotator using subwavelength grating claddings

We present a broadband integrated photonic polarization splitter and rotator (PSR) using adiabatically tapered coupled waveguides with subwavelength grating (SWG) claddings. The PSR adiabatically rotates and splits the fundamental transverse-magnetic (TM 0 ) input to the fundamental transverse-electric (TE 0 ) mode in the coupler waveguide, while passing the TE 0 input through the same waveguide. The SWGs work as an anisotropic metamaterial and facilitate modal conversions, making the PSR efficient and broadband. We rigorously present our design approaches in each section and show the SWG effect by comparing with and without the SWG claddings. The coupling coefficients in each segment explicitly show a stronger coupling effect when the SWGs are included, confirmed by the coupled-mode theory simulations. The full numerical simulation shows that the SWG-PSR operates at 1500–1750 nm (≈250 nm) wavelengths with an extinction ratio larger than 20 dB, confirmed by the experiment for the 1490–1590 nm range. The insertion losses are below 1.3 dB. Since our PSR is designed based on adiabatical mode evolution, the proposed PSR is expected to be tolerant to fabrication variations and should be broadly applicable to polarization management in photonic integrated circuits.

47 OTHER INSTRUMENTATION↗

IER-539 CED-1: Preliminary Design of a New Horizontal Split Table

This report presents the preliminary design (CED-1) of IER-539, focused on a new Horizontal Split Table (HST) critical assembly machine. The Department of Energy (DOE) National Nuclear Criticality Safety Program (NCSP) funds critical experiment R&D and its Mission and Vision document identifies a long-standing technical gap of a general purpose HST. HSTs have historically been used for experiments that are impractical to assemble on a vertical lift machine (VLM), usually due to their large physical size or weight. While many HSTs were historically employed in the US (including at LLNL), there are currently no surviving HSTs in the US. While there are two critical experiments facilities remaining in the US, the National Critical Experiments Research Center (NCERC) operated by LANL and the Critical Experiments Facility (CX) operated by SNL, these facilities are not equipped with a general purpose HST. There are a number of experiment designs that require an HST, including needs for Defense Programs and Nuclear Energy. Lack of an HST constraints the national ability to perform large and heavy experiments. Large footprints are needed to successfully create mockup experiments for advanced reactors (mock-up cold critical reactor cores, which is a urgent need for new reactor designs using novel materials and longer cycle lengths), to create solution experiments in lattices, or to test nuclear data of heavy metals. A new HST is a vital piece of criticality experiment infrastructure that would allow for the conduct of integral critical experiments to address needs in the nuclear criticality safety, nuclear data, and other DOE stakeholder communities.

42 ENGINEERING↗

Modelling detector-specific reconstruction uncertainties in LAr-TPC

The Short-Baseline Neutrino (SBN) program features three Liquid Argon Time Projection Chamber (LAr-TPC) detectors positioned along the Booster Neutrino Beam (BNB) axis: the Short Baseline Neutrino Near Detector, MicroBooNE, and the ICARUS T600. As the largest operational LAr-TPC, ICARUS T600 serves as the far detector, located 600 m from the BNB target. While its primary goal is to record neutrino events, it also detects other ionizing events, including cosmic rays. This work focuses on analyzing and modeling detector-specific reconstruction uncertainties in LAr-TPC. These inefficiencies, identified during the Pattern Recognition phase handled by the PANDORA algorithm, impact subsequent Particle Fits and Offline Analysis. Specifically, inaccuracies in track reconstruction can lead to significant physical consequences, such as erroneous particle energy estimates and poor Particle Identification (PID), reducing the efficiency of neutrino event characterization. A key issue addressed is split tracks, caused by missing hits or incomplete track stitching by PANDORA. The aim of this internship is to characterize, model, and quantify the impact of split tracks on track reconstruction.

43 PARTICLE ACCELERATORS↗

High accuracy computation of fluid-structure interaction in transonic cascades

A coupling strategy for simulating fluid-structure interaction phenomena is formulated and applied to the prediction of flutter in transonic cascades. The flow is governed by the Euler equations and discretized using a finite volume flux-splitting scheme. The structure is modeled using an isoparametric finite element formulation. The coupling strategy successfully reconciles these two formulations at the fluid-structure interface by enforcing both kinematic and kinetic boundary conditions. In particular, the conservation laws applicable to the combined fluid-structure system are preserved across the interface. Since the primary mechanism driving aeroelastic phenomena involves energy exchange occurring at the interface, this highly accurate coupling mechanism is believed to improve the predictive capability of the scheme. The coupled equations are advanced simultaneously in time using an implicit time integration method. Results obtained using the coupling method are presented for cascade geometries operating in transonic flow.

Boschitsch, Alexander↗

Fourth order difference methods for hyperbolic IBVP's

Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.

Gustafsson, Bertil↗

Energy Levels and Intensity Parameters of Ho3(+) Ions in Y3Al5O12 and Lu3Al5O12

The energy levels of the trivalent lanthanide Ho(sup 3+) in Y3Al5O12 (YAG) and Lu3Al5O12 (LuAG) have been measured. The Stark split levels for the first nine Ho manifolds in these materials have been measured, and the results have been fit to a free ion plus crystal field Hamiltonian to generate a theoretical set of energy levels. Crystal field parameters were varied to determine the best fit between experimental and theoretical energy levels. The energy levels of Ho:LuAG are seen to be very similar to those in Ho:YAG. However, subtle changes resulting from replacing Y(sup 3+) with Lu(sup 3+) in the garnet crystal Y3Al5O12 result in different transition wavelengths in LuAG. This has implications for Ho (sup 5)I7yields (sup 5)I8 lasers operating at approximately 2.1 micrometers. Although the energy levels have been measured previously in Ho:YAG, they have not been measured in Ho:LuAG. A comparison of the energy levels in Ho:YAG measured here show some discrepancies with previous measurements. The consistency of the energy level placement between Ho:LuAG and Ho:YAG indicate that the earlier studies may have some errors in the assignments. Finally, a Judd-Ofelt analysis is performed on Ho:YAG and Ho:LuAG to determine the intensity parameters, and thus, the transition probabilities and branching ratios of the first eight excited manifolds.

Walsh, Brian M.↗

Failure Modes of Platinized pn+-GaInP Photocathodes for Solar-Driven H2 Evolution

The long-term stability for the hydrogen-evolution reaction (HER) of homojunction pn+-Ga0.52In0.48P photocathodes (band gap = 1.8 eV) with an electrodeposited Pt catalyst (pn+-GaInP/Pt) has been systematically evaluated in both acidic and alkaline electrolytes. Electrode dissolution during chronoamperometry was correlated with changes over time in the current density-potential (J-E) behavior to reveal the underlying failure mechanism. Pristine pn+-GaInP/Pt photocathodes yielded an open-circuit photopotential (Eoc) as positive as >1.0 V vs the potential of the reversible hydrogen electrode (RHE) and a light-limited current density (Jph) of >12 mA cm-2 (1-sun illumination). However, Eoc and Jph gradually degraded at either pH 0 or pH 14. The performance degradation was attributed to three different failure modes: (1) gradual thinning of the n+-emitter layer due to GaInP dissolution in acid; (2) active corrosion of the underlying GaAs substrate at positive potentials causing delamination of the upper GaInP epilayers; and (3) direct GaAs/electrolyte contact compromising the operational stability of the device. This work reveals the importance of both substrate stability and structural integrity of integrated photoelectrodes toward stable solar fuel generation.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,S↗

SiC Composite Turbine Vanes

Turbine inlet guide vanes have been fabricated from composites of silicon carbide fibers in silicon carbide matrices. A unique design for a cloth made from SiC fibers makes it possible to realize the geometric features necessary to form these vanes in the same airfoil shapes as those of prior metal vanes. The fiber component of each of these vanes was made from SiC-fiber cloth coated with boron nitride. The matrix was formed by chemical-vapor infiltration with SiC, then slurry-casting of SiC, followed by melt infiltration with silicon. These SiC/SiC vanes were found to be capable of withstanding temperatures 400 F (222 C) greater than those that can be withstood by nickel-base-superalloy turbine airfoils now in common use in gas turbine engines. The higher temperature capability of SiC/SiC parts is expected to make it possible to use them with significantly less cooling than is used for metallic parts, thereby enabling engines to operate more efficiently while emitting smaller amounts of NOx and CO. The SiC/SiC composite vanes were fabricated in two different configurations. Each vane of one of the configurations has two internal cavities formed by a web between the suction and the pressure sides of the vane. Each vane of the other configuration has no web (see Figure 1). It is difficult to fabricate components having small radii, like those of the trailing edges of these vanes, by use of stiff stoichiometric SiC fibers currently preferred for SiC/SiC composites. To satisfy the severe geometric and structural requirements for these vanes, the aforementioned unique cloth design, denoted by the term Y-cloth, was conceived (see Figure 2). In the regions away from the trailing edge, the Y-cloth features a fiber architecture that had been well characterized and successfully demonstrated in combustor liners. To form a sharp trailing edge (having a radius of 0.3 mm), the cloth was split into two planes during the weaving process. The fiber tows forming the trailing-edge section were interlocked, thereby enhancing through-thickness strength of the resulting composite material. For vanes of the webless configuration, each made from a layup of six plies of Ycloth, the length of each Y-cloth layer was cut so that the two strips corresponding to the aforementioned two planes would wrap around the perimeter of a graphite vane preform tool with a 10-mm overlap. The overlap was used to join the two strips in a fringe splice. To make the external sixth ply, a standard woven cloth was cut to the required final length and a fringe splice joined the two ends of the cloth at the trailing edge. The cloth was then prepregged. The entire assembly was then placed into an aluminum compaction tool designed to form the outer net shape of the vane. After the prepreg material was allowed to dry, the preform was removed from the aluminum tooling and placed into an external graphite tool before being shipped to a vendor for matrix infiltration. To make the SiC fiber preform for a vane having an internal web, a slightly different initial approach was followed. Each of two sections forming the internal cavities (and ultimately the web) was created by first slipping two concentric layers of a two-dimensional, 2-by-2, 45 - braided tube around a net-shape graphite mandrel. The tubes on both mandrels were prepregged and allowed to dry. The resulting two subassemblies were put together, then four additional plies were wrapped around them in the same fashion as that described above for the six plies of the vaneless configuration. The consolidation of the SiC fiber preforms into SiC/SiC composite parts was performed by commercial vendors using their standard processes. The capability of two of the webless SiC/SiC turbine vanes was demonstrated in tests in a turbine environment. The tests included 50 hours of steady-state operation and 102 two-minute thermal cycles. A surface temperature of 1,320 C was reached during the tests.

Calomino, Anthony M.↗

TEA of a Unique Two-Pathways Process for Post-Combustion CO 2 Capture

A unique two-Pathways process using aqueous sodium glycinate for CO 2 capture from a split flue gas stream emitted from 600 MWe post-combustion coal power plant was developed in Aspen Plus v.10. The split gas flow rate used was 44.75 ton/h and contained 0.0023 mol% SO 2 and 13.33 mol% CO 2 . The process includes a washing unit, a CO 2 absorption unit, a reverse osmosis unit, and a solvent regeneration unit or an ultrafiltration unit. The washing unit uses deionized water to completely remove SO 2 and the CO 2 absorption unit uses SGS to capture at least 90 mol% of the CO 2 in the split flue gas stream. Upon CO 2 and SGS reactions, the resulting liquid products exhibit phase-separation into CO 2 -lean phase and CO 2 -rich phase, allow two distinct pathways. Pathway (i) is to regenerate mostly the CO 2 -rich phase, collect the released CO 2 , and compress it for sequestration purposes. Pathway (ii) is to send the liquid stream from the CO 2 absorption unit to the ultrafiltration unit to separate the solid nanomaterials. The hydraulics and mass transfer characteristics in the washing and CO 2 absorption units were obtained; and techno-economic analysis (TEA) for Pathways (i) and (ii), including Capital Expenditure (CAPEX), Operating Expenditure (OPEX), and Levelized Cost of CO 2 Captured (LCOC), were calculated and compared. The simulation results revealed that the CAPEX, OPEX, and LCOC for Pathway (i) were ($\$12,039,251$), (261 dollar/h), and (54.01 dollar/ton-CO 2 captured), respectively, and those for Pathway (ii) were ($\$5,908,000$), (237.2 dollar/h), and (39.90 dollar/ton-CO 2 captured), respectively. Moreover, in Pathway (ii), 8.19 ton/h of CO 2 were captured to produce 15.62 ton/h NaHCO 3 nanomaterials, which were sold to offset the overall process cost. The LCOC values indicate that Pathway (ii) is more cost-effective than Pathway (i) because LCOC values for Pathway (ii) are much lower than those for Pathway (i).COC values for Pathway (ii) are much lower than those for Pathway (i).

20 FOSSIL-FUELED POWER PLANTS↗

The Integrated Airport Competition Model, 1998

This paper addresses recent model development by the Directorate General of Civil Aviation (DGCA) and Hague Consulting Group (HCG) concerning long-distance travel. Long-distance travel demand is growing very quickly and raising a great deal of economic and policy issues. There is increasing competition among the main Western European airports, and smaller, regional airports are fighting for market share. New modes of transport, such as high speed rail, are also coming into the picture and affect the mode split for medium distance transport within Europe. Developments such as these are demanding the attention of policy makers and a tool is required for their analysis. For DGCA, Hague Consulting Group has developed a model system to provide answers to the policy questions posed by these expected trends, and to identify areas where policy makers can influence the traveller choices. The development of this model system, the Integrated Airport Competition Model/integraal Luchthaven Competitie Model (ILCM), began in 1992. Since that time the sub-models, input data and user interface have been expanded, updated and improved. HCG and DGCA have transformed the ILCM from a prototype into an operational forecasting tool.

Veldhuis, J.↗