Environmental simulation

Environmental Simulation for AEC Computational Design

This skill provides a comprehensive reference for environmental performance simulation in the Architecture, Engineering, and Construction industry. It covers daylight analysis, solar radiation studies, energy modeling, computational fluid dynamics for wind, thermal comfort assessment, acoustic simulation, and the Ladybug Tools ecosystem that ties these workflows together within parametric design environments.

---

1. Performance-Driven Design Philosophy

Why Simulate Early

The single most consequential decision in building performance is made in the first five percent of the design timeline. Orientation, massing, window-to-wall ratio, and floor plate depth are locked in during concept design, yet their impact on energy consumption, daylight quality, and occupant comfort persists for the entire operational life of the building — typically 50 to 100 years.

Late-stage simulation is a diagnostic exercise. Early-stage simulation is a generative tool. The distinction matters because retrofitting performance into an already-resolved form is orders of magnitude more expensive than shaping the form around performance from the beginning.

Cost of late performance analysis:

Stage where issue discovered	Relative cost to fix
Concept design	1x (baseline)
Schematic design	5x
Design development	15x
Construction documents	50x
Construction	150x
Post-occupancy	500x+

These multipliers are well-documented across the AEC industry and reflect the reality that changing a window-to-wall ratio on a sketch costs nothing, but changing it after curtain wall shop drawings are issued can cost hundreds of thousands of dollars.

Integration of Simulation into the Design Loop

Performance simulation must not live in a separate silo from the design model. The traditional workflow — architect designs, sends model to energy consultant, waits two weeks, receives PDF report, ignores half the recommendations because they conflict with the architectural intent — is fundamentally broken.

The correct workflow is a tight feedback loop:

1. Parametric model defines geometry with variable parameters (orientation, WWR, floor depth, shading depth, envelope composition).
2. Simulation engine evaluates one or more performance metrics against that geometry.
3. Results feed back into the model as color-mapped surfaces, numeric dashboards, or fitness values for optimization.
4. Designer adjusts parameters or launches an automated optimization run.
5. Repeat at increasing fidelity as the design matures.

This loop should complete in seconds to minutes during early design and in minutes to hours during detailed design. Any simulation that takes days to return results is, by definition, excluded from the design loop and relegated to post-rationalization.

Performance as a Design Driver

Performance-driven design does not mean surrendering architectural intent to a spreadsheet. It means treating environmental performance as a first-class design variable alongside spatial quality, structural logic, and aesthetic expression. The best performing buildings in the world — the Bullitt Center, One Angel Court, the Manitoba Hydro Place — are also among the most architecturally compelling because their designers used performance constraints as creative catalysts.

Metrics that can drive design decisions:

• Daylight autonomy drives floor plate depth, section profile, and facade design.
• Solar radiation drives orientation, massing articulation, and shading strategy.
• Energy use intensity drives envelope specification, HVAC selection, and renewable capacity.
• Wind comfort drives podium design, canopy placement, and landscape strategy.
• Thermal comfort drives material selection, ventilation strategy, and public realm design.
• Acoustic performance drives room proportions, material specification, and partition layout.

Iterative vs. Single-Pass Simulation Strategies

Single-pass simulation runs one analysis on a fixed design. It answers the question "how does this design perform?" This is the minimum viable use of simulation and is appropriate only for compliance checking at the end of a design phase.

Iterative simulation runs multiple analyses across a parameter space. It answers the question "which design performs best?" This is the correct use of simulation in a computational design workflow.

Iterative strategies include:

• Manual parameter sweeps: Designer changes one variable at a time and observes the effect on performance. Simple but slow and unable to capture multi-variable interactions.
• Parametric studies: Systematic variation of two or three parameters across defined ranges, producing a response surface that reveals optimal regions.
• Evolutionary optimization: Genetic algorithms (Galapagos, Wallacei) explore a high-dimensional parameter space, evolving toward Pareto-optimal solutions across multiple conflicting objectives.
• Machine learning surrogates: A neural network is trained on a sample of simulation results and then used to predict performance across the full parameter space at near zero computational cost. This enables real-time performance feedback during design.

---

2. Daylight Analysis

Core Metrics

Daylight Factor (DF) The ratio of indoor illuminance to outdoor diffuse horizontal illuminance under a CIE overcast sky, expressed as a percentage. DF is climate-independent and therefore useful for comparing designs but not for predicting actual illuminance.

• Minimum DF of 2% is generally required for a space to be considered "daylit."
• Average DF of 5% or more indicates strong daylight availability.
• DF does not account for direct sunlight, building orientation, or climate.

Spatial Daylight Autonomy (sDA) The percentage of analysis area that achieves at least 300 lux for at least 50% of occupied hours annually. sDA is the primary metric in LEED v4.1 and IES LM-83.

• sDA300/50% >= 55% is the LEED threshold for 2 points.
• sDA300/50% >= 75% earns 3 points.
• Based on annual simulation with climate-specific weather data.
• Accounts for orientation, context, and dynamic shading.

Annual Sunlight Exposure (ASE) The percentage of analysis area that receives more than 1000 lux of direct sunlight for more than 250 occupied hours per year. ASE measures visual discomfort risk from excessive direct sun.

• ASE1000/250h must be <= 10% for LEED compliance.
• High ASE indicates glare risk and potential overheating.
• ASE is often in tension with sDA — maximizing daylight can increase glare risk.

Useful Daylight Illuminance (UDI) The percentage of occupied hours that illuminance falls within a useful range:

UDI category	Illuminance range	Interpretation
UDI-fell-short	< 100 lux	Too dark, electric light needed
UDI-supplementary	100–300 lux	Useful with supplemental light
UDI-autonomous	300–3000 lux	Ideal daylight range
UDI-exceeded	> 3000 lux	Too bright, glare likely

UDI is more nuanced than sDA because it penalizes both under-lit and over-lit conditions.

Daylight Glare Probability (DGP) A metric for evaluating glare from a specific viewpoint, based on luminance distribution in the field of view.

DGP value	Glare perception
< 0.35	Imperceptible
0.35–0.40	Perceptible
0.40–0.45	Disturbing
> 0.45	Intolerable

DGP is evaluated using Evalglare from the Radiance suite and requires a rendered fisheye luminance image from the occupant's viewpoint.

Standards Reference

LEED v4.1 IEQ Credit: Daylight

• Option 1 (Simulation): sDA300/50% >= 55% in 55% of regularly occupied area (2 pts),

  = 75% (3 pts). ASE1000/250h <= 10%. Analysis grid at workplane height (0.76 m).

• Option 2 (Measurement): Illuminance between 300–3000 lux at 9 AM and 3 PM on equinox.

EN 17037: Daylight in Buildings

• Target illuminance of 300 lux for 50% of daylight hours at the worst point.
• Minimum illuminance of 100 lux for 50% of daylight hours across 95% of area.
• Recommends evaluation of view out, sunlight exposure, and glare protection.

BREEAM Hea 01: Visual Comfort

• Average daylight factor >= 2% or point daylight factor >= 0.8% at the worst point.
• Uniformity ratio (minimum DF / average DF) >= 0.4.
• View of sky from desk height (requirement for higher credits).

Simulation Engines

Radiance The gold standard for physically-based lighting simulation. Uses backward ray-tracing to compute illuminance and luminance with high accuracy. All major daylight standards reference Radiance as a validated engine.

Key Radiance parameters for annual simulation:

Parameter	Low quality	Medium quality	High quality
-ab (bounces)	2	3	5
-ad (divisions)	512	2048	4096
-as (super-samples)	128	1024	2048
-ar (resolution)	64	128	256
-aa (accuracy)	0.25	0.15	0.1

3-Phase Method Decomposes daylight transport into three matrices:

1. View matrix (V): sensor points to interior surface patches
2. Transmission matrix (T): window group behavior
3. Daylight matrix (D): sky patches to exterior window surfaces

Annual daylight = V * T * D * sky_vector. This allows rapid recalculation when only the window system changes (swap T matrix) without recomputing V or D.

5-Phase Method Adds direct solar contribution more accurately by separating the direct sun component from the diffuse sky component and computing it with higher resolution. Essential for accurate ASE calculations and glare analysis.

Material Properties

Common material reflectances for Radiance modeling:

Material	Reflectance	Specularity	Roughness
White painted wall	0.70	0.0	0.0
Light grey painted wall	0.50	0.0	0.0
Dark grey painted wall	0.20	0.0	0.0
Concrete (raw)	0.30	0.0	0.05
Wood flooring (light)	0.40	0.02	0.02
Carpet (medium)	0.20	0.0	0.0
Ceiling tile (white)	0.80	0.0	0.0
Glass (clear single)	Tvis: 0.88	—	—
Glass (clear double)	Tvis: 0.78	—	—
Glass (low-e double)	Tvis: 0.65	—	—
Ground (grass)	0.20	0.0	0.0
Ground (asphalt)	0.10	0.0	0.0
Ground (concrete paving)	0.30	0.0	0.0
Red brick	0.25	0.0	0.03
Aluminum (brushed)	0.60	0.50	0.05

Grid-Based vs. Point-in-Time Analysis

Grid-based analysis places sensors at regular intervals across a horizontal workplane (typically 0.76 m above floor for offices, 0.85 m for standing desks). Grid spacing should be no larger than 0.5 m for final analysis and no larger than 1.0 m for early design studies. Results are spatial maps showing illuminance distribution.

Point-in-time analysis evaluates illuminance or luminance at a specific moment (date, time, sky condition). Useful for worst-case glare checks (e.g., December 21 at 3 PM with low sun angle) or for visualizing light distribution at critical moments.

Parametric Daylight Optimization Workflow

1. Define variable parameters: window-to-wall ratio (20–80%), window sill height, head height, external shading depth (0–1.5 m), light shelf depth, room depth.
2. Set up Honeybee model with Radiance modifiers for all surfaces.
3. Run annual daylight simulation (Honeybee-Radiance recipe).
4. Extract sDA, ASE, and UDI from results.
5. Feed metrics into multi-objective optimizer (Wallacei recommended).
6. Objectives: maximize sDA, minimize ASE, minimize facade cost.
7. Analyze Pareto front to select design variants that balance performance and cost.

Common Modeling Mistakes

• Missing context buildings: Omitting neighboring buildings that cast shadows leads to dramatically overstated daylight predictions. Always include context within at least 200 m radius for urban sites.
• Wrong material reflectance: Using default white (0.50) for all surfaces when actual materials are much darker (0.20–0.30) can overpredict daylight by 30–50%.
• Insufficient grid resolution: A 2 m grid misses spatial variation entirely. Use 0.5 m maximum for compliance analysis.
• Ignoring furniture: Large obstructions like tall storage units significantly affect daylight distribution in deep plans.
• Wrong analysis period: Using calendar year instead of occupied hours misrepresents sDA values. Always define occupancy schedule correctly.
• Not modeling window frames: Frames reduce glazed area by 15–30%. Omitting them overpredicts daylight proportionally.

---

3. Solar Radiation Analysis

Solar Geometry

The position of the sun is defined by two angles:

• Altitude (elevation): Angle above the horizon (0° at horizon, 90° at zenith).
• Azimuth: Angle measured clockwise from true north (0° = N, 90° = E, 180° = S, 270° = W).

Solar position depends on latitude, longitude, date, and time. Key concepts:

• Solar noon: The moment the sun crosses the local meridian (highest altitude).
• Equation of time: Correction for the difference between solar time and clock time due to Earth's elliptical orbit and axial tilt. Varies by +/- 16 minutes.
• Declination: The angle between the sun's rays and the equatorial plane. Ranges from +23.45° (summer solstice) to -23.45° (winter solstice).
• Hour angle: Angular displacement of the sun from solar noon. 15° per hour.

Solar Radiation Components

Component	Source	Behavior	Typical share (annual)
Direct beam	Sun disk	Casts sharp shadows	50–70%
Diffuse	Sky vault (scattered by atmosphere)	Uniform, no shadows	25–40%
Reflected	Ground and surrounding surfaces	Depends on albedo	5–15%

Total solar radiation (global) = direct + diffuse + reflected.

Radiation Analysis Types

Cumulative annual radiation: Total solar energy received by a surface over an entire year, measured in kWh/m². This is the most common analysis for facade studies and PV potential assessment. Typical values for a horizontal surface range from 900 kWh/m²/yr (northern Europe) to 2200 kWh/m²/yr (desert regions).

Monthly radiation: Breakdown of annual radiation by month, revealing seasonal patterns. Critical for understanding heating vs. cooling season dynamics.

Hourly radiation: Instantaneous radiation values for specific hours, used for peak load calculations and shadow studies.

Solar Access Hours

Solar access analysis counts the number of hours per year (or per day for a specific date) that a point receives direct sunlight. Used for:

• Right-to-light assessments (UK planning requirement: 2+ hours on March 21).
• BRE Site Layout Planning guidelines (25+ hours on March 21 for outdoor amenity).
• Passive solar design (minimum 4 hours winter solstice for south-facing facades).

Shadow Range Analysis

Shadow range diagrams overlay all shadow positions for a specific date or period, showing which areas are always in shade, always in sun, or intermittently shaded. Useful for:

• Public space design (ensuring benches get winter sun).
• Urban planning (assessing impact of new tall buildings on surrounding neighborhoods).
• PV array siting (avoiding partially shaded zones).

Radiation on Tilted Surfaces

Radiation on a tilted surface differs from horizontal radiation. The optimal tilt angle for annual energy collection is approximately equal to the site latitude. For winter-optimized collection, add 15° to latitude. For summer-optimized, subtract 15°.

Isotropic sky model: Diffuse radiation is uniform across the sky dome. Simple but inaccurate — underestimates radiation near the horizon and circumsolar region.

Perez all-weather model: Divides the sky into circumsolar, horizon brightening, and isotropic components. Most accurate model for tilted surface calculations.

PV Potential Assessment

1. Calculate annual radiation on candidate surfaces (roof, facade, canopy).
2. Subtract shading losses from context geometry and self-shading.
3. Apply panel efficiency (monocrystalline: 20–22%, polycrystalline: 16–18%, thin-film: 10–13%).
4. Apply system losses (inverter: 3–5%, wiring: 2–3%, soiling: 2–5%, temperature: 5–10%).
5. Result: annual energy yield in kWh per panel or per m².

Typical yield: 150–200 kWh/m²/yr in northern Europe, 250–350 kWh/m²/yr in southern US and Middle East.

Sky Models

Sky model	Use case	Characteristics
CIE overcast	Daylight factor calculation	Luminance varies only with altitude
CIE clear	Sunny day analysis	Includes circumsolar and horizon zones
CIE intermediate	Partly cloudy conditions	Blend of clear and overcast
Perez all-weather	Annual simulation	Uses weather data, most accurate
Uniform	Simple diffuse analysis	Equal luminance everywhere (unrealistic)

EPW Weather Files

EPW (EnergyPlus Weather) files contain hourly data for a typical meteorological year:

• Dry bulb temperature, dew point temperature, relative humidity
• Direct normal radiation, diffuse horizontal radiation, global horizontal radiation
• Wind speed, wind direction
• Atmospheric pressure, sky cover, precipitation

Sources: EnergyPlus.net, Climate.OneBuilding.org (most comprehensive), ASHRAE IWEC2, Meteonorm (commercial, can generate for any location).

TMY methodology: Typical Meteorological Year data is constructed by selecting the most representative month from a 15–30 year record for each calendar month. TMY data represents typical conditions, not extreme conditions. For resilience analysis, use actual year data or future climate projections.

---

4. Energy Simulation

Building Energy Model Components

A complete building energy model requires the following layers:

1. Geometry: Thermal zones defined by enclosed volumes. Each zone has a single air temperature assumption. Zones should be separated where significantly different thermal conditions exist (perimeter vs. core, different orientations, different uses).

2. Constructions: Material layers for each opaque surface (walls, roofs, floors) and glazing properties for windows. Each layer defined by thickness, conductivity, density, and specific heat.

3. Schedules: Hourly profiles for occupancy, lighting, equipment, thermostat setpoints, HVAC availability, and ventilation rates. Schedules are the single most influential input in energy modeling after geometry.

4. Internal loads: Heat gains from people (sensible + latent), lighting, and equipment. Specified as peak values modulated by schedules.

5. HVAC systems: Heating, cooling, and ventilation equipment. Ranges from ideal air (unlimited capacity, perfect efficiency — for early design) to fully detailed systems with specific equipment curves.

6. Infiltration: Uncontrolled air leakage through the envelope. Specified as ACH (air changes per hour) or flow per unit envelope area.

7. Natural ventilation: Window opening behavior, stack effect, wind-driven ventilation. Can be modeled as scheduled or as airflow network.

EnergyPlus Engine Fundamentals

EnergyPlus is a whole-building energy simulation engine developed by the US Department of Energy. It uses a heat-balance method that simultaneously solves:

• Conduction through opaque surfaces (conduction transfer functions)
• Solar radiation through windows (detailed optical calculations)
• Internal convection (TARP algorithm)
• Longwave radiation exchange between surfaces
• Air heat balance (zone energy balance)
• HVAC system response

EnergyPlus runs sub-hourly timesteps (typically 6 per hour = 10-minute intervals) with weather data interpolated from hourly EPW values.

Honeybee-Energy Workflow

1. Create rooms: Define thermal zone geometry from Rhino surfaces or procedural generation. Each room = one thermal zone.
2. Assign program types: Predefined combinations of loads, schedules, and setpoints for common space types (office, residential, retail, etc.).
3. Set construction sets: Climate-appropriate envelope specifications. Honeybee includes ASHRAE baseline construction sets by climate zone.
4. Define HVAC: Start with IdealAirSystem for early design. Progress to detailed systems (VAV, VRF, DOAS, radiant) for later stages.
5. Set simulation parameters: Run period, timestep, terrain, solar distribution.
6. Run simulation: Honeybee generates IDF file and calls EnergyPlus.
7. Parse results: Monthly/annual energy by end use, zone temperatures, comfort hours.

Key Outputs

Energy Use Intensity (EUI) Total annual energy consumption divided by gross floor area, measured in kWh/m²/yr (or kBtu/ft²/yr in US practice). EUI is the primary benchmark for building energy performance.

Building type	Good EUI (kWh/m²/yr)	Average EUI	Poor EUI
Office	80–120	150–200	250+
Residential (apt)	50–80	100–140	180+
Retail	100–150	200–280	350+
Education	70–110	130–170	220+
Hospital	200–300	350–450	550+
Hotel	120–180	220–300	400+
Laboratory	250–400	500–700	900+
Warehouse	30–50	60–100	150+

Heating/Cooling Loads Annual energy required for space heating and cooling, broken down monthly. The ratio of heating to cooling reveals the building's dominant load and guides passive strategy selection.

Peak Loads Maximum instantaneous heating and cooling demand, measured in kW or W/m². Peak loads size the HVAC equipment. Reducing peak loads through passive design (thermal mass, shading, high-performance envelope) allows smaller, less expensive mechanical systems.

Early-Stage Energy Estimation

For concept design, full EnergyPlus simulation may be premature. Simplified methods:

• Degree-day method: Estimate heating/cooling energy from HDD/CDD and envelope UA-value. Fast but ignores solar gains, internal gains, and thermal mass.
• Lookup tables: Reference EUI values from benchmarking databases (CBECS, CIBSE TM46, EU building typologies) adjusted for climate and building features.
• Shoebox models: Single-zone energy models representing a typical floor section. Capture the essential physics in 1% of the modeling effort.
• Regression surrogates: Train a model on thousands of EnergyPlus runs, then predict EUI from a handful of input parameters in milliseconds.

Parametric Energy Optimization

Key parameters for early-stage energy optimization:

Parameter	Typical range	Impact on EUI
Orientation	0–360°	5–15% (climate-dependent)
WWR (north)	15–60%	5–20%
WWR (south)	15–60%	10–30%
WWR (east/west)	15–40%	10–25%
Wall U-value	0.15–0.50 W/m²K	5–15%
Roof U-value	0.10–0.30 W/m²K	3–10%
Glazing U-value	0.8–3.0 W/m²K	5–20%
Glazing SHGC	0.20–0.60	10–25%
Shading depth	0–2.0 m	5–15%
Infiltration rate	0.1–1.0 ACH	5–15%
Lighting power	4–12 W/m²	10–20%

Passive Design Strategies Validated by Simulation

• Orientation: Long axis east-west maximizes southern exposure (northern hemisphere), improving winter solar gain while allowing effective summer shading.
• Thermal mass: Exposed concrete soffits absorb daytime gains and release heat at night, reducing peak cooling loads by 10–20%.
• Night purge ventilation: Flushing the building with cool night air pre-cools thermal mass, shifting cooling loads and reducing HVAC energy by 15–30%.
• External shading: Horizontal overhangs on south, vertical fins on east/west. Fixed shading design derived from sun path analysis.
• High-performance glazing: Triple glazing with low-e coatings (U < 1.0 W/m²K, SHGC 0.25–0.40) dramatically reduces both heating and cooling loads.
• Daylighting with dimming: Photosensor-controlled dimming of electric lights in daylit zones reduces lighting energy by 40–60%.

---

5. CFD Wind Analysis

Outdoor Wind Comfort — Lawson Criteria

The Lawson comfort criteria classify wind conditions by acceptable activity:

Category	Mean wind speed threshold	Gust equivalent	Acceptable activity
Sitting (long)	< 2.5 m/s	< 4.0 m/s	Outdoor dining, reading
Sitting (short)	< 4.0 m/s	< 6.0 m/s	Bus stops, café terraces
Standing	< 6.0 m/s	< 8.0 m/s	Window shopping, waiting
Walking (leisure)	< 8.0 m/s	< 10.0 m/s	Strolling, walking routes
Business walking	< 10.0 m/s	< 12.0 m/s	Walking to work, commuting
Uncomfortable	10–15 m/s	12–18 m/s	Hair disturbed, clothing flaps
Dangerous	> 15 m/s	> 20 m/s	Structural damage risk

Wind comfort is typically assessed at pedestrian level (1.5 m above ground) for the annual wind climate, reporting the worst-season or annual probability of exceedance.

Computational Wind Tunnel Setup

Domain sizing rules (H = tallest building height):

Boundary	Distance from buildings	Rationale
Inlet	5H upstream	Allow boundary layer to develop
Outlet	15H downstream	Allow wake to dissipate
Lateral walls	5H from edge of model	Prevent blockage effects (< 3%)
Top	5H above tallest point	Prevent artificial acceleration

Blockage ratio (frontal area of buildings / domain cross-section) must be < 3%. Higher blockage artificially accelerates flow around buildings.

Mesh resolution guidelines:

Region	Cell size
Near building surfaces	0.5–1.0 m
Pedestrian level (0–3 m)	0.5–1.0 m
Near-field (within 2H)	1.0–3.0 m
Far-field (beyond 2H)	3.0–10.0 m
Boundary layer refinement	First cell < 0.5 m

Total cell count for a typical urban study: 2–10 million cells.

Turbulence Models

RANS (Reynolds-Averaged Navier-Stokes): Time-averaged equations with turbulence closure models. Solves for mean flow quantities. Computationally affordable.

Model	Strengths	Weaknesses	Use case
Standard k-epsilon	Robust, well-validated	Poor for separation, wakes	Initial screening
Realizable k-epsilon	Better for recirculation	Still struggles with anisotropy	General urban studies
k-omega SST	Excellent near-wall, good separation	Higher cost than k-epsilon	Detailed pedestrian comfort
Spalart-Allmaras	Low cost, good for attached flow	Poor for complex urban geometry	Aerospace, not urban

LES (Large Eddy Simulation): Resolves large-scale turbulent structures directly, models only the smallest scales. Much more accurate for urban flows but 100–1000x more expensive than RANS. Reserved for research and critical safety assessments.

Butterfly (OpenFOAM for Grasshopper) Workflow

1. Create geometry in Rhino/Grasshopper (buildings as closed Breps).
2. Define wind tunnel (domain) using Butterfly components.
3. Set inlet boundary condition: atmospheric boundary layer profile with reference wind speed, direction, roughness length (z0), and reference height.
4. Set mesh parameters: base cell size, refinement levels near buildings.
5. Select turbulence model (k-epsilon or k-omega SST).
6. Run snappyHexMesh for mesh generation.
7. Run simpleFoam (steady-state RANS solver).
8. Post-process: extract wind speed at 1.5 m height plane, map to Lawson categories.

Simplified Wind Analysis for Early Design

Full CFD is expensive and time-consuming. For early design, use:

• Desktop wind assessment: Classify building form using standard typologies (slab, tower, podium-tower, courtyard) and apply known wind behavior patterns.
• Wind comfort rules of thumb: Tall buildings create downwash proportional to their height. Corner acceleration increases with building width. Through-building passages create venturi effects.
• Lawson screening tool: Quick assessment based on building height, width, and surroundings without running CFD.
• Reduced-order models: Pre-computed wind pressure databases for simple shapes.

Natural Ventilation Potential

CFD can assess whether natural ventilation is viable:

• Calculate pressure coefficients (Cp) on building facades from wind simulation.
• Derive pressure difference between windward and leeward openings.
• Estimate ventilation flow rate: Q = Cd * A * sqrt(2 * dP / rho).
• Compare achieved air change rate against minimum ventilation requirements.

Wind-Driven Rain

Wind-driven rain analysis predicts wetting patterns on building facades:

• Combines rainfall intensity with wind speed and direction.
• Identifies facade zones at risk of water penetration.
• Informs material selection, detailing, and drainage design.
• Quantified as catch ratio: rain on vertical surface / rain on horizontal surface.

---

6. Thermal Comfort

Indoor Thermal Comfort

PMV/PPD (Fanger Model) Predicted Mean Vote (PMV) is a steady-state thermal comfort index on a 7-point scale:

PMV value	Thermal sensation
-3	Cold
-2	Cool
-1	Slightly cool
0	Neutral
+1	Slightly warm
+2	Warm
+3	Hot

Predicted Percentage Dissatisfied (PPD) is derived from PMV. At PMV = 0, PPD = 5% (some people are always dissatisfied). ASHRAE 55 requires -0.5 < PMV < +0.5 (PPD < 10%).

PMV inputs:

Parameter	Symbol	Typical indoor range
Air temperature	Ta	18–28°C
Mean radiant temperature	Tr	16–35°C
Air speed	Va	0.05–0.5 m/s
Relative humidity	RH	30–70%
Metabolic rate	Met	1.0–2.0 met (office: 1.1)
Clothing insulation	Clo	0.5–1.5 clo

Adaptive Model For naturally ventilated buildings, the adaptive model (ASHRAE 55 Section 5.4, EN 15251) defines acceptable indoor temperature as a function of outdoor running mean temperature.

Acceptable operative temperature = 17.8 + 0.31 * outdoor running mean temperature (80% acceptability band: +/- 3.5°C; 90% band: +/- 2.5°C).

The adaptive model acknowledges that occupants in naturally ventilated buildings tolerate wider temperature ranges because they have more control (opening windows, adjusting clothing).

Outdoor Thermal Comfort

Universal Thermal Climate Index (UTCI) An equivalent temperature that represents the physiological response of the human body to the outdoor thermal environment.

UTCI range (°C)	Stress category	Thermal perception
> 46	Extreme heat stress	Unbearable
38–46	Very strong heat stress	Very hot
32–38	Strong heat stress	Hot
26–32	Moderate heat stress	Warm
9–26	No thermal stress	Comfortable
0–9	Slight cold stress	Slightly cool
-13–0	Moderate cold stress	Cool
-27–(-13)	Strong cold stress	Cold
-40–(-27)	Very strong cold stress	Very cold
< -40	Extreme cold stress	Extreme cold

Physiological Equivalent Temperature (PET) The air temperature at which the body's heat balance would be the same in a standard indoor reference environment. More intuitive for non-specialists.

Mean Radiant Temperature (MRT) The uniform temperature of an imaginary black enclosure that would result in the same net radiation heat exchange as the actual environment. MRT is often the dominant factor in outdoor comfort and is strongly influenced by:

• Direct solar radiation (sun exposure vs. shade)
• Longwave radiation from surrounding surfaces (hot pavement vs. vegetation)
• Sky view factor (open sky vs. enclosed courtyard)

MRT can be 20–30°C higher in direct sun than in shade. This is why shade trees and canopies are the single most effective microclimate intervention in hot climates.

Microclimate Simulation Workflow

1. Model urban geometry (buildings, ground surfaces, vegetation).
2. Assign surface materials (albedo, emissivity, thermal admittance).
3. Set meteorological boundary conditions from EPW data.
4. Run coupled simulation: radiation (solar + longwave) + CFD (wind field) + energy balance (surface temperatures).
5. Calculate MRT at pedestrian height from radiation field.
6. Compute UTCI or PET at grid of points.
7. Map comfort categories spatially and temporally.

Tools: ENVI-met (most comprehensive urban microclimate tool), Ladybug Tools (MRT and UTCI from weather data and simplified radiation), RayMan, SOLWEIG.

---

7. Acoustic Simulation

Room Acoustics Metrics

Metric	Full name	Target (speech)	Target (music)
RT60	Reverberation time (60 dB decay)	0.4–0.8 s (office)	1.5–2.5 s (concert)
EDT	Early Decay Time	≈ RT60	≈ RT60 (uniform)
C80	Clarity (ratio early/late)	—	-2 to +4 dB
D50	Definition (early/total ratio)	> 0.50	—
STI	Speech Transmission Index	> 0.60 (good)	—
G	Strength (dB re: 10 m free field)	—	0 to +10 dB

RT60 is the most commonly specified acoustic metric. It depends on room volume and total absorption:

Sabine equation: RT60 = 0.161 * V / A

Where V = room volume (m³) and A = total absorption (m² Sabins).

Eyring equation (more accurate for highly absorptive rooms): RT60 = 0.161 * V / (-S * ln(1 - alpha_avg))

Where S = total surface area and alpha_avg = average absorption coefficient.

Material Absorption Coefficients

Material	125 Hz	250 Hz	500 Hz	1 kHz	2 kHz	4 kHz
Concrete (painted)	0.01	0.01	0.02	0.02	0.02	0.03
Brick (unglazed)	0.03	0.03	0.03	0.04	0.05	0.07
Plasterboard on studs	0.29	0.10	0.06	0.05	0.04	0.04
Glass (window)	0.35	0.25	0.18	0.12	0.07	0.04
Timber floor	0.15	0.11	0.10	0.07	0.06	0.07
Carpet (heavy on pad)	0.08	0.24	0.57	0.69	0.71	0.73
Acoustic ceiling tile	0.50	0.70	0.60	0.70	0.70	0.50
Curtain (heavy, draped)	0.07	0.31	0.49	0.75	0.70	0.60
Mineral wool (50 mm)	0.15	0.45	0.70	0.80	0.80	0.80
Perforated metal + absorber	0.40	0.70	0.80	0.85	0.75	0.65
Upholstered seat (occupied)	0.60	0.75	0.85	0.90	0.90	0.85
Open doorway	1.00	1.00	1.00	1.00	1.00	1.00

Ray-Tracing Acoustic Simulation

Pachyderm Acoustics is a Grasshopper plugin for room acoustics simulation:

1. Define room geometry as closed meshes in Rhino.
2. Assign absorption and scattering coefficients to each surface.
3. Place source and receiver points.
4. Run ray-tracing simulation (10,000–100,000 rays).
5. Extract impulse response at each receiver.
6. Compute RT60, EDT, C80, D50, STI from impulse response.
7. Visualize sound pressure level distribution as color map.

Sound Insulation

STC (Sound Transmission Class): Single-number rating for airborne sound insulation of partitions (North American standard, ASTM E413).

STC rating	Performance
25	Normal speech easily understood
30	Loud speech understood, normal speech audible
35	Loud speech audible but not easily understood
40	Loud speech audible as murmur
45	Loud speech not audible
50	Very loud sounds barely heard
55+	Most sounds inaudible

Rw (Weighted Sound Reduction Index): ISO equivalent of STC (ISO 717-1).

Typical constructions:

• Single plasterboard on studs: STC 33–38
• Double plasterboard on studs with insulation: STC 45–50
• Concrete block (200 mm): STC 45–50
• Concrete slab (200 mm): STC 50–55
• Double stud wall with resilient channels: STC 55–60

Outdoor Noise Propagation

Environmental noise from roads, railways, and aircraft is assessed using:

• ISO 9613-2: Attenuation of sound during outdoor propagation (geometric spreading, atmospheric absorption, ground effect, screening by barriers).
• CNOSSOS-EU: Harmonized European noise calculation method.
• Distance attenuation: Point source: -6 dB per doubling of distance. Line source (road): -3 dB per doubling of distance.
• Barrier effect: A solid barrier provides 5–15 dB insertion loss depending on path length difference.
• Noise mapping: Color-coded maps of facade or free-field noise levels, typically at 4 m height. Required for Environmental Impact Assessments.

---

8. Ladybug Tools Ecosystem

Overview

Ladybug Tools is an open-source collection of plugins for Grasshopper (Rhino) that provides comprehensive environmental analysis capabilities for building and urban design. It connects parametric geometry to validated simulation engines.

Ladybug (Weather Data and Outdoor Analysis)

Core capabilities:

• EPW import and visualization: Dry bulb temperature, radiation, wind, humidity as hourly heatmaps, bar charts, and statistical summaries.
• Sun path diagram: 3D stereographic or orthographic sun path with hourly sun positions, analemmas, and sun vectors for any location.
• Wind rose: Frequency and speed distribution by direction for any analysis period.
• Radiation rose: Directional radiation distribution for facade orientation studies.
• Shadow study: Calculate sunlight hours on test surfaces with context geometry.
• Outdoor comfort: UTCI, PET, and other outdoor comfort models from weather data.
• Sky dome: Visualize sky luminance/radiance distribution (Tregenza sky patches).
• Direct sun hours: Mesh-based analysis of solar access hours.
• View analysis: Assess visual exposure from points to targets.

Honeybee (Building Simulation)

Honeybee-Radiance (daylight):

• Create Radiance models from Grasshopper geometry.
• Assign material modifiers (plastic, glass, trans, BSDF).
• Define sensor grids and views.
• Run point-in-time illuminance, daylight factor, annual daylight (sDA/ASE/UDI).
• Glare analysis with Evalglare.
• Parametric blind/shade studies.

Honeybee-Energy (thermal/energy):

• Create thermal zones (rooms) from geometry.
• Assign program types (loads and schedules by space type).
• Define construction sets (opaque + glazing assemblies).
• Set HVAC systems (ideal air, detailed systems).
• Run EnergyPlus simulations.
• Parse results: energy by end use, zone temperatures, comfort.

Key Honeybee concepts:

• Room: A closed volume representing one thermal zone.
• Face: A planar surface bounding a room (wall, floor, roof/ceiling).
• Aperture: A window or skylight within a face.
• Door: An opaque or glass door within a face.
• Shade: An external shading surface (overhang, fin, context building).
• Modifier: Radiance material properties (reflectance, transmittance).
• Construction: Energy material layers (conductivity, density, specific heat).
• ProgramType: Combination of people, lighting, equipment, ventilation, setpoints.
• ConstructionSet: Collection of constructions for all face types.

Butterfly (CFD)

• Wraps OpenFOAM for use in Grasshopper.
• Creates computational domain (wind tunnel) around building geometry.
• Generates block-structured mesh with refinement regions.
• Sets boundary conditions (inlet velocity profile, outlet, walls).
• Runs steady-state RANS solver (simpleFoam).
• Extracts velocity and pressure fields for visualization.
• Current status: Less actively maintained than Ladybug/Honeybee. For production CFD work, consider standalone OpenFOAM or commercial tools.

Dragonfly (Urban Scale)

• Urban weather generator (UWG): Modifies rural EPW data to account for urban heat island effect, producing urban-specific weather data.
• District energy modeling: Aggregates building energy models for neighborhood-scale analysis.
• Urban geometry: Creates building footprint extrusions from GIS data.
• REopt integration: Optimizes distributed energy resources (PV, storage, CHP).
• Urban microclimate: Couples with other tools for outdoor comfort mapping.

Version Compatibility

Component	Current version	Rhino compatibility	Python
Ladybug	1.8.x	Rhino 7, Rhino 8	IronPython / CPython
Honeybee	1.8.x	Rhino 7, Rhino 8	IronPython / CPython
Honeybee-Radiance	1.66.x	Rhino 7, Rhino 8	CPython 3.7+
Honeybee-Energy	1.104.x	Rhino 7, Rhino 8	CPython 3.7+
Butterfly	0.0.x	Rhino 6, Rhino 7	IronPython
Dragonfly	1.8.x	Rhino 7, Rhino 8	CPython 3.7+

Note: Ladybug Tools 1.x (LBT) uses the Pollination installer, which manages Radiance, EnergyPlus, and OpenStudio dependencies automatically.

Installation

1. Install Rhino 7 or 8.
2. Download Pollination Grasshopper installer from pollination.cloud.
3. Run installer — it sets up Ladybug Tools, Radiance 5.4+, EnergyPlus 23.1+, and OpenStudio 3.7+ automatically.
4. Alternatively, install from Food4Rhino and manually configure engine paths.
5. Verify installation: drop LB Versioner component on Grasshopper canvas.

Common Workflows

Workflow	Components used	Engine	Output
Daylight factor	HB Model, HB Modifier, HB Grid, HB DF	Radiance	DF spatial map
Annual daylight (sDA/ASE)	HB Model, HB Annual Daylight	Radiance	sDA, ASE percentages
Point-in-time illuminance	HB Model, HB Point-in-Time	Radiance	Illuminance grid
Glare analysis	HB Model, HB Glare, Evalglare	Radiance	DGP images
Energy model (annual)	HB Room, HB Program, HB Construction, HB Energy	EnergyPlus	EUI, loads, temps
Outdoor solar radiation	LB Direct Sun Hours, LB Radiation	Built-in	kWh/m² map
Sun path + shadow study	LB Sun Path, LB Shadow Study	Built-in	Hours of sunlight
Wind rose	LB Wind Rose	Built-in	Directional wind plot
Outdoor comfort (UTCI)	LB UTCI, LB MRT	Built-in	UTCI spatial map
Parametric optimization	Above + Galapagos / Wallacei	Various	Pareto-optimal designs

Integration with Optimization

Galapagos (single-objective): Built into Grasshopper. Uses genetic algorithm. Connect performance metric to fitness input. Good for single-metric optimization (e.g., minimize EUI). Limited to one objective.

Wallacei (multi-objective): NSGA-2 algorithm for Grasshopper. Handles 2–10+ objectives simultaneously. Produces Pareto front of non-dominated solutions. Includes built-in analytics for exploring the solution space. Recommended for performance-driven design where multiple conflicting objectives exist.

Colibri (design space exploration): Records every iteration of a parametric study (inputs and outputs) to a structured data file. Pairs with Design Explorer web app for parallel coordinates visualization. Useful for understanding parameter sensitivity before launching optimization.

---

9. Simulation Quality Assurance

Validation Strategies

Every simulation result should be viewed with appropriate skepticism. Validation builds confidence in results:

• Analytical validation: Compare simulation output against known analytical solutions for simple cases (e.g., steady-state heat flow through a wall, illuminance from a point source). If the engine cannot reproduce analytical results, something is fundamentally wrong.
• Empirical validation: Compare simulation output against measured data from real buildings or controlled experiments. ASHRAE Standard 140 (BESTEST) provides standardized test cases for energy simulation engines.
• Comparative validation: Compare results from two different simulation engines on the same model. Agreement builds confidence; disagreement demands investigation.
• Sensitivity analysis: Vary uncertain inputs (material properties, schedules, weather data) and observe the effect on outputs. If results are highly sensitive to an uncertain input, invest effort in getting that input right.

Mesh Independence Studies

For any spatially discretized simulation (daylight grids, CFD meshes, FEM meshes):

1. Run the simulation with a coarse mesh.
2. Refine the mesh (halve cell size) and rerun.
3. Compare results at key locations.
4. If results change by more than 5%, refine further.
5. Repeat until results converge (< 2% change between successive refinements).
6. Report the final mesh resolution and the convergence behavior.

Mesh independence is non-negotiable for CFD and strongly recommended for daylight grids. A simulation on an inadequate mesh is not worth the computation time.

Result Interpretation Guidelines

• Order of magnitude: Do results fall within the expected range for this building type and climate? An EUI of 500 kWh/m²/yr for a well-insulated office is wrong.
• Spatial distribution: Do daylight, temperature, and wind maps show physically plausible patterns? Symmetric buildings should produce symmetric results.
• Temporal patterns: Do energy loads follow expected seasonal patterns? Cooling should peak in summer, heating in winter (for most climates).
• Comparative ranking: Even if absolute values are uncertain, the relative ranking of design options is usually reliable. Use simulation for comparison, not prediction.

Common Errors and Diagnostic Checklist

Error	Symptom	Fix
Missing surfaces	Unrealistic heat loss / daylight	Check for gaps in geometry
Wrong boundary conditions	Adiabatic exterior walls	Verify face types (wall/floor/roof)
Interior walls as exterior	Excessive heating/cooling loads	Check adjacencies between rooms
Wrong weather file	Results don't match climate	Verify EPW location matches project
Schedule errors	Overnight loads in unoccupied building	Review occupancy and HVAC schedules
Unit confusion	Values 10x too high or low	Check W vs kW, m² vs ft², °C vs °F
Insufficient Radiance bounces	Dark interiors, low daylight values	Increase -ab parameter
Coarse CFD mesh	Smoothed-out wind patterns	Refine mesh, check independence
Wrong material properties	Unrealistic surface temperatures	Verify conductivity, reflectance, etc.
Unconverged CFD	Oscillating residuals	Improve mesh quality, check BC setup

When to Trust Results vs. Use Engineering Judgment

Trust simulation results when:

• The model has been validated against analytical or empirical benchmarks.
• Mesh independence has been demonstrated.
• Material properties and boundary conditions are well-characterized.
• Results are within expected physical ranges.
• Sensitivity analysis shows results are robust to uncertain inputs.

Use engineering judgment when:

• Input data is highly uncertain (occupant behavior, future climate).
• The simulation engine has known limitations for the specific case.
• Results contradict well-established physical principles.
• The model uses significant simplifications (shoebox model, ideal HVAC).
• Computational constraints prevented adequate mesh resolution.

Reporting Simulation Results

Professional simulation reports should include:

1. Executive summary: Key findings and recommendations in non-technical language.
2. Model description: Geometry, materials, boundary conditions, assumptions.
3. Simulation setup: Engine, version, settings, mesh details, weather data.
4. Results: Spatial maps, charts, tables with clear units and legends.
5. Interpretation: What the results mean for design decisions.
6. Limitations: What the simulation does not capture, uncertainty bounds.
7. Recommendations: Specific design actions supported by the analysis.
8. Appendices: Detailed input data, convergence plots, sensitivity results.

Never present simulation results without stating the assumptions and limitations. A number without context is more dangerous than no number at all.

---

Quick Reference: Simulation Engine Selection

Analysis need	Recommended engine	Ladybug tool	Accuracy level
Daylight factor	Radiance	Honeybee	High
Annual daylight (sDA/ASE)	Radiance (3/5-phase)	Honeybee	High
Glare analysis	Radiance + Evalglare	Honeybee	High
Solar radiation	Built-in (Tregenza method)	Ladybug	Medium-High
Sun hours / shadow	Built-in (ray intersection)	Ladybug	High
Building energy (annual)	EnergyPlus	Honeybee	High
Building energy (early stage)	Degree-day / lookup	Manual	Low-Medium
Outdoor wind comfort	OpenFOAM	Butterfly	Medium-High
Indoor thermal comfort (PMV)	EnergyPlus	Honeybee	High
Outdoor thermal comfort (UTCI)	Built-in	Ladybug	Medium
Urban microclimate	ENVI-met	External	High
Room acoustics	Pachyderm	External (GH)	Medium-High
Outdoor noise	ISO 9613 / CadnaA	External	Medium-High

---

Key Formulas Reference

Daylight Factor: DF = (Ei / Eo) * 100% Where Ei = indoor illuminance, Eo = unobstructed outdoor diffuse horizontal illuminance.

Sabine reverberation time: RT60 = 0.161 * V / A Where V = volume (m³), A = total absorption (m² Sabins).

Natural ventilation flow rate: Q = Cd * A * sqrt(2 * deltaP / rho) Where Cd = discharge coefficient (~0.6), A = opening area, deltaP = pressure difference, rho = air density.

Solar heat gain through glazing: Qsolar = A_glazing * SHGC * I_solar Where SHGC = solar heat gain coefficient, I_solar = incident solar radiation (W/m²).

Thermal transmittance (U-value): U = 1 / R_total R_total = R_si + sum(d_i / k_i) + R_se Where R_si/R_se = surface resistances, d = thickness, k = conductivity.

UTCI (simplified approximation): UTCI ≈ f(Ta, Tr, va, RH) — computed via 6th-order polynomial regression.

Reynolds number: Re = rho * v * L / mu Where rho = air density, v = velocity, L = characteristic length, mu = dynamic viscosity. Re > 4000 indicates turbulent flow (for external aerodynamics, Re >> 10^6 always turbulent).