Chapter 13
Computer-Aided, Multiphysics and Optimization Methods
Bachelors Program in Electrical Engineering
Electrical Machine Design
•
Bachelor
Computer-Aided, Multiphysics and Optimization Methods
Introduction and Context
Electrical-machine design has traditionally relied on analytical magnetic-circuit models, empirical coefficients, prototype testing, and iterative engineering judgment. These methods remain essential because they are fast, transparent, and useful during preliminary sizing. However, modern machines operate under increasingly demanding conditions: high power density, rapid transients, acoustic constraints, thermal limits, reduced material usage, variable-speed operation, and stringent efficiency standards. Analytical calculations alone may not adequately represent magnetic saturation, leakage flux, slotting, eddy-current loss, three-dimensional end effects, rotor stress, or temperature-dependent material properties.
Computer-aided engineering extends classical design by combining analytical models, numerical field solutions, thermal networks, structural calculations, optimization algorithms, and experimental validation. The central principle is not that simulation replaces engineering judgment. Rather, simulation provides a hierarchy of models with different levels of fidelity, while engineering judgment determines which model is appropriate, how its assumptions should be checked, and whether its predictions agree with physical evidence.
A modern machine-design workflow is therefore an evidence-based loop:
Define requirements and constraints.
Build a low-order analytical model.
Select relevant electromagnetic, thermal, mechanical, and acoustic physics.
Generate and solve a numerical model.
Verify numerical accuracy and validate physical accuracy.
Perform sensitivity, uncertainty, and tolerance studies.
Optimize the design under realistic operating conditions.
Manufacture or prototype critical designs and compare measurements with predictions.
perform mesh, time-step, sensitivity, uncertainty, and tolerance studies;
formulate single-objective and multi-objective machine optimization problems;
interpret Pareto-optimal solutions and engineering trade-offs;
identify sources of modelling error and validate predictions against measurements;
apply good practices for reproducibility, data governance, and responsible use of automated design tools.
1. Model Hierarchy in Electrical-Machine Design
1.1 Analytical and lumped-parameter models
An analytical model expresses machine behaviour using idealized geometry and physical approximations. Examples include magnetic-circuit calculations, winding-function models, dq-axis equivalent circuits, thermal-resistance networks, and lumped rotor-dynamic models. Such models are computationally inexpensive and reveal parameter dependence.
For a simple magnetic path, the reluctance is
ℜ = l/(μA)
where ℜ is reluctance in A/Wb, l is magnetic-path length in metres, μ is permeability in H/m, and A is cross-sectional area in m². If several sections carry the same flux in series, their reluctances are added. The corresponding flux is approximately
Φ = NI/ℜtotal
where N is the number of turns and I is current. This relation assumes a one-dimensional path, negligible leakage and fringing, and a known or approximately constant permeability. Ferromagnetic saturation violates the constant-μ assumption; a nonlinear B–H curve must then be used.
Analytical models are especially valuable for design-space exploration, initial dimensions, parameter identification, and checking whether a finite-element result is physically plausible. They should not be discarded after a finite-element model is created.
1.2 Numerical field models
A numerical model divides the machine domain into elements and approximates field variables within each element. It can represent complicated geometry, nonlinear materials, motion, local saturation, and distributed losses. The price is increased computational cost and the need for careful verification.
The principal model classes are:
Model
Typical strength
Typical limitation
2-D magnetostatic
Fast flux-density and torque estimation
Cannot represent important axial effects
2-D transient electromagnetic
Ripple, harmonics, dynamic currents
Requires time stepping and motion treatment
3-D electromagnetic
End leakage, skew, axial current paths
High computational cost
Thermal network
Rapid temperature estimates
Local hot spots may be missed
Thermal finite element
Detailed temperature distribution
Requires credible loss and boundary data
Structural finite element
Stress, deformation, modal response
Sensitive to contact, preload, and material data
Coupled multiphysics
Interaction among domains
Complex convergence and calibration
1.3 Verification and validation
Verification asks, “Was the mathematical model solved correctly?” It includes checking element quality, mesh convergence, solver residuals, energy balance, symmetry, and limiting cases. Validation asks, “Does the model represent the real machine adequately?” It requires comparison with independent measurements, such as no-load back-emf, load torque, efficiency, winding resistance, temperature, vibration, or natural frequency.
A small numerical residual does not prove physical validity. A highly converged model with incorrect material data or boundary conditions can still be wrong.
2. Electromagnetic Finite-Element Analysis
2.1 Governing equations and formulation
Finite-element electromagnetic analysis begins with Maxwell’s equations. In magnetoquasistatic machine analysis, displacement current is often neglected because the characteristic dimensions and frequencies make conduction and magnetic effects dominant. A common two-dimensional formulation uses the out-of-plane magnetic vector potential A_z. For a conducting region, a representative equation is
∇·(ν∇Az) + σ(∂Az/∂t + v·∇Az) = Js
Here A_z is magnetic vector potential in Wb/m, ν = 1/μ is reluctivity in m/H, σ is electrical conductivity in S/m, v is conductor velocity in m/s, and J_s is imposed source-current density in A/m². The exact sign and motion terms depend on coordinate convention and formulation. In a stationary, nonconducting, magnetostatic region, the equation reduces to
∇·(ν∇Az) = Js
The magnetic flux density is obtained from the vector potential. In Cartesian two-dimensional geometry,
Bx = ∂Az/∂y, By = −∂Az/∂x
where B is in tesla. The field solution is converted into flux linkage, induced voltage, force, and torque. Torque may be calculated by Maxwell stress, virtual work, co-energy, or an appropriate volume-force method. These methods should agree within numerical error when applied consistently.
For a linear magnetic system, magnetic co-energy is
W′m = 1/2 LI²
and torque at fixed current may be estimated as
T = ∂W′m/∂θ
where L is inductance in H, I is current in A, and θ is rotor position in radians. In saturated machines, co-energy must be computed from the nonlinear constitutive relation rather than by the simple 1/2LI² expression.
2.2 Geometry, materials, and boundary conditions
A credible model requires more than a drawing. The analyst must specify:
actual air gaps, slot openings, bridges, fillets, laminations, magnets, conductors, and end regions;
nonlinear steel B–H data, including rolling-direction effects where relevant;
electrical conductivity and temperature dependence;
permanent-magnet recoil permeability, remanence, and temperature coefficient;
lamination stacking factor and loss model;
excitation, circuit connection, phase sequence, and rotor position;
periodic, anti-periodic, or external boundary conditions.
Symmetry reduces computation, but only if the excitation and geometry possess the assumed symmetry. Incorrect symmetry can suppress real harmonics and produce misleading torque or force results.
2.3 Loss calculation
Copper loss at a specified temperature is commonly estimated as
PCu = m Irms² R(T)
where m is the number of identical phase windings, I_rms is phase RMS current in A, and R(T) is phase resistance in Ω at temperature T. For a conductor with reference resistance R₀ at T₀,
R(T) = R₀[1 + α(T − T₀)]
where α is the temperature coefficient in K⁻¹. This linear relation is suitable over moderate temperature ranges; material-specific data should be used at high temperature.
Iron loss often contains hysteresis and eddy-current contributions. A simplified empirical form is
pfe = kh f B̂n + ke f² B̂²
where p_fe is specific loss, f is frequency, B̂ is peak flux density, k_h and k_e are fitted coefficients, and n is an empirical exponent. This expression is not universal: waveform, minor loops, stress, manufacturing, and rotational flux can require improved loss models. Loss separation should therefore be calibrated against appropriate material-test data.
The physics are interdependent. Current produces electromagnetic torque and copper loss. Loss raises temperature. Temperature changes resistance, magnet strength, insulation life, lubricant viscosity, and mechanical dimensions. Electromagnetic forces create stress, vibration, and acoustic noise. Deformation changes the air gap and therefore changes the electromagnetic field.
A one-way sequential coupling is often adequate for preliminary design:
solve electromagnetic fields;
map averaged or distributed losses to the thermal model;
calculate temperatures;
update temperature-dependent material properties;
calculate mechanical stress and deformation.
A fully coupled iterative solution is required when feedback is strong, such as high-speed rotors, large thermal drift, narrow air gaps, or temperature-sensitive permanent magnets.
3.2 Thermal modelling
For a lumped thermal network, a node energy balance is
Cth dT/dt = Ploss − (T − Tamb)/Rth
where C_th is thermal capacitance in J/K, T is node temperature in °C or K, T_amb is ambient temperature, P_loss is heat input in W, and R_th is thermal resistance in K/W. At steady state,
T = Tamb + Ploss Rth
This result assumes constant loss and thermal resistance. In reality, convection coefficients, radiation, contact resistance, and losses may vary with speed and temperature.
Thermal finite-element analysis solves the heat equation. In a stationary isotropic region,
−∇·(k∇T) = q̇
where k is thermal conductivity in W/(m·K) and q̇ is volumetric heat generation in W/m³. Boundary conditions may include prescribed temperature, convection q = h(T − T_amb), radiation, or contact conductance. The convection coefficient h must be justified by correlation, experiment, or fluid simulation; assigning an arbitrary large value can conceal thermal problems.
3.3 Mechanical and rotor-dynamic analysis
Electromagnetic forces can be calculated from the Maxwell stress tensor. In vacuum or air, the traction associated with magnetic field B is represented by
Tm = (1/μ₀)(BB − 1/2 B²I)
where μ₀ is the permeability of free space, B is the magnetic-flux-density vector, I is the identity tensor, and BB denotes the dyadic product. Integrating traction over a suitable surface gives force. The surface must be placed in a region where the field is accurately resolved and should not cross poorly modelled interfaces.
For a rotating component, centrifugal stress depends on speed, density, geometry, and restraint. Structural analysis must include interference fits, adhesive layers, retaining sleeves, keyways, thermal expansion, and manufacturing tolerances where relevant. A rotor may satisfy static stress limits yet fail because of resonance or fatigue. Modal and harmonic-response analyses are therefore important for high-speed machines.
4. Numerical Quality and Model Credibility
4.1 Mesh and time-step convergence
A mesh-convergence study repeats the calculation with progressively finer meshes and monitors quantities of interest such as average torque, peak tooth flux density, force, or loss. Local refinement is required near air gaps, slot openings, magnet edges, bridges, and narrow laminations. Peak values at mathematically sharp corners may be mesh-dependent; averaged or physically integrated quantities are often more meaningful.
For transient analysis, the time step must resolve electrical frequency, switching events, rotor motion, and the fastest relevant thermal or mechanical dynamics. A practical study halves the time step and checks whether the target output changes by an acceptable amount. “Acceptable” must be defined before the study; for example, a designer may require torque change below 1%, but the threshold depends on the application.
4.2 Verification checklist
A professional verification record includes:
dimensional and unit checks;
comparison with analytical limiting cases;
mesh and time-step refinement;
residual and convergence histories;
energy and power balance;
independence from irrelevant boundary distance;
correct symmetry and periodicity;
sensitivity to solver tolerances;
conservation checks for current, flux, and torque.
4.3 Validation hierarchy
Validation should use quantities not employed to tune the model. A useful sequence is:
Stage
Evidence
Typical purpose
Material
B–H curve, loss test, thermal conductivity
Calibrate constitutive data
Component
Coil resistance, inductance, magnet flux
Check subassemblies
Machine
Back-emf, torque, efficiency
Check integrated electromagnetic model
Operating test
Temperature, vibration, acoustic spectrum
Check multiphysics behaviour
Qualification
Endurance and overload tests
Assess lifetime and robustness
Agreement should be expressed with uncertainty. If predicted temperature is 120 °C and measured temperature is 115 °C, the difference is 5 °C; its significance depends on sensor accuracy, repeatability, operating-point uncertainty, and spatial mismatch between sensor and predicted hot spot.
5. Sensitivity, Uncertainty, and Tolerance Studies
5.1 Sensitivity analysis
Sensitivity measures how an output changes when an input changes. A local normalized sensitivity of output y to parameter x_i is
Si = (xi/y)(∂y/∂xi)
It is dimensionless and applicable near a nominal design where the derivative is meaningful. A finite-difference estimate is
∂y/∂xi ≈ [y(xi + Δxi) − y(xi − Δxi)]/(2Δxi)
provided Δx_i is sufficiently small and numerical noise is controlled. Large parameter interactions make one-at-a-time studies inadequate; design-of-experiments or global methods such as variance-based sensitivity analysis are then preferable.
Typical influential variables include air-gap length, magnet remanence, slot fill factor, conductor resistivity, lamination loss coefficients, cooling coefficient, eccentricity, and rotor position. Sensitivity results guide measurement effort and tolerance allocation.
5.2 Uncertainty propagation
If independent input variables x_i have small uncertainties with standard deviations σ_xi, the first-order variance approximation is
σy² ≈ ∑(∂y/∂xi)² σxi²
For correlated inputs, covariance terms must be added:
σy² ≈ ∑∑(∂y/∂xi)(∂y/∂xj)Cov(xi,xj)
These approximations assume local linearity and appropriately characterized input distributions. Monte Carlo simulation is more flexible: sample inputs from their distributions, solve the model for each sample, and estimate the output distribution. The number of samples should be justified by convergence of statistics or tail probabilities, not selected arbitrarily.
5.3 Manufacturing tolerances
A nominal optimum may be unsuitable if small dimensional errors cause large performance variation. Tolerance studies should include air-gap variation, magnet placement, skew angle, eccentricity, winding resistance, lamination thickness, shaft runout, and cooling variation. Robust design seeks acceptable performance over the tolerance distribution rather than only at nominal dimensions.
6. Optimization of Electrical Machines
6.1 Problem formulation
An optimization problem has design variables x, objective functions f(x), equality constraints h(x) = 0, and inequality constraints g(x) ≤ 0. A general statement is
minimize F(x) = [f₁(x), f₂(x), …, fk(x)] subject to gj(x) ≤ 0, hl(x) = 0, xL ≤ x ≤ xU
Design variables may include slot dimensions, pole arc, magnet thickness, air gap, turns, conductor size, lamination grade, cooling parameters, and control settings. Objectives may include mass, cost, average torque, torque ripple, loss, temperature, acoustic noise, and material usage.
The objectives and constraints must be evaluated at all required operating points. Optimizing only rated-point efficiency can produce poor field-weakening performance, overload heating, or unacceptable starting behaviour.
6.2 Multi-objective optimization and Pareto fronts
In multi-objective design, improving one objective may worsen another. A design x dominates design y if it is no worse in every objective and strictly better in at least one. A Pareto-optimal design is not dominated by any feasible alternative. The collection of such designs is the Pareto front.
A weighted-sum method forms
J(x) = ∑ wi fi^*(x)
where w_i are nonnegative weights and f_i* are normalized objectives. Normalization is essential because kilograms, watts, degrees Celsius, and monetary cost cannot be added meaningfully without scaling. Weighted sums may miss non-convex portions of a Pareto front. Evolutionary algorithms, ε-constraint methods, and surrogate-assisted optimization are alternatives.
6.3 Surrogate models and computational efficiency
A full three-dimensional transient multiphysics solve may take minutes or hours, making thousands of evaluations impractical. A surrogate approximates the expensive model using sampled designs. Polynomial response surfaces, Gaussian-process models, radial-basis functions, and neural networks are possible choices. Surrogates require training data, cross-validation, and a defined domain of validity. They must not be extrapolated confidently beyond the sampled design space.
An effective workflow is to use analytical models for broad screening, two-dimensional finite elements for refinement, and three-dimensional coupled simulations for a small number of finalists. Final candidates are then tested against tolerances and operating cycles.
7. Worked Example: Thermal Constraint in a Copper Winding
Consider a phase winding with resistance 0.18 Ω at 20 °C. The temperature coefficient of copper is approximated as α = 0.00393 K⁻¹. The phase RMS current is 12 A, and there are three identical phases. The estimated thermal resistance from winding to ambient is 0.42 K/W. Ambient temperature is 30 °C. Estimate copper loss and steady winding temperature, assuming the thermal resistance is constant.
At temperature T, phase resistance is
R(T) = 0.18[1 + 0.00393(T − 20)]
The loss is
PCu = 3(12²)R(T) = 432R(T)
Because resistance depends on temperature, solve the steady-state relation iteratively:
T = 30 + 0.42PCu
First estimate at 20 °C gives P_Cu = 432 × 0.18 = 77.76 W and T = 62.66 °C. At 62.66 °C,
R = 0.18[1 + 0.00393(62.66 − 20)] ≈ 0.2102 Ω
Thus P_Cu ≈ 432 × 0.2102 = 90.8 W, giving T ≈ 30 + 0.42 × 90.8 = 68.1 °C. Repeating once gives R ≈ 0.2141 Ω, P_Cu ≈ 92.5 W, and T ≈ 68.9 °C. The converged estimate is approximately 69 °C.
This example illustrates three important points. First, copper loss must be evaluated at operating temperature. Second, a thermal network is useful for rapid screening but does not identify local hot spots. Third, the result is only as reliable as the assumed thermal resistance. A validation test should measure winding temperature under the same current, speed, cooling, and ambient conditions.
8. Illustrative Case Study: Traction Permanent-Magnet Motor
Background
This illustrative case concerns a 75 kW interior permanent-magnet traction motor. The design team must meet average torque, efficiency, maximum winding temperature, rotor stress, and acoustic-noise requirements while limiting active material cost.
Problem
An initial analytical design meets rated torque but exhibits high torque ripple and a predicted rotor temperature close to the magnet demagnetization limit. The machine also has a narrow air gap, so manufacturing eccentricity may produce significant unbalanced magnetic force.
Analysis
A two-dimensional transient electromagnetic model is first created using measured steel B–H data and temperature-dependent magnet properties. The model is verified through mesh refinement and comparison of no-load back-emf with an analytical winding-function model. Distributed copper and iron losses are transferred to a thermal finite-element model. The thermal model includes estimated housing convection and measured contact conductance from a similar machine.
A structural model evaluates centrifugal stress, electromagnetic radial force, thermal expansion, and modal frequencies. A tolerance study varies air gap, magnet remanence, rotor eccentricity, and cooling coefficient. The results show that magnet remanence and cooling coefficient strongly affect temperature, while eccentricity dominates the second spatial harmonic of radial force.
Approach and solution
The optimization variables are magnet-arc ratio, barrier dimensions, slot-opening width, and selected skew angle. Objectives are normalized loss, torque ripple, mass, and peak temperature. Constraints require rated torque, allowable stress, maximum temperature, and minimum first critical speed margin. A design-of-experiments study generates training points, followed by a surrogate-assisted multi-objective search. The Pareto set reveals that a small increase in active material reduces temperature substantially but provides diminishing returns beyond a particular magnet thickness.
The selected design is not the nominal minimum-loss design. It is a robust Pareto solution with slightly higher nominal loss but lower sensitivity to eccentricity and cooling variation. Three prototypes are manufactured. Measurements include back-emf, torque ripple, efficiency map, winding temperature, rotor vibration, and order-tracked acoustic noise.
Outcome
The prototype agrees with simulation within the project’s pre-established uncertainty bands for average torque and winding temperature. Torque ripple is higher than predicted because the manufactured slot opening and magnet placement differ from nominal geometry. Updating the tolerance model improves prediction of the measured harmonic content. The final design uses a manufacturing-control limit on magnet placement and an inspection procedure for rotor eccentricity.
Learning points
A nominal finite-element optimum is not necessarily a robust production design.
Electromagnetic, thermal, and mechanical analyses must exchange physically consistent data.
Validation must include manufacturing variation and operating conditions.
The Pareto front supports engineering decisions; it does not automatically select the final product.
Measurement should be planned early enough to expose model deficiencies before design release.
9. Professional Practice, Ethics, and Contemporary Tools
Computer-aided design tools can conceal assumptions behind attractive colour plots. Engineers must retain traceable records of geometry versions, material sources, solver settings, mesh statistics, convergence criteria, scripts, and post-processing definitions. Results should be reproducible by another competent engineer.
AI and machine learning can assist with surrogate modelling, anomaly detection, automated meshing, design-space exploration, and extraction of patterns from test data. However, an AI-generated design or prediction is not evidence of correctness. Training-data bias, extrapolation, hidden correlations, intellectual-property restrictions, and undocumented software changes must be controlled. Safety-critical decisions require human review, independent checks, and physical validation.
Ethically, designers must not report only favourable operating points, hide uncertainty, or present a converged numerical result as a certified performance guarantee. Regulatory requirements may govern insulation systems, temperature rise, rotating-part containment, electromagnetic compatibility, efficiency, and environmental compliance. The exact standard depends on machine type, market, and application; compliance must be demonstrated using the applicable current edition of the relevant standard.
10. Best-Practice Design Procedure
A disciplined procedure is as follows:
Define requirements, duty cycle, ambient conditions, tolerances, and acceptance criteria.
Build and check an analytical baseline.
Select the lowest-complexity model that can answer the design question.
Document materials, excitations, boundaries, and assumptions.
Verify mesh, time step, solver convergence, and conservation quantities.
Couple losses, temperatures, forces, and deformation using consistent coordinates and units.
Evaluate all important operating points, not only rated operation.
Quantify sensitivity and uncertainty before allocating design margins.
Optimize with constraints and realistic manufacturing tolerances.
Validate finalists using independent measurements and update the model when discrepancies arise.
Preserve a version-controlled computational record and design rationale.
Chapter Summary
Computer-aided electrical-machine design combines analytical models, finite-element field solutions, thermal and structural analysis, optimization, and experimental evidence. Analytical models provide speed and physical insight; finite-element methods represent nonlinear and geometrically complex fields; multiphysics coupling captures interactions among loss, temperature, force, deformation, and vibration. Verification establishes numerical credibility, whereas validation establishes correspondence with the physical machine. Sensitivity and uncertainty studies identify influential variables and quantify risk. Multi-objective optimization exposes trade-offs through Pareto-optimal designs, while robust optimization accounts for manufacturing and operating variation. Professional practice requires traceability, appropriate model fidelity, honest reporting of uncertainty, and physical validation. AI and surrogate methods are useful accelerators but do not replace engineering accountability.
Key Terms
Analytical model: A mathematical approximation based on simplified physical assumptions.
Finite-element analysis: Numerical solution of field or structural equations over a discretized domain.
Magnetic vector potential: A field variable from which magnetic flux density is derived.
Magnetostatic analysis: Electromagnetic analysis in which fields and excitations are treated as time-invariant.
Transient analysis: Simulation that resolves time-dependent electrical, magnetic, thermal, or mechanical behaviour.
Multiphysics: Simultaneous consideration of interacting physical domains.
Verification: Demonstration that numerical equations have been solved correctly.
Validation: Demonstration that a model represents physical behaviour adequately.
Mesh convergence: Reduction of discretization error demonstrated through mesh refinement.
Thermal network: Lumped representation of heat capacities and thermal resistances.
Sensitivity: Measure of output response to variation of an input parameter.
Uncertainty propagation: Estimation of output uncertainty from uncertain inputs.
Tolerance study: Evaluation of performance under manufacturing or assembly variation.
Design variable: Quantity deliberately changed during a design study.
Constraint: Requirement that a candidate design must satisfy.
Pareto-optimal design: Feasible design not dominated by another design in all objectives.
Surrogate model: Approximation of an expensive computational model.
Robust design: Design whose performance remains acceptable under variation.
Co-energy: Energy-based electromagnetic quantity useful for calculating force or torque.
Model calibration: Adjustment of model parameters using suitable experimental evidence.
Review Questions
Short Answer Questions
Distinguish verification from validation in machine simulation.
Why are analytical models retained after a finite-element model is developed?
State two conditions under which a two-dimensional model may be inadequate.
What is meant by a Pareto-optimal design?
Why must copper resistance be evaluated at operating temperature?
Long Answer Questions
Explain the complete workflow for coupled electromagnetic–thermal–mechanical machine design.
Describe the main sources of error in finite-element electromagnetic analysis and methods for controlling them.
Explain how mesh convergence, time-step convergence, and experimental validation differ.
Formulate a multi-objective optimization problem for a permanent-magnet motor and identify suitable constraints.
Discuss the advantages and limitations of surrogate models in computational machine design.
Analytical / Application Questions
A three-phase winding has 0.20 Ω phase resistance at 20 °C, carries 10 A RMS, and has a thermal resistance of 0.50 K/W to a 25 °C ambient. Using α = 0.00393 K⁻¹, estimate the steady temperature by iteration.
A design output depends on two independent inputs. The local derivatives are ∂y/∂x₁ = 4 and ∂y/∂x₂ = −2, while σ_x₁ = 0.10 and σ_x₂ = 0.20. Use first-order variance propagation to estimate σ_y.
Propose a verification and validation plan for a high-speed interior permanent-magnet motor whose principal concerns are magnet temperature, rotor stress, and acoustic noise.
Multiple Choice Questions with Answers
Question 1
Which statement best distinguishes verification from validation?
A. Verification compares a model with a prototype, while validation checks mesh quality.
B. Verification checks whether the equations are solved correctly, while validation checks agreement with physical behaviour.
C. Verification is performed only after manufacturing, while validation is performed before simulation.
D. Verification and validation are identical terms.
Correct Answer: B. Verification checks whether the equations are solved correctly, while validation checks agreement with physical behaviour.
Question 2
For a linear magnetic system, which expression represents magnetic co-energy?
A. W′_m = LI
B. W′_m = L/I
C. W′_m = 1/2 LI²
D. W′_m = 2LI²
Correct Answer: C. W′_m = 1/2 LI²
Question 3
Which region usually requires particularly careful mesh refinement in a rotating electrical machine?
A. A distant air region with nearly uniform field
B. The air gap and narrow geometric features
C. A perfectly uniform outer boundary only
D. A region containing no electromagnetic field
Correct Answer: B. The air gap and narrow geometric features
Question 4
If copper resistance is 0.20 Ω at 20 °C, α = 0.004 K⁻¹, and the temperature is 70 °C, what is the approximate resistance?
A. 0.160 Ω
B. 0.200 Ω
C. 0.240 Ω
D. 0.280 Ω
Correct Answer: C. 0.240 Ω
Question 5
Which quantity is most directly associated with uncertainty propagation?
A. Input standard deviations and model sensitivities
B. Only the nominal output value
C. The colour scale of a field plot
D. The number of geometric dimensions alone
Correct Answer: A. Input standard deviations and model sensitivities
Question 6
A design is Pareto-optimal when it is
A. the heaviest feasible design
B. not dominated by another feasible design in all objectives
C. the design with every objective equal to zero
D. always the design with minimum cost
Correct Answer: B. not dominated by another feasible design in all objectives
Question 7
What is a major limitation of a weighted-sum multi-objective method?
A. It cannot use numerical objectives
B. It always violates constraints
C. It may miss non-convex portions of the Pareto front
D. It requires no objective normalization
Correct Answer: C. It may miss non-convex portions of the Pareto front
Question 8
In a steady lumped thermal model with constant loss, the temperature rise is
A. P_loss/R_th
B. P_loss R_th
C. R_th/P_loss
D. P_loss + R_th
Correct Answer: B. P_loss R_th
Question 9
Which practice most improves reproducibility of a computational design study?
A. Saving only the final colour plot
B. Changing solver settings without recording them
C. Recording geometry, materials, mesh, solver settings, and version-controlled scripts
D. Reporting only the best operating point
Correct Answer: C. Recording geometry, materials, mesh, solver settings, and version-controlled scripts
Question 10
Why should a nominal optimization result be subjected to a tolerance study?
A. To remove all physical constraints
B. To determine whether manufacturing variation causes unacceptable performance changes
C. To guarantee that simulation is unnecessary
D. To replace experimental validation completely
Correct Answer: B. To determine whether manufacturing variation causes unacceptable performance changes