Criteria

Text:
Affiliation:
Display:

Results

Viewing 1 to 30 of 121
2011-05-17
Technical Paper
2011-01-1512
Jaegon Yoo, Klaus Pfeiffer, Koo-Tae Kang
Advanced powertrain test, which is simulating real road load condition, was performed on the dynamic test bed. This cutting edge system can reproduce real road resistance based upon the vehicle dynamic model and wheel slip model. This wheel slip function is simulating the real behavior of the powertrain wheel as close as possible at each wheel independently. Additionally, low inertia of dynamometer motor themselves is another advantage for this purpose. This test bed is capable of testing all kinds of 2WD and 4WD powertrain configuration regardless of transmission type. Also, vehicle configuration can be mounted and tested on this test bed with small addition of supporting system alternatively. For the application, a four wheel drive powertrain was mounted on the test bed and driveline noise and vibration behavior such as transfer rattling noise and tip in/out shock were reproduced on this test bed.
2011-05-17
Technical Paper
2011-01-1526
Gianluca Montenegro, Augusto Della Torre, Angelo Onorati, Robert Fairbrother, Andreas Dolinar
The acoustic simulation of internal combustion engine exhaust systems is an important aspect to meet customer expectations and legislation targets. One dimensional gas dynamic simulation tools are used for the calculation of the exhaust orifice noise in the early stages of the engine development process. This includes the prediction of the acoustic performance of individual components in the exhaust line. One common element used in exhaust systems to increase the acoustic damping is the plug flow muffler. This study looks at the prediction of acoustic performance of various plug mufflers at different flow velocities. These include a single plug muffler, a double plug muffler and an eccentric plug muffler with different porosities for the perforated sections. To this purpose a generic 3D cell approach was developed and applied.
2011-04-12
Technical Paper
2011-01-1438
Johann C. Wurzenberger, Sophie Bardubitzki, Peter Bartsch, Tomaz Katrasnik
Modern Diesel engines require an integrated development of combustion strategies, air management and exhaust aftertreatment. This study presents a comprehensive simulation approach with the aim to support engine development activities in the virtual environment. A real-time capable engine, vehicle and control model is extended by three key features. First, a pollutant production model is embedded in a two-zone cylinder model. Second, a framework for catalytic pollutant conversion is built focusing on modern diesel exhaust aftertreatment systems. Third, an extended species transport model is introduced considering the transport of pollutants through the air path. The entire plant model is validated on the example of a passenger car Diesel engine. The predicted engine behavior is compared with steady-state measurements. The NO formation model is investigated for a series of steady-state and transient operating conditions.
2011-04-12
Technical Paper
2011-01-1439
Christoph Poetsch, Herwig Ofner, Eberhard Schutting
The paper describes a universally structured simulation platform which is used for the analysis and prediction of combustion in compression ignition (CI) engines. The models are on a zero-dimensional crank angle resolved basis as commonly used for engine cycle simulations. This platform represents a kind of thermodynamic framework which can be linked to single and multi zone combustion models. It is mainly used as work environment for the development and testing of new models which thereafter are implemented to other codes. One recent development task focused on a multi zone combustion model which corresponds to the approach of Hiroyasu. This model was taken from literature, extended with additional features described in this paper, and implemented into the thermodynamic simulation platform.
2011-04-12
Journal Article
2011-01-0395
Damien Le Guen, Thomas Weck, Adrien Balihe, Benoit Verbeke
Today's powertrains are becoming more and more complex due to the increasing number of gear box types requiring gearshift patterns like conventional (equipped with GSI) and automatic-manual transmissions (AT, AMT), double clutch and continuous variable transmissions (DCT, CVT). This increasing variety of gear boxes requires a higher effort for the overall optimization of the powertrain. At the same time, it is necessary to assess the impact of different powertrains and control strategies on CO₂ emissions very early in the development process. The optimization of Gear Shift Patterns (G.S.P.) has to fulfill multiple constraints in terms of objective customers' requirements, like driveability, NVH, performance, emissions and fuel consumption. For these reasons, RENAULT and AVL entered an engineering collaboration in order to develop a dedicated simulation tool: CRUISE GSP.
2011-04-12
Journal Article
2011-01-1219
Paul Whitaker, Paul Kapus, Martin Ogris, Peter Hollerer
Particulate emission reduction has long been a challenge for diesel engines as the diesel diffusion combustion process can generate high levels of soot which is one of the main constituents of particulate matter. Gasoline engines use a pre-mixed combustion process which produces negligible levels of soot, so particulate emissions have not been an issue for gasoline engines, particularly with modern port fuel injected (PFI) engines which provide excellent mixture quality. Future European and US emissions standards will include more stringent particulate limits for gasoline engines to protect against increases in airborne particulate levels due to the more widespread use of gasoline direct injection (GDI). While GDI engines are typically more efficient than PFI engines, they emit higher particulate levels, but still meet the current particulate standards.
2013-04-08
Journal Article
2013-01-1570
Jacob Swanson, David Kittelson, Barouch Giechaskiel, Alexander Bergmann, Martyn Twigg
The European Emissions Stage 5b standard for diesel passenger cars regulates particulate matter to 0.0045 g/km and non-volatile part/km greater than 23 nm size to 6.0x10₁₁ as determined by the PMP procedure that uses a heated evaporation tube to remove semi-volatile material. Measurement artifacts associated with the evaporation tube technique prevents reliable extension of the method to a lower size range. Catalytic stripper (CS) technology removes possible sources of these artifacts by effectively removing all hydrocarbons and sulfuric acid in the gas phase in order to avoid any chemical reactions or re-nucleation that may cause measurement complications. The performance of a miniature CS was evaluated and experimental results showed solid particle penetration was 50% at 10.5 nm. The sulfate storage capacity integrated into the CS enabled it to chemically remove sulfuric acid vapor rather than rely on dilution to prevent nucleation.
2013-04-08
Journal Article
2013-01-1563
Leonidas Ntziachristos, Stavros Amanatidis, Zissis Samaras, Barouch Giechaskiel, Alexander Bergmann
The Particle Measurement Programme (PMP) developed an exhaust particle number measurement protocol that has been adopted by current light duty vehicle emission regulations in Europe. This includes thermal treatment of the exhaust aerosol to isolate solid particles only and a number counting device with a lower cutpoint of 23 nm to avoid measurement of smaller particles that may affect the repeatability of the measurement. In this paper, we examine a potential alternative to the PMP system, where the thermal treatment is replaced by a catalytic stripper (CS). This offers oxidation and not just evaporation of the volatile components. Alternative sampling systems, either fulfilling the PMP recommendations or utilizing a CS, have been explored in terms of their volatile particle removal efficiency. Tests have been conducted on diesel exhaust, diesel equipped with DPF and gasoline direct injection emissions.
2014-04-01
Technical Paper
2014-01-1171
Aleš Kolar, Ralf Cerna, Werner Hofegger, Christoph Pichler, Markus Riener, Nathan Murphy, Georg Zembacher
The correct information about legal demands of the On-Board-Diagnostic (OBD) system in a vehicle project is required throughout the entire development process. Usually, the main obstacle in succeeding is to provide the company's expertise of some few experts for all employees who work in OBD related projects. The paper describes the AVL solution for knowledge management and tool supported control system design and calibration: OBD System Development Database. The software enables the user to access the regulatory requirements for a specific application and legislation from past, present and future (proposed rule-making) point of view. Information concerning already available and stored monitoring concepts is linked to the requirements in order to re-use potentially suitable concepts and to enable an efficient knowledge exchange within the company.
2014-04-01
Journal Article
2014-01-1094
Johann C. Wurzenberger, Tomaz Katrasnik
This works presents a real-time capable simulation model for dual fuel operated engines. The computational performance is reached by an optimized filling and emptying modeling approach applying tailored models for in-cylinder combustion and species transport in the gas path. The highly complex phenomena taking place during Diesel and gasoline type combustion are covered by explicit approaches supported by testbed data. The impact of the thermodynamic characteristics induced by the different fuels is described by an appropriate set of transport equations in combination with specifically prepared property databases. A thermodynamic highly accurate 6-species approach is presented. Additionally, a 3-species and a 1-species transport approach relying on the assumption of a lumped fuel are investigated regarding accuracy and computational performance. The comparison of measured and simulated pressure and temperature traces shows very good agreement.
2014-04-01
Technical Paper
2014-01-1095
Christoph Poetsch
Abstract The present works presents a real-time capable engine model with physical based description of the fuel injection and the combustion process. The model uses a crank-angle resolved cylinder model and a filling and emptying approach for cylinder and gas-path interaction. A common rail injection system model is developed and implemented into the real-time engine framework. The injection model calculates injection quantity and injection rate profile from the input of the ECU signals target injection pressure and injection timing. The model accounts for pressure oscillations in the injection system. A phenomenological combustion model for Diesel engines is implemented, which is based on the mixing controlled combustion modeling approach. The combustion model calculates the rate of heat release from the injection rate given by the injection model.
2013-09-08
Technical Paper
2013-24-0030
Christoph Poetsch, Henrik Schuemie, Herwig Ofner, Reinhard Tatschl, Oldrich Vitek
In spark-ignition engines, fluctuations of the in-cylinder pressure trace and the apparent rate of heat release are usually observed from one cycle to another. These Cycle-to-Cycle Variations (CCV) are affected by the early flame development and the subsequent flame front propagation. The CCV are responsible for engine performance (e.g. fuel consumption) and the knock behavior. The occurrence of the phenomena is unpredictable and the stochastic nature offers challenges in the optimization of engine control strategies. In the present work, CCV are analyzed in terms of their impact on the engine knock behavior and the related efficiency. Target is to estimate the possible fuel consumption savings in steady-state operation and in the drivecycle, when CCV are reduced. Since CCV are immanent on real engines, such a study can only be done by means of simulation.
2014-04-01
Technical Paper
2014-01-0668
Armin Traussnig, Heinz Petutschnig, Andreas Ennemoser, Michael Stolz, Mauro Tizianel
Abstract In order to meet current and future emission and CO2 targets, an efficient vehicle thermal management system is one of the key factors in conventional as well as in electrified powertrains. Furthermore the increasing number of vehicle configurations leads to a high variability and degrees of freedom in possible system designs and the control thereof, which can only be handled by a comprehensive tool chain of vehicle system simulation and a generic control system architecture. The required model must comprise all relevant systems of the vehicle (control functionality, cooling system, lubrication system, engine, drive train, HV components etc.). For proper prediction with respect to energy consumption all interactions and interdependencies of those systems have to be taken into consideration, i.e. all energy fluxes (mechanical, hydraulically, electrical, thermal) have to be exchanged among the system boundaries accordingly.
2014-04-01
Journal Article
2014-01-0641
Rok Kopun, Dongsheng Zhang, Wilfried Edelbauer, Bernhard Stauder, Branislav Basara, David Greif
In this paper, a recently improved Computational Fluid Dynamics (CFD) methodology for virtual prototyping of the heat treatment of cast aluminum parts, above most of cylinder heads of internal combustion engines (ICE), is presented. The comparison between measurement data and numerical results has been carried out to simulate the real time immersion quenching cooling process of realistic cylinder head structure using the commercial CFD code AVL FIRE®. The Eulerian multi-fluid modeling approach is used to handle the boiling flow and the heat transfer between the heated structure and the sub-cooled liquid. While for the fluid region governing equations are solved for each phase separately, only the energy equation is solved in the solid region. Heat transfer coefficients depend on the boiling regimes which are separated by the Leidenfrost temperature.
2015-04-14
Technical Paper
2015-01-0932
Satoru Sasaki, Masaaki Kato, Takamasa Yokota, Mitsuru Konno, Denis Gill
Abstract DiMethyl Ether (DME) has been known to be an outstanding fuel for combustion in diesel cycle engines for nearly twenty years. DME has a vapour pressure of approximately 0.5MPa at ambient temperature (293K), thus it requires pressurized fuel systems to keep it in liquid state which are similar to those for Liquefied Petroleum Gas (mixtures of propane and butane). The high vapour pressure of DME permits the possibility to optimize the fuel injection characteristic of direct injection diesel engines in order to achieve a fast evaporation and mixing with the charged gas in the combustion chamber, even at moderate fuel injection pressures. To understand the interrelation between the fuel flow inside the nozzle spray holes tests were carried out using 2D optically accessed nozzles coupled with modelling approaches for the fuel flow, cavitation, evaporation and the gas dynamics of 2-phase (liquid and gas) flows.
2015-04-14
Technical Paper
2015-01-1197
Chao Chen, Franz Diwoky, Zoran Pavlovic, Johann Wurzenberger
This paper presents a system-level thermal model of a fluid-cooled Li-Ion battery module. The model is a reduced order model (ROM) identified by results from finite element analysis (FEA)/computational fluid dynamic (CFD) coupling simulation using the linear and time-invariant (LTI) method. The ROM consists of two LTI sub-systems: one of which describes the battery temperature response to a transient battery current, and the other of which takes into account of the battery temperature variation due to a heat flux induced by a varied inlet temperature of the battery cooling circuit. The thermal LTI model can be coupled to an electrical model to build a complete system-level battery ROM. Test examples show that the ROM is able to provide as accurate results as those from FEA/CFD coupling simulations.
2015-04-14
Technical Paper
2015-01-1275
Wolfgang Johann Schoeffmann, Caroline Truffinet, Michael Howlett, Norbert Ausserhofer, Andreas Zurk
Abstract When designing a new internal combustion engine, the choice of technology for the timing drive system is one of the key decisions that determines the overall characteristics of the engine with far reaching implications on the remaining architecture and overall packaging of the engine. For Passenger car engines there are two mainstream technologies: toothed belts and chains. Each of these offers several sub-variants, such as dry vs. wet belt, or toothed vs. roller chain. This paper examines the differences between these technologies in relation to the key engine attributes including package, cost, weight, durability, NVH and frictional losses. A quantitative evaluation is made where possible, based on data collected from recent engine development programs, backed up by literature study and data from the component supply industry.
2014-10-13
Journal Article
2014-01-2848
Matthew McAllister, Stephen Smith, Paul Kapus, Khai Vidmar, Alexander Hochnetz
This paper describes the findings of a design, simulation and test study into how to reduce particulate number (Pn) emissions in order to meet EU6c legislative limits. The objective of the study was to evaluate the Pn potential of a modern 6-cylinder engine with respect to hardware and calibration when fitted to a full size SUV. Having understood this capability, to redesign the combustion system and optimise the calibration in order to meet an engineering target value of 3×1011 Pn #/km using the NEDC drive cycle. The design and simulation tasks were conducted by JLR with support from AVL. The calibration and all of the vehicle testing was conducted by AVL, in Graz. Extensive design and CFD work was conducted to refine the inlet port, piston crown and injector spray pattern in order to reduce surface wetting and improve air to fuel mixing homogeneity. The design and CFD steps are detailed along with the results compared to target.
2015-04-14
Technical Paper
2015-01-0148
Georg Macher, Harald Sporer, Eric Armengaud, Christian Kreiner
Abstract Increasing demands for safety, security, and certifiability of embedded automotive systems require additional development effort to generate the required evidences that the developed system can be trusted for the application and environment it is intended for. Safety standards such as ISO 26262 for road vehicles have been established to provide guidance during the development of safety-critical systems. The challenge in this context is to provide evidence of consistency, correctness, and completeness of system specifications over different work-products. One of these required work-products is the hardware-software interface (HSI) definition. This work-product is especially important since it defines the interfaces between different technologies. Model-based development (MBD) is a promising approach to support the description of the system under development in a more structured way, thus improving resulting consistency.
2010-04-12
Technical Paper
2010-01-1202
Sasa Bukovnik, Tigran Parikyan, Thomas Resch, Wolfgang Hellinger
Mechanical variable valve systems are being increasingly used for modern combustion engines. It is typical for such systems that the cam and valve are connected via intermediate levers. Different maximum valve lifts and duration can be achieved with the same cam profile. The intermediate levers increase the system inertia and reduce the overall stiffness. Such systems offer more flexibility, but it is more complex to create optimal design compared to the conventional systems. In this paper a new kinematic design methodology for a CVVL (Continuously Variable Valve Lift) system is presented. Additionally, dynamic analysis of the valve train system is performed. The investigated valve train is completely developed and patented by OEM. The main characteristic of the CVVL system is a set of intermediate levers between the cam and the finger follower like ( 1 , 2 ). One cam drives two intake valves over a set of levers.
2010-04-12
Technical Paper
2010-01-0784
Johann C. Wurzenberger, Peter Bartsch, Tomaz Katrasnik
The present work introduces a fully integrated real-time (RT) capable engine and vehicle model. The gas path and drive line are described in the time domain of seconds whereas the reciprocating characteristics of an IC engine are reflected by a crank angle resolved cylinder model. The RT engine model is derived from a high fidelity 1D cycle simulation and gas exchange model to support an efficient and consistent transfer of model data like geometries, heat transfer or combustion. The workflow of model calibration and application is outlined and base ECU functionalities for boost pressure, EGR, smoke and idle speed control are applied for transient engine operation. Steady state results of the RT engine model are compared to experimental data and 1D high fidelity simulations for 19 different engine load points. In addition an NEDC (New European Drive Cycle) is simulated and results are evaluated with data from chassis dynamometer measurements.
2010-04-12
Technical Paper
2010-01-0329
Manuel Nebel, Marie-Sophie Vogels, Timo Combe, Thomas Winsel, Horst Pfluegl, Christoph Hametner
The modern power train calibration process is characterized by shorter development cycles and a reduced number of prototypes. However, simultaneously exhaust aftertreatment and emission testing is becoming increasingly more sophisticated. The introduction of predictive simulation tools that represent the complete power train can likely contribute to improving the efficiency of the calibration process using an integral model based workflow. Engine models, which are purely based on complex physical principles, are usually not capable of real-time applications, especially if the simulation is focused on transient emission optimization. Methods, structures and the realization of a global dynamic real-time model are presented in this paper, combining physical knowledge and experimental models and also static and dynamic sub-structures. Such a model, with physical a priori information embedded in the model structure, provides excellent generalization capability.
2010-04-12
Technical Paper
2010-01-0331
Henry Joswa Rajan, Joseph Kelly, Harald Hoetzendorfer, Nikolaus Keuth, Horst Pfluegl, Thomas Winsel, Siegfried Roeck
Today's calibration process for ECU functions is often based on a wide variety of proprietary tools and individual expert knowledge of calibration engineers. Automatic calibration with an industrialized tool chain provides high potential to reduce testbed time, calibration time and project costs. Based on an efficient measurement procedure in combination with an offline calibration methodology the capability is validated, e.g. for calibrating the ECU function “Air Charge Determination” for SI engines. In this article the implementation, in a series production project of a major OEM, is shown. The whole workflow - which can also be applied to other calibration tasks - will be described in detail. Presented here will be how General Motors Corporation (GM) is able to speed up the calibration of the ECU functions, whilst maintaining at least the same quality of calibration as before, by the use of this tool chain.
2014-10-13
Technical Paper
2014-01-2658
Denis W. Gill, Herwig Ofner, Carsten Stoewe, Karl Wieser, Ernst Winklhofer, Masaaki Kato, Takamasa Yokota, Jost Weber
Abstract For nearly twenty years, DiMethyl Ether has been known to be an outstanding fuel for combustion in diesel cycle engines. Not only does it have a high Cetane number, it burns absolutely soot free and produces lower NOx exhaust emissions than the equivalent diesel. However, the physical properties of DME such as its low viscosity, lubricity and bulk modulus have negative effects for the fuel injection system, which have both limited the achievable injection pressures to about 500 bar and DME's introduction into the market. To overcome some of these effects, a common rail fuel injection system was adapted to operate with DME and produce injection pressures of up to 1000 bar. To understand the effect of the high injection pressure, tests were carried out using 2D optically accessed nozzles. This allowed the impact of the high vapour pressure of DME on the onset of cavitation in the nozzle hole to be assessed and improve the flow characteristics.
2014-10-13
Technical Paper
2014-01-2690
Werner E. Holly, Thomas Lauer, Henrik Schuemie, Shinsuke Murakami
Abstract The combustion efficiency of large gas engines is limited by knocking combustion. Due to fact that the quality of the fuel gas has a high impact on the self-ignition of the mixture, it is the aim of this work to model the knocking combustion for fuel gases with different composition using detailed chemistry. A cycle-resolved knock simulation of the fast burning cycles was carried out in order to assume realistic temperatures and pressures in the unburned mixture Therefore, an empirical model that predicts the cyclic variations on the basis of turbulent and chemical time scales was derived from measured burn rates and implemented in a 1D simulation model. Based on the simulation of the fast burning engine cycles the self-ignition process of the unburned zone was calculated with a stochastic reactor model and correlated to measurements from the engines test bench. A good agreement of the knock onset could be achieved with this approach.
2014-10-13
Technical Paper
2014-01-2626
Timothy J. Jacobs, Louis Camilli, Matthias Neubauer
Abstract A key element to achieving vehicle emission certification for most light-duty vehicles using spark-ignition engine technology is prompt catalyst warming. Emission mitigation largely does not occur while the catalyst is below its “light-off temperature”, which takes a certain time to achieve when the engine starts from a cold condition. If the catalyst takes too long to light-off, the vehicle could fail its emission certification; it is necessary to minimize the catalyst warm up period to mitigate emissions as quickly as possible. One technique used to minimize catalyst warm up is to calibrate the engine in such a way that it delivers high temperature exhaust. At idle or low speed/low-load conditions, this can be done by retarding spark timing with a corresponding increase in fuel flow rate and / or leaning the mixture. Both approaches, however, encounter limits as combustion stability degrades and / or nitrogen oxide emissions rise excessively.
2014-06-30
Technical Paper
2014-01-2067
Michael Klanner, Mathias Mair, Franz Diwoky, Oszkar Biro, Katrin Ellermann
Abstract The noise vibration and harshness (NVH) simulation of electric machines becomes increasingly important due to the use of electric machines in vehicles. This paper describes a method to reduce the calculation time and required memory of the finite element NVH simulation of electrical machines. The stator of a synchronous electrical machine is modeled as a two-dimensional problem to reduce investigation effort. The electromagnetic forces acting on the stator are determined by FE-simulation in advance. Since these forces need to be transferred from the electromagnetic model to the structural model, a coupling algorithm is necessary. In order to reduce the number of nodes, which are involved in the coupling between the electromagnetic and structural model, multipoint constraints (MPC) are used to connect several coupling nodes to one new coupling node. For the definition of the new coupling nodes, the acting load is analyzed with a 2D-FFT.
2014-06-30
Technical Paper
2014-01-2093
Vishal Parmar, Daniele Di Rocco, Martin Sopouch, Philippe Albertini
Abstract Over the past 30 years, simulation of the N&V (Noise and Vibration) behaviour of automotive drivelines became an integral part of the powertrain development process. With current and future HEVs (Hybrid-Electrical Vehicles), additional phenomena and effects have entered the scene and need to be taken into account during layout/design as well as optimization phase. Beside effects directly associated with the e-components (namely electric whistle and whine), torque changes caused by activation/deactivation of the e-machine give rise to vibration issues (e.g. driveline shuffle or clonk) as well. This is in particular true for transient operation conditions like boosting and recuperation. Moreover, aspects of starting the Internal Combustion Engine (ICE) using the built-in e-machine in conjunction with the dynamic behaviour of torsional decoupling devices become increasingly important. In order to cope with above-mentioned effects a multi-physics simulation approach is required.
2014-04-01
Technical Paper
2014-01-1421
David Greif, Wilfried Edelbauer, Jure Strucl
Abstract The paper addresses aspects of modeling cavitating flows within high pressure injection equipment while considering the effects of liquid compressibility. The presented numerical study, performed using the commercial CFD code AVL FIRE®, mimics common rail conditions, where the variation in liquid density as a function of pressure may be relevant owing to very high pressure injection scenarios. The flow through the injector has been calculated and the conditions at the outlet of the nozzle orifice have been applied as inlet condition for subsequent Euler-Lagrangian spray calculations to investigate the effects of liquid compressibility treatment on spray propagation. Flows of such nature are of interest within automotive and other internal combustion (IC) related industries to obtain good spray and emissions characteristics.
2014-06-30
Technical Paper
2014-01-2041
Alexander Ulz, Alfred Rust, Bernhard Graf, Alois Sontacchi
Abstract Due to future directives of the European Union regarding fuel consumption and CO2 emissions the automotive industry is forced to develop new and unconventional technologies. These include for example stop-start-systems, cylinder deactivation or even reduction of the number of cylinders which however lead to unusual acoustical perceptions and customer complaints. Therefore, it is necessary to evaluate the sound character of engines with low numbers of cylinders (2 and 3 cylinders) and also the differences to the character of the more common 4-cylinder engines. Psychoacoustic parameters are used to describe and understand the differences. Based on the gained knowledge possible potentials for improvement can be derived in the future. The used data base consists of artificial head recordings of car interior noise according to defined driving conditions measured on the AVL test track.
Viewing 1 to 30 of 121

Filter

  • Range:
    to:
  • Year: